Controlling Electronic Product Branching at Conical Intersections in the UV Photolysis of para-Substituted Thiophenols

Bond analysis in meta- and para-substituted thiophenols: overlap descriptors, local mode analysis, and QTAIM

CONTEXT

This study delves into the chemical nuances of thiophenols and their derivatives through a comprehensive computational analysis, moving beyond traditional energetic perspectives such as bond dissociation enthalpy and S-H dissociation dynamics. By employing the overlap model along with its topological descriptors (OP/TOP), quantum theory of atoms in molecules (QTAIM), and local vibrational mode (LVM) theories, the research provides a deeper understanding of the S-H and C-S bonding scenarios in substituted thiophenols. The investigation follows the electron-donating capacity of S-H substituent variation with the nature and positioning of other ring substituents. Energy profile analyses indicate distinct stability differences in the cis and trans conformations of meta- and para-PhSH systems, influenced by the electron-donating strength of these substituents. The study also uncovers significant variations in S-H bond distances and descriptor values, particularly in para-substituted PhSH, reflecting the influence of electron-donating or withdrawing substituents. In contrast, alterations at the meta-position show minimal effects on C-S bond descriptors, while para-substitutions markedly influence C-S bond characteristics, demonstrating a clear correlation with the electron-donating or withdrawing capabilities of the substituents. This research sheds light on the intricate bond dynamics in aromatic systems with diverse substituents, highlighting the complex interaction between electronic effects and molecular conformation.

METHODS

The study employs the ω B97X-D/Def2TZVP level of theory for molecular geometries, ensuring accurate characterization of structures as true minima via analytical harmonic frequency determination. The electronic properties of S-H and C-S bonds in variously substituted thiophenols were analyzed using OP/TOP, QTAIM, and LVM methodologies. Computational processes, including conformational scans, geometry optimizations, and vibrational frequency calculations, were conducted using Gaussian 09, with ultra-fine integration grids and tight convergence criteria for the SCF procedure. Bond descriptors were computed utilizing ChemBOS, Multiwfn, and LModeA software, providing a robust and detailed examination of bond properties.

Photodegradation of Flutamide and Halogen Derivatives of Nitrobenzotrifluoride: The NO Release Channel

The photodegradation of the nonsteroidal antiandrogen drug flutamide has been long linked to the photoisomerization involving the nitro group. In this work, the dynamics of NO photoelimination upon photolysis at 266 nm of flutamide, nitrobenzotrifluoride, and its halogen derivatives were investigated. Similar to nitrobenzene and its derivatives, a bimodal translational energy distribution was observed for the NO photofragment indicating the presence of two distinct elimination channels resulting in slow and fast components. The trends in the slow/fast branching ratio show that halogen substitution at the para position increases the triplet state yield due to the internal heavy-atom effect leading to enhancement of the fast component. Furthermore, the topology of the triplet state potential energy surface showed that the minimum energy path favors the oxaziridine ring-type intermediate over the NO2 roaming mechanism in all five molecules investigated. The steric interaction between the NO2 group and the CF3 group, which are placed in the ortho position, lowers the barrier for the formation of the oxaziridine transition state compared to that of nitrobenzene.

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.

Anharmonic excited state frequencies of para-difluorobenzene, toluene and catechol using analytic RI-CC2 second derivatives

Analytic second nuclear derivatives for excited electronic state energies have been implemented for the resolution-of-the-identity accelerated CC2, CIS(D_∞) and ADC(2) models. Our efficient implementation with O(N²) memory demands enables the treatment of medium sized molecules with large basis sets and high numerical precision and thereby paves the way for semi-numerical evaluation of the higher-order derivatives required for anharmonic corrections to excited state vibrational frequencies. We compare CC2 harmonic and anharmonic excited state frequencies with experimental values for para-difluorobenzene, toluene and catechol. Basis set problems occur for out-of-plane bending vibrations due to intramolecular basis set superposition error. For non-planar molecules and in plane modes of planar molecules, the agreement between theory and experiment is better than 30 cm^-1 on average and we reassign a number of experimental bands on the basis of the ab initio predictions.

CASPT2, CASSCF and non-adiabatic molecular dynamics (NAMD) studies on the low-lying electronic states of 1H-1,2,3-triazole photolysis

The photolysis mechanisms of 1H-1,2,3-triazole and 1H-1,2,3-benzotriazole were elucidated by employing multiconfigurational methods (CASSCF and CASPT2) and non-adiabatic molecular dynamics.

Taking the plunge: chemical reaction dynamics in liquids

The dynamics of chemical reactions in liquid solutions are now amenable to direct study using ultrafast laser spectroscopy techniques and advances in computer simulation methods. The surrounding solvent affects the chemical reaction dynamics in numerous ways, which include: (i) formation of complexes between reactants and solvent molecules; (ii) modifications to transition state energies and structures relative to the reactants and products; (iii) coupling between the motions of the reacting molecules and the solvent modes, and exchange of energy; (iv) solvent caging of reactants and products; and (v) structural changes to the solvation shells in response to the changing chemical identity of the solutes, on timescales which may be slower than the reactive events. This article reviews progress in the study of bimolecular chemical reaction dynamics in solution, concentrating on reactions which occur on ground electronic states. It illustrates this progress with reference to recent experimental and computational studies, and considers how the various ways in which a solvent affects the chemical reaction dynamics can be unravelled. Implications are considered for research in fields such as mechanistic synthetic chemistry.

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.

Exploring the Dynamics of the Photoinduced Ring-Opening of Heterocyclic Molecules

Excited states formed by electron promotion to an antibonding σ* orbital are now recognized as key to understanding the photofragmentation dynamics of a broad range of heteroatom containing small molecules: alcohols, thiols, amines, and many of their aromatic analogues. Such excited states may be populated by direct photoexcitation, or indirectly by nonadiabatic transfer of population from some other optically excited state (e.g., a ππ* state). This Perspective explores the extent to which the fast-growing literature pertaining to such (n/π)σ*-state mediated bond fissions can inform and enhance our mechanistic understanding of photoinduced ring-opening in heterocyclic molecules.

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.

Photochemical reaction dynamics of 2,2'-dithiobis(benzothiazole): direct observation of the addition product of an aromatic thiyl radical to an alkene with time-resolved vibrational and electronic absorption spectroscopy

The photochemical reaction dynamics of the benzothiazole-2-thiyl (BS) radical, produced by 330 nm ultraviolet photolysis of 2,2'-dithiobis(benzothiazole) (BSSB), are examined on the picosecond time scale. The initial addition product of a thiol-ene reaction between the BS radical and styrene is directly observed by transient vibrational absorption spectroscopy (TVAS). Transient electronic absorption spectroscopy (TEAS) in the ultraviolet and visible spectral regions reveals rapid formation of the ground state BS radical with a time constant of ∼200 fs. The photolytically generated BS radical decays through geminate recombination to the parent molecule BSSB and competitive formation of a BS radical dimer with a rate coefficient of (3.7 ± 0.2) × 10(10) M(-1) s(-1) in methanol, and thereafter (36 ± 1)% of the initially formed BS radicals survive at the longest time delay (1.3 ns). In styrene solution, in contrast to methanol and toluene solutions, kinetic traces of the BS radical show an additional decay with a time constant of 305 ± 13 ps, and a broad band at 345-500 nm grows with the same time constant, suggesting a bimolecular reaction of the BS radical with styrene. The TVAS measurements reveal an absorption band of the ground state BS radical at 1301 cm(-1) in toluene solution, and the band decays with a time constant of 294 ± 32 ps in styrene solution. Two product bands grow at 1239 cm(-1) and 1429 cm(-1) with respective time constants of 312 ± 68 ps and 325 ± 33 ps, and are attributed to the addition product BS-St radical formed from the BS radical and styrene. A bimolecular reaction rate coefficient of kreact = (3.8 ± 0.2) × 10(8) M(-1) s(-1) is deduced and 22 ± 1% of the initially formed BS radicals are converted to the BS-St radical in neat styrene solution.

Femtosecond to microsecond observation of the photochemical reaction of 1,2-di(quinolin-2-yl)disulfide with methyl methacrylate

Multiple radical reaction steps have been observed in a continuous sequence with sub-picosecond to microsecond transient absorption spectroscopy.

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.

Exploring Autoionization and Photoinduced Proton-Coupled Electron Transfer Pathways of Phenol in Aqueous Solution

The excited state dynamics of phenol in water have been investigated using transient absorption spectroscopy. Solvated electrons and vibrationally cold phenoxyl radicals are observed upon 200 and 267 nm excitation, but with formation time scales that differ by more than 4 orders of magnitude. The impact of these findings is assessed in terms of the relative importance of autoionization versus proton-coupled electron transfer mechanisms in this computationally tractable model system.

Tuning photochemistry: substituent effects on πσ* state mediated bond fission in thioanisoles

The electronic branching in the thiophenoxyl radicals formed by UV photolysis of thioanisole can be tuned by placing electron withdrawing/donating substituents at the 4-position.

Spin-polarized hydrogen Rydberg time-of-flight: experimental measurement of the velocity-dependent H atom spin-polarization

We have developed a new experimental method allowing direct detection of the velocity dependent spin-polarization of hydrogen atoms produced in photodissociation. The technique, which is a variation on the H atom Rydberg time-of-flight method, employs a double-resonance excitation scheme and experimental geometry that yields the two coherent orientation parameters as a function of recoil speed for scattering perpendicular to the laser propagation direction. The approach, apparatus, and optical layout we employ are described here in detail and demonstrated in application to HBr and DBr photolysis at 213 nm. We also discuss the theoretical foundation for the approach, as well as the resolution and sensitivity we achieve.

Photodissociation dynamics of the thiophenoxy radical at 248, 193, and 157 nm

The photodissociation dynamics of the thiophenoxy radical (C6H5S) have been investigated using fast beam coincidence translational spectroscopy. Thiophenoxy radicals were produced by photodetachment of the thiophenoxide anion followed by photodissociation at 248 nm (5.0 eV), 193 nm (6.4 eV), and 157 nm (7.9 eV). Experimental results indicate two major competing dissociation channels leading to SH + C6H4 (o-benzyne) and CS + C5H5 (cyclopentadienyl) with a minor contribution of S + C6H5 (phenyl). Photofragment mass distributions and translational energy distributions were measured at each dissociation wavelength. Transition states and minima for each reaction pathway were calculated using density functional theory to facilitate experimental interpretation. The proposed dissociation mechanism involves internal conversion from the initially prepared electronic excited state to the ground electronic state followed by statistical dissociation. Calculations show that SH loss involves a single isomerization step followed by simple bond fission. For both SH and S loss, C-S bond cleavage proceeds without an exit barrier. By contrast, the CS loss pathway entails multiple transition states and minima as it undergoes five membered ring formation and presents a small barrier with respect to products. The calculated reaction pathway is consistent with the experimental translational energy distributions in which the CS loss channel has a broader distribution peaking farther away from zero than the corresponding distributions for SH loss.

Comparing molecular photofragmentation dynamics in the gas and liquid phases

This article explores the extent to which insights gleaned from detailed studies of molecular photodissociations in the gas phase (i.e. under isolated molecule conditions) can inform our understanding of the corresponding photofragmentation processes in solution. Systems selected for comparison include a thiophenol (p-methylthiophenol), a thioanisole (p-methylthioanisole) and phenol, in vacuum and in cyclohexane solution. UV excitation in the gas phase results in RX-Y (X = O, S; Y = H, CH3) bond fission in all cases, but over timescales that vary by ~4 orders of magnitude - all of which behaviours can be rationalised on the basis of the relevant bound and dissociative excited state potential energy surfaces (PESs) accessed by UV photoexcitation, and of the conical intersections that facilitate radiationless transfer between these PESs. Time-resolved UV pump-broadband UV/visible probe and/or UV pump-broadband IR probe studies of the corresponding systems in cyclohexane solution reveal additional processes that are unique to the condensed phase. Thus, for example, the data clearly reveal evidence of (i) vibrational relaxation of the photoexcited molecules prior to their dissociation and of the radical fragments formed upon X-Y bond fission, and (ii) geminate recombination of the RX and Y products (leading to reformation of the ground state parent and/or isomeric adducts). Nonetheless, the data also show that, in each case, the characteristics (and the timescale) of the initial bond fission process that occurs under isolated molecule conditions are barely changed by the presence of a weakly interacting solvent like cyclohexane. These condensed phase studies are then extended to an ether analogue of phenol (allyl phenyl ether), wherein UV photo-induced RO-allyl bond fission constitutes the first step of a photo-Claisen rearrangement.

Photofragmentation Dynamics in Solution Probed by Transient IR Absorption Spectroscopy: πσ*-Mediated Bond Cleavage in p-Methylthiophenol and p-Methylthioanisole

The 267 nm photodissociation dynamics of p-methylthiophenol (p-MePhSH) and p-methylthioanisole (p-MePhSMe) dissolved in CD3CN have been probed by subpicosecond time-resolved broadband infrared spectroscopy. Prompt (τ < 1 ps) S-H bond fission in p-MePhSH is confirmed by monitoring the time-evolution of the parent (S0) bleach and the transient absorption of the p-MePhS products. Vibrational relaxation of the latter occurs on a ∼8.5 ps time scale, and ∼40% of the total radical population undergoes geminate recombination over a ∼150 ps time scale, yielding (mainly) the p-MePhSH(S0) parent. S-Me bond fission following photoexcitation to the S1 state of p-MePhSMe occurs over a much longer timescale, with a rate that is very dependent on the degree of vibrational excitation within S1. The various findings are compared and contrasted with results from complementary gas-phase photofragmentation studies of both molecules, which are shown to provide a valuable starting point for describing the solution-phase dynamics.