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 trajectories sans isochrones: Quantum barrier scattering with rectilinear constant velocity trajectories

Brad A. Rowland, Róbert E. Wyatt

Published Oct 23, 2007• 28 citations
Publisher

Abstract

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.

DOI
10.1063/1.2790006

Related papers

Complex quantum trajectories for barrier scattering

Brad A. Rowland

We have directed much attention towards developing quantum trajectory methods which can accurately predict the transmission probabilities for a variety of quantum mechanical barrier scattering processes. One promising method involves solving the complex quantum Hamilton-Jacobi equation with the Derivative Propagation Method (DPM). We present this method, termed complex valued DPM (CVDPM(n)). CVDPM(n) has been successfully employed in the Lagrangian frame to accurately compute transmission probabilities on 'thick' one dimensional Eckart and Gaussian potential surfaces. CVDPM(n) is able to reproduce accurate results with a much lower order of approximation than is required by real valued quantum trajectory methods, from initial wave packet energies ranging from the tunneling case (E[subscript o]=0) to high energy cases (twice the barrier height). We successfully extended CVDPM(n) to two-dimensional problems (one translational degree of freedom representing an Eckart or Gaussian barrier coupled to a vibrational degree of freedom) in the Lagrangian framework with great success. CVDPM helps to explain why barrier scattering from "thick" barriers is a much more well posed problem than barrier scattering from "thin" barriers. Though results in these two cases are in very good agreement with grid methods, the search for an appropriate set of initial conditions (termed an 'isochrone) from which to launch the trajectories leads to a time-consuming search problem that is reminiscent of the rootsearching problem from semi-classical dynamics. In order to circumvent the isochrone problem, we present CVDPM(n) equations of motion which are derived and implemented in the arbitrary Lagrangian-Eulerian frame for a metastable potential as well as the Eckart and Gaussian surfaces. In this way, the isochrone problem can be circumvented but at the cost of introducing other computational difficulties. In order to understand why CVDPM may give better transmission probabilities than real valued counterparts, much attention we have been studying and applying numerical analytic continuation techniques to visualize complex-extended wave packets as well as the complex-extended quantum potential. Numerical analytic continuation techniques have also been used to analytically continue a discrete real-valued potential into the complex plane for CVDPM with very promising results.

•Dec 1, 2007•0 citations

Barrier scattering with complex-valued quantum trajectories: Taxonomy and analysis of isochrones

Julianne K. David, Róbert E. Wyatt

To facilitate the search for isochrones when using complex-valued trajectory methods for quantum barrier scattering calculations, the structure and shape of isochrones in the complex plane were studied. Isochrone segments were categorized based on their distinguishing features, which are shared by each situation studied: High and low energy wave packets, scattering from both thick and thin Gaussian and Eckart barriers of varying height. The characteristic shape of the isochrone is a trifurcated system: Trajectories that transmit the barrier are launched from the lower branch (T), while the middle and upper branches form the segments for reflected trajectories (F and B). In addition, a model is presented for the curved section of the lower branch (from which transmitted trajectories are launched), and important features of the complex extension of the initial wave packet are identified.

The Journal of Chemical Physics•Mar 4, 2008•16 citations

Analysis of Barrier Scattering with Real and Complex Quantum Trajectories

Brad A. Rowland, Róbert E. Wyatt

In this study, an analysis of the one-dimensional Eckart and Gaussian barrier scattering problems is undertaken using approximate quantum trajectories. Individual quantum trajectories are computed using the derivative propagation method (DPM). Both real-valued and complex-valued DPM quantum trajectories are employed. Of interest are the deep tunneling and the higher energy barrier scattering problems in cases in which the scattering barrier is "thick" by comparison to the width of the initial wave packet. For higher energy scattering problems, it is found that real-valued DPM trajectories very accurately reproduce the transmitted probability densities at low orders when compared to large fixed-grid calculations. However, higher orders must be introduced to obtain good probabilities for deep tunneling problems. Complex-valued DPM is found to accurately reproduce transmitted probability densities at low order for both the deep tunneling and the higher energy scattering problems. Of particular note, complex-classical trajectories are found to very nearly give the exact result for the deep barrier tunneling scattering problem, and the complex DPM converges well at high orders for these thick barrier scattering problems. A variety of analyses are performed to elucidate the dynamics of complex-valued DPM trajectories. The complex-extended barrier potentials are examined in detail, including an analysis of the complex force. Of particular interest are initial conditions for complex-valued DPM trajectories known as isochrones. All trajectories launched from an isochrone arrive on the real axis on the transmitted side of the barrier at the same time. The computation and properties of isochrones as well as the behavior of the initial wave packet in the complex plane are also examined.

The Journal of Physical Chemistry A•Jul 24, 2007•39 citations

Quantum trajectories in complex phase space: Multidimensional barrier transmission

Róbert E. Wyatt, Brad A. Rowland

The quantum Hamilton-Jacobi equation for the action function is approximately solved by propagating individual Lagrangian quantum trajectories in complex-valued phase space. Equations of motion for these trajectories are derived through use of the derivative propagation method (DPM), which leads to a hierarchy of coupled differential equations for the action function and its spatial derivatives along each trajectory. In this study, complex-valued classical trajectories (second order DPM), along which is transported quantum phase information, are used to study low energy barrier transmission for a model two-dimensional system involving either an Eckart or Gaussian barrier along the reaction coordinate coupled to a harmonic oscillator. The arrival time for trajectories to reach the transmitted (product) region is studied. Trajectories launched from an "equal arrival time surface," defined as an isochrone, all reach the real-valued subspace in the transmitted region at the same time. The Rutherford-type diffraction of trajectories around poles in the complex extended Eckart potential energy surface is described. For thin barriers, these poles are close to the real axis and present problems for computing the transmitted density. In contrast, for the Gaussian barrier or the thick Eckart barrier where the poles are further from the real axis, smooth transmitted densities are obtained. Results obtained using higher-order quantum trajectories (third order DPM) are described for both thick and thin barriers, and some issues that arise for thin barriers are examined.

The Journal of Chemical Physics•Jul 24, 2007•45 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

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
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.
We have directed much attention towards developing quantum trajectory methods which can accurately predict the transmission probabilities for a variety of quantum mechanical barrier scattering processes. One promising method involves solving the complex quantum Hamilton-Jacobi equation with the Derivative Propagation Method (DPM). We present this method, termed complex valued DPM (CVDPM(n)). CVDPM(n) has been successfully employed in the Lagrangian frame to accurately compute transmission probabilities on 'thick' one dimensional Eckart and Gaussian potential surfaces. CVDPM(n) is able to reproduce accurate results with a much lower order of approximation than is required by real valued quantum trajectory methods, from initial wave packet energies ranging from the tunneling case (E[subscript o]=0) to high energy cases (twice the barrier height). We successfully extended CVDPM(n) to two-dimensional problems (one translational degree of freedom representing an Eckart or Gaussian barrier coupled to a vibrational degree of freedom) in the Lagrangian framework with great success. CVDPM helps to explain why barrier scattering from "thick" barriers is a much more well posed problem than barrier scattering from "thin" barriers. Though results in these two cases are in very good agreement with grid methods, the search for an appropriate set of initial conditions (termed an 'isochrone) from which to launch the trajectories leads to a time-consuming search problem that is reminiscent of the rootsearching problem from semi-classical dynamics. In order to circumvent the isochrone problem, we present CVDPM(n) equations of motion which are derived and implemented in the arbitrary Lagrangian-Eulerian frame for a metastable potential as well as the Eckart and Gaussian surfaces. In this way, the isochrone problem can be circumvented but at the cost of introducing other computational difficulties. In order to understand why CVDPM may give better transmission probabilities than real valued counterparts, much attention we have been studying and applying numerical analytic continuation techniques to visualize complex-extended wave packets as well as the complex-extended quantum potential. Numerical analytic continuation techniques have also been used to analytically continue a discrete real-valued potential into the complex plane for CVDPM with very promising results.
To facilitate the search for isochrones when using complex-valued trajectory methods for quantum barrier scattering calculations, the structure and shape of isochrones in the complex plane were studied. Isochrone segments were categorized based on their distinguishing features, which are shared by each situation studied: High and low energy wave packets, scattering from both thick and thin Gaussian and Eckart barriers of varying height. The characteristic shape of the isochrone is a trifurcated system: Trajectories that transmit the barrier are launched from the lower branch (T), while the middle and upper branches form the segments for reflected trajectories (F and B). In addition, a model is presented for the curved section of the lower branch (from which transmitted trajectories are launched), and important features of the complex extension of the initial wave packet are identified.
In this study, an analysis of the one-dimensional Eckart and Gaussian barrier scattering problems is undertaken using approximate quantum trajectories. Individual quantum trajectories are computed using the derivative propagation method (DPM). Both real-valued and complex-valued DPM quantum trajectories are employed. Of interest are the deep tunneling and the higher energy barrier scattering problems in cases in which the scattering barrier is "thick" by comparison to the width of the initial wave packet. For higher energy scattering problems, it is found that real-valued DPM trajectories very accurately reproduce the transmitted probability densities at low orders when compared to large fixed-grid calculations. However, higher orders must be introduced to obtain good probabilities for deep tunneling problems. Complex-valued DPM is found to accurately reproduce transmitted probability densities at low order for both the deep tunneling and the higher energy scattering problems. Of particular note, complex-classical trajectories are found to very nearly give the exact result for the deep barrier tunneling scattering problem, and the complex DPM converges well at high orders for these thick barrier scattering problems. A variety of analyses are performed to elucidate the dynamics of complex-valued DPM trajectories. The complex-extended barrier potentials are examined in detail, including an analysis of the complex force. Of particular interest are initial conditions for complex-valued DPM trajectories known as isochrones. All trajectories launched from an isochrone arrive on the real axis on the transmitted side of the barrier at the same time. The computation and properties of isochrones as well as the behavior of the initial wave packet in the complex plane are also examined.
The quantum Hamilton-Jacobi equation for the action function is approximately solved by propagating individual Lagrangian quantum trajectories in complex-valued phase space. Equations of motion for these trajectories are derived through use of the derivative propagation method (DPM), which leads to a hierarchy of coupled differential equations for the action function and its spatial derivatives along each trajectory. In this study, complex-valued classical trajectories (second order DPM), along which is transported quantum phase information, are used to study low energy barrier transmission for a model two-dimensional system involving either an Eckart or Gaussian barrier along the reaction coordinate coupled to a harmonic oscillator. The arrival time for trajectories to reach the transmitted (product) region is studied. Trajectories launched from an "equal arrival time surface," defined as an isochrone, all reach the real-valued subspace in the transmitted region at the same time. The Rutherford-type diffraction of trajectories around poles in the complex extended Eckart potential energy surface is described. For thin barriers, these poles are close to the real axis and present problems for computing the transmitted density. In contrast, for the Gaussian barrier or the thick Eckart barrier where the poles are further from the real axis, smooth transmitted densities are obtained. Results obtained using higher-order quantum trajectories (third order DPM) are described for both thick and thin barriers, and some issues that arise for thin barriers are examined.
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.
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.