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

Complex time paths for semiclassical wave packet propagation with complex trajectories

Jakob Petersen, Kenneth G. Kay

Published Aug 7, 2014• 9 citations
Publisher

Abstract

The use of complex-valued trajectories in semiclassical wave packet methods can lead to problems that prevent calculation of the wave function in certain regions of the configuration space. We investigate this so-called bald spot problem in the context of generalized Gaussian wave packet dynamics. The analysis shows that the bald spot phenomenon is essentially due to the complex nature of the initial conditions for the trajectories. It is, therefore, expected to be a general feature of several semiclassical methods that rely on trajectories with such initial conditions. A bald region is created when a trajectory, needed to calculate the wave function at a given time, reaches a singularity of the potential energy function in the complex plane at an earlier, real time. This corresponds to passage of a branch point singularity across the real axis of the complex time plane. The missing portions of the wave function can be obtained by deforming the time path for the integration of the equations of motion into the complex plane so that the singularity is circumvented. We present examples of bald spots, singularity times, and suitable complex time paths for one-dimensional barrier transmission in the Eckart and Gaussian systems. Although the bald regions for the Eckart system are often localized, they are found to be semi-infinite for the Gaussian system. For the case of deep tunneling, the bald regions for both systems may encompass the entire portion of space occupied by the transmitted wave packet. Thus, the use of complex time paths becomes essential for a treatment of barrier tunneling.

DOI
10.1063/1.4891918

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

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

Quantum trajectories in complex space

Chia‐Chun Chou, Róbert E. Wyatt

Quantum trajectories in complex space in the framework of the quantum Hamilton-Jacobi formalism are investigated. For time-dependent problems, the complex quantum trajectories determined from the exact analytical wave function for the free Gaussian wave packet and the coherent state in the harmonic potential are used to demonstrate that the information transported by the particles in the complex space can be used to synthesize the time-dependent wave function on the real axis. For time-independent problems, the exact complex quantum trajectories for the Eckart potential are obtained by numerically integrating the equations of motion. The unusual structure of the total potential (the sum of the classical and the quantum potentials) for the stationary states for the Eckart potential is pointed out. The variations of the complex-valued kinetic energy, classical potential, and quantum potential along the complex quantum trajectories are analyzed. This paper not only analyzes complex quantum trajectories for time-dependent and time-independent systems but also provides a unified description for complex quantum trajectories for one-dimensional problems in the quantum Hamilton-Jacobi formalism.

Physical Review A•Jul 20, 2007•45 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

Computational Investigation of Wave Packet Scattering in the Complex Plane: Dynamics of Exact Quantum Trajectories

Róbert E. Wyatt, Brad A. Rowland

In two previous studies, the time-dependent scattering of a wave packet from a Gaussian barrier was investigated computationally in the complex z-plane. One of these involved the 'direct' propagation of the wave packet in the complex space, and the other used numerical analytic continuation techniques to generate the dynamics in the complex plane from the wave function computed on the real-axis. In the current study, the dynamics of exact quantum trajectories are analyzed for the same barrier scattering problem. Thousands of quantum trajectories were launched from positions near the center of the initial wave packet. These trajectories were computed by integrating equations-of-motion involving the quantum momentum function, which was obtained from the time-dependent wave function and its derivative. In order to analyze the dynamics, many trajectories were plotted on space-time diagrams. Particular emphasis was placed upon trajectories undergoing reflection in the barrier region. Some groups of strongly correlated trajectories form long-lived highly organized patterns, including helical wrappings around a series of stagnation filaments. These curves alternate with quasi-nodes where the amplitude of the wave function reaches low values. In addition, other trajectories for short times follow hyperbolic paths as they propagate near vorticity tubes surrounding these quasi-nodes.

Journal of Chemical Theory and Computation•Jan 29, 2009•4 citations

Complex trajectories sans isochrones: Quantum barrier scattering with rectilinear constant velocity trajectories

Brad A. Rowland, Róbert E. Wyatt

One of the major obstacles in employing complex-valued trajectory methods for quantum barrier scattering calculations is the search for isochrones. In this study, complex-valued derivative propagation method trajectories in the arbitrary Lagrangian-Eulerian frame are employed to solve the complex Hamilton-Jacobi equation for quantum barrier scattering problems employing constant velocity trajectories moving along rectilinear paths whose initial points can be in the complex plane or even along the real axis. It is shown that this effectively removes the need for isochrones for barrier transmission problems. Model problems tested include the Eckart, Gaussian, and metastable quadratic+cubic potentials over a variety of wave packet energies. For comparison, the "exact" solution is computed from the time-dependent Schrodinger equation via pseudospectral methods.

The Journal of Chemical Physics•Oct 23, 2007•28 citations
The use of complex-valued trajectories in semiclassical wave packet methods can lead to problems that prevent calculation of the wave function in certain regions of the configuration space. We investigate this so-called bald spot problem in the context of generalized Gaussian wave packet dynamics. The analysis shows that the bald spot phenomenon is essentially due to the complex nature of the initial conditions for the trajectories. It is, therefore, expected to be a general feature of several semiclassical methods that rely on trajectories with such initial conditions. A bald region is created when a trajectory, needed to calculate the wave function at a given time, reaches a singularity of the potential energy function in the complex plane at an earlier, real time. This corresponds to passage of a branch point singularity across the real axis of the complex time plane. The missing portions of the wave function can be obtained by deforming the time path for the integration of the equations of motion into the complex plane so that the singularity is circumvented. We present examples of bald spots, singularity times, and suitable complex time paths for one-dimensional barrier transmission in the Eckart and Gaussian systems. Although the bald regions for the Eckart system are often localized, they are found to be semi-infinite for the Gaussian system. For the case of deep tunneling, the bald regions for both systems may encompass the entire portion of space occupied by the transmitted wave packet. Thus, the use of complex time paths becomes essential for a treatment of barrier tunneling.
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 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.
Quantum trajectories in complex space in the framework of the quantum Hamilton-Jacobi formalism are investigated. For time-dependent problems, the complex quantum trajectories determined from the exact analytical wave function for the free Gaussian wave packet and the coherent state in the harmonic potential are used to demonstrate that the information transported by the particles in the complex space can be used to synthesize the time-dependent wave function on the real axis. For time-independent problems, the exact complex quantum trajectories for the Eckart potential are obtained by numerically integrating the equations of motion. The unusual structure of the total potential (the sum of the classical and the quantum potentials) for the stationary states for the Eckart potential is pointed out. The variations of the complex-valued kinetic energy, classical potential, and quantum potential along the complex quantum trajectories are analyzed. This paper not only analyzes complex quantum trajectories for time-dependent and time-independent systems but also provides a unified description for complex quantum trajectories for one-dimensional problems in the quantum Hamilton-Jacobi formalism.
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 two previous studies, the time-dependent scattering of a wave packet from a Gaussian barrier was investigated computationally in the complex z-plane. One of these involved the 'direct' propagation of the wave packet in the complex space, and the other used numerical analytic continuation techniques to generate the dynamics in the complex plane from the wave function computed on the real-axis. In the current study, the dynamics of exact quantum trajectories are analyzed for the same barrier scattering problem. Thousands of quantum trajectories were launched from positions near the center of the initial wave packet. These trajectories were computed by integrating equations-of-motion involving the quantum momentum function, which was obtained from the time-dependent wave function and its derivative. In order to analyze the dynamics, many trajectories were plotted on space-time diagrams. Particular emphasis was placed upon trajectories undergoing reflection in the barrier region. Some groups of strongly correlated trajectories form long-lived highly organized patterns, including helical wrappings around a series of stagnation filaments. These curves alternate with quasi-nodes where the amplitude of the wave function reaches low values. In addition, other trajectories for short times follow hyperbolic paths as they propagate near vorticity tubes surrounding these quasi-nodes.
One of the major obstacles in employing complex-valued trajectory methods for quantum barrier scattering calculations is the search for isochrones. In this study, complex-valued derivative propagation method trajectories in the arbitrary Lagrangian-Eulerian frame are employed to solve the complex Hamilton-Jacobi equation for quantum barrier scattering problems employing constant velocity trajectories moving along rectilinear paths whose initial points can be in the complex plane or even along the real axis. It is shown that this effectively removes the need for isochrones for barrier transmission problems. Model problems tested include the Eckart, Gaussian, and metastable quadratic+cubic potentials over a variety of wave packet energies. For comparison, the "exact" solution is computed from the time-dependent Schrodinger equation via pseudospectral methods.