Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Semiclassical tunneling in the initial value representation

Kenneth G. Kay

Published Aug 15, 1997• 105 citations
Publisher

Abstract

Tunneling in the one-dimensional Eckart system is treated by a semiclassical method that describes the S-matrix in terms of an integral over the initial momenta of real-valued classical trajectories. The results are found to be sensitive to a certain parameter γ which is expected to be essentially arbitrary for classically allowed processes. Analysis of the semiclassical error allows formulation of conditions for the validity of the tunneling treatment. This, in turn, leads to an explanation for the sensitivity of the results to γ and an understanding of how this parameter should be chosen. With an optimized choice, the semiclassical method is found to yield very accurate tunneling results even for probabilities as small as 10−10. The relationship between the present method and the conventional uniform semiclassical treatment of barrier tunneling is discussed.

DOI
10.1063/1.475143

Related papers

Wave packet propagation across barriers by semiclassical initial value methods

Jakob Petersen, Kenneth G. Kay

Semiclassical initial value representation (IVR) formulas for the propagator have difficulty describing tunneling through barriers. A key reason is that these formulas do not automatically reduce, in the classical limit, to the version of the Van Vleck-Gutzwiller (VVG) propagator required to treat barrier tunneling, which involves trajectories that have complex initial conditions and that follow paths in complex time. In this work, a simple IVR expression, that has the correct tunneling form in the classical limit, is derived for the propagator in the case of one-dimensional barrier transmission. Similarly, an IVR formula, that reduces to the Generalized Gaussian Wave Packet Dynamics (GGWPD) expression [D. Huber, E. J. Heller, and R. Littlejohn, J. Chem. Phys. 89, 2003 (1988)] in the classical limit, is derived for the transmitted wave packet. Uniform semiclassical versions of the IVR formulas are presented and simplified expressions in terms of real trajectories and WKB penetration factors are described. Numerical tests show that the uniform IVR treatment gives good results for wave packet transmission through the Eckart and Gaussian barriers in all cases examined. In contrast, even when applied with the proper complex trajectories, the VVG and GGWPD treatments are inaccurate when the mean energy of the wave packet is near the classical transmission threshold. The IVR expressions for the propagator and wave packet are cast as contour integrals in the complex space of initial conditions and these are generalized to potentially allow treatment of a larger variety of systems. A steepest descent analysis of the contour integral formula for the wave packet in the present cases confirms its relationship to the GGWPD method, verifies its semiclassical validity, and explains results of numerical calculations.

The Journal of Chemical Physics•Jul 2, 2015•8 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 Theoretical•Aug 22, 2008•11 citations

An analytically soluble model for the penetration of a two-dimensional barrier: Quantal and semiclassical treatment

Herbert Massmann, K. Möhring, Peter A. Ring

The reflection and transmission of a wave on a two-dimensional barrier is studied quantum mechanically and semiclassically. The simple model consists of two degrees of freedom. One of them contains a barrier (of the Eckart type); the second allows only a bound motion in discrete quantum states. These two degrees of freedom are coupled. The transition probabilities for transitions between the different quantum states can, in the model considered here, be evaluated analytically, both quantum mechanically and semiclassically. Several of the so called uniform semiclassical approximations and also the initial value integral representation for the semiclassical transition amplitudes are considered and compared to the exact quantum mechanical results. Some new insight on the semiclassical treatment of the tunneling process is gained, specially concerning the proper choice of the integration path for the initial value integral representation of the transition amplitudes.

The Journal of Chemical Physics•Jun 15, 1977•9 citations

Time dependent semiclassical tunneling through one dimensional barriers using only real valued trajectories

Michael F. Herman

The time independent semiclassical treatment of barrier tunneling has been understood for a very long time. Several semiclassical approaches to time dependent tunneling through barriers have also been presented. These typically involve trajectories for which the position variable is a complex function of time. In this paper, a method is presented that uses only real valued trajectories, thus avoiding the complications that can arise when complex trajectories are employed. This is accomplished by expressing the time dependent wave packet as an integration over momentum. The action function in the exponent in this expression is expanded to second order in the momentum. The expansion is around the momentum, p0(*) , at which the derivative of the real part of the action is zero. The resulting Gaussian integral is then taken. The stationary phase approximation requires that the derivative of the full action is zero at the expansion point, and this leads to a complex initial momentum and complex tunneling trajectories. The "pseudo-stationary phase" approximation employed in this work results in real values for the initial momentum and real valued trajectories. The transmission probabilities obtained are found to be in good agreement with exact quantum results.

The Journal of Chemical Physics•Oct 27, 2015•0 citations

Some New Classical and Semiclassical Models for Describing Tunneling Processes with Real-Valued Classical Trajectories

Jianhua Xing, Eduardo A. Coronado, William Hughes Miller

A model for describing barrier tunneling (or other classically forbidden processes) using purely real-valued classical trajectories is presented. The basic idea is to introduce an auxiliary degree of freedom that allows for fluctuations in the potential-energy surface for the original classical degrees of freedom. The model can be applied purely classically or better semiclassically (e . g ., via the initial value representation). Numerical results for 1D barrier tunneling are presented to illustrate the model.

The Journal of Physical Chemistry B•Apr 11, 2001•15 citations

Coherent Classical-Path Description of Deep Tunneling

Dong Hui Zhang, Eli Pollak

A central challenge to the semiclassical description of quantum mechanics is the quantum phenomenon of "deep" tunneling. Here we show that real time classical trajectories suffice to account correctly even for deep quantum tunneling, using a recently formulated semiclassical initial value representation series of the quantum propagator and a prefactor free semiclassical propagator. Deep quantum tunneling is effected through what we term as coherent classical paths which are composed of one or more classical trajectories that lead from reactant to product but are discontinuous along the way. The end and initial phase space points of consecutive classical trajectories contributing to the coherent path are close to each other in the sense that the distance between them is weighted by a coherent state overlap matrix element. Results are presented for thermal and energy dependent tunneling through a symmetric Eckart barrier.

Physical Review Letters•Sep 27, 2004•31 citations
Tunneling in the one-dimensional Eckart system is treated by a semiclassical method that describes the S-matrix in terms of an integral over the initial momenta of real-valued classical trajectories. The results are found to be sensitive to a certain parameter γ which is expected to be essentially arbitrary for classically allowed processes. Analysis of the semiclassical error allows formulation of conditions for the validity of the tunneling treatment. This, in turn, leads to an explanation for the sensitivity of the results to γ and an understanding of how this parameter should be chosen. With an optimized choice, the semiclassical method is found to yield very accurate tunneling results even for probabilities as small as 10−10. The relationship between the present method and the conventional uniform semiclassical treatment of barrier tunneling is discussed.
Semiclassical initial value representation (IVR) formulas for the propagator have difficulty describing tunneling through barriers. A key reason is that these formulas do not automatically reduce, in the classical limit, to the version of the Van Vleck-Gutzwiller (VVG) propagator required to treat barrier tunneling, which involves trajectories that have complex initial conditions and that follow paths in complex time. In this work, a simple IVR expression, that has the correct tunneling form in the classical limit, is derived for the propagator in the case of one-dimensional barrier transmission. Similarly, an IVR formula, that reduces to the Generalized Gaussian Wave Packet Dynamics (GGWPD) expression [D. Huber, E. J. Heller, and R. Littlejohn, J. Chem. Phys. 89, 2003 (1988)] in the classical limit, is derived for the transmitted wave packet. Uniform semiclassical versions of the IVR formulas are presented and simplified expressions in terms of real trajectories and WKB penetration factors are described. Numerical tests show that the uniform IVR treatment gives good results for wave packet transmission through the Eckart and Gaussian barriers in all cases examined. In contrast, even when applied with the proper complex trajectories, the VVG and GGWPD treatments are inaccurate when the mean energy of the wave packet is near the classical transmission threshold. The IVR expressions for the propagator and wave packet are cast as contour integrals in the complex space of initial conditions and these are generalized to potentially allow treatment of a larger variety of systems. A steepest descent analysis of the contour integral formula for the wave packet in the present cases confirms its relationship to the GGWPD method, verifies its semiclassical validity, and explains results of numerical calculations.
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.
The reflection and transmission of a wave on a two-dimensional barrier is studied quantum mechanically and semiclassically. The simple model consists of two degrees of freedom. One of them contains a barrier (of the Eckart type); the second allows only a bound motion in discrete quantum states. These two degrees of freedom are coupled. The transition probabilities for transitions between the different quantum states can, in the model considered here, be evaluated analytically, both quantum mechanically and semiclassically. Several of the so called uniform semiclassical approximations and also the initial value integral representation for the semiclassical transition amplitudes are considered and compared to the exact quantum mechanical results. Some new insight on the semiclassical treatment of the tunneling process is gained, specially concerning the proper choice of the integration path for the initial value integral representation of the transition amplitudes.
The time independent semiclassical treatment of barrier tunneling has been understood for a very long time. Several semiclassical approaches to time dependent tunneling through barriers have also been presented. These typically involve trajectories for which the position variable is a complex function of time. In this paper, a method is presented that uses only real valued trajectories, thus avoiding the complications that can arise when complex trajectories are employed. This is accomplished by expressing the time dependent wave packet as an integration over momentum. The action function in the exponent in this expression is expanded to second order in the momentum. The expansion is around the momentum, p0(*) , at which the derivative of the real part of the action is zero. The resulting Gaussian integral is then taken. The stationary phase approximation requires that the derivative of the full action is zero at the expansion point, and this leads to a complex initial momentum and complex tunneling trajectories. The "pseudo-stationary phase" approximation employed in this work results in real values for the initial momentum and real valued trajectories. The transmission probabilities obtained are found to be in good agreement with exact quantum results.
A model for describing barrier tunneling (or other classically forbidden processes) using purely real-valued classical trajectories is presented. The basic idea is to introduce an auxiliary degree of freedom that allows for fluctuations in the potential-energy surface for the original classical degrees of freedom. The model can be applied purely classically or better semiclassically (e . g ., via the initial value representation). Numerical results for 1D barrier tunneling are presented to illustrate the model.
A central challenge to the semiclassical description of quantum mechanics is the quantum phenomenon of "deep" tunneling. Here we show that real time classical trajectories suffice to account correctly even for deep quantum tunneling, using a recently formulated semiclassical initial value representation series of the quantum propagator and a prefactor free semiclassical propagator. Deep quantum tunneling is effected through what we term as coherent classical paths which are composed of one or more classical trajectories that lead from reactant to product but are discontinuous along the way. The end and initial phase space points of consecutive classical trajectories contributing to the coherent path are close to each other in the sense that the distance between them is weighted by a coherent state overlap matrix element. Results are presented for thermal and energy dependent tunneling through a symmetric Eckart barrier.