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

Transactions of Beijing Institute of Technology

Numerical Solution of Quantum Tunneling

Qian‐Shu Li

Published Jan 1, 2010• 0 citations
OpenAlex

Abstract

The quantum hydrodynamics dynamics (QHD) is introduced to study the probability of electron through the tunneling junction. The complex form of the wave function was transformed into polar form, the time evolution equations of C and S were given in the Lagrange framework, and the explicit central differencing scheme was adopted for time integration of the QHD equations. In the Lagrangian framework, the grid points of the system will inevitable become scattered and unstructured during the evolution, the spatially smooth functions C and S are especially suitable for rational function approximation method to enable the implementation of a robust numerical scheme. As an example, the dynamics of a Gaussian wave packet is calculated and results are presented in the figures. In these figures, wave packet penetration of the barrier can be clearly observed, showing the quantum effect obviously.

Related papers

The Quantum Hydrodynamic Description of Tunneling

Brian K. Kendrick

The quantum hydrodynamic approach is based on the de Broglie-Bohm formulation of quantum mechanics. The resulting fluid-like equations of motion describe the flow of probability and an accurate solution to these equations is equivalent to solving the time-dependent Schroedinger equation. Furthermore, the hydrodynamic approach provides new insight into the mechanisms as well as an alternative computational approach for treating tunneling phenomena. New concepts include well-defined 'quantum trajectories', 'quantum potential', and 'quantum force' all of which have classical analogues. The quantum potential and its associated force give rise to all quantum mechanical effects such as zero point energy, tunneling, and interference. A new numerical approach called the Iterative Finite Difference Method (IFDM) will be discussed. The IFDM is used to solve the set of non-linear coupled hydrodynamic equations. It is 2nd-order accurate in both space and time and exhibits exponential convergence with respect to the iteration count. The stability and computational efficiency of the IFDM is significantly improved by using a 'smart' Eulerian grid which has the same computational advantages as a Lagrangian or Arbitrary Lagrangian Eulerian (ALE) grid. The IFDM is also capable of treating anharmonic potentials. Example calculations using the IFDM will be presented which include: a one-dimensional Gaussian wave packet tunneling through an Eckart barrier, a one-dimensional bound-state Morse oscillator, and a two-dimensional (2D) model collinear reaction using an anharmonic potential energy surface. Approximate treatments of the quantum hydrodynamic equations will also be discussed which could allow scaling of the calculations to hundreds of degrees of freedom which is important for treating tunneling phenomena in condensed phase systems.

•Jun 15, 2012•0 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

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

Courtney L. Lopreore, Róbert E. Wyatt

The de Broglie--Bohm causual (hydrodynamic) formulation of quantum mechanics is computationally implemented in the Lagrangian (moving with the fluid) viewpoint. The quantum potential and force are accurately evaluated with a moving weighted least squares algorithm. The quantum trajectory method is then applied to barrier tunneling on smooth potential surfaces. Analysis of the tunneling mechanism leads to a novel and accurate approximation: shortly after the wave packet is launched, completely neglect all quantum terms in the dynamical equations for motion along the tunneling coordinate.

Physical Review Letters•Jun 28, 1999•441 citations

Quantum tomography as a new approach to simulating quantum processes

A. S. Arkhipov, Yu. E. Lozovik

A new method is proposed for ab initio calculations of nonstationary quantum processes on the basis of a probability representation of quantum mechanics with the help of a positive definite function (quantum tomogram). The essence of the method is that an ensemble of trajectories associated with the characteristics of the evolution equation for the quantum tomogram is considered in the space where the quantum tomogram is defined. The method is applied for detailed analysis of transient tunneling of a wave packet. The results are in good agreement with the exact numerical solution to the Schrödinger equation for this system. The probability density distributions are obtained in the coordinate and momentum spaces at consecutive instances. For transient tunneling of a wave packet, the probability of penetration behind the barrier and the time of tunneling are calculated as functions of the initial energy.

Journal of Experimental and Theoretical Physics•Feb 1, 2004•7 citations

Quantum tunneling dynamics using hydrodynamic trajectories

Eric R. Bittner

In this paper we compute quantum trajectories arising from Bohm’s causal description of quantum mechanics. Our computational methodology is based upon a finite-element moving least-squares method (MWLS) presented recently by Wyatt and co-workers [Lopreore and Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. This method treats the “particles” in the quantum Hamilton–Jacobi equation as Lagrangian fluid elements that carry the phase, S, and density, ρ, required to reconstruct the quantum wave function. Here, we compare results obtained via the MWLS procedure to exact results obtained either analytically or by numerical solution of the time-dependent Schrödinger equation. Two systems are considered: first, dynamics in a harmonic well and second, tunneling dynamics in a double well potential. In the case of tunneling in the double well potential, the quantum potential acts to lower the barrier, separating the right- and left-hand sides of the well, permitting trajectories to pass from one side to another. However, as probability density passes from one side to the other, the effective barrier begins to rise and eventually will segregate trajectories in one side from the other. We note that the MWLS trajectories exhibited long time stability in the purely harmonic cases. However, this stability was not evident in the barrier crossing dynamics. Comparisons to exact trajectories obtained via wave packet calculations indicate that the MWLS trajectories tend to underestimate the effects of constructive and destructive interference effects.

arXiv:quant-ph/0001119•The Journal of Chemical Physics•Jun 8, 2000•92 citations
The quantum hydrodynamics dynamics (QHD) is introduced to study the probability of electron through the tunneling junction. The complex form of the wave function was transformed into polar form, the time evolution equations of C and S were given in the Lagrange framework, and the explicit central differencing scheme was adopted for time integration of the QHD equations. In the Lagrangian framework, the grid points of the system will inevitable become scattered and unstructured during the evolution, the spatially smooth functions C and S are especially suitable for rational function approximation method to enable the implementation of a robust numerical scheme. As an example, the dynamics of a Gaussian wave packet is calculated and results are presented in the figures. In these figures, wave packet penetration of the barrier can be clearly observed, showing the quantum effect obviously.
The quantum hydrodynamic approach is based on the de Broglie-Bohm formulation of quantum mechanics. The resulting fluid-like equations of motion describe the flow of probability and an accurate solution to these equations is equivalent to solving the time-dependent Schroedinger equation. Furthermore, the hydrodynamic approach provides new insight into the mechanisms as well as an alternative computational approach for treating tunneling phenomena. New concepts include well-defined 'quantum trajectories', 'quantum potential', and 'quantum force' all of which have classical analogues. The quantum potential and its associated force give rise to all quantum mechanical effects such as zero point energy, tunneling, and interference. A new numerical approach called the Iterative Finite Difference Method (IFDM) will be discussed. The IFDM is used to solve the set of non-linear coupled hydrodynamic equations. It is 2nd-order accurate in both space and time and exhibits exponential convergence with respect to the iteration count. The stability and computational efficiency of the IFDM is significantly improved by using a 'smart' Eulerian grid which has the same computational advantages as a Lagrangian or Arbitrary Lagrangian Eulerian (ALE) grid. The IFDM is also capable of treating anharmonic potentials. Example calculations using the IFDM will be presented which include: a one-dimensional Gaussian wave packet tunneling through an Eckart barrier, a one-dimensional bound-state Morse oscillator, and a two-dimensional (2D) model collinear reaction using an anharmonic potential energy surface. Approximate treatments of the quantum hydrodynamic equations will also be discussed which could allow scaling of the calculations to hundreds of degrees of freedom which is important for treating tunneling phenomena in condensed phase systems.
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 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 de Broglie--Bohm causual (hydrodynamic) formulation of quantum mechanics is computationally implemented in the Lagrangian (moving with the fluid) viewpoint. The quantum potential and force are accurately evaluated with a moving weighted least squares algorithm. The quantum trajectory method is then applied to barrier tunneling on smooth potential surfaces. Analysis of the tunneling mechanism leads to a novel and accurate approximation: shortly after the wave packet is launched, completely neglect all quantum terms in the dynamical equations for motion along the tunneling coordinate.
A new method is proposed for ab initio calculations of nonstationary quantum processes on the basis of a probability representation of quantum mechanics with the help of a positive definite function (quantum tomogram). The essence of the method is that an ensemble of trajectories associated with the characteristics of the evolution equation for the quantum tomogram is considered in the space where the quantum tomogram is defined. The method is applied for detailed analysis of transient tunneling of a wave packet. The results are in good agreement with the exact numerical solution to the Schrödinger equation for this system. The probability density distributions are obtained in the coordinate and momentum spaces at consecutive instances. For transient tunneling of a wave packet, the probability of penetration behind the barrier and the time of tunneling are calculated as functions of the initial energy.
In this paper we compute quantum trajectories arising from Bohm’s causal description of quantum mechanics. Our computational methodology is based upon a finite-element moving least-squares method (MWLS) presented recently by Wyatt and co-workers [Lopreore and Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. This method treats the “particles” in the quantum Hamilton–Jacobi equation as Lagrangian fluid elements that carry the phase, S, and density, ρ, required to reconstruct the quantum wave function. Here, we compare results obtained via the MWLS procedure to exact results obtained either analytically or by numerical solution of the time-dependent Schrödinger equation. Two systems are considered: first, dynamics in a harmonic well and second, tunneling dynamics in a double well potential. In the case of tunneling in the double well potential, the quantum potential acts to lower the barrier, separating the right- and left-hand sides of the well, permitting trajectories to pass from one side to another. However, as probability density passes from one side to the other, the effective barrier begins to rise and eventually will segregate trajectories in one side from the other. We note that the MWLS trajectories exhibited long time stability in the purely harmonic cases. However, this stability was not evident in the barrier crossing dynamics. Comparisons to exact trajectories obtained via wave packet calculations indicate that the MWLS trajectories tend to underestimate the effects of constructive and destructive interference effects.