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:1607.04854 · Physical Review A

Bohmian mechanics in the exact factorization of electron-nuclear wave functions

Yasumitsu Suzuki, Kazuyuki Watanabe

Published Sep 29, 2016• 32 citations
arXivHTMLPDF

Abstract

The exact factorization of an electron-nuclear wave function [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)] allows us to define the rigorous nuclear time-dependent Schr\"odinger equation (TDSE) with a time-dependent potential-energy surface (TDPES) that fully accounts for the coupling to the electronic motion and drives the nuclear wave-packet dynamics. Here, we study whether the propagation of multiple classical trajectories can reproduce the quantum nuclear motion in strong-field processes when their motions are governed by the quantum Hamilton-Jacobi equation derived by applying Bohmian mechanics to this exact nuclear TDSE. We demonstrate that multiple classical trajectories propagated by the force from the gradient of the exact TDPES plus the Bohmian quantum potential can reproduce the strong-field dissociation dynamics of a one-dimensional model of the ${{\mathrm{H}}_{2}}^{+}$ molecule. Our results show that the force from the Bohmian quantum potential plays a non-negligible role in yielding quantum nuclear dynamics in the strong-field process studied here, where ionization and/or splitting of nuclear probability density occurs.

DOI
10.1103/physreva.94.032517

Related papers

Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test

Francesco Talotta, Federica Agostini, Giovanni Ciccotti

In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.

The Journal of Physical Chemistry A•Jul 30, 2020•24 citations

Coupled-Trajectory Quantum-Classical Approach to Electronic Decoherence in Nonadiabatic Processes

Seung Kyu Min, Federica Agostini, Eberhard K. U. Gross

We present a novel quantum-classical approach to nonadiabatic dynamics, deduced from the coupled electronic and nuclear equations in the framework of the exact factorization of the electron-nuclear wave function. The method is based on the quasiclassical interpretation of the nuclear wave function, whose phase is related to the classical momentum and whose density is represented in terms of classical trajectories. In this approximation, electronic decoherence is naturally induced as an effect of the coupling to the nuclei and correctly reproduces the expected quantum behavior. Moreover, the splitting of the nuclear wave packet is captured as a consequence of the correct approximation of the time-dependent potential of the theory. This new approach offers a clear improvement over Ehrenfest-like dynamics. The theoretical derivation presented in this Letter is supported by numerical results that are compared to quantum mechanical calculations.

arXiv:1504.00250•Physical Review Letters•Aug 10, 2015•180 citations

Bohmian dynamics on subspaces using linearized quantum force

Vitaly A. Rassolov, Sophya V. Garashchuk

In the de Broglie-Bohm formulation of quantum mechanics the time-dependent Schrodinger equation is solved in terms of quantum trajectories evolving under the influence of quantum and classical potentials. For a practical implementation that scales favorably with system size and is accurate for semiclassical systems, we use approximate quantum potentials. Recently, we have shown that optimization of the nonclassical component of the momentum operator in terms of fitting functions leads to the energy-conserving approximate quantum potential. In particular, linear fitting functions give the exact time evolution of a Gaussian wave packet in a locally quadratic potential and can describe the dominant quantum-mechanical effects in the semiclassical scattering problems of nuclear dynamics. In this paper we formulate the Bohmian dynamics on subspaces and define the energy-conserving approximate quantum potential in terms of optimized nonclassical momentum, extended to include the domain boundary functions. This generalization allows a better description of the non-Gaussian wave packets and general potentials in terms of simple fitting functions. The optimization is performed independently for each domain and each dimension. For linear fitting functions optimal parameters are expressed in terms of the first and second moments of the trajectory distribution. Examples are given for one-dimensional anharmonic systems and for the collinear hydrogen exchange reaction.

The Journal of Chemical Physics•Mar 29, 2004•49 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

Dissociation dynamics from a de Broglie–Bohm perspective

Zhen Wang, George R. Darling, S. Holloway

Within the framework of Bohm’s reformulation of quantum physics we revisit the activated dissociation of hydrogen molecules at metal surfaces. The quantum-mechanical force, which accounts for most of quantum effects in the method, and is caused by nonlocal, topographical properties of the wave function, is computed using time-dependent wave packets obtained using conventional, spectral methods. Driven by a combination of the classical force together with the quantum force, trajectories carrying probability density either succeed in overcoming the barrier for dissociation or are scattered back into the gas phase. The Bohmian picture for the dissociation process has enabled us to develop a novel mechanism to account for vibrationally enhanced molecular dissociation. This is relevant to the recently observed promotion of dissociation of very highly vibrationally excited NO molecules at Cu surfaces.

The Journal of Chemical Physics•Dec 8, 2001•28 citations

Complex quantum Hamilton-Jacobi equation with Bohmian trajectories: Application to the photodissociation dynamics of NOCl

Chia‐Chun Chou

The complex quantum Hamilton-Jacobi equation-Bohmian trajectories (CQHJE-BT) method is introduced as a synthetic trajectory method for integrating the complex quantum Hamilton-Jacobi equation for the complex action function by propagating an ensemble of real-valued correlated Bohmian trajectories. Substituting the wave function expressed in exponential form in terms of the complex action into the time-dependent Schrödinger equation yields the complex quantum Hamilton-Jacobi equation. We transform this equation into the arbitrary Lagrangian-Eulerian version with the grid velocity matching the flow velocity of the probability fluid. The resulting equation describing the rate of change in the complex action transported along Bohmian trajectories is simultaneously integrated with the guidance equation for Bohmian trajectories, and the time-dependent wave function is readily synthesized. The spatial derivatives of the complex action required for the integration scheme are obtained by solving one moving least squares matrix equation. In addition, the method is applied to the photodissociation of NOCl. The photodissociation dynamics of NOCl can be accurately described by propagating a small ensemble of trajectories. This study demonstrates that the CQHJE-BT method combines the considerable advantages of both the real and the complex quantum trajectory methods previously developed for wave packet dynamics.

The Journal of Chemical Physics•Mar 12, 2014•21 citations
The exact factorization of an electron-nuclear wave function [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)] allows us to define the rigorous nuclear time-dependent Schr"odinger equation (TDSE) with a time-dependent potential-energy surface (TDPES) that fully accounts for the coupling to the electronic motion and drives the nuclear wave-packet dynamics. Here, we study whether the propagation of multiple classical trajectories can reproduce the quantum nuclear motion in strong-field processes when their motions are governed by the quantum Hamilton-Jacobi equation derived by applying Bohmian mechanics to this exact nuclear TDSE. We demonstrate that multiple classical trajectories propagated by the force from the gradient of the exact TDPES plus the Bohmian quantum potential can reproduce the strong-field dissociation dynamics of a one-dimensional model of the H2+{{\mathrm{H}}_{2}}^{+}H2​+ molecule. Our results show that the force from the Bohmian quantum potential plays a non-negligible role in yielding quantum nuclear dynamics in the strong-field process studied here, where ionization and/or splitting of nuclear probability density occurs.
In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.
We present a novel quantum-classical approach to nonadiabatic dynamics, deduced from the coupled electronic and nuclear equations in the framework of the exact factorization of the electron-nuclear wave function. The method is based on the quasiclassical interpretation of the nuclear wave function, whose phase is related to the classical momentum and whose density is represented in terms of classical trajectories. In this approximation, electronic decoherence is naturally induced as an effect of the coupling to the nuclei and correctly reproduces the expected quantum behavior. Moreover, the splitting of the nuclear wave packet is captured as a consequence of the correct approximation of the time-dependent potential of the theory. This new approach offers a clear improvement over Ehrenfest-like dynamics. The theoretical derivation presented in this Letter is supported by numerical results that are compared to quantum mechanical calculations.
In the de Broglie-Bohm formulation of quantum mechanics the time-dependent Schrodinger equation is solved in terms of quantum trajectories evolving under the influence of quantum and classical potentials. For a practical implementation that scales favorably with system size and is accurate for semiclassical systems, we use approximate quantum potentials. Recently, we have shown that optimization of the nonclassical component of the momentum operator in terms of fitting functions leads to the energy-conserving approximate quantum potential. In particular, linear fitting functions give the exact time evolution of a Gaussian wave packet in a locally quadratic potential and can describe the dominant quantum-mechanical effects in the semiclassical scattering problems of nuclear dynamics. In this paper we formulate the Bohmian dynamics on subspaces and define the energy-conserving approximate quantum potential in terms of optimized nonclassical momentum, extended to include the domain boundary functions. This generalization allows a better description of the non-Gaussian wave packets and general potentials in terms of simple fitting functions. The optimization is performed independently for each domain and each dimension. For linear fitting functions optimal parameters are expressed in terms of the first and second moments of the trajectory distribution. Examples are given for one-dimensional anharmonic systems and for the collinear hydrogen exchange reaction.
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.
Within the framework of Bohm’s reformulation of quantum physics we revisit the activated dissociation of hydrogen molecules at metal surfaces. The quantum-mechanical force, which accounts for most of quantum effects in the method, and is caused by nonlocal, topographical properties of the wave function, is computed using time-dependent wave packets obtained using conventional, spectral methods. Driven by a combination of the classical force together with the quantum force, trajectories carrying probability density either succeed in overcoming the barrier for dissociation or are scattered back into the gas phase. The Bohmian picture for the dissociation process has enabled us to develop a novel mechanism to account for vibrationally enhanced molecular dissociation. This is relevant to the recently observed promotion of dissociation of very highly vibrationally excited NO molecules at Cu surfaces.
The complex quantum Hamilton-Jacobi equation-Bohmian trajectories (CQHJE-BT) method is introduced as a synthetic trajectory method for integrating the complex quantum Hamilton-Jacobi equation for the complex action function by propagating an ensemble of real-valued correlated Bohmian trajectories. Substituting the wave function expressed in exponential form in terms of the complex action into the time-dependent Schrödinger equation yields the complex quantum Hamilton-Jacobi equation. We transform this equation into the arbitrary Lagrangian-Eulerian version with the grid velocity matching the flow velocity of the probability fluid. The resulting equation describing the rate of change in the complex action transported along Bohmian trajectories is simultaneously integrated with the guidance equation for Bohmian trajectories, and the time-dependent wave function is readily synthesized. The spatial derivatives of the complex action required for the integration scheme are obtained by solving one moving least squares matrix equation. In addition, the method is applied to the photodissociation of NOCl. The photodissociation dynamics of NOCl can be accurately described by propagating a small ensemble of trajectories. This study demonstrates that the CQHJE-BT method combines the considerable advantages of both the real and the complex quantum trajectory methods previously developed for wave packet dynamics.