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

Tunneling dynamics and spawning with adaptive semiclassical wave packets

Vasile Catrinel Grădinaru, George A. Hagedorn, Alain Joye

Published May 13, 2010• 18 citations
Publisher

Abstract

Tunneling through a one-dimensional Eckart barrier is investigated using a recently developed propagation scheme based on semiclassical wave packets. This version of the time-dependent discrete variable representation method yields linear equations for the parameters, is fully adaptive, and does not require a frozen ansatz in order to approximate the exact solution of the Schrödinger equation accurately. We rely on an analytical result to derive a new algorithm to spawn a second family of semiclassical wave packets after the tunneling has occurred. Numerical results for a benchmark problem demonstrate the accuracy of the new method.

DOI
10.1063/1.3429607

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 Real-Time Tunneling by Multiple Spawning of Classical Trajectories

Frank Großmann

We systematically extend semiclassical real-time propagation methods based on Gaussian wave packets by using the composition property of the time-dependent quantum-mechanical Green function. For the first time, in an application to the problem of barrier tunneling, a semiclassical time-domain calculation of the transmission probability exhibits good agreement with exact quantum mechanical values at energies below the barrier top. Only two insertions of unity are needed for these results in a benchmark model of reactive scattering.

Physical Review Letters•Jul 31, 2000•42 citations

Semiclassical tunneling in real time: Wave-packet dynamics in static and driven barrier potentials

Markus Saltzer, Joachim Ankerhold

We study to what extent propagation of wave packets through potential barriers can be described semiclassically in real time. By analyzing the semiclassical approximation to the exact path-integral representation of the quantum-mechanical time evolution operator, we show that in the long-time sector it is dominated by quasiclassical fluctuations associated with caustics. Contrary to the description of tunneling in imaginary times, simple tunneling paths even for complex energies do not exist. Our findings lead to an extended semiclassical propagation scheme for wave-packet dynamics which is numerical efficient and accurately describes deep incoherent tunneling through static and externally driven barrier potentials.

Physical Review A•Oct 13, 2003•20 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

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

Description of tunneling with the help of coupled frozen Gaussians

Dmitrii V. Shalashilin, Mark Sheard Child

A previous initial value coherent state representation is applied to Gaussian wave packet propagation on one-dimensional and two-dimensional double well potentials. The method uses a basis of trajectory guided frozen Gaussian coherent states, sampled from a Monte Carlo ensemble, for numerical solutions of the quantum time dependent integro-differential Schrödinger equation. Accurate tunneling splittings are obtained. Comparisons are also made with corresponding results obtained by the Herman–Kluk semiclassical initial value representation wave packet method.

The Journal of Chemical Physics•Jun 1, 2001•53 citations
Tunneling through a one-dimensional Eckart barrier is investigated using a recently developed propagation scheme based on semiclassical wave packets. This version of the time-dependent discrete variable representation method yields linear equations for the parameters, is fully adaptive, and does not require a frozen ansatz in order to approximate the exact solution of the Schrödinger equation accurately. We rely on an analytical result to derive a new algorithm to spawn a second family of semiclassical wave packets after the tunneling has occurred. Numerical results for a benchmark problem demonstrate the accuracy of the new method.
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.
We systematically extend semiclassical real-time propagation methods based on Gaussian wave packets by using the composition property of the time-dependent quantum-mechanical Green function. For the first time, in an application to the problem of barrier tunneling, a semiclassical time-domain calculation of the transmission probability exhibits good agreement with exact quantum mechanical values at energies below the barrier top. Only two insertions of unity are needed for these results in a benchmark model of reactive scattering.
We study to what extent propagation of wave packets through potential barriers can be described semiclassically in real time. By analyzing the semiclassical approximation to the exact path-integral representation of the quantum-mechanical time evolution operator, we show that in the long-time sector it is dominated by quasiclassical fluctuations associated with caustics. Contrary to the description of tunneling in imaginary times, simple tunneling paths even for complex energies do not exist. Our findings lead to an extended semiclassical propagation scheme for wave-packet dynamics which is numerical efficient and accurately describes deep incoherent tunneling through static and externally driven barrier potentials.
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.
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.
A previous initial value coherent state representation is applied to Gaussian wave packet propagation on one-dimensional and two-dimensional double well potentials. The method uses a basis of trajectory guided frozen Gaussian coherent states, sampled from a Monte Carlo ensemble, for numerical solutions of the quantum time dependent integro-differential Schrödinger equation. Accurate tunneling splittings are obtained. Comparisons are also made with corresponding results obtained by the Herman–Kluk semiclassical initial value representation wave packet method.