Quantum Entanglement from Classical Trajectories

Coherence in Chemistry: Foundations and Frontiers

Coherence refers to correlations in waves. Because matter has a wave-particle nature, it is unsurprising that coherence has deep connections with the most contemporary issues in chemistry research (e.g., energy harvesting, femtosecond spectroscopy, molecular qubits and more). But what does the word "coherence" really mean in the context of molecules and other quantum systems? We provide a review of key concepts, definitions, and methodologies, surrounding coherence phenomena in chemistry, and we describe how the terms "coherence" and "quantum coherence" refer to many different phenomena in chemistry. Moreover, we show how these notions are related to the concept of an interference pattern. Coherence phenomena are indeed complex, and ambiguous definitions may spawn confusion. By describing the many definitions and contexts for coherence in the molecular sciences, we aim to enhance understanding and communication in this broad and active area of chemistry.

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.

A mapping approach to surface hopping

We present a nonadiabatic classical-trajectory approach that offers the best of both worlds between fewest-switches surface hopping (FSSH) and quasiclassical mapping dynamics. This mapping approach to surface hopping (MASH) propagates the nuclei on the active adiabatic potential-energy surface, such as in FSSH. However, unlike in FSSH, transitions between active surfaces are deterministic and occur when the electronic mapping variables evolve between specified regions of the electronic phase space. This guarantees internal consistency between the active surface and the electronic degrees of freedom throughout the dynamics. MASH is rigorously derivable from exact quantum mechanics as a limit of the quantum-classical Liouville equation (QCLE), leading to a unique prescription for momentum rescaling and frustrated hops. Hence, a quantum-jump procedure can, in principle, be used to systematically converge the accuracy of the results to that of the QCLE. This jump procedure also provides a rigorous framework for deriving approximate decoherence corrections similar to those proposed for FSSH. We apply MASH to simulate the nonadiabatic dynamics in various model systems and show that it consistently produces more accurate results than FSSH at a comparable computational cost.

On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids

Trajectory-based methods that propagate classical nuclei on multiple quantum electronic states are often used to simulate nonadiabatic processes in the condensed phase. A long-standing problem of these methods is their lack of detailed balance, meaning that they do not conserve the equilibrium distribution. In this article, we investigate ideas for restoring detailed balance in mixed quantum-classical systems by tailoring the previously proposed spin-mapping approach to thermal equilibrium. We find that adapting the spin magnitude can recover the correct long-time populations but is insufficient to conserve the full equilibrium distribution. The latter can however be achieved by a more flexible mapping of the spin onto an ellipsoid, which is constructed to fulfill detailed balance for arbitrary potentials. This ellipsoid approach solves the problem of negative populations that has plagued previous mapping approaches and can therefore be applied also to strongly asymmetric and anharmonic systems. Because it conserves the thermal distribution, the method can also exploit efficient sampling schemes used in standard molecular dynamics, which drastically reduces the number of trajectories needed for convergence. The dynamics does however still have mean-field character, as is observed most clearly by evaluating reaction rates in the golden-rule limit. This implies that although the ellipsoid mapping provides a rigorous framework, further work is required to find an accurate classical-trajectory approximation that captures more properties of the true quantum dynamics.

Non-adiabatic ring polymer molecular dynamics in the phase space of the SU(N) Lie group

We derive the non-adiabatic ring polymer molecular dynamics (RPMD) approach in the phase space of the SU(N) Lie Group. This method, which we refer to as the spin mapping non-adiabatic RPMD (SM-NRPMD), is based on the spin-mapping formalism for the electronic degrees of freedom (DOFs) and ring polymer path-integral description for the nuclear DOFs. Using the Stratonovich-Weyl transform for the electronic DOFs and the Wigner transform for the nuclear DOFs, we derived an exact expression of the Kubo-transformed time-correlation function (TCF). We further derive the spin mapping non-adiabatic Matsubara dynamics using the Matsubara approximation that removes the high frequency nuclear normal modes in the TCF and derive the SM-NRPMD approach from the non-adiabatic Matsubara dynamics by discarding the imaginary part of the Liouvillian. The SM-NRPMD method has numerical advantages compared to the original NRPMD method based on the Meyer-Miller-Stock-Thoss (MMST) mapping formalism due to a more natural mapping using the SU(N) Lie Group that preserves the symmetry of the original system. We numerically compute the Kubo-transformed position auto-correlation function and electronic population correlation function for three-state model systems. The numerical results demonstrate the accuracy of the SM-NRPMD method, which outperforms the original MMST-based NRPMD. We envision that the SM-NRPMD method will be a powerful approach to simulate electronic non-adiabatic dynamics and nuclear quantum effects accurately.

Coherent and Dephasing Spectroscopy for Single-Impurity Probing of an Ultracold Bath

We report Ramsey spectroscopy on the clock states of individual Cs impurities immersed in an ultracold Rb bath. We record both the interaction-driven phase evolution and the decay of fringe contrast of the Ramsey interference signal to obtain information about bath density or temperature nondestructively. The Ramsey fringe is modified by a differential shift of the collisional energy when the two Cs states superposed interact with the Rb bath. This differential shift is directly affected by the mean gas density and the details of the Rb-Cs interspecies scattering length, affecting the phase evolution and the contrast of the Ramsey signal. Additionally, we enhance the temperature dependence of the phase shift preparing the system close to a low-magnetic-field Feshbach resonance where the s-wave scattering length is significantly affected by the collisional (kinetic) energy. Analyzing coherent phase evolution and decay of the Ramsey fringe contrast, we probe the Rb cloud's density and temperature. Our results point at using individual impurity atoms as nondestructive quantum probes in complex quantum systems.

Explaining the Efficiency of Photosynthesis: Quantum Uncertainty or Classical Vibrations?

Photosynthetic organisms are known to use a mechanism of vibrationally assisted exciton energy transfer to efficiently harvest energy from light. The importance of quantum effects in this mechanism is a long-standing topic of debate, which has traditionally focused on the role of excitonic coherences. Here, we address another recent claim: that the efficient energy transfer in the Fenna-Matthews-Olson complex relies on nuclear quantum uncertainty and would not function if the vibrations were classical. We present a counter-example to this claim, showing by trajectory-based simulations that a description in terms of quantum electrons and classical nuclei is indeed sufficient to describe the funneling of energy to the reaction center. We analyze and compare these findings to previous classical-nuclear approximations that predicted the absence of an energy funnel and conclude that the key difference and the reason for the discrepancy is the ability of the trajectories to properly account for Newton's third law.