Ensemble Of Classical Trajectories Research Articles

Trajectory-based Simulation of Far-infrared Collision-induced Absorption Profiles of CH4–N2 for Modeling Titan’s Atmosphere

Abstract We report the results of the trajectory-based simulation of far-infrared collision-induced absorption (CIA) due to CH4–N2 pairs at temperatures between 70 and 400 K. Our analysis utilizes recently calculated high-level potential energy and induced dipole surfaces. Treating collision partners as rigid rotors, the time evolution of interaction-induced dipole is accumulated over a vast ensemble of classical trajectories and subsequently transformed into a CIA spectrum via Fourier transform. In our calculations, both bound and unbound states are properly accounted for, and the rigorous theory of lower-order spectral moments is addressed to check the accuracy of simulated profiles. Classically derived trajectory-based profiles are subject to two approximate desymmetrization procedures so that resulting profiles conform to the quantum principle of detailed balance. The simulated profiles are compared to laboratory measurements and employed for modeling Titan’s spectra in the 50–500 cm−1 range. Based on the desymmetrized simulated profiles, a new semiempirical model for CH4–N2 CIA is proposed for modeling Titan’s infrared spectra. Synthetic spectra derived using this model yield an excellent agreement with the data recorded by the Composite Infrared Spectrometer aboard the Cassini spacecraft at low and high emission angles.

A partially linearized spin-mapping approach for simulating nonlinear optical spectra.

We present a partially linearized method based on spin-mapping for computing both linear and nonlinear optical spectra. As observables are obtained from ensembles of classical trajectories, the approach can be applied to the large condensed-phase systems that undergo photosynthetic light-harvesting processes. In particular, the recently derived spin partially linearized density matrix method has been shown to exhibit superior accuracy in computing population dynamics compared to other related classical-trajectory methods. Such a method should also be ideally suited to describing the quantum coherences generated by interaction with light. We demonstrate that this is, indeed, the case by calculating the nonlinear optical response functions relevant for the pump-probe and 2D photon-echo spectra for a Frenkel biexciton model and the Fenna-Matthews-Olsen light-harvesting complex. One especially desirable feature of our approach is that the full spectrum can be decomposed into its constituent components associated with the various Liouville-space pathways, offering a greater insight beyond what can be directly obtained from experiments.

Floquet Phases of Matter via Classical Prethermalization.

We demonstrate that the prethermal regime of periodically driven (Floquet), classical many-body systems can host nonequilibrium phases of matter. In particular, we show that there exists an effective Hamiltonian that captures the dynamics of ensembles of classical trajectories despite the breakdown of this description at the single trajectory level. In addition, we prove that the effective Hamiltonian can host emergent symmetries protected by the discrete time-translation symmetry of the drive. The spontaneous breaking of such an emergent symmetry leads to a subharmonic response, characteristic of time crystalline order, that survives to exponentially late times in the frequency of the drive. To this end, we numerically demonstrate the existence of classical prethermal time crystals in systems with different dimensionalities and ranges of interaction. Extensions to higher order and fractional time crystals are also discussed.

Simulation of collision-induced absorption spectra based on classical trajectories and ab initio potential and induced dipole surfaces. II. CO2-Ar rototranslational band including true dimer contribution.

This paper presents further development of the new semi-classical trajectory-based formalism described in Paper I [Chistikov et al., J. Chem. Phys. 151, 194106 (2019)]. We report the results of simulation and analysis of the low-frequency collision-induced absorption (CIA) in CO2-Ar, including its true dimer component. Our consideration relies on the use of ab initio intermolecular potential energy and induced dipole surfaces for CO2-Ar calculated in an assumption of a rigid CO2 structure using the CCSD(T) method. The theory, the details of which are reported in Paper I [Chistikov et al., J. Chem. Phys. 151, 194106 (2019)], permits taking into account the effect of unbound and quasi-bound classical trajectories on the CIA in the range of a rototranslational band. This theory is largely extended by trajectory-based simulation of the true bound dimer absorption in the present paper. The spectra are obtained from a statistical average over a vast ensemble of classical trajectories restricted by properly chosen domains in the phase space. Rigorous classical theory is developed for two low-order spectral moments interpreted as the Boltzmann-weighted average of the respective dipole functions. These spectral moments were then used to check the accuracy of our trajectory-based spectra, for which both spectral moments can be evaluated independently in terms of specific integrals over the trajectory-based calculated spectral profiles. Good agreement between the spectral moments calculated as integrals over the frequency domain or the phase space largely supports the reliability of our simulated CIA spectra, which conform with the available microwave and far-infrared observations.

Two-dimensional vibrational-electronic spectra with semiclassical mechanics.

Two-dimensional vibrational-electronic (2DVE) spectra probe the effects on vibronic spectra of initial vibrational excitation in an electronic ground state. The optimized mean trajectory (OMT) approximation is a semiclassical method for computing nonlinear spectra from response functions. Ensembles of classical trajectories are subject to semiclassical quantization conditions, with the radiation-matter interaction inducing discontinuous transitions. This approach has been previously applied to two-dimensional infrared and electronic spectra and is extended here to 2DVE spectra. For a system including excitonic coupling, vibronic coupling, and interaction of a chromophore vibration with a resonant environment, the OMT method is shown to well approximate exact quantum dynamics.

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD computer code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Potential applications and caveats to the usage of LD are included as a summary.

Case studies of the time-dependent potential energy surface for dynamics in cavities.

The exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.

Challenges in constructing accurate methods for hydrogen transfer reactions in large biological assemblies: rare events sampling for mechanistic discovery and tensor networks for quantum nuclear effects.

We present two methods that address the computational complexities arising in hydrogen transfer reactions in enzyme active sites. To address the challenge of reactive rare events, we begin with an ab initio molecular dynamics adaptation of the Caldeira-Leggett system-bath Hamiltonian and apply this approach to the study of the hydrogen transfer rate-determining step in soybean lipoxygenase-1. Through direct application of this method to compute an ensemble of classical trajectories, we discuss the critical role of isoleucine-839 in modulating the primary hydrogen transfer event in SLO-1. Notably, the formation of the hydrogen bond between isoleucine-839 and the acceptor-OH group regulates the electronegativity of the donor and acceptor groups to affect the hydrogen transfer process. Curtailing the formation of this hydrogen bond adversely affects the probability of hydrogen transfer. The second part of this paper deals with complementing the rare event sampled reaction pathways obtained from the aforementioned development through quantum nuclear wavepacket dynamics. Essentially the idea is to construct quantum nuclear dynamics on the potential surfaces obtained along the biased trajectories created as noted above. Here, while we are able to obtain critical insights on the quantum nuclear effects from wavepacket dynamics, we primarily engage in providing an improved computational approach for efficient representation of quantum dynamics data such as potential surfaces and transmission probabilities using tensor networks. We find that utilizing tensor networks yields an accurate and efficient description of time-dependent wavepackets, reduced dimensional nuclear eigenstates and associated potential energy surfaces at much reduced cost.

On-the-fly dynamics simulations of transient anions.

A novel theoretical framework for describing the dynamics of transient anions is presented. An ensemble of classical trajectories is propagated on-the-fly, where resonance energies are computed with bound state techniques, and resonance widths are modeled with a combination of bound state and scattering calculations. The methodology was benchmarked against quantum dynamics results for model potential energy curves, and excellent agreement was attained. As a first application, we considered the electron induced dissociation of chloroethane. We found that electron attachment readily stretches the C-Cl bond, which stabilizes the transient anion within ∼10 fs and leads to the release of fast chloride ions. Both magnitude and shape of the computed dissociative electron attachment cross sections are very similar to the available experimental data, even though we found the results to be very sensitive on the accuracy of the underlying methods. These encouraging results place the proposed methodology as a promising approach for studies on transient anions' dynamics of medium sized molecules.

A classical ride through a conical intersection

Regarding the correlated electron-nuclear motion in a model system, we investigate the dynamics in the vicinity of a conical intersection (CoIn) between two excited state potential surfaces. It is documented that an ensemble of classical trajectories which move in the complete electronic-nuclear phase space tracks the quantum wave-packet motion through the CoIn which is accompanied by a non-adiabatic population transfer. On the contrary, for an adiabatic circular motion around the position of the CoIn, the quantum mechanical and classical densities deviate substantially. In the latter case, the Born-Oppenheimer classical nuclear motion on a single potential surface is able to track the quantum dynamics.

Fractional dynamics using an ensemble of classical trajectories.

A trajectory-based formulation for fractional dynamics is presented and the trajectories are generated deterministically. In this theoretical framework, we derive a new class of estimators in terms of confluent hypergeometric function (_{1}F_{1}) to represent the Riesz fractional derivative. Using this method, the simulation of free and confined Lévy flight are in excellent agreement with the exact numerical and analytical results. In addition, the barrier crossing in a bistable potential driven by Lévy noise of index α is investigated. In phase space, the behavior of trajectories reveal the feature of Lévy flight in a better perspective.

Mixed Quantum/Classical Method for Nonadiabatic Quantum Dynamics in Explicit Solvent Models: The ππ*/nπ* Decay of Thymine in Water as a Test Case.

We present a novel mixed quantum classical dynamical method to include solvent effects on internal conversion (IC) processes. All the solute degrees of freedom are represented by a wavepacket moving according to nonadiabatic quantum dynamics, while the motion of an explicit solvent model is described by an ensemble of classical trajectories. The mutual coupling of the solute and solvent dynamics is included within a mean-field framework and the quantum and classical equations of motions are solved simultaneously. As a test case we apply our method to the ultrafast ππ* → nπ* decay of thymine in water. Solvent dynamical response modifies IC yield already on the 50 fs time scale. This effect is due to water librational motions that stabilize the most populated state. Pure static disorder, that is, the existence of different solvent configurations when photoexcitation takes place, also has a remarkable impact on the dynamics.

A mapping variable ring polymer molecular dynamics study of condensed phase proton-coupled electron transfer.

We investigate the mechanisms of condensed phase proton-coupled electron transfer (PCET) using Mapping-Variable Ring Polymer Molecular Dynamics (MV-RPMD), a recently developed method that employs an ensemble of classical trajectories to simulate nonadiabatic excited state dynamics. Here, we construct a series of system-bath model Hamiltonians for the PCET, where four localized electron-proton states are coupled to a thermal bath via a single solvent mode, and we employ MV-RPMD to simulate state population dynamics. Specifically, for each model, we identify the dominant PCET mechanism, and by comparing against rate theory calculations, we verify that our simulations correctly distinguish between concerted PCET, where the electron and proton transfer together, and sequential PCET, where either the electron or the proton transfers first. This work represents a first application of MV-RPMD to multi-level condensed phase systems; we introduce a modified MV-RPMD expression that is derived using a symmetric rather than asymmetric Trotter discretization scheme and an initialization protocol that uses a recently derived population estimator to constrain trajectories to a dividing surface. We also demonstrate that, as expected, the PCET mechanisms predicted by our simulations are robust to an arbitrary choice of the initial dividing surface.

Communication: Adiabatic and non-adiabatic electron-nuclear motion: Quantum and classical dynamics

Using a model for coupled electronic-nuclear motion we investigate the range from negligible to strong non-adiabatic coupling. In the adiabatic case, the quantum dynamics proceeds in a single electronic state, whereas for strong coupling a complete transition between two adiabatic electronic states takes place. It is shown that in all coupling regimes the short-time wave-packet dynamics can be described using ensembles of classical trajectories in the phase space spanned by electronic and nuclear degrees of freedom. We thus provide an example which documents that the quantum concept of non-adiabatic transitions is not necessarily needed if electronic and nuclear motion is treated on the same footing.

State-specific tunneling lifetimes from classical trajectories: H-atom dissociation in electronically excited pyrrole.

A trajectory method of calculating tunneling probabilities from phase integrals along straight line tunneling paths, originally suggested by Makri and Miller [J. Chem. Phys. 91, 4026 (1989)] and recently implemented by Truhlar and co-workers [Chem. Sci. 5, 2091 (2014)], is tested for one- and two-dimensional ab initio based potentials describing hydrogen dissociation in the (1)B1 excited electronic state of pyrrole. The primary observables are the tunneling rates in a progression of bending vibrational states lying below the dissociation barrier and their isotope dependences. Several initial ensembles of classical trajectories have been considered, corresponding to the quasiclassical and the quantum mechanical samplings of the initial conditions. It is found that the sampling based on the fixed energy Wigner density gives the best agreement with the quantum mechanical dissociation rates.

A New Phase Space Density for Quantum Expectations

We introduce a new density for the representation of quantum states on phase space. It is constructed as a weighted difference of two smooth probability densities using the Husimi function and first-order Hermite spectrograms. In contrast to the Wigner function, it is accessible by sampling strategies for positive densities. In the semiclassical regime, the new density allows to approximate expectation values to second order with respect to the high frequency parameter and is thus more accurate than the uncorrected Husimi function. As an application, we combine the new phase space density with Egorov's theorem for the numerical simulation of time-evolved quantum expectations by an ensemble of classical trajectories. We present supporting numerical experiments in different settings and dimensions.

Prethermalization and glassiness in the bosonic Hubbard model

We investigate the nonequilibrium dynamics of the bosonic Hubbard model starting from inhomogeneous superfluid or Mott insulator initial states using the truncated Wigner approximation (TWA). We find that the relaxation of the system develops in two steps for sufficiently large interaction strengths: after a fast relaxation the system gets caught in metastable prethermalized states that precede the true equilibrium state. We find that the lifetime of these prethermalized states increases by several orders of magnitude as we increase the on-site interaction strength beyond a threshold value. We show that the emergence of long-lived metastable states in the quantum dynamics is associated with an ergodic (active) to nonergodic (inactive) dynamical phase transition in the ensemble of classical trajectories that contribute to the semiclassical limit. This dynamical phase transition, which is very similar to that found in different classical models of glasses, is closely related to the dynamic heterogeneity of the classical relaxation.

Calculating time-resolved differential absorbance spectra for ultrafast pump-probe experiments with surface hopping trajectories.

We report a surface hopping approach for modeling the full time- and frequency-resolved differential absorbance spectra (beyond the inhomogenous limit) obtained in ultrafast pump-probe experiments. In our approach, we combine dynamical information obtained from ensembles of classical trajectories propagated on the ground and on the excited potential energy surfaces to directly calculate optical response functions and hence spectral lineshapes. We demonstrate that our method is exact for the model problem of two shifted harmonic potentials with identical harmonic frequencies in the absence of electronic relaxation. We then consider a model three state system with electronic relaxation and show that our method is able to capture the effects of nonadiabatic excited state dynamics on the time-dependent differential absorbance spectra. Furthermore, by comparing our spectra against those spectra calculated with either an (1) inhomogenous expression, (2) ground-state Kubo theory, or (3) excited-state Kubo theory, we show that including dynamical information from both the ground and excited potential energy surfaces significantly improves the reliability of the semiclassical approximations. As such, our surface hopping method should find immediate use in modeling the time-dependent differential abosrbance spectra of ultrafast pump-probe experiments.

Quantum and semiclassical phase-space dynamics of a wave packet in strong fields using initial-value representations

We assess the suitability of quantum and semiclassical initial-value representations (IVRs), exemplified by the coupled coherent states (CCS) method and the Herman–Kluk (HK) propagator, respectively, for modeling the dynamics of an electronic wave packet in a strong laser field, if this wave packet is initially bound. Using Wigner quasiprobability distributions and ensembles of classical trajectories, we identify signatures of over-the-barrier and tunnel ionization in phase space for static and time-dependent fields and the relevant sets of phase-space trajectories to model such features. Overall, we find good agreement with the full solution of the time-dependent Schrödinger equation (TDSE) for Wigner distributions constructed with both IVRs. Our results indicate that the HK propagator does not fully account for tunneling and over-the-barrier reflections. This leads to a dephasing in the time-dependent wave function, which becomes more pronounced for longer times. However, it is able to partly reproduce features associated with the wave packet crossing classically forbidden regions, although the trajectories employed in its construction always obey classical phase-space constraints. We also show that the CCS method represents a fully quantum initial value representation and accurately reproduces the results of a standard TDSE solver. Finally, we show that the HK propagator may be successfully employed to compute the time-dependent dipole acceleration and high-harmonic spectra. Nevertheless, the outcome of the semiclassical computation exhibits disagreements with the TDSE, as a consequence of the previously mentioned dephasing.

How to calculate linear absorption spectra with lifetime broadening using fewest switches surface hopping trajectories: a simple generalization of ground-state Kubo theory.

In this paper, we develop a surface hopping approach for calculating linear absorption spectra using ensembles of classical trajectories propagated on both the ground and excited potential energy surfaces. We demonstrate that our method allows the dipole-dipole correlation function to be determined exactly for the model problem of two shifted, uncoupled harmonic potentials with the same harmonic frequency. For systems where nonadiabatic dynamics and electronic relaxation are present, preliminary results show that our method produces spectra in better agreement with the results of exact quantum dynamics calculations than spectra obtained using the standard ground-state Kubo formalism. As such, our proposed surface hopping approach should find immediate use for modeling condensed phase spectra, especially for expensive calculations using ab initio potential energy surfaces.