An official website of the United States government
We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.
more »
« less
Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation
Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.
more »
« less
- Award ID(s):
- 1856707
- PAR ID:
- 10358117
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 24
- Issue:
- 8
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 4820 to 4831
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Accurate quantum dynamics simulations of nonadiabatic processes are important for studies of electron transfer, energy transfer, and photochemical reactions in complex systems. In this comparative study, we benchmark various approximate nonadiabatic dynamics methods with mapping variables against numerically exact calculations based on the tensor-train (TT) representation of high-dimensional arrays, including TT-KSL for zero-temperature dynamics and TT-thermofield dynamics for finite-temperature dynamics. The approximate nonadiabatic dynamics methods investigated include mixed quantum–classical Ehrenfest mean-field and fewest-switches surface hopping, linearized semiclassical mapping dynamics, symmetrized quasiclassical dynamics, the spin-mapping method, and extended classical mapping models. Different model systems were evaluated, including the spin-boson model for nonadiabatic dynamics in the condensed phase, the linear vibronic coupling model for electronic transition through conical intersections, the photoisomerization model of retinal, and Tully’s one-dimensional scattering models. Our calculations show that the optimal choice of approximate dynamical method is system-specific, and the accuracy is sensitively dependent on the zero-point-energy parameter and the initial sampling strategy for the mapping variables.more » « less
-
In the recently discovered proton-coupled energy transfer (PCEnT) mechanism, the transfer of electronic excitation energy between donor and acceptor chromophores is coupled to a proton transfer reaction. Herein, we develop a general theory for PCEnT and derive an analytical expression for the nonadiabatic PCEnT rate constant. This theory treats the transferring hydrogen nucleus quantum mechanically and describes the PCEnT process in terms of nonadiabatic transitions between reactant and product electron–proton vibronic states. The rate constant is expressed as a summation over these vibronic states, and the contribution of each pair of vibronic states depends on the square of the vibronic coupling as well as the spectral convolution integral, which can be viewed as a generalization of the Förster-type spectral overlap integral for vibronic rather than electronic states. The convolution integral also accounts for the common vibrational modes shared by the donor and acceptor chromophores for intramolecular PCEnT. We apply this theory to model systems to investigate the key features of PCEnT processes. The excited vibronic states can contribute significantly to the total PCEnT rate constant, and the common modes can either slow down or speed up the process. Because the pairs of vibronic states that contribute the most to the PCEnT rate constant may correspond to spectroscopically dark states, PCEnT could occur even when there is no apparent overlap between the donor emission and acceptor absorption spectra. This theory will assist in the interpretation of experimental data and will guide the design of additional PCEnT systems.more » « less
-
Total angular momentum conservation in Ehrenfest dynamics with a truncated basis of adiabatic statesWe show that standard Ehrenfest dynamics does not conserve linear and angular momentum when using a basis of truncated adiabatic states. However, we also show that previously proposed effective Ehrenfest equations of motion [M. Amano and K. Takatsuka, “Quantum fluctuation of electronic wave-packet dynamics coupled with classical nuclear motions,” J. Chem. Phys. 122, 084113 (2005) and V. Krishna, “Path integral formulation for quantum nonadiabatic dynamics and the mixed quantum classical limit,” J. Chem. Phys. 126, 134107 (2007)] involving the non-Abelian Berry force do maintain momentum conservation. As a numerical example, we investigate the Kramers doublet of the methoxy radical using generalized Hartree–Fock with spin–orbit coupling and confirm that angular momentum is conserved with the proper equations of motion. Our work makes clear some of the limitations of the Born–Oppenheimer approximation when using ab initio electronic structure theory to treat systems with unpaired electronic spin degrees of freedom, and we demonstrate that Ehrenfest dynamics can offer much improved, qualitatively correct results.more » « less
-
Polaritons have gained significant attention for the tantalizing possibility of modifying chemical properties and dynamics by coupling molecules to resonant cavity modes to create hybrid light–matter quantum states. Herein, we implement the semiclassical nuclear–electronic orbital time-dependent configuration interaction (NEO-TDCI) approach, which treats electrons and specified nuclei on the same quantum mechanical level, while treating the cavity mode classically. This ab initio dynamics approach can describe both the electronic strong coupling and the vibrational strong coupling regimes at the same level of theory without invoking the Born–Oppenheimer separation between the quantum nuclei and the electrons. This approach is used to simulate resonant and off-resonant vibronic strong coupling, where the cavity mode couples to one or many vibronic transitions associated with joint electronic–nuclear excitations within a vibronic progression. In this case, the cavity mode couples to nuclear motions even for cavity frequencies typically associated with electronic strong coupling. This approach is also used to illustrate that coupling a molecule to a cavity mode can alter hydrogen tunneling dynamics. The semiclassical NEO-TDCI approach provides the foundation for investigating how polaritons may be able to influence chemical reactions involving tunneling and nonadiabatic effects.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1039/d1cp04090a
