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:2012.03569 · Journal of Physics A Mathematical and Theoretical

The bohmion method in nonadiabatic quantum hydrodynamics

Darryl D. Holm, Jonathan Rawlinson, Cesare Tronci

Published Sep 29, 2021• 12 citations
arXivHTMLPDF

Abstract

Abstract Starting with the exact factorization of the molecular wavefunction, this paper presents the results from the numerical implementation in nonadiabatic molecular dynamics of the recently proposed bohmion method. Within the context of quantum hydrodynamics, we introduce a regularized nuclear Bohm potential admitting solutions comprising a train of δ-functions which provide a finite-dimensional sampling of the hydrodynamic flow paths. The bohmion method inherits all the basic conservation laws from its underlying variational structure and captures electronic decoherence. After reviewing the general theory, the method is applied to the well-known Tully models, which are used here as benchmark problems. In the present case of study, we show that the new method accurately reproduces both electronic decoherence and nuclear population dynamics.

DOI
10.1088/1751-8121/ac2ae8

Related papers

The bohmion method in nonadiabatic molecular dynamics

Darryl D. Holm, Jonathan Rawlinson, Cesare Tronci

Starting with the exact factorization of the molecular wavefunction, this paper presents the results from the numerical implementation in ab-initio nonadiabatic dynamics of the recently proposed bohmion method. Within the context of quantum hydrodynamics, the regularized nuclear Bohm potential in the bohmion method admits solutions comprising a train of $\delta$-functions which serve as a finite-dimensional sampling of the hydrodynamic flow paths. In addition, the nonlocal structure of the regularized Bohm potential admits nuclear quantum tunneling events. After reviewing the general theory, the bohmion method is applied to the well-known Tully models, which are used here as benchmark problems for the comparison of the new bohmion method with previous mixed quantum-classical methods.

•Dec 7, 2020•0 citations

Geometry of Nonadiabatic Quantum Hydrodynamics

Michael S. Foskett, Darryl D. Holm, Cesare Tronci

The Hamiltonian action of a Lie group on a symplectic manifold induces a momentum map generalizing Noether’s conserved quantity occurring in the case of a symmetry group. Then, when a Hamiltonian function can be written in terms of this momentum map, the Hamiltonian is called ‘collective’. Here, we derive collective Hamiltonians for a series of models in quantum molecular dynamics for which the Lie group is the composition of smooth invertible maps and unitary transformations. In this process, different fluid descriptions emerge from different factorization schemes for either the wavefunction or the density operator. After deriving this series of quantum fluid models, we regularize their Hamiltonians for finite ħ $\hbar $ by introducing local spatial smoothing. In the case of standard quantum hydrodynamics, the ħ ≠ 0 $\hbar \ne 0$ dynamics of the Lagrangian path can be derived as a finite-dimensional canonical Hamiltonian system for the evolution of singular solutions called ‘Bohmions’, which follow Bohmian trajectories in configuration space. For molecular dynamics models, application of the smoothing process to a new factorization of the density operator leads to a finite-dimensional Hamiltonian system for the interaction of multiple (nuclear) Bohmions and a sequence of electronic quantum states.

arXiv:1807.01031•Acta Applicandae Mathematicae•May 7, 2019•27 citations

Koopmon Trajectories in Nonadiabatic Quantum-Classical Dynamics

Werner Bauer, Paul Bergold, François Gay–Balmaz +1 more

Abstract. In order to alleviate the computational costs of fully quantum nonadiabatic dynamics, we present a mixed quantum-classical (MQC) particle method based on the theory of Koopman wavefunctions. Although conventional MQC models often suffer from consistency issues such as the violation of Heisenberg’s principle, we overcame these difficulties by blending Koopman’s classical mechanics on Hilbert spaces with methods in symplectic geometry. The resulting continuum model enjoys both a variational and a Hamiltonian structure, while its nonlinear character calls for suitable closures. Benefiting from the underlying action principle, here we apply a regularization technique previously developed within our team. This step allows for a singular solution ansatz that introduces the trajectories of computational particles—the koopmons—sampling the Lagrangian classical paths in phase space. In the case of Tully’s nonadiabatic problems, the method reproduces the results of fully quantum simulations with levels of accuracy that are not achieved by standard MQC Ehrenfest simulations. In addition, the koopmon method is computationally advantageous over similar fully quantum approaches, which are also considered in our study. As a further step, we probe the limits of the method by considering the Rabi problem in both the ultrastrong and the deep strong coupling regimes, where MQC treatments appear hardly applicable. In this case, the method succeeds in reproducing parts of the fully quantum results.

arXiv:2312.13878•Multiscale Modeling and Simulation•Oct 14, 2024•4 citations

Electronic quantum trajectories with quantum nuclei

Axel Schild

Quantum trajectory calculations for electrons are a useful tool in the field of molecular dynamics, e.g. to understand processes in ultrafast spectroscopy. They have, however, two limitation: On the one hand, such calculations are typically based on the Born-Oppenheimer approximation (BOA) and the electron dynamics for stationary nuclei is considered, thus neglecting quantum effects of the nuclei. On the other hand, even if the quantum nuclear motion would be taken into account, a BOA dynamics on a single potential energy surface would not provide any electron trajectories because the electronic part is treated as a stationary problem. By using the exact factorization method, we overcome these limitations and generalize the theory of electronic quantum trajectories to a fully quantum-mechanical treatment of the nuclei. After reviewing the time-dependent theory of quantum hydrodynamics and quantum trajectories, we show that the nuclei can be viewed as a quantum clock for the electronic motion and we develop a fully quantum-mechanical clock-dependent version of quantum hydrodynamics. This theory is used to obtain electronic trajectories for quantum nuclei, as is exemplified for a model system of a proton-coupled electron transfer dynamics. Our work generalizes the concept of quantum trajectories and lays the foundations for the development of trajectory-based simulation methods of electron dynamics beyond the BOA.

arXiv:2109.13632•Sep 28, 2021•0 citations

Geometry of quantum hydrodynamics in theoretical chemistry

Michael S. Foskett

This thesis investigates geometric approaches to quantum hydrodynamics (QHD) in order to develop applications in theoretical quantum chemistry. Based upon the momentum map geometric structure of QHD and the associated Lie-Poisson and Euler-Poincaré equations, alternative geometric approaches to the classical limit in QHD are presented. These include a new regularised Lagrangian which allows for singular solutions called 'Bohmions' as well as a 'cold fluid' classical closure quantum mixed states. The momentum map approach to QHD is then applied to the nuclear dynamics in a chemistry model known as exact factorization. The geometric treatment extends existing approaches to include unitary electronic evolution in the frame of the nuclear flow, with the resulting dynamics carrying both Euler-Poincaré and Lie-Poisson structures. A new mixed quantum-classical model is then derived by considering a generalised factorisation ansatz at the level of the molecular density matrix. A new alternative geometric formulation of QHD is then constructed. Introducing a $\mathfrak{u}(1)$ connection as the new fundamental variable provides a new method for incorporating holonomy in QHD, which follows from its constant non-zero curvature. The fluid flow is no longer irrotational and carries a non-trivial circulation theorem, allowing for vortex filament solutions. Finally, non-Abelian connections are then considered in quantum mechanics. The dynamics of the spin vector in the Pauli equation allows for the introduction of an $\mathfrak{so}(3)$ connection whilst a more general $\mathfrak{u}(\mathscr{H})$ connection is introduced from the unitary evolution of a quantum system. This is used to provide a new geometric picture for the Berry connection and quantum geometric tensor, whilst relevant applications to quantum chemistry are then considered.

arXiv:2009.13601•Sep 28, 2020•1 citation

Adiabatic and Nonadiabatic Dynamics with Interacting Quantum Trajectories

Lucien Dupuy, Francesco Talotta, Federica Agostini +3 more

We present a quantum dynamics method based on the propagation of interacting quantum trajectories to describe both adiabatic and nonadiabatic processes within the same formalism. The idea originates from the work of Poirier [ Chem. Phys. 2010, 370, 4–14] and Schiff and Poirier [ J. Chem. Phys. 2012, 136, 031102] on quantum dynamics without wavefunctions. It consists of determining the quantum force arising in the Bohmian hydrodynamic formulation of quantum dynamics using only information about quantum trajectories. The particular time-dependent propagation scheme proposed here results in very stable dynamics. Its performance is discussed by applying the method to analytical potentials in the adiabatic regime, and by combining it with the exact factorization method in the nonadiabatic regime.

Journal of Chemical Theory and Computation•Oct 17, 2022•19 citations
Abstract Starting with the exact factorization of the molecular wavefunction, this paper presents the results from the numerical implementation in nonadiabatic molecular dynamics of the recently proposed bohmion method. Within the context of quantum hydrodynamics, we introduce a regularized nuclear Bohm potential admitting solutions comprising a train of δ-functions which provide a finite-dimensional sampling of the hydrodynamic flow paths. The bohmion method inherits all the basic conservation laws from its underlying variational structure and captures electronic decoherence. After reviewing the general theory, the method is applied to the well-known Tully models, which are used here as benchmark problems. In the present case of study, we show that the new method accurately reproduces both electronic decoherence and nuclear population dynamics.
Starting with the exact factorization of the molecular wavefunction, this paper presents the results from the numerical implementation in ab-initio nonadiabatic dynamics of the recently proposed bohmion method. Within the context of quantum hydrodynamics, the regularized nuclear Bohm potential in the bohmion method admits solutions comprising a train of δ\deltaδ-functions which serve as a finite-dimensional sampling of the hydrodynamic flow paths. In addition, the nonlocal structure of the regularized Bohm potential admits nuclear quantum tunneling events. After reviewing the general theory, the bohmion method is applied to the well-known Tully models, which are used here as benchmark problems for the comparison of the new bohmion method with previous mixed quantum-classical methods.
The Hamiltonian action of a Lie group on a symplectic manifold induces a momentum map generalizing Noether’s conserved quantity occurring in the case of a symmetry group. Then, when a Hamiltonian function can be written in terms of this momentum map, the Hamiltonian is called ‘collective’. Here, we derive collective Hamiltonians for a series of models in quantum molecular dynamics for which the Lie group is the composition of smooth invertible maps and unitary transformations. In this process, different fluid descriptions emerge from different factorization schemes for either the wavefunction or the density operator. After deriving this series of quantum fluid models, we regularize their Hamiltonians for finite ħ $\hbar $ by introducing local spatial smoothing. In the case of standard quantum hydrodynamics, the ħ ≠ 0 ℏ≠0\hbar \ne 0ℏ=0 dynamics of the Lagrangian path can be derived as a finite-dimensional canonical Hamiltonian system for the evolution of singular solutions called ‘Bohmions’, which follow Bohmian trajectories in configuration space. For molecular dynamics models, application of the smoothing process to a new factorization of the density operator leads to a finite-dimensional Hamiltonian system for the interaction of multiple (nuclear) Bohmions and a sequence of electronic quantum states.
Abstract. In order to alleviate the computational costs of fully quantum nonadiabatic dynamics, we present a mixed quantum-classical (MQC) particle method based on the theory of Koopman wavefunctions. Although conventional MQC models often suffer from consistency issues such as the violation of Heisenberg’s principle, we overcame these difficulties by blending Koopman’s classical mechanics on Hilbert spaces with methods in symplectic geometry. The resulting continuum model enjoys both a variational and a Hamiltonian structure, while its nonlinear character calls for suitable closures. Benefiting from the underlying action principle, here we apply a regularization technique previously developed within our team. This step allows for a singular solution ansatz that introduces the trajectories of computational particles—the koopmons—sampling the Lagrangian classical paths in phase space. In the case of Tully’s nonadiabatic problems, the method reproduces the results of fully quantum simulations with levels of accuracy that are not achieved by standard MQC Ehrenfest simulations. In addition, the koopmon method is computationally advantageous over similar fully quantum approaches, which are also considered in our study. As a further step, we probe the limits of the method by considering the Rabi problem in both the ultrastrong and the deep strong coupling regimes, where MQC treatments appear hardly applicable. In this case, the method succeeds in reproducing parts of the fully quantum results.
Quantum trajectory calculations for electrons are a useful tool in the field of molecular dynamics, e.g. to understand processes in ultrafast spectroscopy. They have, however, two limitation: On the one hand, such calculations are typically based on the Born-Oppenheimer approximation (BOA) and the electron dynamics for stationary nuclei is considered, thus neglecting quantum effects of the nuclei. On the other hand, even if the quantum nuclear motion would be taken into account, a BOA dynamics on a single potential energy surface would not provide any electron trajectories because the electronic part is treated as a stationary problem. By using the exact factorization method, we overcome these limitations and generalize the theory of electronic quantum trajectories to a fully quantum-mechanical treatment of the nuclei. After reviewing the time-dependent theory of quantum hydrodynamics and quantum trajectories, we show that the nuclei can be viewed as a quantum clock for the electronic motion and we develop a fully quantum-mechanical clock-dependent version of quantum hydrodynamics. This theory is used to obtain electronic trajectories for quantum nuclei, as is exemplified for a model system of a proton-coupled electron transfer dynamics. Our work generalizes the concept of quantum trajectories and lays the foundations for the development of trajectory-based simulation methods of electron dynamics beyond the BOA.
This thesis investigates geometric approaches to quantum hydrodynamics (QHD) in order to develop applications in theoretical quantum chemistry. Based upon the momentum map geometric structure of QHD and the associated Lie-Poisson and Euler-Poincaré equations, alternative geometric approaches to the classical limit in QHD are presented. These include a new regularised Lagrangian which allows for singular solutions called 'Bohmions' as well as a 'cold fluid' classical closure quantum mixed states. The momentum map approach to QHD is then applied to the nuclear dynamics in a chemistry model known as exact factorization. The geometric treatment extends existing approaches to include unitary electronic evolution in the frame of the nuclear flow, with the resulting dynamics carrying both Euler-Poincaré and Lie-Poisson structures. A new mixed quantum-classical model is then derived by considering a generalised factorisation ansatz at the level of the molecular density matrix. A new alternative geometric formulation of QHD is then constructed. Introducing a u(1)\mathfrak{u}(1)u(1) connection as the new fundamental variable provides a new method for incorporating holonomy in QHD, which follows from its constant non-zero curvature. The fluid flow is no longer irrotational and carries a non-trivial circulation theorem, allowing for vortex filament solutions. Finally, non-Abelian connections are then considered in quantum mechanics. The dynamics of the spin vector in the Pauli equation allows for the introduction of an so(3)\mathfrak{so}(3)so(3) connection whilst a more general u(H)\mathfrak{u}(\mathscr{H})u(H) connection is introduced from the unitary evolution of a quantum system. This is used to provide a new geometric picture for the Berry connection and quantum geometric tensor, whilst relevant applications to quantum chemistry are then considered.
We present a quantum dynamics method based on the propagation of interacting quantum trajectories to describe both adiabatic and nonadiabatic processes within the same formalism. The idea originates from the work of Poirier [ Chem. Phys. 2010, 370, 4–14] and Schiff and Poirier [ J. Chem. Phys. 2012, 136, 031102] on quantum dynamics without wavefunctions. It consists of determining the quantum force arising in the Bohmian hydrodynamic formulation of quantum dynamics using only information about quantum trajectories. The particular time-dependent propagation scheme proposed here results in very stable dynamics. Its performance is discussed by applying the method to analytical potentials in the adiabatic regime, and by combining it with the exact factorization method in the nonadiabatic regime.