© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

The Journal of Chemical Physics

Further analysis of solutions to the time-independent wave packet equations for quantum dynamics: General initial wave packets

David K. Hoffman, Youhong Huang, Wei Kang Zhu, Donald Jack Kouri

Published Jul 15, 1994• 30 citations

Abstract

In this paper we reexamine, and analyze solutions to, the recently derived time-independent wave packet-Schrödinger (TIWS) and time-independent wave packet-Lippmann–Schwinger (TIWLS) equations. These equations are so named because they are inhomogeneous, with the inhomogeneity being the initial L2 wave packet from an underlying time-dependent treatment of the dynamics. We explicitly show that a particular solution of the homogeneous Schrödinger equation can be constructed out of two particular solutions of the inhomogeneous TIW equation satisfying causal and anticausal boundary conditions. The structure of this solution of the homogeneous equation is shown to depend sensitively on the nature of the initial wave packet inhomogeneity, but, as we demonstrate, correct scattering information can be obtained even when the initial wave packet is nonzero only in the target region. It thus becomes possible to carry out quantum scattering calculations in which one need not propagate the wave packet from the noninteracting to the interacting region. The method is illustrated by calculations for two 1-D scattering problems, namely the transmission of an electron through a single barrier and through a double barrier. The latter is especially challenging because of the occurrence of long-lived resonances with the electron trapped inside the double barrier. In addition, we show that the method can still make use of absorbing potentials to decrease grid size, and we show how the formalism can be used to treat inelastic and reactive scattering, and radial scattering variables, thereby constituting a general approach to time-independent wave packet quantum scattering.

DOI: 10.1063/1.468439

Related papers

Time-dependent wave-packet forms of Schrödinger and Lippmann-Schwinger equations

Wei Kang Zhu, Youhong Huang, Donald Jack Kouri +2 more

Time-independent wave-packet forms of the Schr\"odinger equation (TIWSE) and Lippmann-Schwinger equation (TIWLSB) have been derived by a partial time-to-energy Fourier transform of ${\mathit{L}}^{2}$ wave-packet solutions to the time-dependent Schr\"odinger equation. The new equations retain the initial wave packet \ensuremath{\chi}(${\mathit{t}}_{0}$) as a ``universal source'' of scattered waves, which applies for all collision energies E contained in the initial wave packet. The relationship between the solution ${\mathrm{\ensuremath{\Psi}}}_{\mathit{t}0}$(E) of the TIWSE or TIWLSE and the scattering solution ${\mathrm{\ensuremath{\Psi}}}^{(+)}$(E) of the standard time-independent Lippmann-Schwinger or Schr\"odinger equation is given and the method illustrated by a computational application.

Physical Review Letters•Feb 28, 1994•41 citations

Scattered-wave-packet formalism with applications to barrier scattering and quantum transistors

Chia‐Chun Chou, Róbert E. Wyatt

The scattered wave formalism developed for a quantum subsystem interacting with reservoirs through open boundaries is applied to one- or two-dimensional barrier scattering and quantum transistors. The total wave function is divided into incident and scattered components. Markovian outgoing wave boundary conditions are imposed on the scattered or total wave function by either the ratio or polynomial methods. For barrier scattering problems, accurate time-dependent transmission probabilities are obtained through the integration of the modified time-dependent Schrödinger equations for the scattered wave function. For quantum transistors, the time-dependent transport is studied for a quantum wave packet propagating through the conduction channel of a field effect transistor. This study shows that the scattered wave formalism significantly reduces computational effort relative to other open boundary methods and demonstrates wide applications to quantum dynamical processes.

Physical Review E•Nov 23, 2011•3 citations

Two-particle Schrödinger equation animations of wave packet–wave packet scattering

Jon J. V. Maestri, Rubin H. Landau, Manuel José Páez

A simple and explicit technique for the numerical solution of the two-particle, time-dependent Schrödinger equation is assembled and tested. The technique can handle interparticle potentials that are arbitrary functions of the coordinates of each particle, arbitrary initial and boundary conditions, and multidimensional equations. Plots and animations are given here and on the World Wide Web of the scattering of two wave packet in one dimension.

arXiv:physics/9909042•American Journal of Physics•Dec 1, 2000•19 citations

Time-dependent treatment of scattering. II. Novel integral equation approach to quantum wave packets

Omar A. Sharafeddin, Richard Judson, Donald Jack Kouri +1 more

The time-dependent form of the Lippmann–Schwinger integral equation is used as the basis for a novel wave-packet propagation scheme. The method has the advantage over a previous integral equation treatment in that it does not require extensive matrix inversions involving the potential. This feature will be important when applications are made to systems where in some degrees of freedom the potential is expressed in a basis expansion. As was the case for the previous treatment, noniterated and iterated versions of the equations are given; the iterated equations, which are much simpler in the present new scheme than in the old, eliminate a matrix inversion that is required for solving the earlier noniterated equations. In the present noniterated equations, the matrix to be inverted is a function of the kinetic energy operator and thus is diagonal in a Bessel function basis set (or a sine basis set, if the centrifugal potential operator is incorporated into an effective potential). Transition amplitudes for various orbital angular momentum quantum numbers can be obtained from: (1) Fourier transform of the amplitude density from the time to the energy domain, and (2) direct analysis of the scattered wave packet. The approach is illustrated by an application to a standard potential scattering model problem.

The Journal of Chemical Physics•Oct 15, 1990•21 citations

Classical Hamiltonian structures in wave packet dynamics

Stephen K. Gray, John M. Verosky

The general, N state matrix representation of the time-dependent Schrödinger equation is equivalent to an N degree of freedom classical Hamiltonian system. We describe how classical mechanical methods and ideas can be applied towards understanding and modeling exact quantum dynamics. Two applications are presented. First, we illustrate how qualitative insights may be gained by treating the two state problem with a time-dependent coupling. In the case of periodic coupling, Poincaré surfaces of section are used to view the quantum dynamics, and features such as the Floquet modes take on interesting interpretations. The second application illustrates computational implications by showing how Liouville’s theorem, or more generally the symplectic nature of classical Hamiltonian dynamics, provides a new perspective for carrying out numerical wave packet propagation. We show how certain simple and explicit symplectic integrators can be used to numerically propagate wave packets. The approach is illustrated with an application to the problem of a diatomic molecule interacting with a laser, although it and related approaches may be useful for describing a variety of problems.

The Journal of Chemical Physics•Apr 1, 1994•78 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

In this paper we reexamine, and analyze solutions to, the recently derived time-independent wave packet-Schrödinger (TIWS) and time-independent wave packet-Lippmann–Schwinger (TIWLS) equations. These equations are so named because they are inhomogeneous, with the inhomogeneity being the initial L2 wave packet from an underlying time-dependent treatment of the dynamics. We explicitly show that a particular solution of the homogeneous Schrödinger equation can be constructed out of two particular solutions of the inhomogeneous TIW equation satisfying causal and anticausal boundary conditions. The structure of this solution of the homogeneous equation is shown to depend sensitively on the nature of the initial wave packet inhomogeneity, but, as we demonstrate, correct scattering information can be obtained even when the initial wave packet is nonzero only in the target region. It thus becomes possible to carry out quantum scattering calculations in which one need not propagate the wave packet from the noninteracting to the interacting region. The method is illustrated by calculations for two 1-D scattering problems, namely the transmission of an electron through a single barrier and through a double barrier. The latter is especially challenging because of the occurrence of long-lived resonances with the electron trapped inside the double barrier. In addition, we show that the method can still make use of absorbing potentials to decrease grid size, and we show how the formalism can be used to treat inelastic and reactive scattering, and radial scattering variables, thereby constituting a general approach to time-independent wave packet quantum scattering.

Time-independent wave-packet forms of the Schr\"odinger equation (TIWSE) and Lippmann-Schwinger equation (TIWLSB) have been derived by a partial time-to-energy Fourier transform of

{\mathit{L}}^{2}

wave-packet solutions to the time-dependent Schr\"odinger equation. The new equations retain the initial wave packet \ensuremath{\chi}(

{\mathit{t}}_{0}

) as a ``universal source'' of scattered waves, which applies for all collision energies E contained in the initial wave packet. The relationship between the solution ${\mathrm{\ensuremath{\Psi}}}_{\mathit{t}0}$(E) of the TIWSE or TIWLSE and the scattering solution ${\mathrm{\ensuremath{\Psi}}}^{(+)}$(E) of the standard time-independent Lippmann-Schwinger or Schr\"odinger equation is given and the method illustrated by a computational application.

The scattered wave formalism developed for a quantum subsystem interacting with reservoirs through open boundaries is applied to one- or two-dimensional barrier scattering and quantum transistors. The total wave function is divided into incident and scattered components. Markovian outgoing wave boundary conditions are imposed on the scattered or total wave function by either the ratio or polynomial methods. For barrier scattering problems, accurate time-dependent transmission probabilities are obtained through the integration of the modified time-dependent Schrödinger equations for the scattered wave function. For quantum transistors, the time-dependent transport is studied for a quantum wave packet propagating through the conduction channel of a field effect transistor. This study shows that the scattered wave formalism significantly reduces computational effort relative to other open boundary methods and demonstrates wide applications to quantum dynamical processes.

The time-dependent form of the Lippmann–Schwinger integral equation is used as the basis for a novel wave-packet propagation scheme. The method has the advantage over a previous integral equation treatment in that it does not require extensive matrix inversions involving the potential. This feature will be important when applications are made to systems where in some degrees of freedom the potential is expressed in a basis expansion. As was the case for the previous treatment, noniterated and iterated versions of the equations are given; the iterated equations, which are much simpler in the present new scheme than in the old, eliminate a matrix inversion that is required for solving the earlier noniterated equations. In the present noniterated equations, the matrix to be inverted is a function of the kinetic energy operator and thus is diagonal in a Bessel function basis set (or a sine basis set, if the centrifugal potential operator is incorporated into an effective potential). Transition amplitudes for various orbital angular momentum quantum numbers can be obtained from: (1) Fourier transform of the amplitude density from the time to the energy domain, and (2) direct analysis of the scattered wave packet. The approach is illustrated by an application to a standard potential scattering model problem.

The general, N state matrix representation of the time-dependent Schrödinger equation is equivalent to an N degree of freedom classical Hamiltonian system. We describe how classical mechanical methods and ideas can be applied towards understanding and modeling exact quantum dynamics. Two applications are presented. First, we illustrate how qualitative insights may be gained by treating the two state problem with a time-dependent coupling. In the case of periodic coupling, Poincaré surfaces of section are used to view the quantum dynamics, and features such as the Floquet modes take on interesting interpretations. The second application illustrates computational implications by showing how Liouville’s theorem, or more generally the symplectic nature of classical Hamiltonian dynamics, provides a new perspective for carrying out numerical wave packet propagation. We show how certain simple and explicit symplectic integrators can be used to numerically propagate wave packets. The approach is illustrated with an application to the problem of a diatomic molecule interacting with a laser, although it and related approaches may be useful for describing a variety of problems.

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.