Attosecond Monitoring of Nonadiabatic Molecular Dynamics by Transient X-ray Transmission Spectroscopy

Tracing the evolution of molecular coherences can provide a direct, unambiguous probe of nonadiabatic molecular processes, such as the passage through conical intersections of electronic states. Two techniques, attosecond transient absorption spectroscopy (ATAS) and Transient Redistribution of Ultrafast Electronic Coherences in Attosecond Raman Signals (TRUECARS), have been used or proposed for monitoring nonadiabatic molecular dynamics. Both techniques employ the transmission of a weak attosecond extreme-ultraviolet or X-ray probe to interrogate the molecule at controllable time delays with respect to an optical pump, thereby extracting dynamical information from transient spectral features. The connection between these techniques has not been firmly established yet. In this theoretical study, we provide a unified description of both transient transmission techniques, establishing their relationship as limits of the same pump-probe spectroscopy technique for different pulse parameter regimes. We demonstrate this by quantum dynamical simulations of thiophenol photodissociation and show how complementary coherence information can be revealed by the two techniques.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Femtosecond Wavepacket Dynamics Reveals the Molecular Structures in the Excited (S1) and Cationic (D0) States

Molecular structures in the electronically excited (S₁) and cationic (D₀) states of 2-fluorothioanisole (2-FTA) have been precisely refined from the real-time dynamics of the femtosecond (fs) wavepacket prepared by the coherent excitation of the Franck-Condon active out-of-plane torsional modes in the S₁ ← S₀ transition at 285 nm. The simulation to reproduce the experiment in terms of the beating frequencies gives the nonplanar geometry of 2-FTA in S₁, where the out-of-plane dihedral angle (φ) of the S-CH₃ moiety is 51° with respect to the molecular plane. The behavior of the fs wavepacket in terms of the amplitudes and phases with the change of the probe (ionization) wavelength (λ_probe = 300-330 nm) provides the otherwise veiled structure of the cationic D₀ state. While the 2-FTA cation adopts the planar geometry (φ = 0°) at the global minimum, it is found to have a vertical minimum at φ ≈ 135° from the perspective of the D₀ ← S₁ vertical transition. Ab initio calculations support the experiment quite well although the simulation using the model potentials could improve the match with the experiment, giving the new interpretation for the previously disputed photoelectron spectroscopic results.

Wavepacket propagations for the early time dynamics of proton-coupled electron transfer in the charge-transfer state of NH3Cl complex

A charge-transfer (CT) excited state of NH₃Cl, generated by photo-detachment of an electron from the anionic NH₃Cl^- precursor, can be represented as H₂N⁺-H-Cl^- and proceeds to two chemical reactions: one reaction generating NH₂ and HCl resulting from a proton transfer (PT) and the other reaction producing NH₃ and a Cl atom resulting from an electron transfer (ET); both are coupled to form a typical proton-coupled electron transfer (PCET) process. The early time dynamics of this CT were studied using time-dependent wavepacket propagation on three nonadiabatically coupled electronic states in a reduced three-dimensional space. The electronic states were treated using the XMS-CASPT2/aug-cc-pVTZ ab initio methodology. The population dynamics of the three coupled electronic states were analyzed in detail to reveal the initial stage of the PCET process up to ∼100 fs, while the branching ratio, χ = PT/(ET+PT), was determined after wavepacket propagations of up to 2000 fs. Another main result is the dependence of χ on the vibration levels of the initial precursor anion and the isotope substitution of the connecting H atom with deuterium and tritium. Our study reveals the detailed microscopic features of the PCET process embedded in the CT state of the NH₃Cl complex and certain systematic dependences of the branching ratio χ on the above factors.

Real-Time Tunneling Dynamics through Adiabatic Potential Energy Surfaces Shaped by a Conical Intersection

Dynamic shaping of the adiabatic tunneling barrier in the S-H bond extension coordinate of several ortho-substituted thiophenols has been found to be mediated by low-frequency out-of-plane vibrational modes, which are parallel to the coupling vector of the branching plane comprising the conical intersection. The S-H predissociation tunneling rate (k) measured when exciting to the S₁ zero-point level of 2-methoxythiophenol (44 ps)^-1 increases abruptly, to k ≈ (22 ps)^-1, at the energy corresponding to excitation of the 152 cm^-1 out-of-plane vibrational mode and then falls back to k ≈ (40 ps)^-1 when the in-plane mode is excited at 282 cm^-1. Similar resonance-like peaks in plots of S₁ tunneling rate versus internal energy are observed when exciting the corresponding low-frequency out-of-plane modes in the S₁ states of 2-fluorothiophenol and 2-chlorothiophenol. This experiment provides clear-cut evidence for dynamical "shaping" of the lower-lying adiabatic potential energy surfaces by the higher-lying conical intersection seam, which dictates the multidimensional tunneling dynamics.

Practical approximation of the non-adiabatic coupling terms for same-symmetry interstate crossings by using adiabatic potential energies only

A very simple equation, F_ij^App=[(∂²(V_i^a-V_j^a)/∂Q²)/(V_i^a-V_j^a)]^1/2/2, giving a reliable magnitude of non-adiabatic coupling terms (NACTs, F_ij's) based on adiabatic potential energies only (V_i^a and V_j^a) was discovered, and its reliability was tested for several prototypes of same-symmetry interstate crossings in LiF, C₂, NH₃Cl, and C₆H₅SH molecules. Our theoretical derivation starts from the analysis of the relationship between the Lorentzian dependence of NACTs along a diabatization coordinate and the well-established linear vibronic coupling scheme. This analysis results in a very simple equation, α=2κ/Δ_c, enabling the evaluation of the Lorentz function α parameter in terms of the coupling constant κ and the energy gap Δ_c (Δ_c=|V_i^a-V_j^a|_Qc ) between adiabatic states at the crossing point Q_C. Subsequently, it was shown that Q_C corresponds to the point where F_ij^App exhibit maximum values if we set the coupling parameter as κ=[(V_i^a-V_j^a)⋅(∂²(V_i^a-V_j^a)/∂Q²)]_Qc^1/2/2. Finally, we conjectured that this relation could give reasonable values of NACTs not only at the crossing point but also at other geometries near Q_C. In this final approximation, the pre-defined crossing point Q_C is not required. The results of our test demonstrate that the approximation works much better than initially expected. The present new method does not depend on the selection of an ab initio method for adiabatic electronic states but is currently limited to local non-adiabatic regions where only two electronic states are dominantly involved within a nuclear degree of freedom.

Electronic spectrum and characterization of diabatic potential energy surfaces for thiophenol

We present an accurate simulation of the UV spectrum and a diabatization of three singlet potential surfaces along four coordinates.

Conformer specific nonadiabatic reaction dynamics in the photodissociation of partially deuterated thioanisoles (C6H5S-CH2D and C6H5S-CHD2)

In this work, we have investigated nonadiabatic dynamics in the vicinity of conical intersections for predissociation reactions of partially deuterated thioanisole molecules: C₆H₅S-CH₂D and C₆H₅S-CHD₂. Each isotopomer has two distinct rotational conformers according to the geometrical position of D or H of the methyl moiety with respect to the molecular plane for C₆H₅S-CH₂D or C₆H₅S-CHD₂, respectively, as spectroscopically characterized in our earlier report [J. Lee, S.-Y. Kim and S. K. Kim, J. Phys. Chem. A, 2014, 118, 1850]. Since identification and separation of two different rotational conformers of each isotopomer have been unambiguously done, we could interrogate nonadiabatic dynamics of thioanisole in terms of both H/D substitutional and conformational structural effects. Nonadiabatic transition probability, estimated by the experimentally measured branching ratio of the nonadiabatically produced ground-state channel giving C₆H₅S·(X[combining tilde]) versus the adiabatic excited-state channel leading to the C₆H₅S·(Ã) radical, shows resonance-like increases at symmetric (ν_s) or asymmetric (7a) S-CH₂D (or S-CHD₂) stretching mode excitation in S₁ for all conformational isomers of two isotopomers. However, absolute probabilistic value of the nonadiabatic transition is found to vary quite drastically depending on different conformers and isotopomers. The experimental finding that nonadiabatic transition dynamics are very sensitive to subtle changes in the nuclear configuration within the Franck-Condon region induced by the H/D substitution indicates that the S₁/S₂ conical intersection seam is quite narrowly defined in the multi-dimensional nuclear configurational space as far as the S-methyl predissociation reaction is concerned. In order to understand the relation between molecular structure and nonadiabaticity of reaction, potential energy surfaces near S₁/S₂ conical intersections have been theoretically calculated along ν_s and 7a normal mode coordinates for all conformational isomers. Slow-electron velocity map imaging (SEVI) spectroscopy is employed to unravel the extent of intramolecular vibrational redistribution (IVR) for particular mode excitations of S₁, providing insights into the dynamic interplay between IVR and nonadiabatic transition probability near the conical intersection seam.

Extreme population inversion in the fragments formed by UV photoinduced S-H bond fission in 2-thiophenethiol

H atom loss following near ultraviolet photoexcitation of gas phase 2-thiophenethiol molecules has been studied experimentally, by photofragment translational spectroscopy (PTS) methods, and computationally, by ab initio electronic structure calculations. The long wavelength (277.5 ≥ λ(phot) ≥ 240 nm) PTS data are consistent with S-H bond fission after population of the first (1)πσ* state. The partner thiophenethiyl (R) radicals are formed predominantly in their first excited Ã(2)A' state, but assignment of a weak signal attributable to H + R(X˜(2)A'') products allows determination of the S-H bond strength, D0 = 27,800 ± 100 cm(-1) and the Ã-X˜ state splitting in the thiophenethiyl radical (ΔE = 3580 ± 100 cm(-1)). The deduced population inversion between the à and X˜ states of the radical reflects the non-planar ground state geometry (wherein the S-H bond is directed near orthogonal to the ring plane) which, post-photoexcitation, is unable to planarise sufficiently prior to bond fission. This dictates that the dissociating molecules follow the adiabatic fragmentation pathway to electronically excited radical products. π* ← π absorption dominates at shorter excitation wavelengths. Coupling to the same (1)πσ* potential energy surface (PES) remains the dominant dissociation route, but a minor yield of H atoms attributable to a rival fragmentation pathway is identified. These products are deduced to arise via unimolecular decay following internal conversion to the ground (S0) state PES via a conical intersection accessed by intra-ring C-S bond extension. The measured translational energy disposal shows a more striking change once λ(phot) ≤ 220 nm. Once again, however, the dominant decay pathway is deduced to be S-H bond fission following coupling to the (1)πσ* PES but, in this case, many of the evolving molecules are deduced to have sufficiently near-planar geometries to allow passage through the conical intersection at extended S-H bond lengths and dissociation to ground (X˜) state radical products. The present data provide no definitive evidence that complete ring opening can compete with fast S-H bond fission following near UV photoexcitation of 2-thiophenethiol.

Dynamic Symmetry Breaking Hidden in Fano Resonance of a Molecule: S1 State of Diazirine Using Quantum Wave Packet Propagation

Fano resonance in the predissociation of the S1 state of diazirine was studied by applying a time-dependent wave packet propagation method, and dynamic symmetry breaking (DSB) around the stationary structure of S1 was disclosed in a detailed analysis of this theoretical result. The DSB was found to originate in coupling between the asymmetric C-N2 stretching and CH2 wagging modes, suggesting that there is a slight time gap between ring opening and the concurrent dragging of two H atoms of the CH2 moiety. Although the depth of the double well due to DSB is just 0.011 eV, its presence noticeably affects the early time dynamics and observed spectrum.