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

arXiv:1601.03827

Time Evolution of Gaussian Wave Packets under Dirac Equation with Fluctuating Mass and Potential

Atis Yosprakob, Sujin Suwanna

Published Jan 15, 2016• 1 citation
arXivHTMLPDF

Abstract

Localization of relativistic particles have been of great research interests over many decades. We investigate the time evolution of the Gaussian wave packets governed by the one dimensional Dirac equation. For the free Dirac equation, we obtain the evolution profiles analytically in many approximation regimes, and numerical simulations consistent with other numerical schemes. Interesting behaviors such as Zitterbewegung and Klein paradox are exhibited. In particular, the dispersion rate as a function of mass is calculated, and it yields an interesting result that super-massive and massless particles both exhibit no dispersion in free space. For the Dirac equation with random potential or mass, we employ the Chebyshev polynomials expansion of the propagator operator to numerically investigate the probability profiles of the displacement distribution when the potential or mass is uniformly distributed. We observe that the widths of the Gaussian wave packets decrease approximately with the power law of order $o(s^{-ν})$ with $\frac{1}{2}

DOI
10.48550/arxiv.1601.03827

Related papers

Space-time evolution of Dirac wave packets

V. Ya. Demikhovskiǐ, Galina M. Maksimova, А. А. Перов +1 more

In this work we study the dynamics of free three-dimensional relativistic Gaussian wave packets with different spin polarizations. We analyze the connection between the symmetry of initial state and the dynamical characteristics of moving particle. The corresponding solutions of the Dirac equation having different types of symmetry were evaluated analytically and numerically and after that the electron probability densities, as well as the spin densities were visualized. The average values of velocity of the packet center and the average spin were calculated analytically, and the parameters of transient Zitterbewegung in different directions were obtained. These results can be useful for the interpretation of future experiments with trapped ions.

arXiv:1007.1566•Physical Review A•Nov 24, 2010•21 citations

Asymptotic Momentum of Dirac Particles in One Space Dimension

Kabir Narayanan, Abigail Perryman, A. Shadi Tahvildar-Zadeh

We analyze the trajectories of a massive particle in one space dimension whose motion is guided by a spin-half wave function that evolves according to the free Dirac equation, with its initial wave function being a Gaussian wave packet with a nonzero expected value of momentum $k$. We prove that at large times, the wave function is approximately equal to the superposition of two wave packets traveling in opposite directions, which results in trajectories with approximately constant asymptotic momentum $k$ and asymptotic energy $\pm c^2\sqrt{m^2+k^2}$, with $m$ the rest mass of the particle and $c$ the speed of light. The sign of the asymptotic energy is determined by the initial position of the particle. Particles with negative energy will have an asymptotic velocity that is in the opposite direction of their momentum. The proof uses the stationary phase approximation method, for which we establish a rigorous error bound.

math-phCross-listedarXiv:2512.21423•Dec 26, 2025•0 citations

Investigation of an Equation with Logarithmic Coupling with Momentum in Dirac Equation

Mash-Huda Rahman Shipra, Syed Badiuzzaman Faruque

The free particle solutions of the relativistic Dirac equation are characterized by plane waves with infinite uncertainty in position. However, many practical implementations of the solutions require a wave packet configuration, which can be utilized to represent a localized Dirac particle. Unlike the traditional wave packet generation method by superposing multiple plane waves, this study of ours presents an alternative approach towards obtaining a wave packet solution of a free particle relativistic Dirac equation. In this paper, we present Dirac’s free particle equation with a modification in the generalized momentum. The modification is achieved by coupling the momentum with a spatially varying logarithmic function, and this alteration does not affect the relativistic dispersion relation of the particle. Moreover, a solution of this modified Dirac equation is provided as well, which has been calculated using a trial wave function. The wave function solution is carried out in one dimension, where it behaves as a wave packet for a given ratio of the envelope parameter to the reduced Planck's constant greater than unity, where the envelope parameter regulates the width of the wave packet. The solution, being subject to this constraint, represents a bound particle with spin and a continuous energy spectrum.

International Journal of High Energy Physics•Jan 1, 2020•0 citations

Approximations for the free evolution of self-gravitating quantum particles

André Großardt

The evolution of the center-of-mass wave function for a mesoscopic particle according to the Schr\"odinger-Newton equation can be approximated by a harmonic potential if the wave function is narrow compared to the size of the mesoscopic particle. It was noticed by Colin et al. [Phys. Rev. A 93, 062102 (2016).] that, in the regime where self-gravitational effects are weak, intermediate and wider wave functions may be approximated by a harmonic potential as well but with a width-dependent coupling, leading to a time evolution that is determined only by a differential equation for the width of a Gaussian wave function as a single parameter. Such an approximation results in considerably less computational effort in order to predict the self-gravitational effects on the wave-function dynamics. Here, we provide an alternative approach to this kind of approximation, including a rigorous derivation of the equations of motion for an initially Gaussian wave packet under the assumption that its shape is conserved. Our result deviates to some degree from the result by Colin et al. [Phys. Rev. A 93, 062102 (2016).], specifically in the limit of wide wave functions.

arXiv:1503.02622•Physical Review A•Aug 1, 2016•12 citations

Quantum diffusion of Dirac fermions in a tight-binding model

M. Cruz-Méndez, Heber R. da Cruz

In this work, we develop a time-dependent mathematical model to investigate the wave packet evolution of both massive and massless Dirac fermions in a tight-binding lattice. Our study focuses on the quantum diffusion of carriers in a one-dimensional crystal, where the system Hamiltonian is represented by both the Dirac and Schrödinger equations. Interestingly, during the quantum diffusion process, we observe oscillations of the wave packet between the spinor-up and spinor-down bands in the Dirac case. In the massive fermion scenario, these oscillations between the sublattices are influenced by inertial mass effects, resulting in longer oscillation periods as the mass increases. In addition, we find a strong dependence of Dirac fermion dynamics on the Fermi velocity in the tight-binding chain.

APL Quantum•Jun 1, 2025•0 citations

Optical simulation of the free Dirac equation

Thais de Lima Silva, E. R. F. Taillebois, Rafael M. Gomes +2 more

We present a classical optics simulation of the one-dimensional Dirac equation for a free particle. Positive- and negative-energy components are represented by orthogonal polarizations of a free propagating beam, while the spatial profile represents the spatial wave function of the particle. Using a completely tunable time parameter, we observed the oscillation of the average value of the Dirac position operator---known as Zitterbewegung. We are also able to measure the so-called mean-position operator which is a single-particle observable and presents no oscillations. Our work opens the way for optical simulation of interesting phenomenon of relativistic systems, as well as condensed-matter physics, without any requirement for specially engineered medium.

arXiv:1902.04012•Physical Review A•Feb 25, 2019•14 citations
Localization of relativistic particles have been of great research interests over many decades. We investigate the time evolution of the Gaussian wave packets governed by the one dimensional Dirac equation. For the free Dirac equation, we obtain the evolution profiles analytically in many approximation regimes, and numerical simulations consistent with other numerical schemes. Interesting behaviors such as Zitterbewegung and Klein paradox are exhibited. In particular, the dispersion rate as a function of mass is calculated, and it yields an interesting result that super-massive and massless particles both exhibit no dispersion in free space. For the Dirac equation with random potential or mass, we employ the Chebyshev polynomials expansion of the propagator operator to numerically investigate the probability profiles of the displacement distribution when the potential or mass is uniformly distributed. We observe that the widths of the Gaussian wave packets decrease approximately with the power law of order o(s−ν)o(s^{-ν})o(s−ν) with $\frac{1}{2}
In this work we study the dynamics of free three-dimensional relativistic Gaussian wave packets with different spin polarizations. We analyze the connection between the symmetry of initial state and the dynamical characteristics of moving particle. The corresponding solutions of the Dirac equation having different types of symmetry were evaluated analytically and numerically and after that the electron probability densities, as well as the spin densities were visualized. The average values of velocity of the packet center and the average spin were calculated analytically, and the parameters of transient Zitterbewegung in different directions were obtained. These results can be useful for the interpretation of future experiments with trapped ions.
We analyze the trajectories of a massive particle in one space dimension whose motion is guided by a spin-half wave function that evolves according to the free Dirac equation, with its initial wave function being a Gaussian wave packet with a nonzero expected value of momentum kkk. We prove that at large times, the wave function is approximately equal to the superposition of two wave packets traveling in opposite directions, which results in trajectories with approximately constant asymptotic momentum kkk and asymptotic energy ±c2m2+k2\pm c^2\sqrt{m^2+k^2}±c2m2+k2​, with mmm the rest mass of the particle and ccc the speed of light. The sign of the asymptotic energy is determined by the initial position of the particle. Particles with negative energy will have an asymptotic velocity that is in the opposite direction of their momentum. The proof uses the stationary phase approximation method, for which we establish a rigorous error bound.
The free particle solutions of the relativistic Dirac equation are characterized by plane waves with infinite uncertainty in position. However, many practical implementations of the solutions require a wave packet configuration, which can be utilized to represent a localized Dirac particle. Unlike the traditional wave packet generation method by superposing multiple plane waves, this study of ours presents an alternative approach towards obtaining a wave packet solution of a free particle relativistic Dirac equation. In this paper, we present Dirac’s free particle equation with a modification in the generalized momentum. The modification is achieved by coupling the momentum with a spatially varying logarithmic function, and this alteration does not affect the relativistic dispersion relation of the particle. Moreover, a solution of this modified Dirac equation is provided as well, which has been calculated using a trial wave function. The wave function solution is carried out in one dimension, where it behaves as a wave packet for a given ratio of the envelope parameter to the reduced Planck's constant greater than unity, where the envelope parameter regulates the width of the wave packet. The solution, being subject to this constraint, represents a bound particle with spin and a continuous energy spectrum.
The evolution of the center-of-mass wave function for a mesoscopic particle according to the Schr\"odinger-Newton equation can be approximated by a harmonic potential if the wave function is narrow compared to the size of the mesoscopic particle. It was noticed by Colin et al. [Phys. Rev. A 93, 062102 (2016).] that, in the regime where self-gravitational effects are weak, intermediate and wider wave functions may be approximated by a harmonic potential as well but with a width-dependent coupling, leading to a time evolution that is determined only by a differential equation for the width of a Gaussian wave function as a single parameter. Such an approximation results in considerably less computational effort in order to predict the self-gravitational effects on the wave-function dynamics. Here, we provide an alternative approach to this kind of approximation, including a rigorous derivation of the equations of motion for an initially Gaussian wave packet under the assumption that its shape is conserved. Our result deviates to some degree from the result by Colin et al. [Phys. Rev. A 93, 062102 (2016).], specifically in the limit of wide wave functions.
In this work, we develop a time-dependent mathematical model to investigate the wave packet evolution of both massive and massless Dirac fermions in a tight-binding lattice. Our study focuses on the quantum diffusion of carriers in a one-dimensional crystal, where the system Hamiltonian is represented by both the Dirac and Schrödinger equations. Interestingly, during the quantum diffusion process, we observe oscillations of the wave packet between the spinor-up and spinor-down bands in the Dirac case. In the massive fermion scenario, these oscillations between the sublattices are influenced by inertial mass effects, resulting in longer oscillation periods as the mass increases. In addition, we find a strong dependence of Dirac fermion dynamics on the Fermi velocity in the tight-binding chain.
We present a classical optics simulation of the one-dimensional Dirac equation for a free particle. Positive- and negative-energy components are represented by orthogonal polarizations of a free propagating beam, while the spatial profile represents the spatial wave function of the particle. Using a completely tunable time parameter, we observed the oscillation of the average value of the Dirac position operator---known as Zitterbewegung. We are also able to measure the so-called mean-position operator which is a single-particle observable and presents no oscillations. Our work opens the way for optical simulation of interesting phenomenon of relativistic systems, as well as condensed-matter physics, without any requirement for specially engineered medium.