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

Physical Review A

Time-dependent semiclassical tunneling through barriers

Kenneth G. Kay

Published Jul 23, 2013• 19 citations
Publisher

Abstract

There has been a surprising degree of uncertainty about how, or even whether, the standard semiclassical theory of Van Vleck and Gutzwiller (VVG) can be used to treat the time-dependent dynamics of tunneling through barriers. At issue is how semiclassical theory describes the time evolution of tunneling. This uncertainty has hindered development of practical techniques for semiclassical calculations of multidimensional systems. It is shown here that the VVG theory can indeed be applied to barrier tunneling, but generally requires consideration of classical trajectories with both complex initial conditions and complex time paths, even for systems with one degree of freedom. The present approach is verified by applying it to the classically forbidden propagation of a wave packet through a one-dimensional barrier---a very basic problem that cannot be treated by the VVG theory without the modifications described here.

DOI
10.1103/physreva.88.012122

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 computations of time-dependent tunneling

M. Kira, Ilkka Juhani Tittonen, W. K. Lai +1 more

In this paper we consider the time evolution of wave-packet tunneling through a potential barrier. Using a path-summation approach, we can describe the quantum evolution in terms of particle trajectories both in the allowed and forbidden regions. Near the barrier edges, where the potential changes rapidly, the wave aspects dominate. These make the trajectories branch into infinite families of paths, which have to be included in the summation. Combining these features, we treat correctly the complementary aspects of the full quantum process, and this allows us to obtain a computational method that gives accurate numerical results with considerably less expenditure of computer time than a direct integration of the Schr\"odinger equation. We test the method on a simple rectangular potential by comparing the calculations with both conventional wave scattering results and computations from the Schr\"odinger equation. A general potential can be approximated by a staircase, which allows the application of our method. A simple adiabatic extension is shown to work excellently for the case of a periodically modulated barrier.

Physical Review A•Apr 1, 1995•23 citations

Dynamical tunneling of bound systems through a potential barrier: complex way to the top

Fedor L. Bezrukov, D. G. Levkov

A semiclassical method of complex trajectories for the calculation of the tunneling exponent in systems with many degrees of freedom is further developed. It is supplemented with an easily implementable technique, which enables one to single out the physically relevant trajectory from the whole set of complex classical trajectories. The method is applied to semiclassical transitions of a bound system through a potential barrier. We find that the properties of physically relevant complex trajectories are qualitatively different in the cases of potential tunneling at low energy and dynamical tunneling at energies exceeding the barrier height. Namely, in the case of high energies, the physically relevant complex trajectories describe tunneling via creation of a state close to the top of the barrier. The method is checked against exact solutions of the Schrodinger equation in a quantum mechanical system of two degrees of freedom.

•Dec 16, 2003•24 citations

Dynamical tunneling of bound systems through a potential barrier: Complex way to the top

Fedor L. Bezrukov, D. Levkov

A semiclassical method of complex trajectories for the calculation of the tunneling exponent n systems with many degrees of freedom is further developed. It is supplemented with an easily implemented technique that enables one to single out the physically relevant trajectory from the whole set of complex classical trajectories. The method is applied to semiclassical transitions of a bound system through a potential barrier. We find that the properties of physically relevant complex trajectories are qualitatively different in the cases of potential tunneling at low energy and dynamical tunneling at energies exceeding the barrier height. Namely, in the case of high energies, the physically relevant complex trajectories describe tunneling via creation of a state close to the top of the barrier. The method is checked against exact solutions of the Schrödinger equation in a quantum-mechanical system of two degrees of freedom.

arXiv:quant-ph/0312144•Journal of Experimental and Theoretical Physics•Apr 1, 2004•29 citations

Complex-Domain Semiclassical Theory and Periodically Perturbed Barrier Tunneling Problems

Kin’ya Takahashi, Kensuke S. Ikeda

In recent papers we demonstrated a new class of barrier tunneling phenomena which can be observed in periodically driven one dimensional barrier systems (i.e., 1.5D barrier tunneling systems). In such systems the tunneling component is in general accompanied by remarkable interference fringe patterns as the perturbation strength exceeds a certain characteristic value. Such a phenomenon is quite general and is not restricted to periodically driven 1D systems, and it may be also observed in intrinsic two dimensional autonomous systems. The complex-domain semiclassical theory developed by ourselves was applied to the problem, and we succeeded in reconstructing every detail of the fringed tunneling wavefunction only by using the complexified classical tunneling trajectories. This fact greatly motivates us to investigate the mechanism of the tunneling in terms of “classical mechanics”. The later part of the present paper is devoted to asserting the generality of the “classical mechanism” of the fringed tunneling. We develop detailed semiclassical analyses in the complex domain and demonstrate that the key trajectories contributing to the fringed tunneling are quite different from the instanton-like trajectories, which works only in unperturbed or weakly-perturbed regimes. Indeed, such trajectories are traveling in the complex phase space being guided by the complexified stable and unstable manifolds which exhibit heteroclinic-like entanglement with the incoming wave manifolds.

Journal of the Physical Society of Japan•Jan 1, 2003•0 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
There has been a surprising degree of uncertainty about how, or even whether, the standard semiclassical theory of Van Vleck and Gutzwiller (VVG) can be used to treat the time-dependent dynamics of tunneling through barriers. At issue is how semiclassical theory describes the time evolution of tunneling. This uncertainty has hindered development of practical techniques for semiclassical calculations of multidimensional systems. It is shown here that the VVG theory can indeed be applied to barrier tunneling, but generally requires consideration of classical trajectories with both complex initial conditions and complex time paths, even for systems with one degree of freedom. The present approach is verified by applying it to the classically forbidden propagation of a wave packet through a one-dimensional barrier---a very basic problem that cannot be treated by the VVG theory without the modifications described here.
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.
In this paper we consider the time evolution of wave-packet tunneling through a potential barrier. Using a path-summation approach, we can describe the quantum evolution in terms of particle trajectories both in the allowed and forbidden regions. Near the barrier edges, where the potential changes rapidly, the wave aspects dominate. These make the trajectories branch into infinite families of paths, which have to be included in the summation. Combining these features, we treat correctly the complementary aspects of the full quantum process, and this allows us to obtain a computational method that gives accurate numerical results with considerably less expenditure of computer time than a direct integration of the Schr\"odinger equation. We test the method on a simple rectangular potential by comparing the calculations with both conventional wave scattering results and computations from the Schr\"odinger equation. A general potential can be approximated by a staircase, which allows the application of our method. A simple adiabatic extension is shown to work excellently for the case of a periodically modulated barrier.
A semiclassical method of complex trajectories for the calculation of the tunneling exponent in systems with many degrees of freedom is further developed. It is supplemented with an easily implementable technique, which enables one to single out the physically relevant trajectory from the whole set of complex classical trajectories. The method is applied to semiclassical transitions of a bound system through a potential barrier. We find that the properties of physically relevant complex trajectories are qualitatively different in the cases of potential tunneling at low energy and dynamical tunneling at energies exceeding the barrier height. Namely, in the case of high energies, the physically relevant complex trajectories describe tunneling via creation of a state close to the top of the barrier. The method is checked against exact solutions of the Schrodinger equation in a quantum mechanical system of two degrees of freedom.
A semiclassical method of complex trajectories for the calculation of the tunneling exponent n systems with many degrees of freedom is further developed. It is supplemented with an easily implemented technique that enables one to single out the physically relevant trajectory from the whole set of complex classical trajectories. The method is applied to semiclassical transitions of a bound system through a potential barrier. We find that the properties of physically relevant complex trajectories are qualitatively different in the cases of potential tunneling at low energy and dynamical tunneling at energies exceeding the barrier height. Namely, in the case of high energies, the physically relevant complex trajectories describe tunneling via creation of a state close to the top of the barrier. The method is checked against exact solutions of the Schrödinger equation in a quantum-mechanical system of two degrees of freedom.
In recent papers we demonstrated a new class of barrier tunneling phenomena which can be observed in periodically driven one dimensional barrier systems (i.e., 1.5D barrier tunneling systems). In such systems the tunneling component is in general accompanied by remarkable interference fringe patterns as the perturbation strength exceeds a certain characteristic value. Such a phenomenon is quite general and is not restricted to periodically driven 1D systems, and it may be also observed in intrinsic two dimensional autonomous systems. The complex-domain semiclassical theory developed by ourselves was applied to the problem, and we succeeded in reconstructing every detail of the fringed tunneling wavefunction only by using the complexified classical tunneling trajectories. This fact greatly motivates us to investigate the mechanism of the tunneling in terms of “classical mechanics”. The later part of the present paper is devoted to asserting the generality of the “classical mechanism” of the fringed tunneling. We develop detailed semiclassical analyses in the complex domain and demonstrate that the key trajectories contributing to the fringed tunneling are quite different from the instanton-like trajectories, which works only in unperturbed or weakly-perturbed regimes. Indeed, such trajectories are traveling in the complex phase space being guided by the complexified stable and unstable manifolds which exhibit heteroclinic-like entanglement with the incoming wave manifolds.
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.