Journal of Physics A Mathematical and General

Numerical evaluation of coherent-state path integrals in quantum dynamics

Bernd Burghardt, Joachim Stolze

Published Jan 1, 1999• 4 citations

Abstract

The numerical evaluation of coherent-state path integrals for quantum dynamical problems is discussed for one-dimensional examples. To propagate an initial state, we use the normal and antinormal ordered coherent-state path integrals combined with a split-operator technique dividing the Hamiltonian into harmonic and anharmonic parts. For numerical purposes integrations must be approximated by quadrature formulae. This leads to a matrix multiplication scheme which is systematically tested for the double-well and Morse potentials. The method is accurate for propagation times much longer than the natural time scale of the system, and it allows for short as well as long time steps without loss of stability.

DOI: 10.1088/0305-4470/32/11/004

Related papers

Numerical evaluation of coherent-state path integrals with applications to time-dependent problems

Bernd Burghardt, Joachim Stolze

We study the application of the coherent-state path integral as a numerical tool for wave-packet propagation. The numerical evaluation of path integrals is reduced to a matrix-vector multiplication scheme. Together with a split-operator technique we apply our method to a time-dependent double-well potential.

arXiv:quant-ph/9901076Jan 28, 19990 citations

Numerical evaluation of coherent-state path integrals with application to time-dependent problems

Joachim Stolze, Bernd Burghardt

Jan 1, 19990 citations

Numerical path integration without Monte Carlo

Christopher C. Gerry, J. Kiefer

A numerical procedure for evaluating Feynman’s path integral in imaginary time without using the Monte Carlo ‘‘heat bath’’ method is discussed. The method is essentially iterative matrix multiplication and bears a close resemblance to the Feynman procedure of summing all paths by time slicing and then integrating. The procedure is used to determine ground-state energy eigenvalues for the harmonic oscillator, double-well anharmonic oscillator, and the Morse oscillator.

American Journal of PhysicsNov 1, 19888 citations

Semiclassical coherent-state path integrals for scattering

Frank Großmann

We study the application of the semiclassical coherent-state path-integral technique to the problem of scattering. A numerical method based on initial-value trajectories and their stability information is proposed to solve the root search for the complex trajectories that determine the semiclassical propagator. Very good agreement between the numerically exact quantum and semiclassical results for the correlation function can be achieved for scattering from one-dimensional barrier as well as well-type potentials by using just one single complex trajectory at each time.

Physical Review AMay 1, 199830 citations

Numerical path integral techniques for long time dynamics of quantum dissipative systems

Nancy Makri

Recent progress in numerical methods for evaluating the real-time path integral in dissipative harmonic environments is reviewed. Quasi-adiabatic propagators constructed numerically allow convergence of the path integral with large time increments. Integration of the harmonic bath leads to path integral expressions that incorporate the exact dynamics of the quantum particle along the adiabatic path, with an influence functional that describes nonadiabatic corrections. The resulting quasi-adiabatic propagator path integral is evaluated by efficient system-specific quadratures in most regimes of parameter space, although some cases are handled by grid Monte Carlo sampling. Exploiting the finite span of nonlocal influence functional interactions characteristic of broad condensed phase spectra leads to an iterative scheme for calculating the path integral over arbitrary time lengths. No uncontrolled approximations are introduced, and the resulting methodology converges to the exact quantum result with modest amounts of computational power. Applications to tunneling dynamics in the condensed phase are described.

Journal of Mathematical PhysicsMay 1, 1995377 citations

A Comparison of Different Numerical Propagation Schemes for Solving the Time-Dependent Schrödinger Equation in the Position Representation in One Dimension for Mixed Quantum-and Molecular Dynamics Simulations

Salomon R. Billeter, Wilfred F. Van Gunsteren A

Various numerical integration schemes to calculate the propagation of a state following the time-dependent Schrodinger equation in the one dimensional position representation are presented and compared to each other. Three potentials have been used: a harmonic, a double-well and a zero potential. Eigenstates and a coherent state have been chosen as initial states. Special attention has been given to the long-time stability of the algorithms. These are: kinetic referenced split operator (KRSO), kinetic referenced Cayley (KRC), distributed approximating functions (DAF), Chebysheff expansion (CH), residuum minimization (RES), second-order differencing (SOD), an eigenstate expansion (EE) and a corrected kinetic referenced split operator (CKRSO). In addition, a speedup of the KRC and KRSO methods is presented which is specially suited when very few grid points are used. Numerical results are compared to analytically calculated values. Mixed classical/quantum mechanical simulations require a representation of the quantum state on a limited number of grid points, classical integration time steps of about one femtosecond and compatibility with methods to solve the time-ordering problem. For the considered potentials which differ quite essentially from the potentials used for scattering problems in particle physics, the EE method has been found to be faster, more accurate and more stable than the other methods if only a few grid points are required. Otherwise, good results have been obtained with KRC, KRSO, CH, DAF and RES. SOD has been found to be too slow, and CKRSO is not stable enough for long simulation times.

Molecular SimulationNov 1, 199517 citations