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

A partially linearized mapping approach to surface hopping (MASH-PLDM) for nonadiabatic dynamics and nonlinear spectroscopy

Jan Obermeier, Kasra Asnaashari, Jeremy O. Richardson

physics.chem-phquant-phSubmitted Sep 28, 2026• Announced Sep 30, 2026
arXivHTMLPDF

Abstract

We derive a partially linearized version of the mapping approach to surface hopping (MASH) that provides a new method for simulating nonadiabatic dynamics in molecular systems. As in the original MASH formalism, trajectories travel on adiabatic surfaces with hops determined by the dynamics of a spin vector on the Bloch sphere. However, because the new approach is based on the partially linearized density matrix (PLDM), it has two spin vectors, which may be in the same or different hemispheres of the Bloch sphere. The former scenario describes an electronic population moving on a single adiabatic surface, whereas the latter describes a coherence moving on the average surface. We show that this new method, called MASH-PLDM, leads to improved results over the original MASH dynamics for a variety of scattering problems and spin-boson models, especially in cases where coherences play an important role. A further benefit of the new methodology is that it can treat multi-time correlation functions such as those required for nonlinear spectroscopy. We demonstrate this capability with a proof-of-principle simulation of time-resolved pump-probe spectroscopy.

Comments
53 pages, 18 figures

Related papers

A partially linearized spin-mapping approach for nonadiabatic dynamics. I. Derivation of the theory

Jonathan R. Mannouch, Jeremy O. Richardson

We present a new partially linearized mapping-based approach for approximating real-time quantum correlation functions in condensed-phase nonadiabatic systems, called the spin partially linearized density matrix (spin-PLDM) approach. Within a classical trajectory picture, partially linearized methods treat the electronic dynamics along forward and backward paths separately by explicitly evolving two sets of mapping variables. Unlike previously derived partially linearized methods based on the Meyer–Miller–Stock–Thoss mapping, spin-PLDM uses the Stratonovich–Weyl transform to describe the electronic dynamics for each path within the spin-mapping space; this automatically restricts the Cartesian mapping variables to lie on a hypersphere and means that the classical equations of motion can no longer propagate the mapping variables out of the physical subspace. The presence of a rigorously derived zero-point energy parameter also distinguishes spin-PLDM from other partially linearized approaches. These new features appear to give the method superior accuracy for computing dynamical observables of interest when compared with other methods within the same class. The superior accuracy of spin-PLDM is demonstrated in this paper through application of the method to a wide range of spin-boson models as well as to the Fenna–Matthews–Olsen complex.

•Jan 1, 2020•41 citations

Open quantum-classical systems: A hybrid MASH master equation

Kasra Asnaashari, Jeremy O. Richardson

We propose a method which combines the quantum-classical mapping approach to surface hopping (MASH) with the dissipative quantum dynamics of the Lindblad master equation. Like conventional surface-hopping methods, our approach is based on classical trajectories coupled to the dynamics of a quantum subsystem. However, instead of evolving the subsystem wavefunction according to the time-dependent Schrödinger equation, we use stochastic quantum trajectories derived from secular Redfield theory. This enables the simulation of open quantum systems coupled simultaneously to Markovian quantum baths and anharmonic non-Markovian classical degrees of freedom. Applications to the spin--boson model and to the cavity-enhanced fluorescence of an electronically nonadiabatic molecule show excellent agreement with fully quantum-mechanical benchmarks.

quant-pharXiv:2511.05282•The Journal of Chemical Physics•Nov 10, 2025•2 citations

A partially linearized spin-mapping approach for nonadiabatic dynamics. II. Analysis and comparison with related approaches

Jonathan R. Mannouch, Jeremy O. Richardson

In a previous paper [J. R. Mannouch and J. O. Richardson, J. Chem. Phys. 153, 194109 (2020)], we derived a new partially linearized mapping-based classical-trajectory technique called the spin partially linearized density matrix (spin-PLDM) approach. This method describes the dynamics associated with the forward and backward electronic path integrals using a Stratonovich–Weyl approach within the spin-mapping space. While this is the first example of a partially linearized spin-mapping method, fully linearized spin-mapping is already known to be capable of reproducing dynamical observables for a range of nonadiabatic model systems reasonably accurately. Here, we present a thorough comparison of the terms in the underlying expressions for the real-time quantum correlation functions for spin-PLDM and fully linearized spin mapping in order to ascertain the relative accuracy of the two methods. In particular, we show that spin-PLDM contains an additional term within the definition of its real-time correlation function, which diminishes many of the known errors that are ubiquitous for fully linearized approaches. One advantage of partially linearized methods over their fully linearized counterparts is that the results can be systematically improved by re-sampling the mapping variables at intermediate times. We derive such a scheme for spin-PLDM and show that for systems for which the approximation of classical nuclei is valid, numerically exact results can be obtained using only a few “jumps.” Additionally, we implement focused initial conditions for the spin-PLDM method, which reduces the number of classical trajectories that are needed in order to reach convergence of dynamical quantities, with seemingly little difference to the accuracy of the result.

•Jan 1, 2020•34 citations

Surface-hopping dynamics of a spin-boson system

Dónal MacKernan, Giovanni Ciccotti, Raymond Kapral

The spin-boson model is solved within the framework of quantum-classical dynamics using our recently-developed surface-hopping scheme. The quantum-classical equation of motion is expressed in an adiabatic basis and its solution is constructed from an ensemble of trajectories which undergo nonadiabatic transitions and evolve coherently on the adiabatic surfaces. Details of the algorithm for the simulation of the dynamics are presented and the method of simple Monte Carlo sampling used to evaluate the expectation values of observables is discussed. The simulation method is applied to a spin-boson system with a harmonic bath composed of ten oscillators with an Ohmic spectral density. For the spin-boson model the present implementation of quantum-classical dynamics is exact and the results of our surface-hopping simulations are in accord with previous numerically exact results for this model.

The Journal of Chemical Physics•Feb 8, 2002•101 citations

LAND-map, a linearized approach to nonadiabatic dynamics using the mapping formalism

Sara Bonella, D. F. Coker

We present a new approach for calculating quantum time correlation functions for systems whose dynamics exhibits relevant nonadiabatic effects. The method involves partial linearization of the full quantum path-integral expression for the time correlation function written in the nonadiabatic mapping Hamiltonian formalism. Our analysis gives an algorithm which is both numerically efficient and accurate as we demonstrate in test calculations on the spin-boson model where we find results in good agreement with exact calculations. The accuracy of our new approach is comparable to that of calculations performed using other approximate methods over a relatively broad range of model parameters. However, our method converges relatively quickly when compared with most alternative schemes. These findings are very encouraging in view of the application of the new method for studying realistic nonadiabatic model problems in the condensed phase.

The Journal of Chemical Physics•May 11, 2005•151 citations

Assessment of trajectory surface hopping methods in long-time nonadiabatic dynamics

Mohammad Shakiba, Daeho Han, Saikat Mukherjee +1 more

We present an assessment of an extensive set of trajectory surface hopping methodologies for modeling long-time population dynamics in a two-level spin-boson model. The considered methodological recipes involve combinations of three surface hopping approaches with five decoherence correction methods, several decoherence-time computation schemes, and two initial-condition sampling methods. In addition to these combinations, the phase correction approach of Shenvi, Subotnik, and Yang (SSY) is considered. By exploring a wide range of meta-parameters controlling decoherence and dephasing times, we determine the optimal performance of such methods and provide a ranking of the best-performing approaches. We find that inclusion of the SSY correction generally over-accelerates the population relaxation dynamics. We also find a strong dependence of the trajectory surface hopping (TSH) calculations on the type of initial-condition sampling: Wigner sampling leads to overestimated population transfer rates and requires decoherence corrections. The simplified decay of mixing approach is found to perform well in this situation. In contrast, Boltzmann sampling leads to a surprisingly remarkable performance of bare TSH schemes without decoherence, which may be a consequence of error cancellation. We provide a critical discussion of the observed trends in different methods' performance and suggest possible avenues for their further improvement.

The Journal of Chemical Physics•Mar 11, 2026•1 citation
We derive a partially linearized version of the mapping approach to surface hopping (MASH) that provides a new method for simulating nonadiabatic dynamics in molecular systems. As in the original MASH formalism, trajectories travel on adiabatic surfaces with hops determined by the dynamics of a spin vector on the Bloch sphere. However, because the new approach is based on the partially linearized density matrix (PLDM), it has two spin vectors, which may be in the same or different hemispheres of the Bloch sphere. The former scenario describes an electronic population moving on a single adiabatic surface, whereas the latter describes a coherence moving on the average surface. We show that this new method, called MASH-PLDM, leads to improved results over the original MASH dynamics for a variety of scattering problems and spin-boson models, especially in cases where coherences play an important role. A further benefit of the new methodology is that it can treat multi-time correlation functions such as those required for nonlinear spectroscopy. We demonstrate this capability with a proof-of-principle simulation of time-resolved pump-probe spectroscopy.
We present a new partially linearized mapping-based approach for approximating real-time quantum correlation functions in condensed-phase nonadiabatic systems, called the spin partially linearized density matrix (spin-PLDM) approach. Within a classical trajectory picture, partially linearized methods treat the electronic dynamics along forward and backward paths separately by explicitly evolving two sets of mapping variables. Unlike previously derived partially linearized methods based on the Meyer–Miller–Stock–Thoss mapping, spin-PLDM uses the Stratonovich–Weyl transform to describe the electronic dynamics for each path within the spin-mapping space; this automatically restricts the Cartesian mapping variables to lie on a hypersphere and means that the classical equations of motion can no longer propagate the mapping variables out of the physical subspace. The presence of a rigorously derived zero-point energy parameter also distinguishes spin-PLDM from other partially linearized approaches. These new features appear to give the method superior accuracy for computing dynamical observables of interest when compared with other methods within the same class. The superior accuracy of spin-PLDM is demonstrated in this paper through application of the method to a wide range of spin-boson models as well as to the Fenna–Matthews–Olsen complex.
We propose a method which combines the quantum-classical mapping approach to surface hopping (MASH) with the dissipative quantum dynamics of the Lindblad master equation. Like conventional surface-hopping methods, our approach is based on classical trajectories coupled to the dynamics of a quantum subsystem. However, instead of evolving the subsystem wavefunction according to the time-dependent Schrödinger equation, we use stochastic quantum trajectories derived from secular Redfield theory. This enables the simulation of open quantum systems coupled simultaneously to Markovian quantum baths and anharmonic non-Markovian classical degrees of freedom. Applications to the spin--boson model and to the cavity-enhanced fluorescence of an electronically nonadiabatic molecule show excellent agreement with fully quantum-mechanical benchmarks.
In a previous paper [J. R. Mannouch and J. O. Richardson, J. Chem. Phys. 153, 194109 (2020)], we derived a new partially linearized mapping-based classical-trajectory technique called the spin partially linearized density matrix (spin-PLDM) approach. This method describes the dynamics associated with the forward and backward electronic path integrals using a Stratonovich–Weyl approach within the spin-mapping space. While this is the first example of a partially linearized spin-mapping method, fully linearized spin-mapping is already known to be capable of reproducing dynamical observables for a range of nonadiabatic model systems reasonably accurately. Here, we present a thorough comparison of the terms in the underlying expressions for the real-time quantum correlation functions for spin-PLDM and fully linearized spin mapping in order to ascertain the relative accuracy of the two methods. In particular, we show that spin-PLDM contains an additional term within the definition of its real-time correlation function, which diminishes many of the known errors that are ubiquitous for fully linearized approaches. One advantage of partially linearized methods over their fully linearized counterparts is that the results can be systematically improved by re-sampling the mapping variables at intermediate times. We derive such a scheme for spin-PLDM and show that for systems for which the approximation of classical nuclei is valid, numerically exact results can be obtained using only a few “jumps.” Additionally, we implement focused initial conditions for the spin-PLDM method, which reduces the number of classical trajectories that are needed in order to reach convergence of dynamical quantities, with seemingly little difference to the accuracy of the result.
The spin-boson model is solved within the framework of quantum-classical dynamics using our recently-developed surface-hopping scheme. The quantum-classical equation of motion is expressed in an adiabatic basis and its solution is constructed from an ensemble of trajectories which undergo nonadiabatic transitions and evolve coherently on the adiabatic surfaces. Details of the algorithm for the simulation of the dynamics are presented and the method of simple Monte Carlo sampling used to evaluate the expectation values of observables is discussed. The simulation method is applied to a spin-boson system with a harmonic bath composed of ten oscillators with an Ohmic spectral density. For the spin-boson model the present implementation of quantum-classical dynamics is exact and the results of our surface-hopping simulations are in accord with previous numerically exact results for this model.
We present a new approach for calculating quantum time correlation functions for systems whose dynamics exhibits relevant nonadiabatic effects. The method involves partial linearization of the full quantum path-integral expression for the time correlation function written in the nonadiabatic mapping Hamiltonian formalism. Our analysis gives an algorithm which is both numerically efficient and accurate as we demonstrate in test calculations on the spin-boson model where we find results in good agreement with exact calculations. The accuracy of our new approach is comparable to that of calculations performed using other approximate methods over a relatively broad range of model parameters. However, our method converges relatively quickly when compared with most alternative schemes. These findings are very encouraging in view of the application of the new method for studying realistic nonadiabatic model problems in the condensed phase.
We present an assessment of an extensive set of trajectory surface hopping methodologies for modeling long-time population dynamics in a two-level spin-boson model. The considered methodological recipes involve combinations of three surface hopping approaches with five decoherence correction methods, several decoherence-time computation schemes, and two initial-condition sampling methods. In addition to these combinations, the phase correction approach of Shenvi, Subotnik, and Yang (SSY) is considered. By exploring a wide range of meta-parameters controlling decoherence and dephasing times, we determine the optimal performance of such methods and provide a ranking of the best-performing approaches. We find that inclusion of the SSY correction generally over-accelerates the population relaxation dynamics. We also find a strong dependence of the trajectory surface hopping (TSH) calculations on the type of initial-condition sampling: Wigner sampling leads to overestimated population transfer rates and requires decoherence corrections. The simplified decay of mixing approach is found to perform well in this situation. In contrast, Boltzmann sampling leads to a surprisingly remarkable performance of bare TSH schemes without decoherence, which may be a consequence of error cancellation. We provide a critical discussion of the observed trends in different methods' performance and suggest possible avenues for their further improvement.