It is shown how the semiclassical theory of path integrals can be implemented in a practical manner for the analysis of a potential that combines the two-state system of a double well potential (DWP) with decay into a continuous spectrum. This potential may correspond to a variety of physical situations in physics and chemistry. The structure of the formalism and of the results is such that it allows computation not only for analytic but also for numerically given potentials. The central theme is the determination of the energy-dependent Green’s function, which is shown to consist of a regular part and a part containing simple and double complex poles. These poles represent the position of the energy levels, as well as the energy widths and shifts due to the interaction with the continuous spectrum. When applied to the bound DWP without tunneling, the theory is shown to reduce in certain limits to known results from the Jeffreys–Wentzel–Kiamers–Bhrillouin approximation. If the system is taken to be prepared in the first well, the interactions with the remaining of the potential lead to two types of transition rates. One represents the transient motion toward a virtual equilibrium state of the DWP. It emerges as a positive imaginary part of the self-energy. The other represents the decay into the continuum and emerges as a negative imaginary part of the pole. Comparison of the two mechanisms of nonstationarity is made for different magnitudes of the second barrier relative to the first one. Since the system decays to the continuum while oscillating, the theory obtains a correction to the frequency of oscillation in the DWP due to the interaction with the continuum. This phenomenon is observable in real two-state systems, if an external perturbation which affects mainly one state converts it into a resonance state.
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.
Nelson's stochastic mechanics formulates quantum dynamics as a real-time conservative diffusion process in which a particle undergoes Brownian-like motion with a fluctuation amplitude fixed by Planck's constant. While being mathematically equivalent to the Schrödinger formulation, this approach provides an alternative dynamical framework that enables the study of time-resolved quantities that are not straightforwardly defined within the standard operator-based approach. In the present work, Nelson's stochastic mechanics is employed to investigate tunneling-time statistics for bound states in double-well potentials. Using first-passage time theory within this framework, both the mean tunneling time, τˉ, and the full probability distribution, p(τ), are computed. The theoretical predictions are validated through extensive numerical simulations of stochastic trajectories for two representative potentials. For the square double-well potential, analytical expressions for τˉ are derived and are shown to be in excellent agreement with simulations. In the high-barrier limit, the results reveal a direct relation between the stochastic-mechanical and quantum-mechanical tunneling times, expressed as τQM=(π/2)τˉ, where τQM corresponds to half the oscillation period of the probability of finding the particle in either well. This relation is further confirmed for generic double-well systems through a WKB analysis. As a concrete application, the inversion dynamics of the ammonia molecule is analyzed, yielding an inversion frequency of approximately 24 GHz, in close agreement with experimental observations. These results highlight the potential of stochastic mechanics as a conceptually coherent and quantitatively consistent framework for analyzing tunneling phenomena in quantum systems.
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 movement of a proton which moves in a one-dimensional effective double-minimum potential V(x) and can pass by means of the tunnel effect from one well into the other was investigated. The lingering times of the proton in the potential wells (a) and (b) were calculated for more than 150 potential models as functions of the height Vmax of the potential barrier separating the two wells, the energy difference ΔV, and separation d of the minima. In the determination of the lingering times, it was assumed that the proton at time t = 0 is to be found with certainty in one of the wells and can be represented by means of a statistical operator related to the eigenstates of a one-well Hamiltonian H(a)orH(b), respectively. The times determined on this basis were compared with those which result with a quasiclassical method from the classical frequency of oscillation of a particle in a well and from the penetration coefficient of the potential barrier. The deviation between the values determined with both methods is largest if an eigenvalue of H(a) becomes approximately equal to an eigenvalue of H(b). For the potential model V(x) a polynomial of fourth degree in the position coordinate x was chosen. The eigenfunctions of H(a)andH(b) as well as those of the effective Hamiltonian H corresponding to V(x) were represented as linear combinations of eigenfunctions of the harmonic oscillator with angular frequency ω. ω was chosen so that the mean square deviations of the approximated eigenvalues of H became as small as possible for the states of lowest energy.
In some molecules such as ammonia, double vibrations are observed at small energies. These small vibrations are thought to arise from double well potentials. To explain these small energy gaps, double potential wells are investigated, and complex calculations are made. Until recently, since the exact solution of the Schrödinger equation was not known, approximate solutions were always obtained. The exact solution of the Schrödinger equation is now known. In previous studies, a simple procedure for the general solution of the radial Schrödinger Equation has been found for spherical symmetric potentials without making any approximation. In this article, the Schrodinger equation was solved with this new solution method. It has been applied to the oscillations through a barrier tunneling potential for a particle in a central quadratic deformed potential well (double potential well). It is seen that there is no need to look for complex potentials to explain the double wells observed in small vibrations in ammonia and similar molecules, they were compared with experimentally measured energies. The results were found to be fully compatible. It has been seen that there can always be double well potentials at small energies for all molecules, atomic nuclei, and similar particles with angular momentum other than zero.