Ab initio surface hopping excited-state molecular dynamics approach on the basis of spin-orbit coupled states: An application to the A-band photodissociation of CH3 I

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-Vibronic Coupling Controls the Intersystem Crossing of Iodine-Substituted BODIPY Triplet Chromophores

4,4-Difluoro-4-borata-3a-azonia-4a-aza-s-indacene (BODIPY) dyes are extensively used in various applications of their triplet states, ranging from photoredox catalysis, through triplet sensitization to photodynamic therapy. However, the rational design of BODIPY triplet chromophores by ab initio modelling is limited by their strong interactions of spin, electronic and vibrational dynamics. In particular, spin-vibronic coupling is often overlooked when estimating intersystem crossing (ISC) rates. In this study, a combined experimental and theoretical approach using spin-vibronic coupling to correctly describe ISC in BODIPY dyes was developed. For this purpose, seven π-extended BODIPY derivatives with iodine atoms in different positions were examined. It was found that the heavy-atom effect of iodine atoms is site specific, causing high triplet yields in only some positions. This site-specific ISC was explained by El-Sayed rules, so both the contribution and character of the molecular orbitals involved in the excitation must be considered when predicting the ISC rates. Overall, the rational design of BODIPY triplet chromophores requires using (i) the high-quality electronic structure theory, including both static and dynamical correlations; and (ii) the two-component wave function Hamiltonian, and rationalizing; and (iii) ISC based on the character of the molecular orbitals of heavy atoms involved in the excitation, expanding El-Sayed rules beyond their traditional applications.

Multi-state Energy Landscape for Photoreaction of Stilbene and Dimethyl-stilbene

We have recently developed the reaction space projector (ReSPer) method, which constructs a reduced-dimensionality reaction space uniquely determined from reference reaction paths for a polyatomic molecular system and projects classical trajectories into the same reaction space. In this paper, we extend ReSPer to the analysis of photoreaction dynamics and relaxation processes of stilbene and present the concept of a "multi-state energy landscape," incorporating the ground- and excited-state reaction subspaces. The multi-state energy landscape successfully explains the previously established photoreaction processes of cis-stilbene, such as the cis-trans photoisomerization and photocyclization. In addition, we discuss the difference in the excited-state reaction dynamics between stilbene and 1,1'-dimethyl stilbene based on a common reaction subspace determined from the framework part of reference structures with different number of atoms. This approach allows us to target any molecule with a common framework, greatly expanding the applicability of the ReSPer analysis. The multi-state energy landscape provides fruitful insight into photochemical reactions, exploring the excited- and ground-state potential energy surfaces, as well as comprehensive reaction processes with nonradiative transitions between adiabatic states, within the stage of a reduced-dimensionality reaction space.

Mapping wave packet bifurcation at a conical intersection in CH3I by attosecond XUV transient absorption spectroscopy

Extreme ultraviolet (XUV) transient absorption spectroscopy has emerged as a sensitive tool for mapping the real-time structural and electronic evolution of molecules. Here, attosecond XUV transient absorption is used to track dynamics in the A-band of methyl iodide (CH₃I). Gaseous CH₃I molecules are excited to the A-band by a UV pump (277 nm, ∼20 fs) and probed by attosecond XUV pulses targeting iodine I(4d) core-to-valence transitions. Owing to the excellent temporal resolution of the technique, passage through a conical intersection is mapped through spectral signatures of nonadiabatic wave packet bifurcation observed to occur at 15 ± 4 fs following UV photoexcitation. The observed XUV signatures and time dynamics are in agreement with previous simulations [H. Wang, M. Odelius, and D. Prendergast, J. Chem. Phys. 151, 124106 (2019)]. Due to the short duration of the UV pump pulse, coherent vibrational motion in the CH₃I ground state along the C-I stretch mode (538 ± 7 cm^-1) launched by resonant impulsive stimulated Raman scattering and dynamics in multiphoton excited states of CH₃I are also detected.

Spectroscopic Signature of Chemical Bond Dissociation Revealed by Calculated Core-Electron Spectra

The advent of ultrashort soft X-ray pulse sources permits the use of established gas-phase spectroscopy methods to investigate ultrafast photochemistry in isolated molecules with element and site specificity. In the present study, we simulate excited-state wavepacket dynamics of a prototypical process, the ultrafast photodissociation of methyl iodide. Using the simulation, we calculate time-dependent excited-state carbon edge photoelectron and Auger electron spectra. We observe distinct signatures in both types of spectra and show their direct connection to C-I bond dissociation and charge rearrangement processes in the molecule. We demonstrate at the CH₃I molecule that the observed signatures allow us to map the time-dependent dynamics of ultrafast photoinduced bond breaking with unprecedented detail.

A combined multi-reference pump-probe simulation method with application to XUV signatures of ultrafast methyl iodide photodissociation

UV pump-XUV/X-ray probe measurements have been successfully applied in the study of photo-induced chemical reactions. Although rich element-specific electronic structure information is accessible within XUV/X-ray (inner-shell) absorption spectra, it can be difficult to interpret the chemistry directly from the spectrum without supporting theoretical simulations. A multireference method to completely simulate UV pump-XUV/X-ray probe measurement has been developed and applied to study the methyl iodide photodissociation process. Multireference, fewest-switches surface hopping (FSSH) trajectories were used to explore the coupled electronic and ionic dynamics upon photoexcitation of methyl iodide. Interpretation of previous measurements is provided by associated multireference, restricted active space, inner-shell spectral simulations. This combination of multireference FSSH trajectories and XUV spectra provides an interpretation of transient features appearing in previous measurements within the first 100 fs after photoexcitation and validates the significant branching ratio in the final excited-state population. This methodology should prove useful for interpretation of the increasing number of inner-shell probe studies of molecular excited states or for directing new experiments toward interesting regions of the potential energy landscape.

Low-Lying Excited States of hqxcH and Zn-hqxc Complex: Toward Understanding Intramolecular Proton Transfer Emission

Excited state intramolecular proton transfer (ESIPT) has been a topic of interest due to its potential to lead to multiple emissions. Although many organic molecules showing ESIPT emission are already known, studies on metal complexes showing ESIPT and their related theoretical understandings are very limited. In this study, we focus on [Zn(hqxc)₂(DMSO)₂] (Zn-hqxc: hqxc = 3-hydroxy-2-quinoxalinecarboxylate, DMSO = dimethyl sulfoxide), which shows ESIPT emission in the solid state, even though the hqxcH ligand does not show ESIPT emission. To gain insights into the role of the zinc atom and the emission mechanisms, we examined excited states of free hqxcH and the Zn-hqxc complex using time-dependent density functional theory calculations. From the results, it was shown that the zinc atom triggers a structural change of the hqxcH ligand from the lactam form (3,4-dihydro-3-oxo-2-quinoxalinecarboxylic acid) to the enol form (3-hydroxy-2-quinoxalinecarboxylic acid), where the latter form has several stable excited states. Several stable geometries were found for singlet and triplet excited states, suggesting that emissions for the Zn-hqxc complex can be both phosphorescence and fluorescence caused by the enol-enol, keto-keto, and keto-enol forms of the two hqcx ligands in the complex. We found that the photophysical properties of the Zn-hqxc complex are dominated by the ligand due to the filled d¹⁰ of Zn(II). The presented results suggest that, to design new ESIPT metal complexes, one possible approach is to combine a metal atom showing ligand centered emission and a ligand that has separate ESIPT and coordination sites.