Nonadiabatic Transition State Theory and Trajectory Surface Hopping Dynamics: Intersystem Crossing Between 3B1 and 1A1 States of SiH2

Ab initio molecular dynamics study of intersystem crossing dynamics for MH2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces

Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.

Spin-Dependent Stereochemistry: A Nonadiabatic Quantum Dynamics Case Study of S + H2 → SH + H Reaction

We study the spin-dependent stereodynamics of the S + H2 → SH + H reaction by using full-dimensional quantum dynamics calculations with zero total nuclear angular momentum along the triplet 3A″ states and singlet 1A' states. We find that the interplay between the electronic spin direction and the molecular geometry has a measurable influence on the singlet-triplet intersystem crossing reaction probabilities. Our results show that for some incident scattering angles in the body-fixed frame, the relative difference in intersystem crossing reaction probabilities (as determined between spin up and spin down initial states) can be as large as 15%. Our findings are an ab initio demonstration of spin-dependent nonadiabatic dynamics, which we hope will shine light as far as understanding the chiral-induced spin selectivity effect.

Far-UV absorption spectra of SiH2 and dibridged Si2H2 isolated in solid argon

We employ electron bombardment during the deposition of an Ar matrix containing a small proportion of SiH4 to generate various silicon hydrides. Subsequently, the irradiation of a matrix sample at 365 nm decomposes SiH2 and dibridged Si2H2 in solid Ar, which we identify through infrared spectroscopy. We further recorded the corresponding ultraviolet absorption spectra at each experimental stage. An intense band observed in the range of 170-203 nm is largely destroyed upon 365-nm photolysis, which is assigned to the C1B2 ← X1A1 transition of SiH2. Moreover, a moderate band observed in the region of 217-236 nm is reduced slightly, which is assigned to the 31B2 ← X1A1 transition of dibridged Si2H2. These assignments are made based on the observed photolytic behavior, and the prediction of the vertical excitation energies with the corresponding oscillator strengths by using time-dependent density functional theory and equation-of-motion coupled cluster theory.

Two state reactivity (TSR) and hydrogen tunneling reaction kinetics measured in the Co+ mediated decomposition of CH3CHO

The electronic structure of a transition metal atom allows it to act as a catalytic active site by providing lower energy alternative pathways in chemical transformations. We have identified and kinetically characterized three such pathways in the title reaction. One is an adiabatic pathway that occurs on a single potential energy surface described within the Born-Oppenheimer approximation. A second pathway opens microseconds into the reaction as a portion of the reacting population competitively transitions from triplet to singlet multiplicity to circumvent energetic barriers on the triplet surface. These pathways are single- and two-state reactive (SSR and TSR) where the Co+ cation mediates an oxidative addition/reductive elimination sequence of the CH3CHO molecule. The third observed reaction pathway is the aldehyde hydrogen tunneling through an Eyring barrier to form high-spin products. First-order rate constants for the adiabatic and nonadiabatic energy lowered pathways, and the hydrogen tunneling pathway, are each measured using the single photon initiated dissociative rearrangement reaction (SPIDRR) experimental technique. We believe that this is the first experimental study where such disparate dynamic features (SSR, TSR, and H-tunneling) are disentangled in a system's chemistry, attributing specific rate constant values to each effect and quantifying the various competitions. Moreover, multi-reference CASSCF/CASPT2 calculations indicate that structures with covalent Co-H bonds are present exclusively along the excited singlet surface. This phenomenon significantly reduces these structures' energy relative to their triplet counterparts, thus enabling the surface crossing and spin inversion that cause the observed two-state reactivity.

New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States

It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.

Predicting Kinetics and Dynamics of Spin-Dependent Processes

ConspectusPredicting mechanisms and rates of nonadiabatic spin-dependent processes including photoinduced intersystem crossings, thermally activated spin-forbidden reactions, and spin crossovers in metal centers is a very active field of research. These processes play critical roles in transition-metal-based and metalloenzymatic catalysis, molecular magnets, light-harvesting materials, organic light-emitting diodes, photosensitizers for photodynamic therapy, and many other applications. Therefore, accurate modeling of spin-dependent processes in complex systems and on different time scales is important for many problems in chemistry, biochemistry, and materials sciences.Nonadiabatic statistical theory (NAST) and nonadiabatic molecular dynamics (NAMD) are two complementary approaches to modeling the kinetics and dynamics of spin-dependent processes. NAST predicts the probabilities and rate constants of nonradiative transitions between electronic states with different spin multiplicities using molecular properties at only few critical points on the potential energy surfaces (PESs), including the reactant minimum and the minimum energy crossing point (MECP) between two spin states. This makes it possible to obtain molecular properties for NAST calculations using accurate but often computationally expensive electronic structure methods, which is critical for predicting the rate constants of spin-dependent processes. Alternatively, NAST can be used to study spin-dependent processes in very large complex molecular systems using less computationally expensive electronic structure methods. The nuclear quantum effects, such as zero-point vibrational energy, tunneling, and interference between reaction paths can be easily incorporated. However, the statistical and local nature of NAST makes it more suitable for large systems and slow kinetics. In contrast, NAMD explores entire PESs of interacting electronic states, making it ideal for modeling fast barrierless spin-dependent processes. Because the knowledge of large portions of PESs is often needed, the simulations require a very large number of electronic structure calculations, which limits the NAMD applicability to relatively small molecular systems and ultrafast kinetics.In this Account, we discuss our contribution to the development of the NAST and NAMD approaches for predicting the rates and mechanism of spin-dependent processes. First, we briefly describe our NAST and NAMD implementations. The NAST implementation is an extension of the transition state theory to the processes involving two crossing potential energy surfaces of different spin multiplicities. The NAMD approach includes the trajectory surface hopping (TSH) and ab initio multiple spawning (AIMS) methods. Second, we discuss several applications of NAST and NAMD to model spin-dependent processes in different systems. The NAST applicability to large complex systems is demonstrated by the studies of the spin-forbidden isomerization of the active sites of metal-sulfur proteins. Our implementation of the MECP search algorithm within the fully ab initio fragment molecular orbital method allows applying NAST to systems with thousands of atoms, such as the solvated protein rubredoxin. Applications of NAMD to ultrafast spin-dependent processes are represented by the generalized AIMS simulations utilizing the fast GPU-based TeraChem electronic structure program to gain insight into the complex photoexcited state relaxation in 2-cyclopentenone.

Gas-phase formation of silicon monoxide via non-adiabatic reaction dynamics and its role as a building block of interstellar silicates

Silicon monoxide (SiO) is classified as a key precursor and fundamental molecular building block to interstellar silicate nanoparticles, which play an essential role in the synthesis of molecular building blocks connected to the Origins of Life. In the cold interstellar medium, silicon monoxide is of critical importance in initiating a series of elementary chemical reactions leading to larger silicon oxides and eventually to silicates. To date, the fundamental formation mechanisms and chemical dynamics leading to gas phase silicon monoxide have remained largely elusive. Here, through a concerted effort between crossed molecular beam experiments and electronic structure calculations, it is revealed that instead of forming highly-stable silicon dioxide (SiO₂), silicon monoxide can be formed via a barrierless, exoergic, single-collision event between ground state molecular oxygen and atomic silicon involving non-adiabatic reaction dynamics through various intersystem crossings. Our research affords persuasive evidence for a likely source of highly rovibrationally excited silicon monoxide in cold molecular clouds thus initiating the complex chain of exoergic reactions leading ultimately to a population of silicates at low temperatures in our Galaxy.

Incorporating Berry Force Effects into the Fewest Switches Surface-Hopping Algorithm: Intersystem Crossing and the Case of Electronic Degeneracy

We present a preliminary surface-hopping approach for modeling intersystem crossing (ISC) dynamics between four electronic states: one singlet and one (triply degenerate) triplet. In order to incorporate all Berry force effects, the algorithm requires that, when moving along an adiabatic surface associated with the triplet manifold, one must also keep track of a quasi-diabatic index (akin to a "m_s" quantum number) for each trajectory. For a simple model problem, we find that a great deal of new physics can be captured by our algorithm, setting the stage for larger, more realistic (or perhaps even ab initio) simulations in the future.

Molecular Main Group Metal Hydrides

This review serves to document advances in the synthesis, versatile bonding, and reactivity of molecular main group metal hydrides within Groups 1, 2, and 12-16. Particular attention will be given to the emerging use of said hydrides in the rapidly expanding field of Main Group element-mediated catalysis. While this review is comprehensive in nature, focus will be given to research appearing in the open literature since 2001.

Selective Vacuum Evaporation by the Control of the Chemistry of Gas Phase in Vacuum Refining of Si

The evaporation of P from liquid Si under vacuum and reduced pressures of H₂, He, and Ar was studied to evaluate the feasibility of effective P removal with insignificant Si loss. It was found that the introduction of Ar and He inert gases at low pressures reduces the rate of P removal, and their pressure decrease will increase the process rate. Moreover, the kinetics of P removal was higher in He than in Ar, with simultaneous lower Si loss. Under reduced pressures of H₂ gas, however, the P removal rate was higher than that under vacuum conditions with the lowest Si loss. Quantum chemistry and dynamics simulations were applied, and the results indicated that P can maintain its momentum for longer distances in H₂ once it is evaporated from the melt surface and then can travel far away from the surface, while Si atoms lose their momentum in closer distances, yielding less net Si flux to the gas phase. Moreover, this distance is significantly increased with decreasing pressure for H₂, He, and Ar gases; however, it is the largest for H₂ and the lowest for Ar for a given pressure, while the temperature effect is insignificant. The rate of P evaporation was accelerated by applying an additional vacuum tube close to the melt surface for taking out the hot gas particles before they lose their temperature and velocity. It was shown that this technique contributes to the rate of process by preventing condensing gas stream back to the melt surface.

The non-adiabatic nanoreactor: towards the automated discovery of photochemistry

The ab initio nanoreactor has previously been introduced to automate reaction discovery for ground state chemistry. In this work, we present the nonadiabatic nanoreactor, an analogous framework for excited state reaction discovery. We automate the study of nonadiabatic decay mechanisms of molecules by probing the intersection seam between adiabatic electronic states with hyper-real metadynamics, sampling the branching plane for relevant conical intersections, and performing seam-constrained path searches. We illustrate the effectiveness of the nonadiabatic nanoreactor by applying it to benzene, a molecule with rich photochemistry and a wide array of photochemical products. Our study confirms the existence of several types of S₀/S₁ and S₁/S₂ conical intersections which mediate access to a variety of ground state stationary points. We elucidate the connections between conical intersection energy/topography and the resulting photoproduct distribution, which changes smoothly along seam space segments. The exploration is performed with minimal user input, and the protocol requires no previous knowledge of the photochemical behavior of a target molecule. We demonstrate that the nonadiabatic nanoreactor is a valuable tool for the automated exploration of photochemical reactions and their mechanisms.

Nonadiabatic Kinetics in the Intermediate Coupling Regime: Comparing Molecular Dynamics to an Energy-Grained Master Equation

We propose and test an extension of the energy-grained master equation (EGME) for treating nonadiabatic (NA) hopping between different potential energy surfaces, which enables us to model the competition between stepwise collisional relaxation and kinetic processes which transfer population between different electronic states of the same spin symmetry. By incorporating Zhu-Nakamura theory into the EGME, we are able to treat NA passages beyond the simple Landau-Zener approximation, along with the corresponding treatments of zero-point energy and tunneling probability. To evaluate the performance of this NA-EGME approach, we carried out detailed studies of the UV photodynamics of the volatile organic compound C₆-hydroperoxy aldehyde (C₆-HPALD) using on-the-fly ab initio molecular dynamics and trajectory surface hopping. For this multichromophore molecule, we show that the EGME is able to capture important aspects of the dynamics, including kinetic timescales, and diabatic trapping. Such an approach provides a promising and efficient strategy for treating the long-time dynamics of photoexcited molecules in regimes which are difficult to capture using atomistic on-the-fly molecular dynamics.

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.

Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence

We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet-triplet spin-orbit couplings quite well for CH₂S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.

Spin-inversion mechanisms in O2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.

Spin controlled surface chemistry: alkyl desorption from Si(100)-2×1 by nonadiabatic hydrogen elimination

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.

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 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 (i.e. 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.

Reduced-Dimensionality Quantum Dynamics Study of the 3Fe(CO)4 + H2 1FeH2(CO)4 Spin-inversion Reaction

Many chemical reactions of transition metal compounds involve a change in spin state via spin inversion, which is induced by relativistic spin-orbit coupling. In this work, we theoretically study the efficiency of a typical spin-inversion reaction, ³Fe(CO)₄ + H₂ → ¹FeH₂(CO)₄. Structural and vibrational information on the spin-inversion point, obtained through the spin-coupled Hamiltonian approach, is used to construct three degree-of-freedom potential energy surfaces and to obtain singlet-triplet spin-orbit couplings. Using the developed spin-diabatic potential energy surfaces in reduced dimensions, we perform quantum nonadiabatic transition state wave packet calculations to obtain the cumulative reaction probability. The calculated cumulative reaction probability is found to be significantly larger than that estimated from the one-dimensional surface-hopping probability. This indicates the importance of both multidimensional and nuclear quantum effects in spin inversion for polyatomic chemical reaction systems.

A Robust and Unified Solution for Choosing the Phases of Adiabatic States as a Function of Geometry: Extending Parallel Transport Concepts to the Cases of Trivial and Near-Trivial Crossings

We investigate a simple and robust scheme for choosing the phases of adiabatic electronic states smoothly (as a function of geometry) so as to maximize the performance of ab initio non-adiabatic dynamics methods. Our approach is based upon consideration of the overlap matrix (U) between basis functions at successive points in time and selecting the phases so as to minimize the matrix norm of log(U). In so doing, one can extend the concept of parallel transport to cases with sharp curve crossings. We demonstrate that this algorithm performs well under extreme situations where dozens of states cross each other either through trivial crossings (where there is zero effective diabatic coupling), or through non-trivial crossings (when there is a non-zero diabatic coupling), or through a combination of both. In all cases, we compute the time-derivative coupling matrix elements (or equivalently non-adiabatic derivative coupling matrix elements) that are as smooth as possible. Our results should be of interest to all who are interested in either non-adiabatic dynamics, or more generally, parallel transport in large systems.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Nonadiabatic dynamics simulations on internal conversion and intersystem crossing processes in gold(i) compounds

The position at which the second gold(i)-phosphine group is attached was experimentally found to play a noticeable role in intersystem crossing rates of gold(i) naphthalene derivatives. However, the physical origin is ambiguous. Herein we have employed generalized trajectory-based surface-hopping dynamics simulations to simulate the excited-state relaxation dynamics of these gold(i) naphthalene compounds including both the intersystem crossing process from the initially populated first excited singlet states S₁ to triplet manifolds and internal conversion processes within these triplet states. Our predicted intersystem crossing rates are consistent with experiments very well. On the basis of the present results, we have found that (1) ultrafast and subpicosecond intersystem crossing processes are mainly caused by small energy gaps and large spin-orbit couplings between S₁ and T_n; (2) adding the second gold(i)-phosphine group does not increase spin-orbit couplings between S₁ and T_n but decrease their values remarkably, which implies that heavy-atom effects are state-specific, not state-universal; (3) the position at which the second gold(i)-phosphine group is attached has a remarkable influence on the electronic structures of S₁ and T_n and their relative energies, which affect energy gaps and spin-orbit couplings between S₁ and T_n and eventually modulate intersystem crossing rates from S₁ to T_n. These new insights are very useful for the design of gold-containing compounds with excellent photoluminescence properties. Finally, this work also exemplifies that different isomers of a compound could have distinct excited-state relaxation dynamics.