More Like this

1. A phase-space semiclassical approach for modeling nonadiabatic nuclear dynamics with electronic spin

   https://doi.org/10.1063/5.0093345  

   Wu, Yanze ; Bian, Xuezhi ; Rawlinson, Jonathan I. ; Littlejohn, Robert G. ; Subotnik, Joseph E. ( July 2022 , 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
2. Constructing spin‐adiabatic states for the modeling of spin‐crossing reactions. I. A shared‐orbital implementation

   https://doi.org/10.1002/qua.26123  

   Tao, Yunwen ; Pei, Zheng ; Bellonzi, Nicole ; Mao, Yuezhi ; Zou, Zhu ; Liang, Wanzhen ; Yang, Zhibo ; Shao, Yihan ( December 2019 , International Journal of Quantum Chemistry)

   In the modeling of spin‐crossing reactions, it has become popular to directly explore the spin‐adiabatic surfaces. Specifically, through constructing spin‐adiabatic states from a two‐state Hamiltonian (with spin‐orbit coupling matrix elements) at each geometry, one can readily employ advanced geometry optimization algorithms to acquire a “transition state” structure, where the spin crossing occurs. In this work, we report the implementation of a fully‐variational spin‐adiabatic approach based on Kohn‐Sham density functional theory spin states (sharing the same set of molecular orbitals) and the Breit‐Pauli one‐electron spin‐orbit operator. For three model spin‐crossing reactions (predissociation of N₂O, singlet‐triplet conversion in CH₂, and CO addition to Fe(CO)₄), the spin‐crossing points were obtained. Our results also indicated the Breit‐Pauli one‐electron spin‐orbit coupling can vary significantly along the reaction pathway on the spin‐adiabatic energy surface. On the other hand, due to the restriction that low‐spin and high‐spin states share the same set of molecular orbitals, the acquired spin‐adiabatic energy surface shows a cusp (ie, a first‐order discontinuity) at the crossing point, which prevents the use of standard geometry optimization algorithms to pinpoint the crossing point. An extension with this restriction removed is being developed to achieve the smoothness of spin‐adiabatic surfaces.

    

   more » « less
3. On modelling cosmic ray sputtering of interstellar grain ices

   https://doi.org/10.1093/mnras/stac2507  

   Paulive, Alec ; Carder, Joshua T. ; Herbst, Eric ( September 2022 , Monthly Notices of the Royal Astronomical Society)

   In the interstellar medium (ISM), the formation of complex organic molecules (COMs) is largely facilitated by surface reactions. However, in cold dark clouds, thermal desorption of COMs is inefficient because of the lack of thermal energy to overcome binding energies to the grain surface. Non-thermal desorption methods are therefore important explanations for the gas-phase detection of many COMs that are primarily formed on grains. Here, we present a new non-thermal desorption process: cosmic ray sputtering of grain ice surfaces based on water, carbon dioxide, and a simple mixed ice. Our model applies estimated rates of sputtering to the three-phase rate equation model nautilus-1.1, where this inclusion results in enhanced gas-phase abundances for molecules produced by grain reactions such as methanol (CH3OH) and methyl formate (HCOOCH3). Notably, species with efficient gas-phase destruction pathways exhibit less of an increase in models with sputtering compared to other molecules. These model results suggest that sputtering is an efficient, non-specific method of non-thermal desorption that should be considered as an important factor in future chemical models.

    

   more » « less
4. Nonadiabatic Hydrogen Tunneling Dynamics for Multiple Proton Transfer Processes with Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory

   https://doi.org/10.1021/acs.jctc.4c00737  

   Dickinson, Joseph A ; Hammes-Schiffer, Sharon ( September 2024 , Journal of Chemical Theory and Computation)

   Proton transfer and hydrogen tunneling play key roles in many processes of chemical and biological importance. The generalized nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method was developed in order to capture hydrogen tunneling effects in systems involving the transfer and tunneling of one or more protons. The generalized NEO-MSDFT method treats the transferring protons quantum mechanically on the same level as the electrons and obtains the delocalized vibronic states associated with hydrogen tunneling by mixing localized NEO-DFT states in a nonorthogonal configuration interaction scheme. Herein, we present the derivation and implementation of analytical gradients for the generalized NEO-MSDFT vibronic state energies and the nonadiabatic coupling vectors between these vibronic states. We use this methodology to perform adiabatic and nonadiabatic dynamics simulations of the double proton transfer reactions in the formic acid dimer and the heterodimer of formamidine and formic acid. The generalized NEO-MSDFT method is shown to capture the strongly coupled synchronous or asynchronous tunneling of the two protons in these processes. Inclusion of vibronically nonadiabatic effects is found to significantly impact the double proton transfer dynamics. This work lays the foundation for a variety of nonadiabatic dynamics simulations of multiple proton transfer systems, such as proton relays and hydrogen-bonding networks. 

   more » « less
5. Entangled spin-polarized excitons from singlet fission in a rigid dimer

   https://doi.org/10.1038/s41467-023-36529-6  

   Dill, Ryan_D ; Smyser, Kori_E ; Rugg, Brandon_K ; Damrauer, Niels_H ; Eaves, Joel_D ( March 2023 , Nature Communications)

   Singlet fission, a process that splits a singlet exciton into a biexciton, has promise in quantum information. We report time-resolved electron paramagnetic resonance measurements on a conformationally well-defined acene dimer molecule, TIPS-BP1', designed to exhibit strongly state-selective relaxation to specific magnetic spin sublevels. The resulting optically pumped spin polarization is a nearly pure initial state from the ensemble. The long-lived spin coherences modulate the signal intrinsically, allowing a measurement scheme that substantially removes noise and uncertainty in the magnetic resonance spectra. A nonadiabatic transition theory with a minimal number of spectroscopic parameters allows the quantitative assignment and interpretation of the spectra. In this work, we show that the rigid dimer TIPS-BP1' supports persistent spin coherences at temperatures far higher than those used in conventional superconducting quantum hardware.

    

   more » « less