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

The calculation of multidimensional semiclassical wave functions in the forbidden region using real valued coordinates

Michael F. Herman

Published Sep 20, 2010• 3 citations
Publisher

Abstract

A method that uses only real valued coordinates is presented for integrating the many dimensional semiclassical wave function into the forbidden region. The procedure first determines a surface of caustic points by running the set of trajectories that define the wave function in the allowed region. In the forbidden region, the momentum and the action integral are both complex functions of position, and their imaginary parts vanish on the caustic surface. The direction of the imaginary part of the momentum p(I) can be chosen to the perpendicular to the caustic surface at all points on that surface. Equations are derived for integrating the values of the real and imaginary parts of the momentum along the curves that follow the direction of p(I). The equations for the change in the action integral and the prefactor for the semiclassical wave function along these curves are also obtained, allowing for the determination of the semiclassical wave function in the forbidden region. Calculations are performed for a two dimensional problem, and the semiclassical wave function is found to be is excellent agreement with the results of exact quantum calculations.

DOI
10.1063/1.3490087

Related papers

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

Semiclassical construction of Green's functions for multidimensional systems with barrier tunneling

Kenneth G. Kay

We develop a treatment that divides phase space into regions bounded by surfaces and forms the energy-dependent semiclassical Green's function from short classical trajectory segments confined to the individual regions. We use this to generalize Bogomolny's quantization method to systems with barrier tunneling. Combined with an appropriate treatment of tunneling amplitudes, this yields a unified semiclassical framework that uses ordinary, real trajectories in classically allowed regions and imaginary-time trajectories in forbidden regions to obtain energies, decay rates, and wave functions for arbitrary states of multidimensional systems undergoing tunneling. Because the trajectory segments in the allowed region are short, components in the Green's function expression can be numerically evaluated using semiclassical initial value methods, yielding a promising approach for practical calculations.

Physical Review A•Oct 27, 2017•5 citations

Quasilinearization method applied to multidimensional quantum tunneling

Mohsen Razavy, Vincent J. Cote

We apply the quasilinearization method of Bellman and Kalaba [Quasilinearization and Nonlinear Boundary-Value Problems (Elsevier, New York, 1965)] to find approximate solutions for the multidimensional quantum tunneling for separable as well as nonseparable wave equations. By introducing the idea of the complex ``semiclassical trajectory'' which is valid for the motion over and under the barrier, and which, in the proper limit, reduces to the real classical trajectory in the allowed region, we obtain an eigenvalue equation for the characteristic wave numbers. This eigenvalue equation is similar to the corresponding equation obtained from the WKB approximation and yields complex eigenvalues with negative imaginary parts. When the barrier changes very rapidly as a function of the radial distance, we can replace the concept of the semiclassical trajectory, which may not be applicable in this case, by the concept of a complex ``quantum trajectory.'' The trajectory defined either way depends on a constant of integration, and by minimizing the action with respect to this constant we can obtain the minimum escape path. The case of two-dimensional tunneling is discussed as an example of this method.

Physical Review A•Apr 1, 1994•4 citations

Semiclassical propagation: Phase indices and the initial-value formalism

Giuseppe Campolieti, Paul W. Brumer

Initial-value phase-space integral representations for a time-dependent propagator are obtained in the coordinate and momentum representations. To do so we first derive nonuniform semiclassical propagators for the various representations, obtaining the global-time semiclassical phase indices (Maslov indices) due to caustics. Results include readily implementable general phase index formulas for any type of caustic, including cases where the Morse index theorem is inapplicable. The method of obtaining the indices is general and based simply on concatenating uniform short-time propagators which also gives rise to alternative path-integral forms. Initial-value integral representations are then derived by introducing a method of extending short-time initial-value propagator formulas to global times via a simple stationary-phase asympotoic-equivalence approach. The integrals reduce to the nonuniform semiclassical propagators within the stationary-phase approximation, are uniform about caustics, and have integrand phases which properly account for the global-time phases in terms of appropriate Maslov indices. The initial-value integrals are also consistently derived via a canonical mapping procedure on the coordinate-space path integral. Initial-value integrals for time-dependent wave-function propgation are also given. Evaluation of the initial-value integral expressions do not require trajectory root searches for propagation.

Physical Review A•Aug 1, 1994•132 citations

Wentzel-Kramers-Brillouin method in multidimensional tunneling

Z.-H. Huang, T. E. Feuchtwang, Paul H. Cutler +1 more

The WKB method is commonly used in semiclassical approximations to the wave function in both the classically allowed and the forbidden regions of a one-dimensional potential. In a multidimensional space, the method can be adapted to construct wave functions in an ``allowed'' region from classical trajectories or wave normals. However, in the ``forbidden'' region the WKB wave function is in general specified by two sets of wave fronts, the equiphase and the equiamplitude surfaces or equivalently by two sets of paths defined to be normal to these surfaces, respectively. We present a Huygens-type construction for obtaining these wave fronts and paths, which reveals that for non-normal incidence the paths are coupled to each other. The analysis enables us to answer some of the basic questions concerning tunneling in multidimensional nonseparable potentials. A special and important case occurs when the incident wave is normal to the turning surface. We show that for normal incidence the path equations are decoupled and are equivalent to Newton's equations of motion for the inverted potential and energy.

Physical Review A•Jan 1, 1990•80 citations

Multidimensional wave functions from classical trajectories

Michael J. Davis, Eric J. Heller

A new technique is developed to generate semiclassical wave functions. The method uses only information already available from a standard semiclassical quantization of a system. Linear superpositions of Gaussian coherent states that lie along quantizing classical trajectories are used, with phases given by the action integrals plus a Maslov-type correction. Wave functions generated in this way suffer from none of the problems with caustics that primitive semiclassical wave functions encounter. The semiclassical wave functions are convenient for subsequent use in applications, e.g., molecular spectra. By generating wave functions for several simple systems, we show that under most circumstances these wave functions are very accurate approximations to the true quantum states.

The Journal of Chemical Physics•Oct 15, 1981•98 citations
A method that uses only real valued coordinates is presented for integrating the many dimensional semiclassical wave function into the forbidden region. The procedure first determines a surface of caustic points by running the set of trajectories that define the wave function in the allowed region. In the forbidden region, the momentum and the action integral are both complex functions of position, and their imaginary parts vanish on the caustic surface. The direction of the imaginary part of the momentum p(I) can be chosen to the perpendicular to the caustic surface at all points on that surface. Equations are derived for integrating the values of the real and imaginary parts of the momentum along the curves that follow the direction of p(I). The equations for the change in the action integral and the prefactor for the semiclassical wave function along these curves are also obtained, allowing for the determination of the semiclassical wave function in the forbidden region. Calculations are performed for a two dimensional problem, and the semiclassical wave function is found to be is excellent agreement with the results of exact quantum 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.
We develop a treatment that divides phase space into regions bounded by surfaces and forms the energy-dependent semiclassical Green's function from short classical trajectory segments confined to the individual regions. We use this to generalize Bogomolny's quantization method to systems with barrier tunneling. Combined with an appropriate treatment of tunneling amplitudes, this yields a unified semiclassical framework that uses ordinary, real trajectories in classically allowed regions and imaginary-time trajectories in forbidden regions to obtain energies, decay rates, and wave functions for arbitrary states of multidimensional systems undergoing tunneling. Because the trajectory segments in the allowed region are short, components in the Green's function expression can be numerically evaluated using semiclassical initial value methods, yielding a promising approach for practical calculations.
We apply the quasilinearization method of Bellman and Kalaba [Quasilinearization and Nonlinear Boundary-Value Problems (Elsevier, New York, 1965)] to find approximate solutions for the multidimensional quantum tunneling for separable as well as nonseparable wave equations. By introducing the idea of the complex ``semiclassical trajectory'' which is valid for the motion over and under the barrier, and which, in the proper limit, reduces to the real classical trajectory in the allowed region, we obtain an eigenvalue equation for the characteristic wave numbers. This eigenvalue equation is similar to the corresponding equation obtained from the WKB approximation and yields complex eigenvalues with negative imaginary parts. When the barrier changes very rapidly as a function of the radial distance, we can replace the concept of the semiclassical trajectory, which may not be applicable in this case, by the concept of a complex ``quantum trajectory.'' The trajectory defined either way depends on a constant of integration, and by minimizing the action with respect to this constant we can obtain the minimum escape path. The case of two-dimensional tunneling is discussed as an example of this method.
Initial-value phase-space integral representations for a time-dependent propagator are obtained in the coordinate and momentum representations. To do so we first derive nonuniform semiclassical propagators for the various representations, obtaining the global-time semiclassical phase indices (Maslov indices) due to caustics. Results include readily implementable general phase index formulas for any type of caustic, including cases where the Morse index theorem is inapplicable. The method of obtaining the indices is general and based simply on concatenating uniform short-time propagators which also gives rise to alternative path-integral forms. Initial-value integral representations are then derived by introducing a method of extending short-time initial-value propagator formulas to global times via a simple stationary-phase asympotoic-equivalence approach. The integrals reduce to the nonuniform semiclassical propagators within the stationary-phase approximation, are uniform about caustics, and have integrand phases which properly account for the global-time phases in terms of appropriate Maslov indices. The initial-value integrals are also consistently derived via a canonical mapping procedure on the coordinate-space path integral. Initial-value integrals for time-dependent wave-function propgation are also given. Evaluation of the initial-value integral expressions do not require trajectory root searches for propagation.
The WKB method is commonly used in semiclassical approximations to the wave function in both the classically allowed and the forbidden regions of a one-dimensional potential. In a multidimensional space, the method can be adapted to construct wave functions in an ``allowed'' region from classical trajectories or wave normals. However, in the ``forbidden'' region the WKB wave function is in general specified by two sets of wave fronts, the equiphase and the equiamplitude surfaces or equivalently by two sets of paths defined to be normal to these surfaces, respectively. We present a Huygens-type construction for obtaining these wave fronts and paths, which reveals that for non-normal incidence the paths are coupled to each other. The analysis enables us to answer some of the basic questions concerning tunneling in multidimensional nonseparable potentials. A special and important case occurs when the incident wave is normal to the turning surface. We show that for normal incidence the path equations are decoupled and are equivalent to Newton's equations of motion for the inverted potential and energy.
A new technique is developed to generate semiclassical wave functions. The method uses only information already available from a standard semiclassical quantization of a system. Linear superpositions of Gaussian coherent states that lie along quantizing classical trajectories are used, with phases given by the action integrals plus a Maslov-type correction. Wave functions generated in this way suffer from none of the problems with caustics that primitive semiclassical wave functions encounter. The semiclassical wave functions are convenient for subsequent use in applications, e.g., molecular spectra. By generating wave functions for several simple systems, we show that under most circumstances these wave functions are very accurate approximations to the true quantum states.