Abstract
Recent advances in alkaline earth (Ae) metal hydrogenation catalysis have broadened the spectrum of potential catalysts to include candidates from the main group, providing a sustainable alternative to the commonly used transition metals. Although Ae-amides have already been demonstrated to catalyze hydrogenation of imines and alkenes, a lucid understanding of how different metal/ligand combinations influence the catalytic activity is yet to be established. In this article, we use linear scaling relationships and molecular volcano plots to assess the potential of the Ae metal-based catalysts for the hydrogenation of alkenes. By analyzing combinations of eight metals (mono-, bi-, tri-, and tetravalent) and seven ligands, we delineate the impact of metal-ligand interplay on the hydrogenation activity. Our findings highlight that the catalytic activity is majorly determined by the charge and the size of the metal ions. While bivalent Ae metal cations delicately regulate the binding and the release of the reactants and the products, respectively, providing the right balance for this reaction, ligands play only a minor role in determining their catalytic activity. We show how volcano plots can be utilized for the rapid screening of prospective Ae catalysts to establish a guideline to achieve maximum activity in facilitating the hydrogenation process.
Similar content being viewed by others
1 Introduction
Hydrogenation of unsaturated bonds is one of the most fundamental transformations in chemistry finding broad applications at every scale of chemical production [1]. With the conventional hydrogenation catalysts utilizing rare, expensive, and often toxic transition metals, there is a great incentive for chemists to find cheaper and environmentally friendly alternatives. In this context, the recent development of the alkaline earth (Ae) metal amides as effective hydrogenation catalysts for various alkenes and imines is especially significant [2,3,4,5,6]. In fact, these catalysts facilitate direct hydrogenation of alkenes under mild conditions while successfully tackling crucial issues pertinent to hydrogenation protocols such as suppressing polymerization and tolerating various functional groups. Strikingly, Ae metal complexes also catalyze highly selective alkene transfer hydrogenation, using 1,4-cyclohexadiene as a reducing agent [7]. Transfer hydrogenation using alternative hydrogen sources is very attractive since it does not require hazardous pressurized H\(_2\) or elaborate experimental set-ups [8]. Very recently, efficient hydrogenation with a wide range of substrates was also achieved with ligand-free metallic barium [9]. These ground-breaking results confirm the versatility of simple Ae metal complexes in catalysis, establishing a sustainable alternative to transition metal catalysis for industrial applications.
Initially, the catalytic activity of the Ae[(NSiMe\(_3\))]\(_2\) catalysts were demonstrated with their ability to hydrogenate a broad range of alkenes and imines [2,3,4, 7]. While Mg is found to be inactive for hydrogenation of alkenes, heavier Ae-amides exhibit high catalytic activity. Interestingly, alkali (Ak) amides are also capable of facilitating this transformation although their activity is much lower than the Ae-amides. Combined experimental/computational studies by Bauer et al. firmly established the mechanism of such Ae-catalysed hydrogenation processes proceeding through a metal hydride mediated route.[2, 3] Interestingly, further improvements on activity and substrate scope were brought upon by introducing bulkier amide ligands, e.g. N(TRIP)\(_2\) and N(TRIP)(DIPP) (where TRIP=SiiPr\(_3\) and DIPP=2,6-(iPr)\(_2\)C\(_6\)H\(_3\)) instead of N(SiMe\(_3\))\(_2\) [4]. Such observations are rationalized in terms of increased concentration of in situ generated catalytically active smaller metal hydride species due to large ligand size.
Despite the well-established mechanistic picture, a clear-cut understanding of how different Ae metal-ligand combinations influence the catalytic activity is currently lacking. Such knowledge is crucial for the systematic development of Ae-hydrogenation catalysis. One particular tool that can particularly aid in such an understanding of the overall trends in catalytic behavior is molecular volcano plots. Over the past few years, our research group has been developing and utilizing these plots to screen homogeneous catalysts for several important chemical processes [10,11,12,13,14,15,16,17,18]. Molecular volcano plots predict the performance of a catalyst (in terms of turn over frequencies or a particular energetic criterion) based on an easily computed descriptor variable. These plots are constructed after post-processing the linear free energy scaling relationships (LFESRs) [19,20,21] obtained between the descriptor and the relative stability of the reaction intermediates and transition states. The most promising catalysts are then easily identified by visually inspecting their location in the plot (appearing near the volcano top or near the plateau). Inspired by this recent development in Ae metal catalysis, here we investigate main group metal-catalyzed hydrogenation of alkenes using molecular volcano plots to uncover the influence of the metal/ligands and their interplay on the energetics of the catalytic cycle.
2 Results and Discussion
Mechanism The first step to construct a molecular volcano plot is to settle the mechanism of the corresponding catalytic cycle (Scheme 1). Bauer et al. already investigated the mechanism of direct hydrogenation of styrene by Ca(N(SiMe\(_3\)))\(_2\) through DFT computations (Fig. 1) [2]. First, the precatalyst 1 exchanges an amide ligand for a hydride ligand to yield 2, which is the entry point into the catalytic cycle. Hydrogenation begins with the coordination of the alkene substrate to the catalyst leading to the formation of 3. The Ae-H bond then inserts into styrene to produce the benzylmetal intermediate 4 via TS1. The precomplex, 5, precedes the heterolytic cleavage of H\(_2\) through TS2 leading to 6 which upon dissociation releases ethylbenzene and regenerates the catalyst, 2. Figure 2 illustrates a representative free energy profile for the catalytic hydrogenation process (along with the key intermediates and transition states) using a model catalyst, Ca(NMe\(_2\))\(_2\).
LFESRs The prerequisite to constructing a volcano plot is to establish LFESRs which ascertain that the free energies of the various catalytic cycle intermediates are correlated. The LFESRs were determined by analyzing a set of 54 catalysts produced from the combinations of 8 metal cations (Na\(^+\), Mg\(^{2+}\), Ca\(^{2+}\), Sr\(^{2+}\), Ba\(^{2+}\), Sc\(^{3+}\), Ti\(^{4+}\), Zr\(^{4+}\)) and 7 small anionic ligands of varying connecting atom identities (H\(^-\), F\(^-\), BF\(_2\) \(^-\), Me\(^-\), NMe\(_2\) \(^-\), PMe\(_2\) \(^-\), OMe\(^-\)) [22]. Apart from the mono/bivalent metals, tri and tetravalent metals are chosen additionally as a part of a systematic investigation of the catalytic cycle to reveal general trends in reactivity [23].
To ensure general trends in the scaling relationships, the entry point to the catalytic cycle is considered to be a monomeric bisligated species. Therefore, the effects of ligand over-coordination, catalyst aggregations, and the possible deactivation channels are not considered here.
We computed the full catalytic cycle for all 54 catalysts using standard DFT computations (see Computational Details) and the free energies associated with the intermediates and transition states (I/TSs) are estimated relative to the reference state (\(\Delta \)G\(_{RRS}\)), 2, which is the entry point to the catalytic cycle. Based on the quality of the linear correlations, \(\Delta \)G\(_{RRS}\)(4) is chosen as the descriptor variable for the volcano analyses (see Scheme 2). The computed LFESRs along with the quality of the linear fits, established with \(\Delta \)G\(_{RRS}\)(4) as a descriptor, are depicted in Fig. 3.
Kinetic volcano plot The LFESRs are post-processed to obtain the kinetic volcano plot [10, 11] in Fig. 4. The descriptor \(\Delta \)G\(_{RRS}\)(4) is plotted along the x-axis and the y-axis corresponds to the free energy required to complete the most difficult reaction step, kinetics determining step (kds). The plot can be divided into three areas, which correspond to different limiting reaction steps. The limiting reaction in region-I is the binding of styrene (2\(\rightarrow \)3), while the region-III is limited by the release of ethylbenzene (6\(\rightarrow \)2) from the catalyst. In line with Sabatier’s principle terminology, regions-I and III can be readily interpreted as the weak and strong binding sides that characterize the right and left slopes of the volcano, respectively. Region-II (the plateau of the volcano), on the other hand, corresponds to the most kinetically balanced situation in which the most energetic step involves the barrier associated with the heterolytic H\(_2\) cleavage(5\(\rightarrow \)TS2). Catalysts occupying this region are especially interesting as they nearly fulfill Sabatier’s definition of an ideal catalyst.
TOF volcano plot Instead of kds, which involves reaction energies between consecutive two intermediates/TSs, the catalytic activity could be expressed in terms of more realistic turnover frequencies (TOFs) [24] by examining the free energy profiles of the catalysts under the framework of Eyring’s transition state theory (TST). We utilized LFESRs in Fig. 5 to determine hypothetical free energy profiles for a range of descriptor values (\(-115\) to 10 kcal/mol) and these were used as input to derive theoretical TOF values corresponding to each descriptor value. When plotted with a logarithmic y-axis, the TOF volcano takes a more or less similar shape as that of the kinetic volcano (see Fig. 4). The TOF volcano in Fig. 5 is best interpreted using the energy span model [25] (microkinetic analysis) which associates different regions of the plot with a certain pair of intermediate (TDI: turnover determining intermediate) and transition states (TDTS: turnover determining TS), or, two intermediates for a thermodynamically driven association or dissociation process leading to the largest barrier in a free energy profile. The overall energy span associated with a catalytic cycle is typically defined by the following equations,
where \(T_i\) and \(I_j\) are the Gibbs free energies of the ith TS and jth intermediates in the profile, respectively, and \(\Delta G_R\) is the Gibbs free energy of the reaction. Importantly, the TDI and TDTS are not necessarily the highest and lowest states, nor do they have to be adjoined as a single step. The energy difference between the TDI and the TDTS is the apparent activation energy of the entire catalytic cycle determining the catalytic efficiency.
Applying the energy span model, the TOF volcano plot can be divided into three regions depending on the change in TDI/TDTS as a function of the descriptor. The right side region (equivalent to region-I in Fig. 4) corresponds to TS1 as TDTS and 2 as TDI for descriptor values \(\Delta \)G\(_{RRS}\)(4) \(> -16\) kcal/mol. For the top region (equivalent to region-II in Fig. 4) with descriptor values lying between \(-16\) to \(-29\) kcal/mol, TDI and TDTS change to 4 and TS2, respectively. Finally, for \(\Delta \)G\(_{RRS}\)(4) \(< -29\) kcal/mol (equivalent to region-III in Fig. 4) the TDTS and TDI become 2 and 6 [26]. Thus, for the candidates appearing on the right region, the catalysts remain strongly bound to the substrate resulting in a high energy penalty for the hydride insertion step. On the other hand, the catalysts falling on the left slope of the volcano have a overly stabilized intermediate, 6 which makes the product release step strongly endergonic. Approaching the top of the volcano from either side corresponds to balancing the energy requirements of these two steps and increasing the catalytic activity.
Understanding the Influence of Metal and Ligands To examine the potential of the catalysts for the hydrogenation of styrene, we plotted each catalyst candidate according to their descriptor variable (\(\Delta \)G\(_{RRS}\)(4)) on the kinetic as well as on the TOF volcano. The position of a catalyst determines its limiting reaction steps (Fig. 4) or a pair of TDTS and TDI (Fig. 5) within an energy span framework. Nearly all the catalysts considered in the work appear either in region-II or III (TDTS/TDI: 6/2 or TS2/4) of the kinetic (TOF) volcano plot. Their relative positioning suggests that they are roughly clustered based on the metal oxidation state and the total charge of the catalyst. Ae/Ak catalysts with a formally bi/mono positive medium to large size metal ions (Na\(^+\), Ca\(^{2+}\), Sr\(^{2+}\), and Ba\(^{2+}\)) appear at the top of the volcano (i.e., region-II) indicative of their higher activity for styrene hydrogenation. They have nearly ideal kinetic profiles and are limited by the molecular hydrogen cleavage step. Note that the experimentally inactive Mg catalysts separate from the clusters of the other Ae candidates and lie closer to the intersecting area I and II. A careful examination of the structure of the intermediate 4 reveals that, unlike other Ae cations, Mg\(^{2+}\) binds to the benzyl anion in more of a monodentate fashion (involving the benzylic carbon) presumably due to the smaller size of the Mg\(^{2+}\) cation. In fact, the other bigger metal cations (Ca, Sr and Ba) in the Ae catalysts family for which the benzyl anion offers a bidentate binding mode (involving both the ipso and the benzylic carbon centers) making the corresponding benzylmetal complexes more stabilized compared to their Mg-counterpart. As shown in Fig. 5, Mg candidates are anticipated to be somewhat limited by the formation of TS1 (TDTS/TDI: TS1/2), the hydride insertion step, but still predicted as fairly active (although less active compared to other Ae catalysts ). Overall, these findings imply a slightly different behavior of Mg compared to other Ae metals.[27]
Highly charged cations (tri- or tetravalent Sc, Ti and Zr) lead to a more strongly stabilized (more negative \(\Delta \)G\(_{RRS}\)(4)) benzylmetal intermediate than the Ae cations. Accordingly, all tri- or tetravalent Sc, Ti, and Zr-containing candidates fall on the region-III being limited by the release of the product. All of these catalysts overstabilize complex 6 and thus show minimal catalytic activity. A crucial structural feature of 6 is the cation-\(\pi \) interaction in which the central metal cation is bound to the phenyl ring of ethylbenzene. Further insights into this interaction can be obtained by estimating the noncovalent interaction (NCI) index of the intermediate 6 (see SI for computational details of NCI analysis). Figure 6 illustrates the NCI isosurfaces computed for a few representative metal-ligand combinations. All complexes exhibit evidence for attractive interaction between the cation and the product. For a given ligand (Me\(^-\)), the cation-\(\pi \) interaction strengthens with an increased charge of the central metal cation, and thus, highly charged cations lead to the strongest interactions rendering the release of the product more difficult. The corresponding interactions for bivalent metals are nearly similar (see Fig. S6 and S7). In fact, these metals delicately regulate the cation-\(\pi \) interaction so that the energetic cost of association of styrene and dissociation of ethylbenzene from the catalyst balance each other, which is key to achieve high TOF values. The reason for the different behavior of Mg might be related to its small size preventing multiple non-bonding contacts essential for a strong cation-\(\pi \) interaction. Overall, the relative activity of the catalyst candidates increases in the order of tetravalent < trivalent < bi or, monovalent metals. Experimentally, there are substantial activity differences among Ca, Sr and Ba catalysts although based on the catalytic cycle in Scheme 1b, these three metals are predicted to exhibit comparable reactivity. These results likely indicate that the observed reactivity trend in experiment stems from several side reactions pertinent to the catalytic cycle.
In contrast to the metals, the effect of the ligands on the catalytic activity is much less significant. The electron-withdrawing (EW) ligands tend to strengthen the metal-benzyl interaction hence stabilizing 4 (more negative value of the descriptor). For instance, F/BF\(_2\) stabilizes 4 more than electron-donating (ED) ligands such as NMe\(_2\)/OMe. Thus, given a particular metal ion, EW ligand bearing candidates progressively appear towards the left-hand side of the volcano. While the manifestation of this effect is only minor for the bi/monovalent cations, this is amplified when the ligands are combined with highly positive cations (Ti\(^{4+}\) and Zr\(^{4+}\)). The latter exhibits fairly different activity, while replacing F with NMe\(_2\) producing a substantial shift toward the weak binding side (i.e. right along the x-axis). As shown in Figs. S4 and S5, the corresponding NCI analyses are consistent with these results. These observations suggest that it might not be an effective strategy to improve the intrinsic activity of the Ae catalysts through ligand tuning, often a preferred option in TM-based catalysis. Nevertheless, it should be noted that the choice of ligands could have an impact on the overall performance of the catalytic cycle by influencing the side reactions or catalyst-aggregation [4]. The metal charge and size are the most relevant factors affecting the Ae catalysts’ activity towards hydrogenation.
3 Conclusions
In summary, the performance of main group metal catalysts for the hydrogenation of styrene was examined using molecular volcano plots. Our findings reveal that the activity of the metal is essentially determined by the formal charge and size of the metal cation. Catalysts with highly charged cations engage in stronger interactions with styrene mainly through cation-\(\pi \) interactions that limits the release of the product. The hard Mg\(^{2+}\) cation tends to disfavor the hydride insertion. Finally, medium to large mono- and bivalent metals provide the ideal balance for alkene hydrogenation. The influence of the ligands is important when they are combined with tri/tetravalent metals. The intrinsic activity of the Ae catalysts is not significantly affected by the choice of the ligands and hence, most of these are predicted to have comparable performance towards styrene hydrogenation.
4 Computational Details
The geometries of all catalytic cycle intermediates and TSs were optimized at the M06 [28, 29]/def2-SVPD [30] level in implicit benzene solvent using Gaussian16, A.03 [31]. The M06 hybrid functional was proven to be quite accurate for main group thermochemistry and kinetics and our previous works demonstrated that this level of theory is adequate to support the experimental findings of imine and alkene hydrogenations by Ae amides.[2, 3] An analysis of the harmonic vibrational frequencies was performed to ascertain the nature of each optimized structures either as a minimum (no imaginary frequency) or a transition state (one imaginary frequency). Single point energies were computed on the M06 geometries at the PBE0-dDsC/TZ2P [32,33,34,35,36,37] level as implemented in ADF [38]. Free energy corrections were obtained at the M06/def2-SVPD level using the rigid-rotor harmonic oscillator model within the Goodvibes [39] program developed by Paton and Funes-Ardoiz. Solvation corrections for the reported free energy values were obtained using COSMO-RS solvation model at PBE0-dDsC/TZ2P level in benzene [40]. The turnover frequencies (TOF) were calculated at 298.15 K and a concentration of 1 M using AUTOF program developed by Uhe, Kozuch, and Shaik [41,42,43]. The NCI plots were computed with the NCIPLOT program, starting from the M06 wave functions of the optimized geometries [44]. Additional computational details for NCI analysis are provided in the SI.
References
de Vries JG, Elsevier CJ et al (2007) Handbook of homogeneous hydrogenation. WeinheimWiley-VCH, New York
Bauer H, Alonso M, Fischer C, Rösch B, Elsen H, Harder S (2018) Angew Chem Int Ed 57(46):15177
Bauer H, Alonso M, Färber C, Elsen H, Pahl J, Causero A, Ballmann G, De Proft F, Harder S (2018) Nat Catal 1(1):40
Martin J, Knüpfer C, Eyselein J, Färber C, Grams S, Langer J, Thum K, Wiesinger M, Harder S (2020) Angew Chem Int Ed 59(23):9102
Stephan DW (2018) Nature 553:160
H Bauer, S Harder (2020) Early main group metal catalysis: concepts and reactions. Wiley, New York. pp 175–199
Bauer H, Thum K, Alonso M, Fischer C, Harder S (2019) Angew Chem Int Ed 58(13):4248
Wang D, Astruc D (2015) Chem Rev 115(13):6621
Stegner P, Färber C, Zenneck U, Knüpfer C, Eyselein J, Wiesinger M, Harder S (2021) Angew Chem 133(8):4298
Wodrich MD, Busch M, Corminboeuf C (2016) Chemical Science 7(9):5723
Wodrich MD, Sawatlon B, Solel E, Kozuch S, Corminboeuf C (2019) ACS Catal 9(6):5716
Busch M, Wodrich MD, Corminboeuf C (2018) ChemCatChem 10(7):1592
Wodrich MD, Busch M, Corminboeuf C (2018) Helvet Chim Acta 101(9):e1800107
Sawatlon B, Wodrich MD, Corminboeuf C (2020) Org Lett 22(20):7936
Busch M, Wodrich MD, Corminboeuf C (2017) ACS Catal 7(9):5643
Meyer B, Sawatlon B, Heinen S, von Lilienfeld OA, Corminboeuf C (2018) Chem Sci 9(35):7069
Wodrich MD, Fabrizio A, Meyer B, Corminboeuf C (2020) Chem Sci 11(44):12070
Wodrich MD, Sawatlon B, Busch M, Corminboeuf C (2021) Acc Chem Res 54(5):1107
Greeley J (2016) Annu Rev Chem Biomol Eng 7:605
Man IC, Su HY, Calle-Vallejo F, Hansen HA, Martínez JI, Inoglu NG, Kitchin J, Jaramillo TF, Nørskov JK, Rossmeisl J (2011) ChemCatChem 3(7):1159
Calle-Vallejo F, Loffreda D, Koper MT, Sautet P (2015) Nat Chem 7(5):403
The metal-ligand combinations Ti-BF\(_2\) and Ti-PMe\(_2\) are removed from the set of the catalysts due to infeasible geometries present in the catalytic cycle
For Ti/Zr-catalyzed hydrogenation of alkenes, see (a) Journal of Organometallic Chemistry 382, 69 (1990) (b) Journal of Organometallic Chemistry, 462, 191 (1993) (c) Journal of Catalysis, 205, 294 (2002) (d) Chemistry-A European Journal, 26, 2803 (2020). Note that, the mechanisms associated with these reactions are likely different compared to Ae metal-catalyzed hydrogenation considered here. Since the volcano plots are constructed based on a particular mechanism (see Figure 1a), we refrain from making a direct comparison between the performances of these complexes with those of Ae metal catalysts
Kozuch S, Martin JM (2012) ACS Catal 2(12):2787
Kozuch S, Shaik S (2011) Acc Chem Res 44(2):101
2 and 6 as TDTS and TDI, respectively, imply a thermodynamically driven product release
The experimentally observed lack of activity for Mg is due to some other reasons, presumably not associated with the catalytic cycle considered here
Zhao Y, Truhlar DG (2008) Theoret Chem Acc 120(1–3):215
Zhao Y, Truhlar DG (2008) Acc Chem Res 41(2):157
Weigend F, Ahlrichs R (2005) Phys Chem Chem Phys 7(18):3297
M.J. Frisch, G.W. Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, G Scalmani, V Barone, G.A. Petersson, H Nakatsuji, X Li, M Caricato, A.V. Marenich, J Bloino, B.G. Janesko, R Gomperts, B Mennucci, H.P. Hratchian, J.V. Ortiz, A.F. Izmaylov, J.L. Sonnenberg, D Williams-Young, F Ding, F Lipparini, F Egidi, J Goings, B Peng, A Petrone, T Henderson, D Ranasinghe, V.G. Zakrzewski, J Gao, N Rega, G Zheng, W Liang, M Hada, M Ehara, K Toyota, R Fukuda, J Hasegawa, M Ishida, T Nakajima, Y Honda, O Kitao, H Nakai, T Vreven, K Throssell, J.A. Montgomery, Jr., J.E. Peralta, F Ogliaro, M.J. Bearpark, J.J. Heyd, E.N. Brothers, K.N. Kudin, V.N. Staroverov, T.A. Keith, R Kobayashi, J Normand, K Raghavachari, A.P. Rendell, J.C. Burant, S.S. Iyengar, J Tomasi, M Cossi, J.M. Millam, M Klene, C Adamo, R Cammi, J.W. Ochterski, R.L. Martin, K Morokuma, O Farkas, J.B. Foresman, D.J. Fox. Gaussian 16 Revision A.03 (2016). Gaussian Inc. Wallingford CT
Van Lenthe E, Baerends EJ (2003) J Comput Chem 24(9):1142
Perdew JP, Burke K, Ernzerhof M (1996) Phys Rev Lett 77(18):3865
Adamo C, Barone V (1999) J Chem Phys 110(13):6158
Steinmann SN, Corminboeuf C (2010) J Chem Theory Comput 6(7):1990
Steinmann SN, Corminboeuf C (2011) J Chem Phys 134(4):44117
Steinmann SN, Corminboeuf C (2011) J Chem Theory Comput 7(11):3567
Te Velde GT, Bickelhaupt, Baerends EJ, Fonseca Guerra C, van Gisbergen SJ, Snijders JG, Ziegler T (2001) J Comput Chem 22(9): 931
Funes-Ardoiz I, R Paton R (2019) https://doi.org/10.5281/zenodo.595246
Klamt A (2011) Wiley Interdiscipl Rev 1(5):699
Kozuch S, Shaik S (2006) J Am Chem Soc 128(10):3355
Kozuch S, Shaik S (2008) J Phys Chem A 112(26):6032
Uhe A, Kozuch S, Shaik S (2011) J Comput Chem 32(5):978
Contreras-García J, Johnson ER, Keinan S, Chaudret R, Piquemal JP, Beratan DN, Yang W (2011) J Chem Theory Comput 7(3):625
Acknowledgements
S.D. and C.C. acknowledge the EPFL and the Swiss National Science Foundation (Grant No. 200020_175496) for financial support. M.A. thanks the Fund for Scientific Research Flanders (FWO-12F4416N) for a postdoctoral fellowship and the Vrije Universiteit Brussel (VUB) for financial support. B.D.T. thanks FWO for the PhD fellowship received (11F4920N). For the NCI analyses, Tier2 computational resources and services were provided by the Shared ICT Services Centre Funded by the Vrije Universiteit Brussel, the Flemish Supercomputer Center (VSC) and FWO. The authors thank Ruben Laplaza for TOC graphic.
Funding
Open Access funding provided by EPFL Lausanne.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing financial interests.
Supporting Information
Additional computational details and figures, NCI analysis, and free energy profiles for all catalysts (PDF) Computed electronic energies, free energy corrections, solvent corrections, and total free energies (XLS) Cartesian coordinates of all computed structures in XYZ format (ZIP)
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Das, S., Tobel, B.D., Alonso, M. et al. Uncovering the Activity of Alkaline Earth Metal Hydrogenation Catalysis Through Molecular Volcano Plots. Top Catal 65, 289–295 (2022). https://doi.org/10.1007/s11244-021-01480-7
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11244-021-01480-7