skip to main content

Title: Constructing spin‐adiabatic states for the modeling of spin‐crossing reactions. I. A shared‐orbital implementation

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 N2O, singlet‐triplet conversion in CH2, and CO addition to Fe(CO)4), 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
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal of Quantum Chemistry
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The f-block ab initio correlation consistent composite approach was used to predict the dissociation energies of lanthanide sulfides and selenides. Geometry optimizations were carried out using density functional theory and coupled cluster singles, doubles, and perturbative triples with one- and two-component Hamiltonians. For the two-component calculations, relativistic effects were accounted for by utilizing a third-order Douglas–Kroll–Hess Hamiltonian. Spin–orbit coupling was addressed with the Breit–Pauli Hamiltonian within a multireference configuration interaction approach. The state averaged complete active space self-consistent field wavefunctions obtained for the spin–orbit coupling energies were used to assign the ground states of diatomics, and several diagnostics were used to ascertain the multireference character of the molecules. 
    more » « less
  2. We present a new implementation for computing spin–orbit couplings (SOCs) within a time-dependent density-functional theory (TD-DFT) framework in the standard spin-conserving formulation as well in the spin–flip variant (SF-TD-DFT). This approach employs the Breit–Pauli Hamiltonian and Wigner–Eckart’s theorem applied to the reduced one-particle transition density matrices, together with the spin–orbit mean-field treatment of the two-electron contributions. We use a state-interaction procedure and compute the SOC matrix elements using zero-order non-relativistic states. Benchmark calculations using several closed-shell organic molecules, diradicals, and a single-molecule magnet illustrate the efficiency of the SOC protocol. The results for organic molecules (described by standard TD-DFT) show that SOCs are insensitive to the choice of the functional or basis sets, as long as the states of the same characters are compared. In contrast, the SF-TD-DFT results for small diradicals (CH 2 , [Formula: see text], SiH 2 , and [Formula: see text]) show strong functional dependence. The spin-reversal energy barrier in a Fe(III) single-molecule magnet computed using non-collinear SF-TD-DFT (PBE0, ωPBEh/cc-pVDZ) agrees well with the experimental estimate. 
    more » « less
  3. 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
  4. Abstract

    Chirality has been a property of central importance in physics, chemistry and biology for more than a century. Recently, electrons were found to become spin polarized after transmitting through chiral molecules, crystals, and their hybrids. This phenomenon, called chirality-induced spin selectivity (CISS), presents broad application potentials and far-reaching fundamental implications involving intricate interplays among structural chirality, topological states, and electronic spin and orbitals. However, the microscopic picture of how chiral geometry influences electronic spin remains elusive, given the negligible spin-orbit coupling (SOC) in organic molecules. In this work, we address this issue via a direct comparison of magnetoconductance (MC) measurements on magnetic semiconductor-based chiral molecular spin valves with normal metal electrodes of contrasting SOC strengths. The experiment reveals that a heavy-metal electrode provides SOC to convert the orbital polarization induced by the chiral molecular structure tospinpolarization. Our results illustrate the essential role of SOC in the metal electrode for the CISS spin valve effect. A tunneling model with a magnetochiral modulation of the potential barrier is shown to quantitatively account for the unusual transport behavior.

    more » « less
  5. Abstract

    Following the discovery of topological insulators, there has been a renewed interest in superconducting systems that have strong spin-orbit (SO) coupling. Here we address the fundamental question of how the spin properties of a otherwise spin-singlet superconducting ground state evolve with increasing SO impurity density. We have mapped out the Zeeman critical field phase diagram of superconducting Al films that were deposited over random Pb cluster arrays of varying density. These phase diagrams give a direct measure of the Fermi liquid spin renormalization, as well as the spin orbit scattering rate. We find that the spin renormalization is a linear function of the average Pb cluster -to- cluster separation and that this dependency can be used to tune the spin susceptibility of the Al over a surprisingly wide range from 0.8χ0to 4.0χ0, whereχ0is the non-interacting Pauli susceptibility.

    more » « less