Abstract:Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum–classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with <i>ab initio</i> multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jpclett.1c00266?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jpclett.1c00266</a>.Transition densities showing the number of excited states included in the nonadiabatic dynamics and the initial state for each system; AIMC simulation parameters used to generate data plotted in <a class="internalNav" href="#fig5">Figure </a><a class="internalNav" href="#fig5">5</a>; plots of relative differences between adiabatic state populations obtained from EHR and AIMC simulations using different parameters compared to our most accurate AIMC results (<a class="ext-link" href="/doi/suppl/10.1021/acs.jpclett.1c00266/suppl_file/jz1c00266_si_001.pdf">PDF</a>)This article has not yet been cited by other publications.

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. III. Comparison with Born–Oppenheimer dynamics

Ab initio molecular dynamics: Propagating the density matrix with gaussian orbitals. IV. Formal analysis of the deviations from born‐oppenheimer dynamics

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. II. Generalizations based on mass-weighting, idempotency, energy …

Second-Order Mass-Weighting Scheme for Atom-Centered Density Matrix Propagation Molecular Dynamics

Next generation extended Lagrangian first principles molecular dynamics

Effect of exchange-correlation functionals and symmetry constraints of electronic structure on the trajectories of reactive molecular collisions

Nonadiabatic Excited-State Molecular Dynamics Methodologies: Comparison and Convergence

Different flavors of nonadiabatic molecular dynamics

Molecular Dynamics with Nonadiabatic Transitions: A Comparison of Methods

An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics

Extended Lagrangian Born-Oppenheimer Molecular Dynamics with DFT+U

Time evolution of natural orbitals in ab initio molecular dynamics

Nonadiabatic Molecular Dynamics with Subsystem Density Functional Theory: Application to Crystalline Pentacene

Trajectory Surface Hopping with Tight Binding Density Functional Theory applied to Molecular Motors

A Quasi Time-Reversible scheme based on density matrix extrapolation on the Grassmann manifold for Born-Oppenheimer Molecular Dynamics

Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths

Gaussians for Electronic and Rovibrational Quantum Dynamics

The method of Gaussian weighted trajectories. V. On the 1GB procedure for polyatomic processes

Analysis of two overlapping fragmentation approaches in density matrix construction: GMBE-DM vs ADMA