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:2605.05658

Quantum-classical solvation hydrodynamics: a Hamiltonian modeling framework

François Gay-Balmaz, Cesare Tronci

physics.chem-phphysics.comp-phphysics.flu-dynquant-phSubmitted May 7, 2026• Announced May 8, 2026• Revised Sep 1, 2026 (v3)
v1 · May 7, 2026v2 · May 21, 2026v3 · Sep 1, 2026v4 · Oct 3, 2026
arXivHTMLPDF

Abstract

We propose a mixed quantum-classical hydrodynamic framework to model short-time inertial effects in the non-adiabatic evolution of a quantum solute coupled to a classical polar solvent. Drawing upon the work of Burghardt and Bagchi [Chem. Phys. 329 (2006), 343], we employ the Hamiltonian approach to incorporate consistent backreaction and preserve quantum decoherence beyond standard Ehrenfest dynamics. The solvent is treated as an ideal polar fluid and the quantum solute state is coupled to both the position and molecular orientation coordinates of the liquid. This approach retains essential solute-solvent correlations while reducing the computational complexity of previous approaches. We further incorporate dissipative terms to capture both inertial effects and polarization relaxation. After establishing the general setting for non-local dielectric continua, the Marcus local approximation is integrated into the model thereby extending traditional solvation theory to account for collective fluid sloshing on fast timescales.

Comments
35 pages, two appendices. Revised manuscript. Comments welcome

Related papers

Dynamics of coupled Bohmian and phase-space variables: A moment approach to mixed quantum-classical dynamics

Irène Burghardt

The theoretical framework of the mixed quantum-classical description given by Burghardt and Parlant [J. Chem. Phys. 120, 3055 (2004)] is detailed. A representation in terms of partial hydrodynamic moments is developed, the dynamics of which is determined by a hierarchy of equations derived from the quantum Liouville equation. Exact equations of motion are obtained, whose quantum-classical approximants are associated with a fluid-dynamical trajectory representation which couples classical variables to quantum hydrodynamic variables. The latter evolve under a generalized hydrodynamic force which also depends upon the classical phase-space variables. The hydrodynamic moment description is shown to be closely connected to mixed quantum-classical phase-space methods.

The Journal of Chemical Physics•Feb 22, 2005•27 citations

A Hybrid Hydrodynamic−Liouvillian Approach to Mixed Quantum−Classical Dynamics: Application to Tunneling in a Double Well

Keith H. Hughes, Steven M. Parry, Gérard Parlant +1 more

The hybrid quantum-classical approach of Burghardt and Parlant [Burghardt, I.; Parlant, G. J. Chem. Phys. 2004, 120, 3055], referred to here as the quantum-classical moment (QCM) approach, is demonstrated for the dynamics of a quantum double well coupled to a classical harmonic coordinate. The approach combines the quantum hydrodynamic and classical Liouvillian representations by the construction of a particular type of moments (that is, partial hydrodynamic moments) whose evolution is determined by a hierarchy of coupled equations. For pure states, which are at the center of the present study, this hierarchy terminates at the first order. In the Lagrangian picture, the deterministic trajectories result in dynamics which is Hamiltonian in the classical subspace, while the projection onto the quantum subspace evolves under a generalized hydrodynamic force. Importantly, this force also depends upon the classical (Q, P) variables. The present application demonstrates the tunneling dynamics in both the Eulerian and Lagrangian representations. The method is exact if the classical subspace is harmonic, as is the case for the systems studied here.

The Journal of Physical Chemistry A•Aug 9, 2007•14 citations

Quantum Dissipation in Unbounded Systems

Eric R. Bittner, Jeremy B. Maddox

In recent years trajectory based methodologies have become increasingly popular for evaluating the time evolution of quantum systems. A revival of the de Broglie Bohm interpretation of quantum mechanics has spawned several such techniques for examining quantum dynamics from a hydrodynamic perspective. Using techniques similar to those found in computational fluid dynamics one can construct the wavefunction of a quantum system at any time from the trajectories of a discrete ensemble of hydrodynamic fluid elements (Bohm particles) which evolve according to non-classical equations of motion. Until very recently these schemes have been limited to closed systems. In this paper we present our methodology for including the effects of a thermal environment into the hydrodynamic formulation of quantum dynamics. We derive hydrodynamic equations of motion from the Caldeira-Leggett master equation for the reduced density matrix and give a brief overview of our computational scheme which incorporates an adaptive Lagrangian mesh. Our applications focus upon the dissipative dynamics of open unbounded quantum systems. Using both the Wigner phase space representation and the linear entropy we probe the breakdown of the Markov approximation of the bath dynamics at low temperatures. We suggest a criteria for rationalizing the validity of the Markov approximation in open unbound systems and discuss decoherence, energy relaxation, and quantum/classical correspondence in the context of the Bohmian paths.

•Sep 1, 2001•0 citations

Decoherent histories and nonadiabatic quantum molecular dynamics simulations

Eric R. Bittner, Peter J. Rossky

The role of quantum coherence loss in mixed quantum-classical dynamical systems is explored in the context of the theory of quantum decoherence introduced recently by Bittner and Rossky [J. Chem. Phys. 103, 8130 (1995)]. This theory, which is based upon the consistent histories interpretation of quantum mechanics, introduces decoherence in the quantum subsystem by carefully considering the relevant time and length scales over which one must consider the effects of phase interference between alternative histories of the classical subsystem. Such alternative histories are an integral part of any quantum-classical computational scheme which employs transitions between discrete quantum states; consequently, the coherences between alternative histories have a profound effect on the transition probability between quantum states. In this paper, we review the Bittner–Rossky theory and detail a computational algorithm suitable for large-scale quantum molecular dynamics simulations which implements this theory. Application of the algorithm towards the relaxation of a photoexcited aqueous electron compare well to previous estimates of the excited state survival time as well as to the experimental measurements.

arXiv:physics/9708021•The Journal of Chemical Physics•Nov 22, 1997•91 citations

Quantum dissipation in unbounded systems

Jeremy B. Maddox, Eric R. Bittner

In recent years trajectory based methodologies have become increasingly popular for evaluating the time evolution of quantum systems. A revival of the de Broglie--Bohm interpretation of quantum mechanics has spawned several such techniques for examining quantum dynamics from a hydrodynamic perspective. Using techniques similar to those found in computational fluid dynamics one can construct the wave function of a quantum system at any time from the trajectories of a discrete ensemble of hydrodynamic fluid elements (Bohm particles) which evolve according to nonclassical equations of motion. Until very recently these schemes have been limited to conservative systems. In this paper, we present our methodology for including the effects of a thermal environment into the hydrodynamic formulation of quantum dynamics. We derive hydrodynamic equations of motion from the Caldeira-Leggett master equation for the reduced density matrix and give a brief overview of our computational scheme that incorporates an adaptive Lagrangian mesh. Our applications focus upon the dissipative dynamics of open unbounded quantum systems. Using both the Wigner phase space representation and the linear entropy, we probe the breakdown of the Markov approximation of the bath dynamics at low temperatures. We suggest a criteria for rationalizing the validity of the Markov approximation in open unbound systems and discuss decoherence, energy relaxation, and quantum/classical correspondence in the context of the Bohmian paths.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 25, 2002•37 citations

Path integral methods for quantum dynamics in condensed phases

Tuseeta Banerjee

Quantum mechanical effects play an important role in dynamics of condensed phases. It is also well known that the difficulty of solving the full Schrödinger equation grows exponentially with time. One approach to cope up with this difficulty is to restrict quantum treatment to a few particles (system) and treat the dynamics of the rest of particles (bath) using classical trajectories. Path integral methods due to their trajectory like nature provide excellent tool to develop quantum-classical methods, but to capture long time dynamics, the number of classical trajectories grow exponentially. A part of this work focuses on improving Quantum-Classical Path Integral (QCPI) treatments that allow for larger path integral time step by building “smarter” propagators. These improved solvent-driven reference propagators is developed by incorporating physically motivated approximations to the solvent. When used in QCPI expression they allow convergence with larger time steps leading to an exponential reduction of the number of trajectories required. Further advantages of these propagators include improvements in path filtering techniques, smaller number of path integral steps for achieving the memory decoherence time, and smoothing of the integrand which leads to convergence with fewer Monte Carlo sample points. These ideas have been validated on the spin boson model - a prototypical model to study condensed phase dynamics, which consists of a two-level system coupled to a harmonic bath. The new approach of building the reference propagators combined with their iterative evaluation and filtering is validated using parameters that mimic the first electron transfer in wild-type photosynthetic reaction centres. A real-world dynamical simulation contains various anharmonic effects, Ferrocene (donor) - Ferrocenium (acceptor) system in liquid Benzene was studied where the anharmonic effects of the bath of liquid Benzene on the dynamics of the system were treated using Ensemble Averaged Classical Path. The process of electron or proton transfer occurs between two species having different charges and fluctuations in the solvent plays a key role in determining the rates of such reaction dynamics. To correctly account for the effects of initial distribution of the solvent on the dynamics of the quantum system, the concept of equilibrating the bath to the donor state of the system is used. Two different approaches have been used within the influence functional framework.

•Dec 2, 2015•0 citations
We propose a mixed quantum-classical hydrodynamic framework to model short-time inertial effects in the non-adiabatic evolution of a quantum solute coupled to a classical polar solvent. Drawing upon the work of Burghardt and Bagchi [Chem. Phys. 329 (2006), 343], we employ the Hamiltonian approach to incorporate consistent backreaction and preserve quantum decoherence beyond standard Ehrenfest dynamics. The solvent is treated as an ideal polar fluid and the quantum solute state is coupled to both the position and molecular orientation coordinates of the liquid. This approach retains essential solute-solvent correlations while reducing the computational complexity of previous approaches. We further incorporate dissipative terms to capture both inertial effects and polarization relaxation. After establishing the general setting for non-local dielectric continua, the Marcus local approximation is integrated into the model thereby extending traditional solvation theory to account for collective fluid sloshing on fast timescales.
The theoretical framework of the mixed quantum-classical description given by Burghardt and Parlant [J. Chem. Phys. 120, 3055 (2004)] is detailed. A representation in terms of partial hydrodynamic moments is developed, the dynamics of which is determined by a hierarchy of equations derived from the quantum Liouville equation. Exact equations of motion are obtained, whose quantum-classical approximants are associated with a fluid-dynamical trajectory representation which couples classical variables to quantum hydrodynamic variables. The latter evolve under a generalized hydrodynamic force which also depends upon the classical phase-space variables. The hydrodynamic moment description is shown to be closely connected to mixed quantum-classical phase-space methods.
The hybrid quantum-classical approach of Burghardt and Parlant [Burghardt, I.; Parlant, G. J. Chem. Phys. 2004, 120, 3055], referred to here as the quantum-classical moment (QCM) approach, is demonstrated for the dynamics of a quantum double well coupled to a classical harmonic coordinate. The approach combines the quantum hydrodynamic and classical Liouvillian representations by the construction of a particular type of moments (that is, partial hydrodynamic moments) whose evolution is determined by a hierarchy of coupled equations. For pure states, which are at the center of the present study, this hierarchy terminates at the first order. In the Lagrangian picture, the deterministic trajectories result in dynamics which is Hamiltonian in the classical subspace, while the projection onto the quantum subspace evolves under a generalized hydrodynamic force. Importantly, this force also depends upon the classical (Q, P) variables. The present application demonstrates the tunneling dynamics in both the Eulerian and Lagrangian representations. The method is exact if the classical subspace is harmonic, as is the case for the systems studied here.
In recent years trajectory based methodologies have become increasingly popular for evaluating the time evolution of quantum systems. A revival of the de Broglie Bohm interpretation of quantum mechanics has spawned several such techniques for examining quantum dynamics from a hydrodynamic perspective. Using techniques similar to those found in computational fluid dynamics one can construct the wavefunction of a quantum system at any time from the trajectories of a discrete ensemble of hydrodynamic fluid elements (Bohm particles) which evolve according to non-classical equations of motion. Until very recently these schemes have been limited to closed systems. In this paper we present our methodology for including the effects of a thermal environment into the hydrodynamic formulation of quantum dynamics. We derive hydrodynamic equations of motion from the Caldeira-Leggett master equation for the reduced density matrix and give a brief overview of our computational scheme which incorporates an adaptive Lagrangian mesh. Our applications focus upon the dissipative dynamics of open unbounded quantum systems. Using both the Wigner phase space representation and the linear entropy we probe the breakdown of the Markov approximation of the bath dynamics at low temperatures. We suggest a criteria for rationalizing the validity of the Markov approximation in open unbound systems and discuss decoherence, energy relaxation, and quantum/classical correspondence in the context of the Bohmian paths.
The role of quantum coherence loss in mixed quantum-classical dynamical systems is explored in the context of the theory of quantum decoherence introduced recently by Bittner and Rossky [J. Chem. Phys. 103, 8130 (1995)]. This theory, which is based upon the consistent histories interpretation of quantum mechanics, introduces decoherence in the quantum subsystem by carefully considering the relevant time and length scales over which one must consider the effects of phase interference between alternative histories of the classical subsystem. Such alternative histories are an integral part of any quantum-classical computational scheme which employs transitions between discrete quantum states; consequently, the coherences between alternative histories have a profound effect on the transition probability between quantum states. In this paper, we review the Bittner–Rossky theory and detail a computational algorithm suitable for large-scale quantum molecular dynamics simulations which implements this theory. Application of the algorithm towards the relaxation of a photoexcited aqueous electron compare well to previous estimates of the excited state survival time as well as to the experimental measurements.
In recent years trajectory based methodologies have become increasingly popular for evaluating the time evolution of quantum systems. A revival of the de Broglie--Bohm interpretation of quantum mechanics has spawned several such techniques for examining quantum dynamics from a hydrodynamic perspective. Using techniques similar to those found in computational fluid dynamics one can construct the wave function of a quantum system at any time from the trajectories of a discrete ensemble of hydrodynamic fluid elements (Bohm particles) which evolve according to nonclassical equations of motion. Until very recently these schemes have been limited to conservative systems. In this paper, we present our methodology for including the effects of a thermal environment into the hydrodynamic formulation of quantum dynamics. We derive hydrodynamic equations of motion from the Caldeira-Leggett master equation for the reduced density matrix and give a brief overview of our computational scheme that incorporates an adaptive Lagrangian mesh. Our applications focus upon the dissipative dynamics of open unbounded quantum systems. Using both the Wigner phase space representation and the linear entropy, we probe the breakdown of the Markov approximation of the bath dynamics at low temperatures. We suggest a criteria for rationalizing the validity of the Markov approximation in open unbound systems and discuss decoherence, energy relaxation, and quantum/classical correspondence in the context of the Bohmian paths.
Quantum mechanical effects play an important role in dynamics of condensed phases. It is also well known that the difficulty of solving the full Schrödinger equation grows exponentially with time. One approach to cope up with this difficulty is to restrict quantum treatment to a few particles (system) and treat the dynamics of the rest of particles (bath) using classical trajectories. Path integral methods due to their trajectory like nature provide excellent tool to develop quantum-classical methods, but to capture long time dynamics, the number of classical trajectories grow exponentially. A part of this work focuses on improving Quantum-Classical Path Integral (QCPI) treatments that allow for larger path integral time step by building “smarter” propagators. These improved solvent-driven reference propagators is developed by incorporating physically motivated approximations to the solvent. When used in QCPI expression they allow convergence with larger time steps leading to an exponential reduction of the number of trajectories required. Further advantages of these propagators include improvements in path filtering techniques, smaller number of path integral steps for achieving the memory decoherence time, and smoothing of the integrand which leads to convergence with fewer Monte Carlo sample points. These ideas have been validated on the spin boson model - a prototypical model to study condensed phase dynamics, which consists of a two-level system coupled to a harmonic bath. The new approach of building the reference propagators combined with their iterative evaluation and filtering is validated using parameters that mimic the first electron transfer in wild-type photosynthetic reaction centres. A real-world dynamical simulation contains various anharmonic effects, Ferrocene (donor) - Ferrocenium (acceptor) system in liquid Benzene was studied where the anharmonic effects of the bath of liquid Benzene on the dynamics of the system were treated using Ensemble Averaged Classical Path. The process of electron or proton transfer occurs between two species having different charges and fluctuations in the solvent plays a key role in determining the rates of such reaction dynamics. To correctly account for the effects of initial distribution of the solvent on the dynamics of the quantum system, the concept of equilibrating the bath to the donor state of the system is used. Two different approaches have been used within the influence functional framework.