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.

On-The-Fly Symmetrical Quasi-Classical Dynamics With Meyer-Miller Mapping Hamiltonian For The Treatment Of Nonadiabatic Dynamics At Conical Intersections

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

The Development of the QM/MM Interface and Its Application for the on-the-fly QM/MM Nonadiabatic Dynamics in JADE Package: Theory, Implementation and Applications

Numerical simulation of dynamics in the XY model.

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm-Dancoff approximation

Ab initio Spin-Mapping Non-adiabatic Dynamics Simulations of Photochemistry

Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian.

On the use of QM/MM Surface Hopping simulations to understand thermally-activated rare event nonadiabatic transitions in the condensed phase

Accelerating ab initio QM/MM Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semi-empirical QM/MM Hamiltonian

Quantum Quality with Classical Cost: Ab Initio Nonadiabatic Dynamics Simulations using the Mapping Approach to Surface Hopping

Semiclassical Monte Carlo: a first principles approach to non-adiabatic molecular dynamics

Symmetrical windowing for quantum states in quasi-classical trajectory simulations: application to electron transfer

INAQS, a generic interface for non-adiabatic QM/MM dynamics: Design, implementation, and validation for GROMACS/Q-CHEM simulations

Nonadiabatic Excited-State Molecular Dynamics Methodologies: Comparison and Convergence

Assessing mixed quantum-classical molecular dynamics methods for nonadiabatic dynamics of molecules on metal surfaces

Spin-mapping approach for nonadiabatic molecular dynamics

Ab initio nonadiabatic molecular dynamics investigations on the excited carriers in condensed matter systems

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

CT-MQC – a coupled-trajectory mixed quantum/classical method including nonadiabatic quantum coherence effects

Semiclassical Monte-Carlo approach for modelling non-adiabatic dynamics in extended molecules