Semiclassical Transition State Theory (SCTST) can incorporate the non-separable coupling between degrees of freedom (DOFs) of reactive systems and include the effects of reaction path curvature and anharmonicity, as well as quantum tunneling contributions in the rate constant. The rate expression is derived by relying on a perturbative expansion for the vibrational energy which makes it possible to express the semiclassical cumulative reaction probability in a convenient way, without any further assumptions about the separability of the DOFs.[1-4] The main goal of the talk is the extension of this semiclassical methodology to systems of increasing dimensionality upon computation of the densities of states for the reactants and the transition state by means of a convenient parallel implementation of the Wang-Landau algorithm.[5-6] The new strategy is implemented into two codes, “paradensum” and “parsctst”, which are currently distributed with the open source MultiWell program suite for chemical kinetics.[7] The needed input information is just the reaction barrier height, the normal mode frequencies, and the anharmonic force constants, which are routinely calculated by suitable electronic structure packages. After describing the codes and demonstrating their computational accuracy and efficiency, the new implementation is applied to estimate the rate constant of the proton transfer isomerization of the 2,4,6-tri-tert-butylphenyl to 3,5-di-tert-butylneophyl, a reaction involving 145 degrees of freedom and showing a clear tunneling regime below 250K.[6] References [1] Miller, W.H. J. Chem. Phys., 1975, 62, 1899 – 1906. [2] Miller, W.H. Faraday Discuss. Chem. Soc. 1977, 62, 40 – 46. [3] Miller, W.H.; Hernandez, R.; Handy, N.C.; Jayatilaka, D.; Willetts, A. Chem. Phys. Letters 1990, 172, 62 – 68. [4] Hernandez, R.; Miller, W. H. Chem. Phys. Lett. 1993, 214 (2), 129 – 136. [5] Aieta, C.; Gabas, F.; Ceotto, M. J. Phys. Chem. A 2016, 120 (27), 4853 – 4862. [6] Aieta, C.; Gabas, F.; Ceotto, M. J. Chem. Theory Comput. 2019, 15, 2142 − 2153. [7] Barker, J.R.; Nguyen, T.L.; Stanton, J.F.; Aieta, C.; Ceotto, M.; Gabas, F.; Kumar, T.J.D.; Li, C.G.L.; Lohr, L.L.; Maranzana, A.; Ortiz, N.F.; Preses, J.M.; Simmie, J.M.; Sonk, J.A.; Stimac, P.J.; MultiWell-2017 Software Suite; J.R. Barker, University of Michigan, Ann Arbor, Michigan, USA, 2019; http://clasp-research.engin.umich.edu/multiwell/.
Parallel Implementation of Semiclassical Transition State Theory and its application to high-dimensional tunneling reactions / C. Aieta, F. Gabas, M. Ceotto. ((Intervento presentato al convegno International Workshop on Quantum Atomic and Molecular Tunneling Systems (QAMTS) tenutosi a Borovets nel 2019.
Parallel Implementation of Semiclassical Transition State Theory and its application to high-dimensional tunneling reactions
C. AietaPrimo
;F. GabasSecondo
;M. CeottoUltimo
2019
Abstract
Semiclassical Transition State Theory (SCTST) can incorporate the non-separable coupling between degrees of freedom (DOFs) of reactive systems and include the effects of reaction path curvature and anharmonicity, as well as quantum tunneling contributions in the rate constant. The rate expression is derived by relying on a perturbative expansion for the vibrational energy which makes it possible to express the semiclassical cumulative reaction probability in a convenient way, without any further assumptions about the separability of the DOFs.[1-4] The main goal of the talk is the extension of this semiclassical methodology to systems of increasing dimensionality upon computation of the densities of states for the reactants and the transition state by means of a convenient parallel implementation of the Wang-Landau algorithm.[5-6] The new strategy is implemented into two codes, “paradensum” and “parsctst”, which are currently distributed with the open source MultiWell program suite for chemical kinetics.[7] The needed input information is just the reaction barrier height, the normal mode frequencies, and the anharmonic force constants, which are routinely calculated by suitable electronic structure packages. After describing the codes and demonstrating their computational accuracy and efficiency, the new implementation is applied to estimate the rate constant of the proton transfer isomerization of the 2,4,6-tri-tert-butylphenyl to 3,5-di-tert-butylneophyl, a reaction involving 145 degrees of freedom and showing a clear tunneling regime below 250K.[6] References [1] Miller, W.H. J. Chem. Phys., 1975, 62, 1899 – 1906. [2] Miller, W.H. Faraday Discuss. Chem. Soc. 1977, 62, 40 – 46. [3] Miller, W.H.; Hernandez, R.; Handy, N.C.; Jayatilaka, D.; Willetts, A. Chem. Phys. Letters 1990, 172, 62 – 68. [4] Hernandez, R.; Miller, W. H. Chem. Phys. Lett. 1993, 214 (2), 129 – 136. [5] Aieta, C.; Gabas, F.; Ceotto, M. J. Phys. Chem. A 2016, 120 (27), 4853 – 4862. [6] Aieta, C.; Gabas, F.; Ceotto, M. J. Chem. Theory Comput. 2019, 15, 2142 − 2153. [7] Barker, J.R.; Nguyen, T.L.; Stanton, J.F.; Aieta, C.; Ceotto, M.; Gabas, F.; Kumar, T.J.D.; Li, C.G.L.; Lohr, L.L.; Maranzana, A.; Ortiz, N.F.; Preses, J.M.; Simmie, J.M.; Sonk, J.A.; Stimac, P.J.; MultiWell-2017 Software Suite; J.R. Barker, University of Michigan, Ann Arbor, Michigan, USA, 2019; http://clasp-research.engin.umich.edu/multiwell/.File | Dimensione | Formato | |
---|---|---|---|
QAMTS2019_AIETA.pdf
accesso aperto
Tipologia:
Pre-print (manoscritto inviato all'editore)
Dimensione
330.18 kB
Formato
Adobe PDF
|
330.18 kB | Adobe PDF | Visualizza/Apri |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.