Nonadiabatic Excited-State Molecular Dynamics Methodologies: Comparison and Convergence

Victor M. Freixas,Alexander J. White,Tammie Nelson,Huajing Song,Dmitry V. Makhov,Dmitrii Shalashilin,Sebastian Fernandez-Alberti,Sergei Tretiak
DOI: https://doi.org/10.1021/acs.jpclett.1c00266
IF: 6.888
2021-03-17
The Journal of Physical Chemistry Letters
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.
What problem does this paper attempt to address?