skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Modeling Spin-Crossover Dynamics
In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.  more » « less
Award ID(s):
1654547
PAR ID:
10323239
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Annual Review of Physical Chemistry
Volume:
72
Issue:
1
ISSN:
0066-426X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; (, Science Advances)
    We develop an ab initio nonadiabatic molecular dynamics (NAMD) method based on GW plus real-time Bethe-Salpeter equation ( GW + rtBSE-NAMD) for the spin-resolved exciton dynamics. From investigations on MoS 2 , we provide a comprehensive picture of spin-valley exciton dynamics where the electron-phonon (e-ph) scattering, spin-orbit interaction (SOI), and electron-hole (e-h) interactions come into play collectively. In particular, we provide a direct evidence that e-h exchange interaction plays a dominant role in the fast valley depolarization within a few picoseconds, which is in excellent agreement with experiments. Moreover, there are bright-to-dark exciton transitions induced by e-ph scattering and SOI. Our study proves that e-h many-body effects are essential to understand the spin-valley exciton dynamics in transition metal dichalcogenides and the newly developed GW + rtBSE-NAMD method provides a powerful tool for exciton dynamics in extended systems with time, space, momentum, energy, and spin resolution. 
    more » « less
  2. ; ; (, The Journal of Chemical Physics)
    The nonadiabatic states and dynamics are investigated for a linear vibronic coupling Hamiltonian with a static electronic splitting and weak off-diagonal Jahn-Teller coupling through a single vibration with a vibrational-electronic resonance. With a transformation of the electronic basis, this Hamiltonian is also applicable to the anti-correlated vibration in a symmetric homodimer with marginally strong constant off-diagonal coupling, where the non-adiabatic states and dynamics model electronic excitation energy transfer or self-exchange electron transfer. For parameters modeling a free-base naphthalocyanine, the nonadiabatic couplings are deeply quantum mechanical and depend on wavepacket width; scalar couplings are as important as the derivative couplings that are usually interpreted to depend on vibrational velocity in semiclassical curve crossing or surface hopping theories. A colored visualization scheme that fully characterizes the non-adiabatic states using the exact factorization is developed. The nonadiabatic states in this nested funnel have nodeless vibrational factors with strongly avoided zeroes in their vibrational probability densities. Vibronic dynamics are visualized through the vibrational coordinate dependent density of the time-dependent dipole moment in free induction decay. Vibrational motion is amplified by the nonadiabatic couplings, with asymmetric and anisotropic motions that depend upon the excitation polarization in the molecular frame and can be reversed by a change in polarization. This generates a vibrational quantum beat anisotropy in excess of 2/5. The amplitude of vibrational motion can be larger than that on the uncoupled potentials, and the electronic population transfer is maximized within one vibrational period. Most of these dynamics are missed by the adiabatic approximation, and some electronic and vibrational motions are completely suppressed by the Condon approximation of a coordinate-independent transition dipole between adiabatic states. For all initial conditions investigated, the initial nonadiabatic electronic motion is driven towards the lower adiabatic state, and criteria for this directed motion are discussed. 
    more » « less
  3. ; ; ; ; ; (, The Journal of Chemical Physics)
    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
  4. ; ; ; ; (, The Journal of Chemical Physics)
    Chemical relaxation phenomena, including photochemistry and electron transfer processes, form a vigorous area of research in which nonadiabatic dynamics plays a fundamental role. However, for electronic systems with spin degrees of freedom, there are few if any applicable and practical quasiclassical methods. Here, we show that for nonadiabatic dynamics with two electronic states and a complex-valued Hamiltonian that does not obey time-reversal symmetry (as relevant to many coupled nuclear-electronic-spin systems), the optimal semiclassical approach is to generalize Tully’s surface hopping dynamics from coordinate space to phase space. In order to generate the relevant phase-space adiabatic surfaces, one isolates a proper set of diabats, applies a phase gauge transformation, and then diagonalizes the total Hamiltonian (which is now parameterized by both R and P). The resulting algorithm is simple and valid in both the adiabatic and nonadiabatic limits, incorporating all Berry curvature effects. Most importantly, the resulting algorithm allows for the study of semiclassical nonadiabatic dynamics in the presence of spin–orbit coupling and/or external magnetic fields. One expects many simulations to follow as far as modeling cutting-edge experiments with entangled nuclear, electronic, and spin degrees of freedom, e.g., experiments displaying chiral-induced spin selectivity. 
    more » « less
  5. ; ; ; (, Physical Chemistry Chemical Physics)
    null (Ed.)
    An understanding of the role that spin states play in semiconductor surface chemical reactions is currently limited. Herein, we provide evidence of a nonadiabatic reaction involving a localized singlet to triplet thermal excitation of the Si(100) surface dimer dangling bond. By comparing the β-hydrogen elimination kinetics of ethyl adsorbates probed by thermal desorption experiments to electronic structure calculation results, we determined that a coverage-dependent change in mechanism occurs. At low coverage, a nonadiabatic, inter-dimer mechanism is dominant, while adiabatic mechanisms become dominant at higher coverage. Computational results indicate that the spin crossover is rapid near room temperature and the nonadiabatic path is accelerated by a barrier that is 40 kJ mol −1 less than the adiabatic path. Simulated thermal desorption reactions using nonadiabatic transition state theory (NA-TST) for the surface dimer intersystem crossing are in close agreement with experimental observations. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email