The Journal of Chemical Physics

ℏ2 expansion of the transmission probability through a barrier

Eli Pollak, Jianshu Cao

Published Jul 25, 2022• 18 citations

Abstract

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.

DOI: 10.1063/5.0106649

Related papers

Thermal rate of transmission through a barrier: Exact expansion of up to and including terms of order ℏ 4

Eli Pollak, Jianshu Cao

Nine decades after Wigner's formulation of quantum rate theory, his celebrated result was recently generalized to the asymmetric barrier using an exact first-order expansion of the transmission probability in terms of ${\ensuremath{\hbar}}^{2}$. This paper extends the first-order quantum correction to second-order correction of order ${\ensuremath{\hbar}}^{4}$ for the thermally averaged transmission probability through an arbitrary barrier. The derivation employs a systematic expansion of the projection operator onto products and the thermal distribution which involves a Taylor expansion of the potential about the barrier up to eighth order. The resulting exact analytical expression is calibrated with numerical calculations of several model potentials and shows excellent agreement when the ${\ensuremath{\hbar}}^{4}$ term is included. In comparison, the semiclassical transition state theory cannot reproduce the correct ${\ensuremath{\hbar}}^{4}$ terms when the anharmonicity is treated on the level of VPT-4 (vibrational perturbation theory---fourth order) and will potentially need a VPT-6 expansion. Further analysis of the quartic barrier reveals suppressed transmission due to the dominant role of quantum reflection above the barrier. These results not only provide a conceptual framework but can also be applied to heavy atom tunneling and machine learning.

Physical Review AFeb 6, 202311 citations

Tunneling by a semiclassical initial value method with higher order corrections

Gili Hochman, Kenneth G. Kay

Tunneling through the one-dimensional Eckart barrier is treated by applying a higher order semiclassical approximation which adds corrections proportional to powers of ℏ to the Herman–Kluk (HK) initial value approximation. Although the usual, zero-order HK treatment is very poor in this case, the first- and second-order corrections substantially improve the accuracy of the computed tunneling probabilities. To investigate how this works, the HK expression is shown to be equivalent, for the present purposes, to a formula involving a single integral over the initial momentum, with an integrand that has a simple analytical form. Similarly, the most important part of the first-order correction term is shown to be expressible in a very simple form. For a particular range of energies, the integral can be analyzed in terms of a steepest descent treatment along a path through a caustic. In this way, it is verified that the zero-order approximation does not approach the correct classical limit as ℏ → 0, but the first-order term, which does not vanish in this limit, improves the accuracy of the result. More generally, the corrections of each order contain terms that are of all orders in ℏ 1/3 , including those that are O (ℏ 0 ) and which survive in the classical limit. The infinite sum over such terms is performed analytically and shown to yield the correct classical limit for the tunneling amplitude.

Journal of Physics A Mathematical and TheoreticalAug 22, 200811 citations

SEMICLASSICAL TUNNELING EFFECT

Johan Öhman

This work is a contribution to the theory of the quantum tunneling effect. In the literature at least two barriers are studied, for which the third-order contribution to the phase - integral asymptotic approximation of the transmission coefficient do not contribute. These are the parabolic barrier and the inverse Morse barrier. In the present work we will show that with a proper choice of the so called base function there is at least one more barrier in this category namely the Eckart-Epstein potential. The fact that the third - order contribution vanishes is a good indication that we have found an optimal choice of the base function, and the treatment to find an optimal base function may be possible to generalize to other classes of potential barriers. For particles of a low energy compared to the energy near the top of the barrier we obtain a vary low transmission coefficient, which means that the probability for tunneling to occur is very low. There exist some cases, for example that with a double barrier which is transparent, even for certain relatively low energies but no evidence for this kind of transparency for a single barrier has been found. The present work does not give any such evidence. At the same time there are still speculations on cold fusion like effects, which would demand a higher probability for tunneling through for a single barrier.

Jan 1, 20140 citations

Tunneling-time probability distribution

Randall S. Dumont, T. L. Marchioro II

Localized wave-packet tunneling through a one-dimensional potential barrier is investigated with respect to the distribution of tunneling times. Specifically, we construct a quantity--- termed the ``tunneled flux''---which has properties of a tunneling-time probability distribution. The tunneled flux is determined completely in terms of the time evolution of an initially localized wave packet. The tunneled flux and the corresponding tunneling-time distributions are investigated analytically via a semiclassical approximation. Additionally, numerical studies are performed with both the semiclassical and quantal versions of the tunneled flux. The semiclassical calculation qualitatively reproduces numerical tunneled fluxes, converging quantitatively as \ensuremath{\Elzxh}\ensuremath{\rightarrow}0. It also provides a wealth of insight into the nature of tunneling processes. Our principle conclusions are as follows: (i) There is an essential disparity between the time- and energy-domain pictures of wave-packet tunneling; (ii) the tunneling-time distribution can be modeled as an exponential of a skewed Lorentzian function, with a width governed primarily by the phase-space asymmetry of the initial wave packet; and (iii) tunneling is faster than simple classical motion from one side of the barrier to another, even assuming an instantaneous transit between the turning points. Effectively, the barrier acts as a filter for the high-momentum components of the initial wave packet.

Physical Review AJan 1, 199380 citations

Instantaneous Tunneling Flight Time for Wavepacket Transmission through Asymmetric Barriers

Jakob Petersen, Eli Pollak

The time it takes a particle to tunnel through the asymmetric Eckart barrier potential is investigated using Gaussian wavepackets, where the barrier serves as a model for the potential along a chemical reaction coordinate. We have previously shown that the, in principle experimentally measurable, tunneling flight time, which determines the time taken by the transmitted particle to traverse the barrier, vanishes for symmetric potentials like the Eckart and square barrier [ Petersen , J. ; Pollak , E. J. Phys. Chem. Lett. 2017 , 9 , 4017 ]. Here we show that the same result is obtained for the asymmetric Eckart barrier potential, and therefore, the zero tunneling flight time seems to be a general result for one-dimensional time-independent potentials. The wavepacket dynamics is simulated using both an exact quantum mechanical method and a classical Wigner prescription. The excellent agreement between the two methods shows that quantum coherences are not important in pure one-dimensional tunneling and reinforces the conclusion that the tunneling flight time vanishes.

The Journal of Physical Chemistry AMar 20, 201814 citations

Symmetry relationship of quantum tunneling and its applications

Peiwei Wen, C. J. Lin, H. M. Jia +3 more

Quantum mechanics serves as the fundamental cornerstone for numerous areas of research. Finding simple and analytical quantum transmission coefficients for potentials is a rare occurrence. Notably, we have discovered a readily provable yet unreported symmetry tunneling law. For an infinite parabolic potential, it can be proved analytically that the sum of the one-dimensional tunneling probabilities at incident energies equal to the barrier height plus or minus any energy deviation is consistently 1. Based on this relationship, the original WKB approximation and the Kemble's formula can be generalized to the above-barrier energy region in a simplified but more accurate manner, referred to as the symmetric WKB approach in this context. For the realistic potential between two nuclei, numerical results demonstrate that the symmetry relationship at the above-barrier energy region and the sub-barrier energy is also well satisfied, as well as in the cases of the multichannel WKB method. The symmetry tunneling law has been effectively extended to the Eckart potential as well. Such symmetry in tunneling may have universal characteristics but its underlying reason remains to be uncovered, which may help us to further understand the quantum mechanisms.

Physical Review CJun 3, 20242 citations