A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction.
The inclusion of Quantum Mechanical (QM) effects such as zero point energy (ZPE) and tunneling in simulations of chemical reactions, especially in the case of light atom transfer, is an important problem in computational chemistry. In this respect, the hydrogen exchange reaction and its isotopic var...
| Main Authors: | , , , , , , |
|---|---|
| Format: | Journal article |
| Language: | English |
| Published: |
2013
|
| _version_ | 1797085014156902400 |
|---|---|
| author | Suleimanov, Y de Tudela, R Jambrina, P Castillo, J Sáez-Rábanos, V Manolopoulos, D Aoiz, F |
| author_facet | Suleimanov, Y de Tudela, R Jambrina, P Castillo, J Sáez-Rábanos, V Manolopoulos, D Aoiz, F |
| author_sort | Suleimanov, Y |
| collection | OXFORD |
| description | The inclusion of Quantum Mechanical (QM) effects such as zero point energy (ZPE) and tunneling in simulations of chemical reactions, especially in the case of light atom transfer, is an important problem in computational chemistry. In this respect, the hydrogen exchange reaction and its isotopic variants constitute an excellent benchmark for the assessment of approximate QM methods. In particular, the recently developed ring polymer molecular dynamics (RPMD) technique has been demonstrated to give very good results for bimolecular chemical reactions in the gas phase. In this work, we have performed a detailed RPMD study of the H + H(2) reaction and its isotopologues Mu + H(2), D + H(2) and Heμ + H(2), at temperatures ranging from 200 to 1000 K. Thermal rate coefficients and kinetic isotope effects have been computed and compared with exact QM calculations as well as with quasiclassical trajectories and experiment. The agreement with the QM results is good for the heaviest isotopologues, with errors ranging from 15% to 45%, and excellent for Mu + H(2), with errors below 15%. We have seen that RPMD is able to capture the ZPE effect very accurately, a desirable feature of any method based on molecular dynamics. We have also verified Richardson and Althorpe's prediction [J. O. Richardson and S. C. Althorpe, J. Chem. Phys., 2009, 131, 214106] that RPMD will overestimate thermal rates for asymmetric reactions and underestimate them for symmetric reactions in the deep tunneling regime. The ZPE effect along the reaction coordinate must be taken into account when assigning the reaction symmetry in the multidimensional case. |
| first_indexed | 2024-03-07T02:03:13Z |
| format | Journal article |
| id | oxford-uuid:9e17a877-be88-4acf-b0a8-8979e4a82820 |
| institution | University of Oxford |
| language | English |
| last_indexed | 2024-03-07T02:03:13Z |
| publishDate | 2013 |
| record_format | dspace |
| spelling | oxford-uuid:9e17a877-be88-4acf-b0a8-8979e4a828202022-03-27T00:47:39ZA ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction.Journal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:9e17a877-be88-4acf-b0a8-8979e4a82820EnglishSymplectic Elements at Oxford2013Suleimanov, Yde Tudela, RJambrina, PCastillo, JSáez-Rábanos, VManolopoulos, DAoiz, FThe inclusion of Quantum Mechanical (QM) effects such as zero point energy (ZPE) and tunneling in simulations of chemical reactions, especially in the case of light atom transfer, is an important problem in computational chemistry. In this respect, the hydrogen exchange reaction and its isotopic variants constitute an excellent benchmark for the assessment of approximate QM methods. In particular, the recently developed ring polymer molecular dynamics (RPMD) technique has been demonstrated to give very good results for bimolecular chemical reactions in the gas phase. In this work, we have performed a detailed RPMD study of the H + H(2) reaction and its isotopologues Mu + H(2), D + H(2) and Heμ + H(2), at temperatures ranging from 200 to 1000 K. Thermal rate coefficients and kinetic isotope effects have been computed and compared with exact QM calculations as well as with quasiclassical trajectories and experiment. The agreement with the QM results is good for the heaviest isotopologues, with errors ranging from 15% to 45%, and excellent for Mu + H(2), with errors below 15%. We have seen that RPMD is able to capture the ZPE effect very accurately, a desirable feature of any method based on molecular dynamics. We have also verified Richardson and Althorpe's prediction [J. O. Richardson and S. C. Althorpe, J. Chem. Phys., 2009, 131, 214106] that RPMD will overestimate thermal rates for asymmetric reactions and underestimate them for symmetric reactions in the deep tunneling regime. The ZPE effect along the reaction coordinate must be taken into account when assigning the reaction symmetry in the multidimensional case. |
| spellingShingle | Suleimanov, Y de Tudela, R Jambrina, P Castillo, J Sáez-Rábanos, V Manolopoulos, D Aoiz, F A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction. |
| title | A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction. |
| title_full | A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction. |
| title_fullStr | A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction. |
| title_full_unstemmed | A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction. |
| title_short | A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction. |
| title_sort | ring polymer molecular dynamics study of the isotopologues of the h h2 reaction |
| work_keys_str_mv | AT suleimanovy aringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT detudelar aringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT jambrinap aringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT castilloj aringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT saezrabanosv aringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT manolopoulosd aringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT aoizf aringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT suleimanovy ringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT detudelar ringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT jambrinap ringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT castilloj ringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT saezrabanosv ringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT manolopoulosd ringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction AT aoizf ringpolymermoleculardynamicsstudyoftheisotopologuesofthehh2reaction |