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.

Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules

The Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules

On the adiabatic representation of Meyer-Miller electronic-nuclear dynamics

Nonadiabatic Dynamics Via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity.

Different flavors of nonadiabatic molecular dynamics

Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation

Molecular Dynamics with Nonadiabatic Transitions: A Comparison of Methods

On the representation of coupled adiabatic potential energy surfaces using quasi-diabatic Hamiltonians: a distributed origins expansion approach

Research Advances On Nonadiabatic Energy Transfer Dynamics For Triatomic Molecules

Fractional Derivative Models for Energy Levels and Molar Entropy of Diatomic Molecules

Evaluating the Anharmonicity Contributions to the Molecular Excited State Internal Conversion Rates with Finite Temperature TD-DMRG

Analytic–numerical Approach to Calculate Electronic Nonadiabatic Coupling Terms: Study of the C2H Molecule and the H2 + H System

An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces

Coupled 3D ( J ≥ 0) Time-Dependent Wave Packet Calculation for the F + H 2 Reaction on Accurate Ab Initio Multi-State Diabatic Potential Energy Surfaces

Nonadiabatic Excited-State Molecular Dynamics Methodologies: Comparison and Convergence

Quasi-Diabatization Based on Minimizing Derivative Couplings in a Limited Configuration Space: Elimination of Boundary Condition Dependence

Nonadiabatic Derivative Couplings Calculated Using Information of Potential Energy Surfaces without Wavefunctions: Ab Initio and Machine Learning Implementations

Diatomic Molecules in deSitter and Anti-deSitter Spaces

Tddft Study on Quantization Behaviors of Nonadiabatic Couplings in Polyatomic Systems

Diabatization Based on the Dipole and Quadrupole: the DQ Method

The Curl–Divergence equations for the electronic non-adiabatic coupling terms: study of the C2H molecule and the H2+H system