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:2603.15395 · Universe

Quantum-classical diagnostics and Bohmian inequivalence for higher time-derivative Hamiltonians

Sanjib Dey, Andreas Fring

quant-phhep-thmath-phSubmitted Mar 16, 2026• Announced Mar 17, 2026• 0 citations
arXivHTMLPDF

Abstract

We develop a Bohmian analysis of a two-dimensional ghost Hamiltonian and its mapping to the degenerate Pais-Uhlenbeck model. Using Gaussian wavepackets, we derive the corresponding guidance equations, the centre and width evolution, and the quantum potential. We use these quantities to characterise bounded, quasi-semiclassical, spiral, and runaway regimes. The Bohmian trajectories provide a direct dynamical diagnostic of coherence, packet deformation, and quantum-classical separation. We then compare a bi-Hamiltonian pair consisting of the ghost Hamiltonian and a classically equivalent alternative formulation. While the two descriptions produce identical classical trajectories, they lead to different Bohmian trajectories and different quantum potentials evaluated along those trajectories. This demonstrates that classical equivalence need not extend to Bohmian quantum dynamics and identifies a concrete quantum ambiguity in the degenerate higher-derivative system.

Journal reference
Universe 2026, 12(7), 200
DOI
10.3390/universe12070200
Comments
17 pages, 6 figures

Related papers

Bohmian versus semiclassical description of interference phenomena

Yi Zhao, Nancy Makri

The origin of quantum interference characteristic of bound nonlinear systems is investigated within the Bohmian formulation of time-dependent quantum mechanics. By contrast to time-dependent semiclassical theory, whereby interference is a consequence of phase mismatch between distinct classical trajectories, the Bohmian, fully quantum mechanical expression for expectation values has a quasiclassical appearance that does not involve phase factors or cross terms. Numerical calculations reveal that quantum interference in the Bohmian formulation manifests itself directly as sharp spatial/temporal variations of the density surrounding kinky trajectories. These effects are most dramatic in regions where the underlying classical motion exhibits focal points or caustics, and crossing of the Bohmian trajectories is prevented through extremely strong and rapidly varying quantum mechanical forces. These features of Bohmian dynamics, which constitute the hallmark of quantum interference and are ubiquitous in bound nonlinear systems, represent a major source of instability, making the integration of the Bohmian equations extremely demanding in such situations.

The Journal of Chemical Physics•Jun 19, 2003•36 citations

Mixing quantum and classical dynamics using Bohmian trajectories

Etienne Gindensperger, C. Meier, Juan Alberto Beswick

A novel time-dependent hybrid quantum/classical propagation scheme based on Bohmian quantum trajectories is presented. The quantum subsystem is described by a wave packet depending on the quantum variables x and, via the total potential energy of the system, parametrically on the classical trajectories X(t). The wave packet is used to calculate de Broglie–Bohm quantum trajectories x(t) which are used to calculate the force acting on the classical variables. Quantum corrections of the classical equation of motion are also included. The method is applied to a simple case of two coupled oscillators. Comparison between exact quantum and approximate results demonstrates that these MQCB (Mixed Quantum/Classical Bohmian) trajectories provide a good description of the energy exchange between the two oscillators.

The Journal of Chemical Physics•Dec 1, 2000•132 citations

Monte Carlo Bohmian Dynamics from Trajectory Stability Properties

Jian Liu, Nancy Makri

It is shown that the quantum force in the Bohmian formulation of quantum mechanics can be related to the stability properties of the given trajectory. In turn, the evolution of the stability properties is governed by higher order derivatives of the quantum potential, leading to an infinite hierarchy of coupled differential equations whose solution specifies completely all aspects of the dynamics. Neglecting derivatives of the quantum potential beyond a certain order allows truncation of the hierarchy, leading to approximate Bohmian trajectories. Use of the method in conjunction with Bohmian initial value formulations [ J. Chem. Phys. 2003, 119, 60] gives rise to simple position-space representations of observables or time correlation functions. These are analogous to approximate quasiclassical expressions based on the Wigner or Husimi phase space density but involve lower dimensional integrals with smoother integrands and avoid the costly evaluation of phase space transforms. The lowest-order version of the truncated hierarchy can capture large corrections to classical mechanical treatments and yields (with fewer trajectories) results that are somewhat more accurate than those based on quasiclassical phase space treatments.

The Journal of Physical Chemistry A•May 28, 2004•40 citations

Applicability criterion for semiclassical Bohmian dynamics

Sophya V. Garashchuk, Vitaly A. Rassolov

Abstract In this paper we develop a criterion of the applicability for the trajectory‐based semiclassical propagation methods that are exact for Gaussian wavepackets in locally quadratic potentials. The error estimate takes into account a deviation of the wavefunction from the Gaussian form due to the anharmonicity of a potential. The derivation is based on the time‐dependent noncoherent solutions of the Schrödinger equation for a general quadratic potential that form a complete orthonormal time‐dependent basis. The anharmonicity is treated as a time‐dependent cubic perturbation, which is further related to the nonlinearity of a force acting on trajectories. The criterion is expressed in terms of scalar products of the average force, its first moments, and dispersion of the wavefunction and is efficient in the context of dynamics with the Bohmian trajectories. We give an application for the Bohmian dynamics with the approximate quantum potential, for the Gaussian wavepacket dynamics and for the semiclassical initial value representation propagator. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004

International Journal of Quantum Chemistry•Aug 30, 2004•2 citations

Bohmian quantum trajectories from coherent states

Sanjib Dey, Andreas Fring

We find that real and complex Bohmian quantum trajectories resulting from well-localized Klauder coherent states in the quasi-Poissonian regime possess qualitatively the same type of trajectories as those obtained from a purely classical analysis of the corresponding Hamilton-Jacobi equation. In the complex cases we treated the quantum potential results to a constant, such that the agreement is exact. For the real cases we provide conjectures for analytical solutions for the trajectories as well as the corresponding quantum potentials. The overall qualitative behavior is governed by the Mandel parameter determining the regime in which the wave functions evolve as solitonlike structures. We demonstrate these features explicitly for the harmonic oscillator and the P\"oschl-Teller potential.

arXiv:1305.4619•Physical Review A•Aug 19, 2013•37 citations

Semiclassical dynamics with quantum trajectories: Formulation and comparison with the semiclassical initial value representation propagator

Sophya V. Garashchuk, Vitaly A. Rassolov

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.

The Journal of Chemical Physics•Jan 30, 2003•54 citations
We develop a Bohmian analysis of a two-dimensional ghost Hamiltonian and its mapping to the degenerate Pais-Uhlenbeck model. Using Gaussian wavepackets, we derive the corresponding guidance equations, the centre and width evolution, and the quantum potential. We use these quantities to characterise bounded, quasi-semiclassical, spiral, and runaway regimes. The Bohmian trajectories provide a direct dynamical diagnostic of coherence, packet deformation, and quantum-classical separation. We then compare a bi-Hamiltonian pair consisting of the ghost Hamiltonian and a classically equivalent alternative formulation. While the two descriptions produce identical classical trajectories, they lead to different Bohmian trajectories and different quantum potentials evaluated along those trajectories. This demonstrates that classical equivalence need not extend to Bohmian quantum dynamics and identifies a concrete quantum ambiguity in the degenerate higher-derivative system.
The origin of quantum interference characteristic of bound nonlinear systems is investigated within the Bohmian formulation of time-dependent quantum mechanics. By contrast to time-dependent semiclassical theory, whereby interference is a consequence of phase mismatch between distinct classical trajectories, the Bohmian, fully quantum mechanical expression for expectation values has a quasiclassical appearance that does not involve phase factors or cross terms. Numerical calculations reveal that quantum interference in the Bohmian formulation manifests itself directly as sharp spatial/temporal variations of the density surrounding kinky trajectories. These effects are most dramatic in regions where the underlying classical motion exhibits focal points or caustics, and crossing of the Bohmian trajectories is prevented through extremely strong and rapidly varying quantum mechanical forces. These features of Bohmian dynamics, which constitute the hallmark of quantum interference and are ubiquitous in bound nonlinear systems, represent a major source of instability, making the integration of the Bohmian equations extremely demanding in such situations.
A novel time-dependent hybrid quantum/classical propagation scheme based on Bohmian quantum trajectories is presented. The quantum subsystem is described by a wave packet depending on the quantum variables x and, via the total potential energy of the system, parametrically on the classical trajectories X(t). The wave packet is used to calculate de Broglie–Bohm quantum trajectories x(t) which are used to calculate the force acting on the classical variables. Quantum corrections of the classical equation of motion are also included. The method is applied to a simple case of two coupled oscillators. Comparison between exact quantum and approximate results demonstrates that these MQCB (Mixed Quantum/Classical Bohmian) trajectories provide a good description of the energy exchange between the two oscillators.
It is shown that the quantum force in the Bohmian formulation of quantum mechanics can be related to the stability properties of the given trajectory. In turn, the evolution of the stability properties is governed by higher order derivatives of the quantum potential, leading to an infinite hierarchy of coupled differential equations whose solution specifies completely all aspects of the dynamics. Neglecting derivatives of the quantum potential beyond a certain order allows truncation of the hierarchy, leading to approximate Bohmian trajectories. Use of the method in conjunction with Bohmian initial value formulations [ J. Chem. Phys. 2003, 119, 60] gives rise to simple position-space representations of observables or time correlation functions. These are analogous to approximate quasiclassical expressions based on the Wigner or Husimi phase space density but involve lower dimensional integrals with smoother integrands and avoid the costly evaluation of phase space transforms. The lowest-order version of the truncated hierarchy can capture large corrections to classical mechanical treatments and yields (with fewer trajectories) results that are somewhat more accurate than those based on quasiclassical phase space treatments.
Abstract In this paper we develop a criterion of the applicability for the trajectory‐based semiclassical propagation methods that are exact for Gaussian wavepackets in locally quadratic potentials. The error estimate takes into account a deviation of the wavefunction from the Gaussian form due to the anharmonicity of a potential. The derivation is based on the time‐dependent noncoherent solutions of the Schrödinger equation for a general quadratic potential that form a complete orthonormal time‐dependent basis. The anharmonicity is treated as a time‐dependent cubic perturbation, which is further related to the nonlinearity of a force acting on trajectories. The criterion is expressed in terms of scalar products of the average force, its first moments, and dispersion of the wavefunction and is efficient in the context of dynamics with the Bohmian trajectories. We give an application for the Bohmian dynamics with the approximate quantum potential, for the Gaussian wavepacket dynamics and for the semiclassical initial value representation propagator. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004
We find that real and complex Bohmian quantum trajectories resulting from well-localized Klauder coherent states in the quasi-Poissonian regime possess qualitatively the same type of trajectories as those obtained from a purely classical analysis of the corresponding Hamilton-Jacobi equation. In the complex cases we treated the quantum potential results to a constant, such that the agreement is exact. For the real cases we provide conjectures for analytical solutions for the trajectories as well as the corresponding quantum potentials. The overall qualitative behavior is governed by the Mandel parameter determining the regime in which the wave functions evolve as solitonlike structures. We demonstrate these features explicitly for the harmonic oscillator and the P\"oschl-Teller potential.
We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.