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

International Journal of Quantum Chemistry

Moving boundary truncated grid method: Multidimensional quantum dynamics

Chun‐Yaung Lu, Tsung‐Yen Lee, Chia‐Chun Chou

Published Nov 1, 2019• 4 citations
Publisher

Abstract

Abstract The moving boundary truncated grid (TG) method is used to study wave packet dynamics of multidimensional quantum systems. As time evolves, appropriate Eulerian grid points required for propagating a wave packet are activated and deactivated with no advance information about the dynamics. This method is applied to the Henon‐Heiles potential and wave packet barrier scattering in two, three, and four dimensions. Computational results demonstrate that the TG method not only leads to a great reduction in the number of grid points needed to perform accurate calculations but also is computationally more efficient than the full grid calculations.

DOI
10.1002/qua.26055

Related papers

Wavepacket Dynamics on Dynamically Adapting Grids: Application of the Equidistribution Principle

Keith H. Hughes, Róbert E. Wyatt

A moving grid approach to wavepacket dynamics is described that enables grid points to be used efficiently in regions where high resolution of the wavepacket is required. The grid movement is based on the principle of equidistribution and by using a grid smoothing technique the grid points trace a path that continuously adapt to reflect the dynamics of the wavepacket. The technique is robust and allows accurate computations to be obtained for long wavepacket propagation times. Results are presented for two systems: tunnelling dynamics in a double well potential and scattering of a wavepacket from a repulsive Eckart barrier.

APS•Mar 1, 2003•0 citations

Quantum wave packet dynamics with trajectories: Implementation with adaptive Lagrangian grids

Róbert E. Wyatt, Eric R. Bittner

The quantum trajectory method was recently developed to solve the hydrodynamic equations of motion in the Lagrangian, moving-with-the-fluid, picture. In this approach, trajectories are integrated for fluid elements (“particles”) moving under the influence of the combined force from the potential surface and the quantum potential. To accurately compute the quantum potential and the quantum force, it is necessary to obtain the derivatives of a function given only the values on the unstructured mesh defined by the particle locations. However, in some regions of space–time, the particle mesh shows compression and inflation associated with regions of large and small density, respectively. Inflation is especially severe near nodes in the wave function. In order to circumvent problems associated with highly nonuniform grids defined by the particle locations, adaptation of moving grids is introduced in this study. By changing the representation of the wave function in these local regions (which can be identified by diagnostic tools), propagation is possible to much longer times. These grid adaptation techniques are applied to the reflected portion of a wave packet scattering from an Eckart potential.

The Journal of Chemical Physics•Nov 22, 2000•122 citations

Quantum wave packet dynamics with trajectories: Application to reactive scattering

Róbert E. Wyatt

The quantum trajectory method (QTM) for time-dependent wave packet dynamics involves integration of the de Broglie–Bohm hydrodynamic equations for the evolving probability fluid [C. Lopreore and R. E. Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. The equations of motion for discretized elements of the probability fluid (particles) are integrated in the Lagrangian, moving with the fluid, picture. These fluid elements move under the influence of both the usual potential energy function and the quantum potential, which involves the curvature of the quantum amplitude. The quantum potential and the quantum force are evaluated using a moving weighted least squares algorithm. As a demonstration of applicability, the QTM is applied to a model collinear reaction with an activation barrier. The reaction probabilities are in good agreement with exact quantum results, even for a relatively small number of particles in the ensemble. The QTM accurately describes tunneling using only real valued trajectories. In addition to the reaction probability, plots are presented to show the probability density and the flux distributions at several time steps during the reaction.

The Journal of Chemical Physics•Sep 8, 1999•135 citations

Time‐dependent schrödinger equation with Markovian outgoing wave boundary conditions: Applications to quantum tunneling dynamics and photoionization

Chia‐Chun Chou, Róbert E. Wyatt

Abstract Markovian outgoing wave boundary conditions are introduced as an approximate method to reduce the size of the computational grid for time integration of the time‐dependent Schrödinger equation. The ratio and polynomial methods developed as open boundary conditions are applied to the wave function at the boundaries of the computational grid. This computational method is used to study the wave packet dynamics for a metastable well, a double well, and strong‐field ionization of a model atom. Accurate results demonstrate that this method can significantly reduce the number of grid points required in a dynamical calculation for quantum dynamical problems. © 2012 Wiley Periodicals, Inc.

International Journal of Quantum Chemistry•Feb 13, 2012•10 citations

A new method for solving the quantum hydrodynamic equations of motion

Brian K. Kendrick

The quantum hydrodynamic equations associated with the de Broglie–Bohm formulation of quantum mechanics are solved using a meshless method based on a moving least squares approach. An arbitrary Lagrangian–Eulerian frame of reference is used which significantly improves the accuracy and stability of the method when compared to an approach based on a purely Lagrangian frame of reference. A regridding algorithm is implemented which adds and deletes points when necessary in order to maintain accurate and stable calculations. It is shown that unitarity in the time evolution of the quantum wave packet is significantly improved by propagating using averaged fields. As nodes in the reflected wave packet start to form, the quantum potential and force become very large and numerical instabilities occur. By introducing artificial viscosity into the equations of motion, these instabilities can be avoided and the stable propagation of the wave packet for very long times becomes possible. Results are presented for the scattering of a wave packet from a repulsive Eckart barrier.

The Journal of Chemical Physics•Sep 8, 2003•100 citations

Computational Investigation of Wave Packet Scattering in the Complex Plane: Propagation on a Grid

Róbert E. Wyatt, Brad A. Rowland

The time-dependent scattering of a wave packet from a Gaussian barrier is investigated computationally in the complex z-plane. The initial wave packet and the potential energy are obtained through analytic continuation from functions specified on the real-axis. The wave packet is then propagated on the two-dimensional grid. For a low initial wave packet energy, the time evolution is followed by plotting the following functions: |ψ(z,t)|, real(ψ(z,t)), and the quantum momentum function (QMF), p(z,t). In the reflected packet, an important role is played by ripples (quasi-nodes) forming above the real axis. As these quasi-nodes move down across the real axis, they are 'detected' as 'interference oscillations' in the density. In contrast, the component of the packet below the real axis makes a significant contribution to the transmitted packet. Vector maps of the QMF show hyperbolic flow around quasi-nodes and counterclockwise circular flow around transient stagnation points, where the QMF vanishes. However, when the Pólya vector field (defined by P(z,t) = p*(z,t)) is plotted, circular counterclockwise flow is obtained near the quasi-nodes. The real and imaginary parts of the quantum action function S(z,t) are plotted and the vorticity, defined by the curl of the Pólya field, is used to pinpoint regions of nonanalyticity in the QMF.

Journal of Chemical Theory and Computation•Jan 29, 2009•12 citations
Abstract The moving boundary truncated grid (TG) method is used to study wave packet dynamics of multidimensional quantum systems. As time evolves, appropriate Eulerian grid points required for propagating a wave packet are activated and deactivated with no advance information about the dynamics. This method is applied to the Henon‐Heiles potential and wave packet barrier scattering in two, three, and four dimensions. Computational results demonstrate that the TG method not only leads to a great reduction in the number of grid points needed to perform accurate calculations but also is computationally more efficient than the full grid calculations.
A moving grid approach to wavepacket dynamics is described that enables grid points to be used efficiently in regions where high resolution of the wavepacket is required. The grid movement is based on the principle of equidistribution and by using a grid smoothing technique the grid points trace a path that continuously adapt to reflect the dynamics of the wavepacket. The technique is robust and allows accurate computations to be obtained for long wavepacket propagation times. Results are presented for two systems: tunnelling dynamics in a double well potential and scattering of a wavepacket from a repulsive Eckart barrier.
The quantum trajectory method was recently developed to solve the hydrodynamic equations of motion in the Lagrangian, moving-with-the-fluid, picture. In this approach, trajectories are integrated for fluid elements (“particles”) moving under the influence of the combined force from the potential surface and the quantum potential. To accurately compute the quantum potential and the quantum force, it is necessary to obtain the derivatives of a function given only the values on the unstructured mesh defined by the particle locations. However, in some regions of space–time, the particle mesh shows compression and inflation associated with regions of large and small density, respectively. Inflation is especially severe near nodes in the wave function. In order to circumvent problems associated with highly nonuniform grids defined by the particle locations, adaptation of moving grids is introduced in this study. By changing the representation of the wave function in these local regions (which can be identified by diagnostic tools), propagation is possible to much longer times. These grid adaptation techniques are applied to the reflected portion of a wave packet scattering from an Eckart potential.
The quantum trajectory method (QTM) for time-dependent wave packet dynamics involves integration of the de Broglie–Bohm hydrodynamic equations for the evolving probability fluid [C. Lopreore and R. E. Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. The equations of motion for discretized elements of the probability fluid (particles) are integrated in the Lagrangian, moving with the fluid, picture. These fluid elements move under the influence of both the usual potential energy function and the quantum potential, which involves the curvature of the quantum amplitude. The quantum potential and the quantum force are evaluated using a moving weighted least squares algorithm. As a demonstration of applicability, the QTM is applied to a model collinear reaction with an activation barrier. The reaction probabilities are in good agreement with exact quantum results, even for a relatively small number of particles in the ensemble. The QTM accurately describes tunneling using only real valued trajectories. In addition to the reaction probability, plots are presented to show the probability density and the flux distributions at several time steps during the reaction.
Abstract Markovian outgoing wave boundary conditions are introduced as an approximate method to reduce the size of the computational grid for time integration of the time‐dependent Schrödinger equation. The ratio and polynomial methods developed as open boundary conditions are applied to the wave function at the boundaries of the computational grid. This computational method is used to study the wave packet dynamics for a metastable well, a double well, and strong‐field ionization of a model atom. Accurate results demonstrate that this method can significantly reduce the number of grid points required in a dynamical calculation for quantum dynamical problems. © 2012 Wiley Periodicals, Inc.
The quantum hydrodynamic equations associated with the de Broglie–Bohm formulation of quantum mechanics are solved using a meshless method based on a moving least squares approach. An arbitrary Lagrangian–Eulerian frame of reference is used which significantly improves the accuracy and stability of the method when compared to an approach based on a purely Lagrangian frame of reference. A regridding algorithm is implemented which adds and deletes points when necessary in order to maintain accurate and stable calculations. It is shown that unitarity in the time evolution of the quantum wave packet is significantly improved by propagating using averaged fields. As nodes in the reflected wave packet start to form, the quantum potential and force become very large and numerical instabilities occur. By introducing artificial viscosity into the equations of motion, these instabilities can be avoided and the stable propagation of the wave packet for very long times becomes possible. Results are presented for the scattering of a wave packet from a repulsive Eckart barrier.
The time-dependent scattering of a wave packet from a Gaussian barrier is investigated computationally in the complex z-plane. The initial wave packet and the potential energy are obtained through analytic continuation from functions specified on the real-axis. The wave packet is then propagated on the two-dimensional grid. For a low initial wave packet energy, the time evolution is followed by plotting the following functions: |ψ(z,t)|, real(ψ(z,t)), and the quantum momentum function (QMF), p(z,t). In the reflected packet, an important role is played by ripples (quasi-nodes) forming above the real axis. As these quasi-nodes move down across the real axis, they are 'detected' as 'interference oscillations' in the density. In contrast, the component of the packet below the real axis makes a significant contribution to the transmitted packet. Vector maps of the QMF show hyperbolic flow around quasi-nodes and counterclockwise circular flow around transient stagnation points, where the QMF vanishes. However, when the Pólya vector field (defined by P(z,t) = p*(z,t)) is plotted, circular counterclockwise flow is obtained near the quasi-nodes. The real and imaginary parts of the quantum action function S(z,t) are plotted and the vorticity, defined by the curl of the Pólya field, is used to pinpoint regions of nonanalyticity in the QMF.