An official website of the United States government
We have measured the merged-beams rate coefficient for the titular isotope exchange reactions as a function of the relative collision energy in the range of ∼3 meV–10 eV. The results appear to scale with the number of available sites for deuteration. We have performed extensive theoretical calculations to characterize the zero-point energy corrected reaction path. Vibrationally adiabatic minimum energy paths were obtained using a combination of unrestricted quadratic configuration interaction of single and double excitations and internally contracted multireference configuration interaction calculations. The resulting barrier height, ranging from 68 meV to 89 meV, together with the various asymptotes that may be reached in the collision, was used in a classical over-the-barrier model. All competing endoergic reaction channels were taken into account using a flux reduction factor. This model reproduces all three experimental sets quite satisfactorily. In order to generate thermal rate coefficients down to 10 K, the internal excitation energy distribution of each H3+ isotopologue is evaluated level by level using available line lists and accurate spectroscopic parameters. Tunneling is accounted for by a direct inclusion of the exact quantum tunneling probability in the evaluation of the cross section. We derive a thermal rate coefficient of <1×10−12 cm3 s−1 for temperatures below 44 K, 86 K, and 139 K for the reaction of D with H3+, H2D+, and D2H+, respectively, with tunneling effects included. The derived thermal rate coefficients exceed the ring polymer molecular dynamics prediction of Bulut et al. [J. Phys. Chem. A 123, 8766 (2019)] at all temperatures.
more »
« less
ℏ2 expansion of the transmission probability through a barrier
Ninety years ago, Wigner derived the leading order expansion term in ℏ2 for the tunneling rate through a symmetric barrier. His derivation included two contributions: one came from the parabolic barrier, but a second term involved the fourth-order derivative of the potential at the barrier top. He left us with a challenge, which is answered in this paper, to derive the same but for an asymmetric barrier. A crucial element of the derivation is obtaining the ℏ2 expansion term for the projection operator, which appears in the flux-side expression for the rate. It is also reassuring that an analytical calculation of semiclassical transition state theory (TST) reproduces the anharmonic corrections to the leading order of ℏ2. The efficacy of the resulting expression is demonstrated for an Eckart barrier, leading to the conclusion that especially when considering heavy atom tunneling, one should use the expansion derived in this paper, rather than the parabolic barrier approximation. The rate expression derived here reveals how the classical TST limit is approached as a function of ℏ and, thus, provides critical insights to understand the validity of popular approximate theories, such as the classical Wigner, centroid molecular dynamics, and ring polymer molecular dynamics methods.
more »
« less
- Award ID(s):
- 1800301
- PAR ID:
- 10593500
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 157
- Issue:
- 7
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Natural enzymes are powerful catalysts, reducing the apparent activation energy for reaction, enabling chemistry to proceed as much as 1015 times faster than the corresponding solution reaction. It has been suggested for some time that in some cases quantum tunneling can contribute to this rate enhancement by offering pathways through a barrier inaccessible to activated events. A central question of interest to both physical chemists and biochemists is the extent to which evolution introduces below the barrier or tunneling mechanisms. In view of the rapidly expanding chemistries for which artificial enzymes have now been created, it is of interest to see how quantum tunneling has been used in these reactions. In this paper, we study the evolution of possible proton tunneling during C-H bond cleavage in enzymes that catalyze the Morita-Baylis-Hillman (MBH) reaction. The enzymes were generated by theoretical design followed by laboratory evolution. We employ classical and centroid molecular dynamics approaches in path sampling computations to determine if there is a quantum contribution to lowering the free energy of the proton transfer for various experimentally generated protein and substrate combinations. This data is compared to experiments reporting on the observed kinetic isotope effect (KIE) for the relevant reactions. Our results indicate modest involvement of tunneling when laboratory evolution has resulted in a system with a higher classical free energy barrier to chemistry (that is when optimization of processes other than chemistry result in a higher chemical barrier.)more » « less
-
The ultimate regularity of quantum mechanics creates a tension with the assumption of classical chaos used in many of our pictures of chemical reaction dynamics. Out-of-time-order correlators (OTOCs) provide a quantum analog to the Lyapunov exponents that characterize classical chaotic motion. Maldacena, Shenker, and Stanford have suggested a fundamental quantum bound for the rate of information scrambling, which resembles a limit suggested by Herzfeld for chemical reaction rates. Here, we use OTOCs to study model reactions based on a double-well reaction coordinate coupled to anharmonic oscillators or to a continuum oscillator bath. Upon cooling, as one enters the tunneling regime where the reaction rate does not strongly depend on temperature, the quantum Lyapunov exponent can approach the scrambling bound and the effective reaction rate obtained from a population correlation function can approach the Herzfeld limit on reaction rates: Tunneling increases scrambling by expanding the state space available to the system. The coupling of a dissipative continuum bath to the reaction coordinate reduces the scrambling rate obtained from the early-time OTOC, thus making the scrambling bound harder to reach, in the same way that friction is known to lower the temperature at which thermally activated barrier crossing goes over to the low-temperature activationless tunneling regime. Thus, chemical reactions entering the tunneling regime can be information scramblers as powerful as the black holes to which the quantum Lyapunov exponent bound has usually been applied.more » « less
-
Subcritical crack growth can occur under a constant applied load below the threshold value for catastrophic failure, also known as static fatigue. Here, we report how a crack grows under a combination of stress-intensity factor (K) and temperature in a model brittle glass using molecular dynamics simulations. The model glass is under dry conditions, thus avoiding the complexity of corrosion chemistry. The crack growth rate is shown to be inconsistent with the commonly used subcritical crack growth model rooted in the transition state theory (TST), in which the applied stress-intensity factor reduces the transition barrier. A new subcritical crack growth model is proposed with a constant barrier and a K-dependent prefactor in TST, representing the size of the region for potential bond breaking. The thermomechanical condition for subcritical crack growth is also mapped in the K-T domain, in between elastic deformation and catastrophic fracture regimes. Finally, we show substantial crack self-healing once the applied load is removed, under the thermodynamic driving force of surface energy reduction. Our findings provide new insights into the mechanochemical coupling during static fatigue and call for experimental investigation of whether the activation energy is K-dependent.more » « less
-
null (Ed.)The T 1 excited state relaxation in thiophosgene has attracted much attention as a relatively simple model for the intersystem crossing (ISC) transitions in polyatomic molecules. The very short (20–40 ps) T 1 lifetime predicted in several theoretical studies strongly disagrees with the experimental values (∼20 ns) indicating that the kinetics of T 1 → S 0 ISC is not well understood. We use the nonadiabatic transition state theory (NA-TST) with the Zhu–Nakamura transition probability and the multireference perturbation theory (CASPT2) to show that the T 1 → S 0 ISC occurs in the quantum tunneling regime. We also introduce a new zero-point vibrational energy correction scheme that improves the accuracy of the predicted ISC rate constants at low internal energies. The predicted lifetimes of the T 1 vibrational states are between one and two orders of magnitude larger than the experimental values. This overestimation is attributed to the multidimensional nature of quantum tunneling that facilitates ISC transitions along the non-minimum energy path and is not accounted for in the one-dimensional NA-TST.more » « less
