© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

Journal of Physics A Mathematical and Theoretical

Tunneling by a semiclassical initial value method with higher order corrections

Gili Hochman, Kenneth G. Kay

Published Aug 22, 2008• 11 citations

Abstract

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.

DOI: 10.1088/1751-8113/41/38/385303

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

Semiclassical corrections to the Herman-Kluk propagator

Gili Hochman, Kenneth G. Kay

The approximate expression for the propagator of Herman and Kluk is systematically improved by application of semiclassical corrections for a double-well system. The results confirm the semiclassical nature of the Herman-Kluk approximation. The influence of caustics on the effectiveness of the corrections is discussed.

Physical Review A•Jun 27, 2006•14 citations

Semiclassical tunneling in the initial value representation

Kenneth G. Kay

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.

The Journal of Chemical Physics•Aug 15, 1997•105 citations

Determination of the border between “shallow” and “deep” tunneling regions for the Herman-Kluk method by the asymptotic approach

Atushi Tanaka

The evaluation of a tunneling tail by the Herman-Kluk method, which is a quasiclassical way to compute quantum dynamics, is examined by asymptotic analysis. In the shallower part of the tail, as well as in the classically allowed region, it is shown that the leading terms of semiclassical evaluations of quantum theory and the Herman-Kluk formula agree, which is known as an asymptotic equivalence. In the deeper part, it is shown that the asymptotic equivalence breaks down, due to the emergence of unusual ``tunneling trajectory,'' which is an artifact of the Herman-Kluk method.

arXiv:quant-ph/0601135•Physical Review A•Feb 21, 2006•6 citations

Semiclassical dynamics with quantum trajectories: Formulation and comparison with the semiclassical initial value representation propagator

Sophya V. Garashchuk, Vitaly A. Rassolov

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.

The Journal of Chemical Physics•Jan 30, 2003•54 citations

Reply to Comment on Semiclassical approximations in phase space with coherent states

Michel Baranger, Marcus A. M. de Aguiar, F Keck +2 more

The Herman–Kluk (HK) formula was shown in (Baranger et al J. Phys. A: Math. Gen. 34 7227 ) not to be a correct semiclassical limit of an exact quantum mechanical formula. Two previous attempts to derive it using semiclassical arguments contain serious errors. These statements are left totally untouched by Herman and Grossmann's comment. They argue that the formula which we found to be at fault is not the one that should be called the HK formula. However, the formula we criticized is definitely one of the steps, in fact the main step, in these two published derivations of the HK formula. Very recently, a new derivation was published by Miller. It is interesting, but it is not semiclassical.

Journal of Physics A Mathematical and General•Oct 21, 2002•29 citations

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.

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 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.

The evaluation of a tunneling tail by the Herman-Kluk method, which is a quasiclassical way to compute quantum dynamics, is examined by asymptotic analysis. In the shallower part of the tail, as well as in the classically allowed region, it is shown that the leading terms of semiclassical evaluations of quantum theory and the Herman-Kluk formula agree, which is known as an asymptotic equivalence. In the deeper part, it is shown that the asymptotic equivalence breaks down, due to the emergence of unusual ``tunneling trajectory,'' which is an artifact of the Herman-Kluk method.

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.

The Herman–Kluk (HK) formula was shown in (Baranger et al J. Phys. A: Math. Gen. 34 7227 ) not to be a correct semiclassical limit of an exact quantum mechanical formula. Two previous attempts to derive it using semiclassical arguments contain serious errors. These statements are left totally untouched by Herman and Grossmann's comment. They argue that the formula which we found to be at fault is not the one that should be called the HK formula. However, the formula we criticized is definitely one of the steps, in fact the main step, in these two published derivations of the HK formula. Very recently, a new derivation was published by Miller. It is interesting, but it is not semiclassical.