skip to main content


Title: 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.  more » « less
Award ID(s):
1654547
NSF-PAR ID:
10207950
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
29
ISSN:
1463-9076
Page Range / eLocation ID:
16641 to 16647
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Energy transfer in organic materials is extensively studied due to many applications in optoelectronics. The electronic and vibrational relaxations within molecular assemblies can be influenced by stacking arrangements or additions of a backbone that unites them. Here, we present the computational study of the photoexcitation dynamics of a perylene diimide monomer, and face-to-face stacked dimer and trimer. By using non-adiabatic excited-state molecular dynamics simulations, we show that the non-radiative relaxation is accelerated with the number of stacked molecules. This effect is explained by differences in the energy splitting between states that impacts their corresponding nonadiabatic couplings. Additionally, our analysis of the vibronic dynamics reveals that the passage through the different conical intersections that participate in the relaxation of the stacked systems, activate a positive feedback mechanism. This effect involves a narrow set of vibrational normal modes that accelerate the process by increasing the efficiency of its vibronic dynamics. In contrast, an addition of a biologically inspired backbone slows down the relaxation rate due to its participation in the vibronic dynamics of the molecular stacking arrangements. Our results suggest the stacking arrangements and common backbones as strategies to modulate the efficiency of electronic and vibrational relaxation of diimide-based systems and other molecular aggregates. 
    more » « less
  2. 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
  3. Intermolecular interactions between radicals and closed-shell molecules are ubiquitous in chemical processes, ranging from the benchtop to the atmosphere and extraterrestrial space. While energy decomposition analysis (EDA) schemes for closed-shell molecules can be generalized for studying radical–molecule interactions, they face challenges arising from the unique characteristics of the electronic structure of open-shell species. In this work, we introduce additional steps that are necessary for the proper treatment of radical–molecule interactions to our previously developed unrestricted Absolutely Localized Molecular Orbital (uALMO)-EDA based on density functional theory calculations. A “polarize-then-depolarize” (PtD) scheme is used to remove arbitrariness in the definition of the frozen wavefunction, rendering the ALMO-EDA results independent of the orientation of the unpaired electron obtained from isolated fragment calculations. The contribution of radical rehybridization to polarization energies is evaluated. It is also valuable to monitor the wavefunction stability of each intermediate state, as well as their associated spin density profiles, to ensure the EDA results correspond to a desired electronic state. These radical extensions are incorporated into the “vertical” and “adiabatic” variants of uALMO-EDA for studies of energy changes and property shifts upon complexation. The EDA is validated on two model complexes, H 2 O⋯˙F and FH⋯˙OH. It is then applied to several chemically interesting radical–molecule complexes, including the sandwiched and T-shaped benzene dimer radical cation, complexes of pyridine with benzene and naphthalene radical cations, binary and ternary complexes of the hydroxyl radical with water (˙OH(H 2 O) and ˙OH(H 2 O) 2 ), and the pre-reactive complexes and transition states in the ˙OH + HCHO and ˙OH + CH 3 CHO reactions. These examples suggest that this second generation uALMO-EDA is a useful tool for furthering one's understanding of both energetic and property changes associated with radical–molecule interactions. 
    more » « less
  4. Abstract

    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
  5. Hydrogen tunneling plays a critical role in many biologically and chemically important processes. The nuclear–electronic orbital multistate density functional theory (NEO-MSDFT) method was developed to describe hydrogen transfer systems. In this approach, the transferring proton is treated quantum mechanically on the same level as the electrons within multicomponent DFT, and a nonorthogonal configuration interaction scheme is used to produce delocalized vibronic states from localized vibronic states. The NEO-MSDFT method has been shown to provide accurate hydrogen tunneling splittings for fixed molecular systems. Herein, the NEO-MSDFT analytical gradients for both ground and excited vibronic states are derived and implemented. The analytical gradients and semi-numerical Hessians are used to optimize and characterize equilibrium and transition state geometries and to generate minimum energy paths (MEPs), for proton transfer in the deprotonated acetylene dimer and malonaldehyde. The barriers along the resulting MEPs are lower when the transferring proton is quantized because the NEO-MSDFT method inherently includes the zero-point energy of the transferring proton. Analysis of the proton densities along the MEPs illustrates that the proton density can exhibit symmetric or asymmetric bilobal character associated with symmetric or slightly asymmetric double-well potential energy surfaces and hydrogen tunneling. Analysis of the contributions to the intrinsic reaction coordinate reveals that changes in the C–O bond lengths drive proton transfer in malonaldehyde. This work provides the foundation for future reaction path studies and direct nonadiabatic dynamics simulations of a wide range of hydrogen transfer reactions.

     
    more » « less