Structure and dynamic role of conical intersections in the πσ*-mediated photodissociation reactions

Dynamic Role of the Intramolecular Hydrogen Bonding in the S1 State Relaxation Dynamics Revealed by the Direct Measurement of the Mode-Dependent Internal Conversion Rate of 2-Chlorophenol and 2-Chlorothiophenol

The dynamic role of the intramolecular hydrogen bond in the S1 relaxation of cis-2-chlorophenol (2-CP) or cis-2-chlorothiophenol (2-CTP) has been investigated in a state-specific manner. Whereas ultrafast internal conversion is dominant for 2-CP, the H-tunneling competes with internal conversion for 2-CTP even at the S1 origin. The S0-S1 internal conversion rate of 2-CTP could be directly measured from the S1 lifetimes of 2-CTP-d1 (Cl-C6H4-SD) as the D-tunneling is kinetically blocked, allowing distinct estimations of tunneling and internal conversion rates with increasing the energy. The internal conversion rate of 2-CTP increases by two times at the out-of-plane torsional mode excitation, suggesting that the internal conversion is facilitated at the nonplanar geometry. It then sharply increases at ∼600 cm-1, indicating that the S1/S0 conical intersection is readily accessible at the extended C-Cl bond length. The strength of the intramolecular hydrogen bond should be responsible for the distinct dynamic behaviors of 2-CP and 2-CTP.

Mode-dependent H atom tunneling dynamics of the S1 phenol is resolved by the simple topographic view of the potential energy surfaces along the conical intersection seam

Mode-dependent H atom tunneling dynamics of the O-H bond predissociation of the S₁ phenol has been theoretically analyzed. As the tunneling is governed by the complicated multi-dimensional potential energy surfaces that are dynamically shaped by the upper-lying S₁(ππ*)/S₂(πσ*) conical intersection, the mode-specific tunneling dynamics of phenol (S₁) has been quite formidable to be understood. Herein, we have examined the topography of the potential energy surface along the particular S₁ vibrational mode of interest at the nuclear configurations of the S₁ minimum and S₁/S₂ conical intersection. The effective adiabatic tunneling barrier experienced by the reactive flux at the particular S₁ vibrational mode excitation is then uniquely determined by the topographic shape of the potential energy surface extended along the conical intersection seam coordinate associated with the particular vibrational mode. The resultant multi-dimensional coupling of the specific vibrational mode to the tunneling coordinate is then reflected in the mode-dependent tunneling rate as well as nonadiabatic transition probability. Remarkably, the mode-specific experimental result of the S₁ phenol tunneling reaction [K. C. Woo and S. K. Kim, J. Phys. Chem. A 123, 1529-1537 (2019)] (in terms of the qualitative and relative mode-dependent dynamic behavior) could be well rationalized by semi-classical calculations based on the mode-specific topography of the effective tunneling barrier, providing the clear conceptual insight that the skewed potential energy surfaces along the conical intersection seam (strongly or weakly coupled to the tunneling reaction coordinate) may dictate the tunneling dynamics in the proximity of the conical intersection.

πσ*-Mediated Nonadiabatic Tunneling Dynamics of Thiophenols in S1: The Semiclassical Approaches

The S-H bond tunneling predissociation dynamics of thiophenol and its ortho-substituted derivatives (2-fluorothiophenol, 2-methoxythiophenol, and 2-chlorothiphenol) in S₁ (ππ*) where the H atom tunneling is mediated by the nearby S₂ (πσ*) state (which is repulsive along the S-H bond extension coordinate) have been investigated in a state-specific way using the picosecond time-resolved pump-probe spectroscopy for the jet-cooled molecules. The effects of the specific vibrational mode excitations and the SH/SD substitutions on the S-H(D) bond rupture tunneling dynamics have been interrogated, giving deep insights into the multidimensional aspects of the S₁/S₂ conical intersection, which also shapes the underlying adiabatic tunneling potential energy surfaces (PESs). The semiclassical tunneling rate calculations based on the Wentzel-Kramers-Brillouin (WKB) approximation or Zhu-Nakamura (ZN) theory have been carried out based on the ab initio PESs calculated in the (one, two, or three) reduced dimensions to be compared with the experiment. Though the quantitative experimental results could not be reproduced satisfactorily by the present calculations, the qualitative trends among different molecules in terms of the behavior of the tunneling rate versus the (adiabatic) barrier height or the number of PES dimensions could be rationalized. Most interestingly, the H/D kinetic isotope effect observed in the tunneling rate could be much better explained by the ZN theory compared to the WKB approximation, indicating that the nonadiabatic coupling matrix elements should be invoked for understanding the tunneling dynamics taking place in the proximity of the conical intersection.

Photochemical carbon-sulfur bond cleavage in thioethers mediated via excited state Rydberg-to-valence evolution

Time-resolved photoelectron imaging and supporting ab initio quantum chemistry calculations were used to investigate non-adiabatic excess energy redistribution dynamics operating in the saturated thioethers diethylsulfide, tetrahydrothiophene and thietane. In all cases, 200 nm excitation leads to molecular fragmentation on an ultrafast (<100 fs) timescale, driven by the evolution of Rydberg-to-valence orbital character along the S-C stretching coordinate. The C-S-C bending angle was also found to be a key coordinate driving initial internal conversion through the excited state Rydberg manifold, although only small angular displacements away from the ground state equilibrium geometry are required. Conformational constraints imposed by the cyclic ring structures of tetrahydrothiophene and thietane do not therefore influence dynamical timescales to any significant extent. Through use of a high-intensity 267 nm probe, we were also able to detect the presence of some transient (bi)radical species. These are extremely short lived, but they appear to confirm the presence of two competing excited state fragmentation channels - one proceeding directly from the initially prepared 4p manifold, and one involving non-adiabatic population of the 4s state. This is in addition to a decay pathway leading back to the S₀ electronic ground state, which shows an enhanced propensity in the 5-membered ring system tetrahydrothiophene over the other two species investigated.

Non-Born-Oppenheimer effects in molecular photochemistry: an experimental perspective

Non-adiabatic couplings between Born-Oppenheimer (BO)-derived potential energy surfaces are now recognized as pivotal in describing the non-radiative decay of electronically excited molecules following photon absorption. This opinion piece illustrates how non-BO effects provide photostability to many biomolecules when exposed to ultraviolet radiation, yet in many other cases are key to facilitating 'reactive' outcomes like isomerization and bond fission. The examples are presented in order of decreasing molecular complexity, spanning studies of organic sunscreen molecules in solution, through two families of heteroatom containing aromatic molecules and culminating with studies of isolated gas phase H₂O molecules that afford some of the most detailed insights yet available into the cascade of non-adiabatic couplings that enable the evolution from photoexcited molecule to eventual products. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Multiphoton-excited dynamics of the trans or cis structural isomer of 1,2-dibromoethylene

Photofragmentation dynamics of cis and trans isomers of 1,2-dibromoethylene (1,2-DBE) have been investigated by multiphoton excitation using a picosecond (ps) laser pulse. It has been found that the Br₂ ⁺ product ion preferentially originates from the cis isomer rather than from trans. The Boltzmann-type isotropic low kinetic energy components of the Br⁺ and Br₂ ⁺ product state distributions seem to be most likely from the unimolecular reactions of the vibrationally hot cationic ground state generated by the three-photon absorption at the photon energy below ∼38 000 cm^-1. The highly anisotropic kinetic energy components of Br⁺ and Br₂ ⁺ start to appear at the photon energy above ∼38 000 cm^-1, where the D_n (n ≥ 1) - D₀ transition is facilitated within the same ps laser pulse as the parent molecule is efficiently ionized by the two-photon absorption. The transition dipole moment of the D₄ - D₀ transition of the strongest oscillator strength has been theoretically predicted to be parallel to the C-Br bond or C=C bond axis for the trans or cis isomer, respectively. The fast anisotropic with the (β ∼ +2) component in the Br⁺ product distribution is thus likely from the trans isomer, whereas that of Br₂ ⁺ (β ∼ -0.5) should be the consequence of the photodissociation of the cis isomer. The isomer-specific reactivity found here in the picosecond multiphoton excitation of 1,2-DBE provides a nice platform for the better understanding of the structure-reactivity relationship under the harsh condition of the strong or ultrashort optical field.

S1-State Decay Dynamics of Benzenediols (Catechol, Resorcinol, and Hydroquinone) and Their 1:1 Water Clusters

The S₁-state decaying rates of the three different benzenediols, catechol, resorcinol, and hydroquinone, and their 1:1 water clusters have been state-specifically measured using the picosecond time-resolved parent ion transients obtained by the pump (excitation) and probe (ionization) scheme. The S₁ lifetime of catechol is found to be short, giving τ ∼ 5.9 ps at the zero-point level. This is ascribed to the H-atom detachment from the free OH moiety of the molecule. Consistent with a previous report (J. Phys. Chem. Lett. 2013, 4, 3819-3823), the S₁ lifetime gets lengthened with low-frequency vibrational mode excitations, giving τ ∼ 9.0 ps for the 116 cm^-1 band. The S₁ lifetimes at the additional vibronic modes of catechol are newly measured, showing the nonnegligible mode-dependent fluctuations of the tunneling rate. When catechol is complexed with water, the S₁ lifetime is enormously increased to τ ∼ 1.80 ns at the zero-point level while it shows an unusual dip at the intermolecular stretching mode excitation (τ ∼ 1.03 ns at 146 cm^-1). Otherwise, it is shortened monotonically with increasing the internal energy, giving τ ∼ 0.67 ns for the 856 cm^-1 band. Two different asymmetric or symmetric conformers of resorcinol give the respective S₁ lifetimes of 4.5 or 6.3 ns at their zero-point levels according to the estimation from our transients taken within the temporal window of 0-2.7 ns. When resorcinol is 1:1 complexed with H₂O, the S₁ decaying rate is slightly accelerated for both conformers. The S₁ lifetimes of trans and cis forms of hydroquinone are measured to be more or less same, giving τ ∼ 2.8 ns at the zero-point level. When H₂O is complexed with hydroquinone, the S₁ decaying process is facilitated for both conformers, slightly more efficiently for the cis conformer.

Conformer-Specific Tunneling Dynamics Dictated by the Seam Coordinate of the Conical Intersection

The dynamic role of the conical intersection "seam" coordinate has been first revealed in the H fragmentation reaction of ortho(o)-cresol conformers. One of the (3N - 8) dimensional seam coordinates of the S₁(ππ*)/S₂(πσ*) conical intersection has been identified as the CH₃ torsional potential function. The tunneling dynamics of the reactive flux is dictated by its nuclear layout with respect to the CH₃ torsional angle, as the multidimensional tunneling barrier is dynamically shaped along the conical intersection seam. The effective tunneling-barrier weight-averaged over the quantum-mechanical probability along the CH₃ torsional angle perfectly explains the experimental finding: the sharp variation of the tunneling rate ((700-400) ps^-1) with the CH₃ torsional mode excitations within the narrow (0-100 cm^-1) energetic window. The much longer S₁ lifetime of cis compared to trans is ascribed to the higher-lying S₁/S₂ conical intersection of the former. With the use of distinct lifetimes, vibronic bands of each conformer could be completely separated.

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.

Conformer Specific Excited-State Structure of 3-Methylthioanisole

Trans and cis conformers of 3-methylthioanisole have been spectroscopically investigated to reveal the conformer specific structural changes upon the S₁(ππ*)-S₀ excitation. The conformational cooling during the supersonic expansion is found to be quite efficient in the Ar carrier gas giving the trans conformational isomer exclusively in the molecular beam, whereas both trans and cis conformers are populated in the jet when the sample is carried in Ne. Using the Stark deflector, trans and cis conformers are unambiguously identified, showing the distinct Stark deflection profiles according to their sufficiently different dipole moments of 1.013 or 1.670 D, respectively. For the trans conformer, the methyl moiety on the meta-position adopting the eclipsed geometry in S₀ transforms into the staggered geometry in S₁ to activate a series of the CH₃ torsional mode. A Hamiltonian with the one-dimensional sinusoidal torsional potential is solved using the free-rotor basis set to explain the experiment, giving the 3-fold torsional barrier of 34 and 304 cm^-1 for S₀ and S₁, respectively. For the cis conformer, on the other hand, the CH₃ torsion is little activated in the S₁-S₀ transition as both S₀ and S₁ adopt the staggered geometry at the minimum energy points. The doublet of each band of the cis conformer is ascribed to tunneling split due to the very low CH₃ torsional barrier of 27 cm^-1 in S₀. It is found that the cis conformer undergoes a planar to pseudoplanar structural change upon the S₁-S₀ transition. Theoretical calculation based on the double-well model potential curve could explain the experiment quite well, suggesting that the SCH₃ moiety of the cis conformer in S₁ becomes out-of-plane with respect to the plane of the phenyl moiety. This implies that excited-state predissociation dynamics of trans and cis conformers of the title molecule might be different.

Multidimensional characterization of the conical intersection seam in the normal mode space

Multidimensional conical intersection seam has been characterized by utilizing the dynamic resonances in the nonadiabatic transition probability experimentally observed in the predissociation of thioanisole isotopomers.

Vibronic structure and predissociation dynamics of 2-methoxythiophenol (S1): The effect of intramolecular hydrogen bonding on nonadiabatic dynamics

Vibronic spectroscopy and the S-H bond predissociation dynamics of 2-methoxythiophenol (2-MTP) in the S₁ (ππ^*) state have been investigated for the first time. Resonant two-photon ionization and slow-electron velocity map imaging (SEVI) spectroscopies have revealed that the S₁-S₀ transition of 2-MTP is accompanied with the planar to the pseudoplanar structural change along the out-of-plane ring distortion and the tilt of the methoxy moiety. The S₁ vibronic bands up to their internal energy of ∼1000 cm^-1 are assigned from the SEVI spectra taken via various S₁ vibronic intermediate states with the aid of ab initio calculations. Intriguingly, Fermi resonances have been identified for some vibronic bands. The S-H bond breakage of 2-MTP occurs via tunneling through an adiabatic barrier under the S₁/S₂ conical intersection seam, and it is followed by the bifurcation into either the adiabatic or nonadiabatic channel at the S₀/S₂ conical intersection where the diabatic S₂ state (πσ^*) is unbound with respect to the S-H bond elongation coordinate, giving the excited (Ã) or ground (X̃) state of the 2-methoxythiophenoxy radical, respectively. Surprisingly, the nonadiabatic transition probability at the S₀/S₂ conical intersection, estimated from the velocity map ion images of the nascent D fragment from 2-MTP-d₁ (2-CH₃O-C₆H₄SD) at the S₁ zero-point energy level, is found to be exceptionally high to give the X̃/Ã product branching ratio of 2.03 ± 0.20, which is much higher than the value of ∼0.8 estimated for the bare thiophenol at the S₁ origin. It even increases to 2.33 ± 0.17 at the ν₄₅ ² mode (101 cm^-1) before it rapidly decays to 0.69 ± 0.05 at the S₁ internal energy of about 2200 cm^-1. This suggests that the strong intramolecular hydrogen bonding of S⋯D⋯OCH₃ in 2-MTP at least in the low S₁ internal energy region should play a significant role in localizing the reactive flux onto the conical intersection seam. The minimum energy pathway calculations (second-order coupled-cluster resolution of the identity or time-dependent-density functional theory) of the adiabatic S₁ state suggest that the intimate dynamic interplay between the S-H bond cleavage and intramolecular hydrogen bonding could be crucial in the nonadiabatic surface hopping dynamics taking place at the conical intersection.

Experimental observation of nonadiabatic bifurcation dynamics at resonances in the continuum

The surface crossing of bound and unbound electronic states in multidimensional space often gives rise to resonances in the continuum. This situation happens in the πσ*-mediated photodissociation reaction of 2-fluorothioanisole; optically-bright bound S1 (ππ*) vibrational states of 2-fluorothioanisole are strongly coupled to the optically-dark S2 (πσ*) state, which is repulsive along the S-CH3 elongation coordinate. It is revealed here that the reactive flux prepared at such resonances in the continuum bifurcates into two distinct reaction pathways with totally different dynamics in terms of energy disposal and nonadiabatic transition probability. This indicates that the reactive flux in the Franck-Condon region may either undergo nonadiabatic transition funneling through the conical intersection from the upper adiabat, or follow a low-lying adiabatic path, along which multiple dynamic saddle points may be located. Since 2-fluorothioanisole adopts a nonplanar geometry in the S1 minimum energy, the quasi-degenerate S1/S2 crossing seam in the nonplanar geometry, which lies well below the planar S1/S2 conical intersection, is likely responsible for the efficient vibronic coupling, especially in the low S1 internal energy region. As the excitation energy increases, bound-to-continuum coupling is facilitated with the aid of intramolecular vibrational redistribution, along many degrees of freedom spanning the large structural volume. This leads to the rapid domination of the continuum character of the reactive flux. This work reports direct and robust experimental observations of the nonadiabatic bifurcation dynamics of the reactive flux occurring at resonances in the continuum of polyatomic molecules.

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.

Molecular Photofragmentation Dynamics in the Gas and Condensed Phases

Exciting a molecule with an ultraviolet photon often leads to bond fission, but the final outcome of the bond cleavage is typically both molecule and phase dependent. The photodissociation of an isolated gas-phase molecule can be viewed as a closed system: Energy and momentum are conserved, and the fragmentation is irreversible. The same is not true in a solution-phase photodissociation process. Solvent interactions may dissipate some of the photoexcitation energy prior to bond fission and will dissipate any excess energy partitioned into the dissociation products. Products that have no analog in the corresponding gas-phase study may arise by, for example, geminate recombination. Here, we illustrate the extent to which dynamical insights from gas-phase studies can inform our understanding of the corresponding solution-phase photochemistry and how, in the specific case of photoinduced ring-opening reactions, solution-phase studies can in some cases reveal dynamical insights more clearly than the corresponding gas-phase study.

Ultrafast X-Ray Crystallography and Liquidography

Time-resolved X-ray diffraction provides direct information on three-dimensional structures of reacting molecules and thus can be used to elucidate structural dynamics of chemical and biological reactions. In this review, we discuss time-resolved X-ray diffraction on small molecules and proteins with particular emphasis on its application to crystalline (crystallography) and liquid-solution (liquidography) samples. Time-resolved X-ray diffraction has been used to study picosecond and slower dynamics at synchrotrons and can now access even femtosecond dynamics with the recent arrival of X-ray free-electron lasers.

Attosecond Time Delay in Photoemission and Electron Scattering near Threshold

We study the time delay in the primary photoemission channel near the opening of an additional channel and compare it with the Wigner time delay in elastic scattering of the photoelectron near the corresponding inelastic threshold. The photoemission time delay near threshold is significantly enhanced, to a measurable 40 as, in comparison to the corresponding elastic scattering delay. The enhancement is due to the different lowest order of interelectron interaction coupling the primary and additional photoemission channels. We illustrate these findings by considering photodetachment from the H^{-} negative ion, and compare it with electron scattering on the hydrogen atom near the first excitation threshold. Other threshold processes of atomic photoionization and molecular photofragmentation, where photoemission time delay is enhanced, are identified. This opens the possibility of studying threshold behavior utilizing attosecond chronoscopy.

Control of Nonadiabatic Passage through a Conical Intersection by a Dynamic Resonance

Nonadiabatic processes, dominated by dynamic passage of reactive fluxes through conical intersections (CIs), are considered to be appealing means for manipulating reaction paths, particularly via initial vibrational preparation. Nevertheless, obtaining direct experimental evidence of whether specific-mode excitation affects the passage at the CI is challenging, requiring well-resolved time- or frequency-domain experiments. Here promotion of methylamine-d2 (CH3ND2) molecules to spectral-resolved rovibronic states on the excited S1 potential energy surface, coupled to sensitive D photofragment probing, allowed us to follow the N-D bond fission dynamics. The branching ratios between slow and fast D photofragments and the internal energies of the CH3ND(X̃) photofragments confirm correlated anomalies for predissociation initiated from specific rovibronic states. These anomalies reflect the existence of a dynamic resonance that strongly depends on the energy of the initially excited rovibronic states, the evolving vibrational mode on the repulsive S1 part during N-D bond elongation, and the manipulated passage through the CI that leads to CH3ND radicals excited with C-N-D bending. This resonance plays an important role in the bifurcation dynamics at the CI and can be foreseen to exist in other photoinitiated processes and to control their outcome.