1
|
Menéndez M, Veselinova A, Zanchet A, Jambrina PG, Aoiz FJ. Rate coefficients for the O + H 2 and O + D 2 reactions: how well ring polymer molecular dynamics accounts for tunelling. Phys Chem Chem Phys 2024; 26:20947-20961. [PMID: 39046374 DOI: 10.1039/d4cp01711k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
We present here extensive calculations of the O(3P) + H2 and O(3P) + D2 reaction dynamics spanning the temperature range from 200 K to 2500 K. The calculations have been carried out using fully converged time-independent quantum mechanics (TI QM), quasiclassical trajectories (QCT) and ring polymer molecular dynamics (RPMD) on the two lowest lying adiabatic potential energy surfaces (PESs), 13A' and 13A'', calculated by Zanchet et al. [J. Chem. Phys., 2019, 151, 094307]. TI QM rate coefficients were determined using the cumulative reaction probability formalism on each PES including all of the total angular momenta and the Coriolis coupling and can be considered to be essentially exact within the Born-Oppenheimer approximation. The agreement between the rate coefficients calculated by using QM and RPMD is excellent for the reaction with D2 in almost the whole temperature range. For the reaction with H2, although the agreement is very good above 500 K, the deviations are significant at lower temperatures. In contrast, the QCT calculations largely underestimate the rate coefficients for the two isotopic variants due to their inability to account for tunelling. The differences found in the disagreements between RPMD and QM rate coefficients for the reactions for both the isotopologues are indicative of the ability of the RPMD method to accurately describe systems where tunelling plays a relevant role. Considering that both reactions are dominated by tunelling below 500 K, the present results show that RPMD is a very powerful tool for determining rate coefficients. The present QM rate coefficients calculated on adiabatic PESs slightly underestimate the best global fits of the experimental measurements, which we attribute to the intersystem crossing with the singlet 11A' PES.
Collapse
Affiliation(s)
- Marta Menéndez
- Departamento de Química Física, Unidad Asociada CSIC, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | - Anzhela Veselinova
- Departamento de Química Física, Universidad de Salamanca, 37008 Salamanca, Spain
| | - Alexandre Zanchet
- Instituto de Física Fundamental, CSIC, C/Serrano 121-123, 28006 Madrid, Spain
| | - Pablo G Jambrina
- Departamento de Química Física, Universidad de Salamanca, 37008 Salamanca, Spain
| | - F Javier Aoiz
- Departamento de Química Física, Unidad Asociada CSIC, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| |
Collapse
|
2
|
Veselinova A, Menéndez M, González-Sánchez L, Zanchet A, Aoiz FJ, Jambrina PG. Dynamical effects on the O( 3P) + D 2 reaction and its impact on the Λ-doublet population. Phys Chem Chem Phys 2024; 26:6752-6762. [PMID: 38323460 DOI: 10.1039/d3cp05510h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2024]
Abstract
The O(3P) + D2 → OD(2Π) + D reaction presents the peculiarity of taking place on two different potential energy surfaces (PESs) of different symmetry, 3A' and 3A'', which become degenerate for collinear configurations where the saddle-point of the reaction is located. The degeneracy is broken for non-collinear approaches with the energy on the 3A' PES rising more abruptly with the bending angle, making the frequency of this mode higher on the 3A' state. Consequently, the 3A' PES should be less reactive than the 3A'' one. Nevertheless, quantum scattering calculations show that the cross section is higher on the 3A' PES for energies close to the classical reaction threshold and rotationless reactant. It is found that the differences between the reactivity on the two PESs are greater for low values of total angular momentum, where the centrifugal barrier is lower and contribute to the higher population of the Π(A') Λ-doublet states of OD at low collision energies. At high collision energies, the Π(A') Λ-doublet state is also preferentially populated. Analysis of the differential cross sections reveals that the preponderance for the Π(A') Λ-doublet at low energies comes from backward scattering, originating from the reaction on the 3A' PES, while at high energies, it proceeds from a different mechanism that leads to sideways scattering on the 3A'' PES and that populates the Π(A') manifold.
Collapse
Affiliation(s)
- A Veselinova
- Departamento de Química-Física, Universidad de Salamanca, Salamanca, 37008, Spain.
| | - M Menéndez
- Departamento de Química-Física, Universidad Complutense de Madrid, Madrid, 28040, Spain
| | - L González-Sánchez
- Departamento de Química-Física, Universidad de Salamanca, Salamanca, 37008, Spain.
| | - A Zanchet
- Instituto de Física Fundamental (CSIC), 28006, Madrid, Spain
| | - F J Aoiz
- Departamento de Química-Física, Universidad Complutense de Madrid, Madrid, 28040, Spain
| | - P G Jambrina
- Departamento de Química-Física, Universidad de Salamanca, Salamanca, 37008, Spain.
| |
Collapse
|
3
|
Veselinova A, Agúndez M, Goicoechea JR, Menéndez M, Zanchet A, Verdasco E, Jambrina PG, Aoiz FJ. Quantum study of reaction O( 3 P) + H 2 ( v, j) → OH + H: OH formation in strongly UV-irradiated gas. ASTRONOMY AND ASTROPHYSICS 2021; 648:A76. [PMID: 34257462 PMCID: PMC7611199 DOI: 10.1051/0004-6361/202140428] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The reaction between atomic oxygen and molecular hydrogen is an important one in astrochemistry as it regulates the abundance of the hydroxyl radical and serves to open the chemistry of oxygen in diverse astronomical environments. However, the existence of a high activation barrier in the reaction with ground state oxygen atoms limits its efficiency in cold gas. In this study we calculate the dependence of the reaction rate coefficient on the rotational and vibrational state of H2 and evaluate the impact on the abundance of OH in interstellar regions strongly irradiated by far-UV photons, where H2 can be efficiently pumped to excited vibrational states. We use a recently calculated potential energy surface and carry out time-independent quantum mechanical scattering calculations to compute rate coefficients for the reaction O(3 P) + H2 (v, j) → OH + H, with H2 in vibrational states v = 0-7 and rotational states j = 0-10. We find that the reaction becomes significantly faster with increasing vibrational quantum number of H2, although even for high vibrational states of H2 (v = 4-5) for which the reaction is barrierless, the rate coefficient does not strictly attain the collision limit and still maintains a positive dependence with temperature. We implemented the calculated state-specific rate coefficients in the Meudon PDR code to model the Orion Bar PDR and evaluate the impact on the abundance of the OH radical. We find the fractional abundance of OH is enhanced by up to one order of magnitude in regions of the cloud corresponding to A V = 1.3-2.3, compared to the use of a thermal rate coefficient for O + H2, although the impact on the column density of OH is modest, of about 60%. The calculated rate coefficients will be useful to model and interpret JWST observations of OH in strongly UV-illuminated environments.
Collapse
Affiliation(s)
- A. Veselinova
- Departamento de Química Física, University of Salamanca, Plaza Caidos S/N, E-37008, Salamanca, Spain
- Departamento de Química Física (Unidad Asociada al CSIC), Universidad Complutense de Madrid, Ciudad Universitaria, S/N, E-20840, Madrid, Spain
| | - M. Agúndez
- Instituto de Física Fundamental, CSIC, Calle Serrano 121-123, E-28006, Madrid, Spain
| | - J. R. Goicoechea
- Instituto de Física Fundamental, CSIC, Calle Serrano 121-123, E-28006, Madrid, Spain
| | - M. Menéndez
- Departamento de Química Física (Unidad Asociada al CSIC), Universidad Complutense de Madrid, Ciudad Universitaria, S/N, E-20840, Madrid, Spain
| | - A. Zanchet
- Instituto de Física Fundamental, CSIC, Calle Serrano 121-123, E-28006, Madrid, Spain
| | - E. Verdasco
- Departamento de Química Física (Unidad Asociada al CSIC), Universidad Complutense de Madrid, Ciudad Universitaria, S/N, E-20840, Madrid, Spain
| | - P. G. Jambrina
- Departamento de Química Física, University of Salamanca, Plaza Caidos S/N, E-37008, Salamanca, Spain
| | - F. J. Aoiz
- Departamento de Química Física (Unidad Asociada al CSIC), Universidad Complutense de Madrid, Ciudad Universitaria, S/N, E-20840, Madrid, Spain
| |
Collapse
|
4
|
Ortega P, Zanchet A, Sanz-Sanz C, Gómez-Carrasco S, González-Sánchez L, Jambrina PG. DpgC-Catalyzed Peroxidation of 3,5-Dihydroxyphenylacetyl-CoA (DPA-CoA): Insights into the Spin-Forbidden Transition and Charge Transfer Mechanisms*. Chemistry 2020; 27:1700-1712. [PMID: 32975323 DOI: 10.1002/chem.202002993] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Indexed: 11/06/2022]
Abstract
Despite being a very strong oxidizing agent, most organic molecules are not oxidized in the presence of O2 at room temperature because O2 is a diradical whereas most organic molecules are closed-shell. Oxidation then requires a change in the spin state of the system, which is forbidden according to non-relativistic quantum theory. To overcome this limitation, oxygenases usually rely on metal or redox cofactors to catalyze the incorporation of, at least, one oxygen atom into an organic substrate. However, some oxygenases do not require any cofactor, and the detailed mechanism followed by these enzymes remains elusive. To fill this gap, here the mechanism for the enzymatic cofactor-independent oxidation of 3,5-dihydroxyphenylacetyl-CoA (DPA-CoA) is studied by combining multireference calculations on a model system with QM/MM calculations. Our results reveal that intersystem crossing takes place without requiring the previous protonation of molecular oxygen. The characterization of the electronic states reveals that electron transfer is concomitant with the triplet-singlet transition. The enzyme plays a passive role in promoting the intersystem crossing, although spontaneous reorganization of the water wire connecting the active site with the bulk presets the substrate for subsequent chemical transformations. The results show that the stabilization of the singlet radical-pair between dioxygen and enolate is enough to promote spin-forbidden reaction without the need for neither metal cofactors nor basic residues in the active site.
Collapse
Affiliation(s)
- Pablo Ortega
- Departamento de Química Física, University of Salamanca, Salamanca, 37008, Spain
| | - Alexandre Zanchet
- Departamento de Química Física, University of Salamanca, Salamanca, 37008, Spain.,Instituto de Física Fundamental (CSIC), Madrid, 28006, Spain
| | - Cristina Sanz-Sanz
- Departamento de Química Física Aplicada, University Autónoma de Madrid, Madrid, 28049, Spain
| | | | | | - Pablo G Jambrina
- Departamento de Química Física, University of Salamanca, Salamanca, 37008, Spain
| |
Collapse
|
5
|
Jambrina PG, Zanchet A, Menéndez M, Herrero VJ, Aoiz FJ. Unexpected dynamical effects change the lambda-doublet propensity in the tunneling region for the O( 3P) + H 2 reaction. Phys Chem Chem Phys 2019; 21:25389-25396. [PMID: 31709441 DOI: 10.1039/c9cp04690a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
One of the most relevant features of the O(3P) + H2 reaction is that it occurs on two different potential energy surfaces (PESs) of symmetries A' and A'' that correlate reactants and products. The respective saddle points, which correspond to a collinear arrangement, are the same for both PESs, whilst the barrier height rises more abruptly on the 3A' PES than on the 3A'' PES. Accordingly, the reactivity on the 3A'' PES should be always higher than on the 3A' PES. In this work, we present accurate quantum-scattering calculations showing that this is not always the case for rotationless reactants, where dynamical factors near the reaction threshold cause the 3A' PES to dominate at energies around the barrier. Further calculation of cross sections and Λ-doublet populations has allowed us to establish how the reaction mechanism changes from the deep tunneling regime to hyperthermal energies.
Collapse
Affiliation(s)
- P G Jambrina
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad de Salamanca, 37003, Salamanca, Spain
| | - A Zanchet
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad de Salamanca, 37003, Salamanca, Spain and Departamento de Química Física, Facultad de Química, Universidad Complutense de Madrid (Unidad Asociada CSIC), 28040 Madrid, Spain.
| | - M Menéndez
- Departamento de Química Física, Facultad de Química, Universidad Complutense de Madrid (Unidad Asociada CSIC), 28040 Madrid, Spain.
| | - V J Herrero
- Instituto de Estructura de la Materia, IEM-CSIC c/Serrano 123, 28006 Madrid, Spain
| | - F J Aoiz
- Departamento de Química Física, Facultad de Química, Universidad Complutense de Madrid (Unidad Asociada CSIC), 28040 Madrid, Spain.
| |
Collapse
|