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

arXiv:2605.04139

Tunneling from an oscillating state in quantum mechanics

Oliver Janssen, Matthew Kleban, Cameron Norton

quant-phhep-thSubmitted May 5, 2026• Announced May 7, 2026• 0 citations• Revised Oct 29, 2026 (v2)
v1 · May 5, 2026v2 · Oct 29, 2026
arXivHTMLPDF

Abstract

We study the decay of general initial states out of a metastable potential well in quantum mechanics. We provide a closed-form expression for the probability current that tunnels through the barrier in terms of the resonant states into which the initial state can be decomposed. All ingredients are computed analytically to first subleading order in the semiclassical limit. Specializing to a coherently-oscillating state, we derive an equivalent semiclassical approximation to the current via the path integral, which circumvents the decomposition into resonant states. We test the accuracy of the approximation by comparing it to a numerical solution of the Schrödinger equation.

Comments
29 pages + 3 appendices. v2: major revisions, including a generalization of the main result to oscillations of arbitrary amplitude in the well, a path integral derivation of the main result, and the discussion of small oscillations where a simple closed-form expression is found

Related papers

Time-dependent tunnelling via path integrals. Connection to results of the quantum mechanics of decaying states

Theodosios Geo Douvropoulos, Cleanthes A. Nicolaides

The probability, P ( t ), of the irreversible dissipation into a continuous spectrum of an initially ( t = 0) localized (Ψ o ) nonstationary state acquires, as time increases, 'memory' due to the lower energy bound of the spectrum, and eventually follows a nonexponential decay (NED). Regardless of the degree of dependence on energy, the magnitude of this deviation from exponential decay depends on the degree of proximity to threshold, and on whether the theory employs a real energy distribution, one form of which is g ( E ) ≡ ⟨Ψ o |δ( H − E )|Ψ o ⟩, or a complex energy distribution, G ( E ) ≡ ⟨Ψ o |( H − E + i0)|Ψ o ⟩. It is the latter that is physically consistent, since it originates from the singularity at t = 0, which breaks the S -matrix unitarity, in accordance with the non-Hermitian character of decaying states. In order to test the quantum mechanical theory, we carried out semiclassical path integral calculations of the P ( t ) for an isolated narrow tunnelling state, whereby the truncated Fourier transform of a semiclassical Green function, G sc ( E ), is obtained. The results are in agreement with the analytic results of quantum mechanics when energy and time asymmetry are taken into account. It is shown that the analytic structure of G sc ( E ) is [ D regular + D pole ], where D pole is a finite sum over complex poles, which are the complex eigenvalues, W n , that the potential can support. The W n are given by E n + Δ n − (i/2)Γ n , where E n are the real eigenvalues of the corresponding bound potential, Γ n are the energy widths and Δ n are the energy shifts, both expressed in terms of computable semiclassical quantities. The spherical harmonic oscillator (SHO) with and without angular momentum, and unstable ground states of diatomic molecules, are treated as particular cases. The exact spectrum of the SHO is recovered only when the Kramers–Langer semiclassical expression for the centrifugal potential is used, thereby bypassing the difficulty of the singularity at r = 0. The spectrum from the use of the quantal form l ( l + 1) reduces to that of l ( l + 1/2) 2 in the limit of large l , i.e., for orbits far from r = 0. Using previously computed energies and widths for the vibrational levels of He 2 2+ 1σ g 2 1 Σ g + , the application of two formulae for P ( t ), one derived from a Lorentzian real energy distribution and the other from the corresponding complex energy distribution, shows that, for the lowest level, in the former case NED starts after about 193 lifetimes, and in the latter after about 102 lifetimes. The fact that this difference is large should have consequences for the deeper understanding of irreversibility at the quantum level.

Journal of Physics B Atomic Molecular and Optical Physics•Oct 24, 2002•10 citations

Tunneling of Bound Systems at Finite Energies: Complex Paths Through Potential Barriers

Gian Franco Bonini, Andrew G. Cohen, Claudio Rebbi B +1 more

We adapt the semiclassical technique, as used in the context of instanton transitions in quantum field theory, to the description of tunneling transmissions at finite energies through potential barriers by complex quantum mechanical systems. Even for systems initially in their ground state, not generally describable in semiclassical terms, the transmission probability has a semiclassical (exponential) form. The calculation of the tunneling exponent uses analytic continuation of degrees of freedom into a complex phase space as well as analytic continuation of the classical equations of motion into the complex time plane. We test this semiclassical technique by comparing its results with those of a computational investigation of the full quantum mechanical system, finding excellent agreement.

arXiv:quant-ph/9901062•Jan 21, 1999•2 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

Time dependence and line shape of spontaneous quantum tunneling

B. U. Felderhof

The time dependence of spontaneous quantum tunneling of a particle from a potential well through a barrier in one or three dimensions is studied on the basis of the solution of the time-dependent Schrödinger equation as an initial-value problem by the method of Laplace transform. The wavefunction in the outer space is expressed in terms of special functions known from the theory of classical diffusion. A concise expression is derived for the probability current into the outer space. At long times the current oscillates about zero before decaying with a t −4 long-time tail. A comparison is made with the solution obtained by an expansion of the initial wavefunction in terms of the complete set of scattering states and bound states.

Journal of Physics A Mathematical and Theoretical•Oct 7, 2008•2 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

Multidimensional tunneling in terms of complex classical mechanics: Wave functions, energy splittings, and decay rates in nonintegrable systems

Shoji Takada

A semiclassical theory of multidimensional tunneling is formulated to calculate the tunneling wave function, energy splitting in a double well and decay rate constant from a metastable state. First, the tunneling wave function is calculated by analytic continuation of a quantized torus prepared in analytic forms using either the Chapman–Garrett–Miller method or the Birkhoff–Gustavson normal form method. For a weakly nonintegrable system, tunnelings are confirmed to be classified into two qualitatively different domains; pure tunneling in the I region and mixed tunneling in the C region. Semiclassical wave functions agree with quantum mechanical ones within a few percent both in classically allowed and tunneling regions. Breakdown of this simple picture is exemplified for cases of relatively strong couplings. Second, expressions of the tunneling energy splitting in a double well potential and the decay rate from a metastable state are derived. The wave function near a well formulated above is connected with the semiclassical Green’s function in the deep tunneling region. The latter is expressed by complex trajectories which start from the complex quantized torus. A preliminary numerical comparison with the quantum mechanical value is also done for the energy splitting.

The Journal of Chemical Physics•Mar 8, 1996•28 citations
We study the decay of general initial states out of a metastable potential well in quantum mechanics. We provide a closed-form expression for the probability current that tunnels through the barrier in terms of the resonant states into which the initial state can be decomposed. All ingredients are computed analytically to first subleading order in the semiclassical limit. Specializing to a coherently-oscillating state, we derive an equivalent semiclassical approximation to the current via the path integral, which circumvents the decomposition into resonant states. We test the accuracy of the approximation by comparing it to a numerical solution of the Schrödinger equation.
The probability, P ( t ), of the irreversible dissipation into a continuous spectrum of an initially ( t = 0) localized (Ψ o ) nonstationary state acquires, as time increases, 'memory' due to the lower energy bound of the spectrum, and eventually follows a nonexponential decay (NED). Regardless of the degree of dependence on energy, the magnitude of this deviation from exponential decay depends on the degree of proximity to threshold, and on whether the theory employs a real energy distribution, one form of which is g ( E ) ≡ ⟨Ψ o |δ( H − E )|Ψ o ⟩, or a complex energy distribution, G ( E ) ≡ ⟨Ψ o |( H − E + i0)|Ψ o ⟩. It is the latter that is physically consistent, since it originates from the singularity at t = 0, which breaks the S -matrix unitarity, in accordance with the non-Hermitian character of decaying states. In order to test the quantum mechanical theory, we carried out semiclassical path integral calculations of the P ( t ) for an isolated narrow tunnelling state, whereby the truncated Fourier transform of a semiclassical Green function, G sc ( E ), is obtained. The results are in agreement with the analytic results of quantum mechanics when energy and time asymmetry are taken into account. It is shown that the analytic structure of G sc ( E ) is [ D regular + D pole ], where D pole is a finite sum over complex poles, which are the complex eigenvalues, W n , that the potential can support. The W n are given by E n + Δ n − (i/2)Γ n , where E n are the real eigenvalues of the corresponding bound potential, Γ n are the energy widths and Δ n are the energy shifts, both expressed in terms of computable semiclassical quantities. The spherical harmonic oscillator (SHO) with and without angular momentum, and unstable ground states of diatomic molecules, are treated as particular cases. The exact spectrum of the SHO is recovered only when the Kramers–Langer semiclassical expression for the centrifugal potential is used, thereby bypassing the difficulty of the singularity at r = 0. The spectrum from the use of the quantal form l ( l + 1) reduces to that of l ( l + 1/2) 2 in the limit of large l , i.e., for orbits far from r = 0. Using previously computed energies and widths for the vibrational levels of He 2 2+ 1σ g 2 1 Σ g + , the application of two formulae for P ( t ), one derived from a Lorentzian real energy distribution and the other from the corresponding complex energy distribution, shows that, for the lowest level, in the former case NED starts after about 193 lifetimes, and in the latter after about 102 lifetimes. The fact that this difference is large should have consequences for the deeper understanding of irreversibility at the quantum level.
We adapt the semiclassical technique, as used in the context of instanton transitions in quantum field theory, to the description of tunneling transmissions at finite energies through potential barriers by complex quantum mechanical systems. Even for systems initially in their ground state, not generally describable in semiclassical terms, the transmission probability has a semiclassical (exponential) form. The calculation of the tunneling exponent uses analytic continuation of degrees of freedom into a complex phase space as well as analytic continuation of the classical equations of motion into the complex time plane. We test this semiclassical technique by comparing its results with those of a computational investigation of the full quantum mechanical system, finding excellent agreement.
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.
The time dependence of spontaneous quantum tunneling of a particle from a potential well through a barrier in one or three dimensions is studied on the basis of the solution of the time-dependent Schrödinger equation as an initial-value problem by the method of Laplace transform. The wavefunction in the outer space is expressed in terms of special functions known from the theory of classical diffusion. A concise expression is derived for the probability current into the outer space. At long times the current oscillates about zero before decaying with a t −4 long-time tail. A comparison is made with the solution obtained by an expansion of the initial wavefunction in terms of the complete set of scattering states and bound states.
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.
A semiclassical theory of multidimensional tunneling is formulated to calculate the tunneling wave function, energy splitting in a double well and decay rate constant from a metastable state. First, the tunneling wave function is calculated by analytic continuation of a quantized torus prepared in analytic forms using either the Chapman–Garrett–Miller method or the Birkhoff–Gustavson normal form method. For a weakly nonintegrable system, tunnelings are confirmed to be classified into two qualitatively different domains; pure tunneling in the I region and mixed tunneling in the C region. Semiclassical wave functions agree with quantum mechanical ones within a few percent both in classically allowed and tunneling regions. Breakdown of this simple picture is exemplified for cases of relatively strong couplings. Second, expressions of the tunneling energy splitting in a double well potential and the decay rate from a metastable state are derived. The wave function near a well formulated above is connected with the semiclassical Green’s function in the deep tunneling region. The latter is expressed by complex trajectories which start from the complex quantized torus. A preliminary numerical comparison with the quantum mechanical value is also done for the energy splitting.