Quantum-classical Dynamics Research Articles

Erratum: Hybrid completely positive Markovian quantum-classical dynamics [Phys. Rev. A 107 , 062206 (2023)

Erratum: Hybrid completely positive Markovian quantum-classical dynamics [Phys. Rev. A <b>107</b> , 062206 (2023)

Markovian master equations for quantum-classical hybrid systems

The problem of constructing a consistent quantum-classical hybrid dynamics is afforded in the case of a quantum component in a separable Hilbert space and a continuous, finite-dimensional classical component. In the Markovian case, the problem is formalized by the notion of hybrid dynamical semigroup. A classical component can be observed without perturbing the system and information on the quantum component can be extracted, thanks to the quantum-classical interaction. This point is formalized by showing how to introduce positive operator valued measures and operations compatible with the hybrid dynamical semigroup; in this way the notion of hybrid dynamics is connected to quantum measurements in continuous time. Then, the case of the most general quasi-free generator is presented and the various quantum-classical interaction terms are discussed. To bee quasi-free means to send, in the Heisenberg description, hybrid Weyl operators into multiples of Weyl operators; the results on the structure of quasi-free semigroups were proved in Reference [12]. Even in the pure quantum case, a quasi-free semigroup is not restricted to have only a Gaussian structure, but also jump-type terms are allowed. An important result is that, to have interactions producing a flow of information from the quantum component to the classical one, suitable dissipative terms must be present in the generator. Finally, some possibilities are discussed to go beyond the quasi-free case.

Mixed quantum-classical dynamics for near term quantum computers

Mixed quantum-classical dynamics is a set of methods often used to understand systems too complex to treat fully quantum mechanically. Many techniques exist for full quantum mechanical evolution on quantum computers, but mixed quantum-classical dynamics are less explored. We present a modular algorithm for general mixed quantum-classical dynamics where the quantum subsystem is coupled with the classical subsystem. We test it on a modified Shin-Metiu model in the first quantization through Ehrenfest propagation. We find that the Time-Dependent Variational Time Propagation algorithm performs well for short-time evolutions and retains qualitative results for longer-time evolutions.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods.

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation.

Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.

Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics?

A common strategy to simulate mixed quantum-classical dynamics is by propagating classical trajectories with mapping variables, often using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian or the related spin-mapping approach. When mapping the quantum subsystem, the coupled dynamics reduce to a set of equations of motion to integrate. Several numerical algorithms have been proposed, but a thorough performance comparison appears to be lacking. Here, we compare three time-propagation algorithms for the MMST Hamiltonian: the Momentum Integral (MInt) (J. Chem. Phys., 2018, 148, 102326), the Split-Liouvillian (SL) (Chem. Phys., 2017, 482, 124-134), and the algorithm in J. Chem. Phys., 2012, 136, 084101 that we refer to as the Degenerate Eigenvalue (DE) algorithm due to the approximation required during derivation. We analyze the accuracy of individual trajectories, correlation functions, energy conservation, symplecticity, Liouville's theorem, and the computational cost. We find that the MInt algorithm is the only rigorously symplectic algorithm. However, comparable accuracy at a lower computational cost can be obtained with the SL algorithm. The approximation implicitly made within the DE algorithm conserves energy poorly, even for small timesteps, and thus leads to slightly different results. These results should guide future mapping-variable simulations.

Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study.

Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. This approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations.

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Mixed Quantum-Classical Liouville Equation Treatment of Electronic Spectroscopy of Condensed Systems: Harmonic and Anharmonic Electron-Phonon Coupling.

This Review integrates the use of electronic optical response function theory and the mixed quantum-classical (MQC) Liouville equation (MQCLE), thereby leading to electronic spectroscopy in MQC media. It further sheds light on the applicability, utility, and efficiency of the mixed quantum-classical dynamics (MQCD) formalism, which starts off with the MQCLE, in probing spectroscopy and dynamics of condensed systems, whereby quantum mechanics and classical mechanics are combined systematically. The author has been exploring and implementing MQCD to investigate electron-phonon coupling effects on electronic dephasing in harmonic and anharmonic systems by calculating linear and nonlinear optical transition analytically and numerically dipole moment time correlation functions in an MQC environment, thereby presenting an in depth spectral profile analysis and their shape and symmetry. The distinctive capability of the MQC time correlation functions is that ergodicity and stationarity properties are inherently satisfied as part of the mixed quantum-classical dynamics (MQCD) framework, unlike classical correlation functions. While some research groups have applied MQCLE to calculate vibrational spectra to study hydrogen-bonded complexes in a MQC environment and other groups calculated Optical response function to probe electron transfer dynamics using the basis mapping technique, the approach, purpose, rigor, applications, and path to the end results reported herein are different. Finally, the same framework is employed to study dissipative systems in the MQC limit, whereby the zero-phonon line adopts the correct width and eliminates its asymmetry. While the full quantum mechanical model, like the multimode Brownian oscillator (MBO) model, yields the correct width and inaccurate shape in the low-temperature limit, the MQCD formalism seems to produce an accurate zero-phonon profile. Nonlinear optical signals are also reviewed in MQC media to show the applicability and utility of this approach. The vibronic optical response functions developed here will account for geometry change, frequency change, and anharmonicity upon electronic excitation to accurately probe electronic dephasing, electron-phonon coupling, shape, and symmetry of profiles and present differences and similarities to the MBO model on pure electronic dephasing. Frequency change and anharmonicity are vitally crucial for accurately assessing electron-phonon coupling upon electronic excitation. This is an additional unique result obtained by the author to further demonstrate the applicability and utility of this approach over other approximation schemes in probing electronic dephasing including that of the MBO model.

A multiscale approach to coupled nuclear and electronic dynamics. II. Exact and approximated evaluation of nonradiative transition rates.

This work follows a companion article, which will be referred to as Paper I [Campeggio et al., J. Chem. Phys. 158, 244104 (2023)] in which a quantum-stochastic Liouville equation for the description of the quantum-classical dynamics of a molecule in a dissipative bath has been formulated in curvilinear internal coordinates. In such an approach, the coordinates of the system are separated into three subsets: the quantum coordinates, the classical relevant nuclear degrees of freedom, and the classical irrelevant (bath) coordinates. The equation has been derived in natural internal coordinates, which are bond lengths, bond angles, and dihedral angles. The resulting equation needs to be parameterized. In particular, one needs to compute the potential energy surfaces, the friction tensor, and the rate constants for the nonradiative jumps among the quantum states. While standard methods exist for the calculation of energy and dissipative properties, an efficient evaluation of the transition rates needs to be developed. In this paper, an approximated treatment is introduced, which leads to a simple explicit formula with a single adjustable parameter. Such an approximated expression is compared with the exact calculation of transition rates obtained via molecular dynamics simulations. To make such a comparison possible, a simple sandbox system has been used, with two quantum states and a single internal coordinate (together with its conjugate momentum). Results show that the adjustable parameter, which is an effective decoherence time, can be parameterized from the effective relaxation times of the autocorrelation functions of the conjugated momenta of the relevant nuclear coordinates.

Capturing electronic decoherence in quantum-classical dynamics using the ring-polymer-surface-hopping–density-matrix approach

Capturing electronic decoherence in quantum-classical dynamics using the ring-polymer-surface-hopping–density-matrix approach

Hybrid completely positive Markovian quantum-classical dynamics

A concise and self-contained derivation of hybrid quantum-classical dynamics is given in terms of Markovian master equations. Many previously known results are rederived and revised and some of them are completed or corrected. Using a method as simple as possible, our goal is a brief introduction to the state of the art of hybrid dynamics, with a limited discussion of the implications for foundations and without discussion of further relevance in the measurement problem, quantum gravity, chemistry, numeric methods, etc. Hybrid dynamics is defined as a special case of composite quantum dynamics where the observables of one of the two subsystems are restricted to the commuting set of diagonal operators in a fixed basis. With this restriction, the derivation of hybrid dynamical equations is clear conceptually and simple technically. Jump and diffusive dynamics follow in the form of hybrid master equations. Their stochastic interpretation (called unravelings) is derived. We discuss gauge-type ambiguities, problems of uniqueness, and covariance of the diffusive master equation. Also conditions of minimum noise and of monitoring the quantum trajectory are derived. We conclude that the hybrid formalism is equivalent to the standard Markovian theory of time-continuous quantum measurement (monitoring) on the one hand and is a motivating alternative formalism on the other.

Energy-conserving coupled trajectory mixed quantum-classical dynamics.

The coupled-trajectory mixed quantum-classical method (CTMQC), derived from the exact factorization approach, has successfully predicted photo-chemical dynamics in a number of interesting molecules, capturing population transfer and decoherence from first principles. However, due to the approximations made, CTMQC does not guarantee energy conservation. We propose a modified algorithm, CTMQC-E, which redefines the integrated force in the coupled-trajectory term so to restore energy conservation, and demonstrate its accuracy on scattering in Tully's extended coupling region model and photoisomerization in a retinal chromophore model.

A Quantum-Classical Model of Brain Dynamics.

The study of the human psyche has elucidated a bipartite structure of logic reflecting the quantum-classical nature of the world. Accordingly, we posited an approach toward studying the brain by means of the quantum-classical dynamics of a mixed Weyl symbol. The mixed Weyl symbol can be used to describe brain processes at the microscopic level and, when averaged over an appropriate ensemble, can provide a link to the results of measurements made at the meso and macro scale. Within this approach, quantum variables (such as, for example, nuclear and electron spins, dipole momenta of particles or molecules, tunneling degrees of freedom, and so on) can be represented by spinors, whereas the electromagnetic fields and phonon modes can be treated either classically or semi-classically in phase space by also considering quantum zero-point fluctuations. Quantum zero-point effects can be incorporated into numerical simulations by controlling the temperature of each field mode via coupling to a dedicated Nosé-Hoover chain thermostat. The temperature of each thermostat was chosen in order to reproduce quantum statistics in the canonical ensemble. In this first paper, we introduce a general quantum-classical Hamiltonian model that can be tailored to study physical processes at the interface between the quantum and the classical world in the brain. While the approach is discussed in detail, numerical calculations are not reported in the present paper, but they are planned for future work. Our theory of brain dynamics subsumes some compatible aspects of three well-known quantum approaches to brain dynamics, namely the electromagnetic field theory approach, the orchestrated objective reduction theory, and the dissipative quantum model of the brain. All three models are reviewed.

Quantum-quasiclassical analysis of center-of-mass nonseparability in hydrogen atom stimulated by strong laser fields *

We have developed a quantum-quasiclassical computational scheme for quantitative treating of the nonseparable quantum–classical dynamics of the 6D hydrogen atom in a strong laser pulse. In this approach, the electron is treated quantum mechanically and the center-of-mass (CM) motion classically. Thus, the Schrödinger equation for the electron and the classical Hamilton equations for the CM variables, nonseparable due to relativistic effects stimulated by strong laser fields, are integrated simultaneously. In this approach, it is natural to investigate the idea of using the CM- velocity spectroscopy as a classical ‘build-up’ set up for detecting the internal electron quantum dynamics. We have performed such an analysis using the hydrogen atom in linearly polarized laser fields as an example and found a strong correlation between the CM kinetic energy distribution after a laser pulse and the spectral density of electron kinetic energy. This shows that it is possible to detect the quantum dynamics of an electron by measuring the distribution of the CM kinetic energy.

Hybrid quantum-classical dynamics of pure-dephasing systems *

We consider the interaction dynamics of a classical oscillator and a quantum two-level system for different pure-dephasing Hamiltonians of the type . This type of systems represents a severe challenge for popular hybrid quantum–classical descriptions. For example, in the case of the common Ehrenfest model, the classical density evolution is shown to decouple entirely from the pure-dephasing quantum dynamics. We focus on a recently proposed hybrid wave equation that is based on Koopman’s wavefunction description of classical mechanics. This model retains quantum–classical correlations whenever a coupling potential is present. Here, several benchmark problems are considered and the results are compared with those arising from fully quantum dynamics. A good agreement is found for a series of study cases involving harmonic oscillators with linear and quadratic coupling, as well as time-varying coupling parameters. In all these cases the classical evolution coincides exactly with the oscillator dynamics resulting from the fully quantum description. In the special case of time-independent coupling involving a classical oscillator with varying frequency, the quantum Bloch rotation exhibits peculiar features that escape from the hybrid description. In addition, nonlinear corrections to the harmonic Hamiltonian lead to an overall growth of decoherence at long times, which is absent in the fully quantum treatment.

Molecular Basis of RNA-Driven ATP Hydrolysis in DExH-Box Helicases.

The spliceosome machinery catalyzes precursor messenger (pre-m)RNA splicing. In each cycle, the spliceosome experiences massive compositional and conformational remodeling fueled by the concerted action of specific RNA-dependent ATPases/helicases. Intriguingly, these enzymes are allosterically activated to perform ATP hydrolysis and trigger helicase activity only upon pre-mRNA binding. Yet, the molecular mechanism underlying the RNA-driven regulation of their ATPase function remains elusive. Here, we focus on the Prp2 ATPase/helicase which contributes to reshaping the spliceosome into its catalytic competent state. By performing classical and quantum-classical molecular dynamics simulations, we unprecedentedly unlock the molecular terms governing the Prp2 ATPase/helicase function. Namely, we dissect the molecular mechanism of ATP hydrolysis, and we disclose that RNA binding allosterically triggers the formation of a set of interactions linking the RNA binding tunnel to the catalytic site. This activates the Prp2's ATPase function by optimally placing the nucleophilic water and the general base of the enzymatic process to perform ATP hydrolysis. The key structural motifs, mechanically coupling RNA gripping and the ATPase/helicase functions, are conserved across all DExH-box helicases. This mechanism could thus be broadly applicable to all DExH-box helicase family.

Dynamics of mixed quantum–classical spin systems *

Mixed quantum–classical spin systems have been proposed in spin chain theory and, more recently, in magnon spintronics. However, current models of quantum–classical dynamics beyond mean-field approximations typically suffer from long-standing consistency issues, and, in some cases, invalidate Heisenberg’s uncertainty principle. Here, we present a fully Hamiltonian theory of quantum–classical spin dynamics that appears to be the first to ensure an entire series of consistency properties, including positivity of both the classical and the quantum density, so that Heisenberg’s principle is satisfied at all times. We show how this theory may connect to recent energy-balance considerations in measurement theory and we present its Poisson bracket structure explicitly. After focusing on the simpler case of a classical Bloch vector interacting with a quantum spin observable, we illustrate the extension of the model to systems with several spins, and restore the presence of orbital degrees of freedom.

Direct Nonadiabatic Dynamics of Ammonia with Curvature-Driven Coherent Switching with Decay of Mixing and with Fewest Switches with Time Uncertainty: An Illustration of Population Leaking in Trajectory Surface Hopping Due to Frustrated Hops

Mixed quantum-classical nonadiabatic dynamics is a widely used approach to simulate molecular dynamics involving multiple electronic states. There are two main categories of mixed quantum-classical nonadiabatic dynamics algorithms, namely, trajectory surface hopping (TSH) in which the trajectory propagates on a single potential energy surface, interrupted by hops, and self-consistent-potential (SCP) methods, such as semiclassical Ehrenfest, in which propagation occurs on a mean-field surface without hops. In this work, we will illustrate an example of severe population leaking in TSH. We emphasize that such leaking is a combined effect of frustrated hops and long-time simulations that drive the final excited-state population toward zero as a function of time. We further show that such leaking can be alleviated-but not eliminated-by the fewest switches with time uncertainty TSH algorithm (here implemented in the SHARC program); the time uncertainty algorithm slows down the leaking process by a factor of 4.1. The population leaking is not present in coherent switching with decay of mixing (CSDM), which is an SCP method with non-Markovian decoherence included. Another result in this paper is that we find very similar results with the original CSDM algorithm, with time-derivative CSDM (tCSDM), and with curvature-driven CSDM (κCSDM). Not only do we find good agreement for electronically nonadiabatic transition probabilities but also we find good agreement of the norms of the effective nonadiabatic couplings (NACs) that are derived from the curvature-driven time-derivative couplings as implemented in κCSDM with the time-dependent norms of the nonadiabatic coupling vectors computed by state-averaged complete-active-space self-consistent field theory.

A Mean-Field Treatment of Vacuum Fluctuations in Strong Light-Matter Coupling.

Mean-field mixed quantum-classical dynamics could provide a much-needed means to inexpensively model quantum electrodynamical phenomena by describing the optical field and its vacuum fluctuations classically. However, this approach is known to suffer from an unphysical transfer of energy out of the vacuum fluctuations when the light-matter coupling becomes strong. We highlight this issue for the case of an atom in an optical cavity and resolve it by introducing an additional set of classical coordinates to specifically represent vacuum fluctuations whose light-matter interaction is scaled by the instantaneous ground-state population of the atom. This not only rigorously prevents the aforementioned unphysical energy transfer but is also shown to yield a radically improved accuracy in terms of the atomic population and the optical field dynamics, generating results in excellent agreement with full quantum calculations. As such, the resulting method emerges as an attractive solution for the affordable modeling of strong light-matter coupling phenomena involving macroscopic numbers of optical modes.