Excited state hydrogen transfer dynamics in substituted phenols and their complexes with ammonia: ππ∗-πσ∗ energy gap propensity and ortho-substitution effect

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear-electronic coupling in addition to accurately predicting Franck-Condon vertical excitations. In this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

Quantum Dynamics of Photodissociation: Recent Advances and Challenges

Recent advances in constructing accurate potential energy surfaces and nonadiabatic couplings from high-level ab initio data have revealed detailed potential landscapes in not only the ground electronic state but also excited ones. They enabled quantitatively accurate characterization of photoexcited reactive systems using quantum mechanical methods. In this Perspective, we survey the recent progress in quantum mechanical studies of adiabatic and nonadiabatic photodissociation dynamics, focusing on initial state control and product energy disposal. These new insights helped to understand quantum effects in small prototypical systems, and the results serve as benchmarks for developing more approximate theoretical methods.

Photodissociation dynamics of m- and p-cresol in the S1 state: Interplay between the mode-randomization and H atom tunneling reaction

The H atom tunneling dissociation dynamics of the S1 state of meta- or para-cresol has been investigated by using the picosecond time-resolved pump-probe spectroscopy in a state-specific manner. The S1 state lifetime (mainly due to the H atom tunneling reaction) is found to be mode-dependent whereas it quickly converges and remains constant as the rapid intramolecular vibrational energy redistribution (IVR) starts to participate in the S1 state relaxation with the increase of the S1 internal energy (Eint). The IVR rate and its change with increasing Eint have been reflected in the parent ion transients taken by tuning the total energy (hνpump + hνprobe) just above the adiabatic ionization threshold (so that the dissipation of the initial mode-character could be monitored as a function of the reaction time), indicating that the mode randomization rate into the S1 isoenergetic manifolds exceeds the tunneling rate quite early in terms of Eint for m-cresol (≤∼1200 cm-1) or p-cresol (≤∼800 cm-1) compared to the case of phenol (≤∼1800 cm-1). Though the H atom tunneling dynamics of phenol (S1) seems to be little influenced by the methyl substitution on the either m- or p-position, the IVR rate has been found to be strongly accelerated due to the sharply-increasing (S1) density of states with increasing Eint due to the pivotal role of the low-frequency CH3 torsional mode.

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.

πσ*-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.

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'.

Excited state dynamics of protonated dopamine: hydration and conformation effects

Electronic and vibrational spectroscopy in a cryogenic ion trap has been applied to protonated dopamine water clusters and assigned with the help of quantum chemistry calculations performed in the ground and electronic excited states. A dramatic hydration effect is observed when dopamine is solvated by three water molecules. The broad electronic spectra recorded for the bare and small water clusters containing protonated dopamine turn to sharp, well-resolved vibronic transitions in the 1-3 complex. This reflects the change induced by hydration in the photodynamics of protonated dopamine which is initially controlled by an excited state proton transfer (ESPT) reaction from the ammonium group toward the catechol ring. Interestingly, conformer selectivity is revealed in the 1-3 complex which shows two low lying energy conformers for which the ESPT reaction is prevented or not depending on the H-bond network formed between the dopamine and water molecules.

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.

Direct nonadiabatic quantum dynamics simulations of the photodissociation of phenol

Gaussian wavepacket methods are becoming popular for the investigation of nonadiabatic molecular dynamics. In the present work, a recently developed efficient algorithm for the Direct Dynamics variational Multi-Configurational Gaussian (DD-vMCG) method has been used to describe the multidimensional photodissociation dynamics of phenol including all degrees of freedom. Full-dimensional quantum dynamic calculations including for the first time six electronic states (¹ππ, 1¹ππ*, 1¹πσ*, 2¹πσ*, 2¹ππ*, 3¹ππ*), along with a comparison to an existing analytical 4-state model for the potential energy surfaces are presented. Including the fifth singlet excited state is shown to have a significant effect on the nonadiabatic photodissociation of phenol to the phenoxyl radical and hydrogen atom. State population and flux analysis from the DD-vMCG simulations of phenol provided further insights into the decay mechanism, confirming the idea of rapid relaxation to the ground state through the ¹ππ/1¹πσ* conical intersection.

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.

Controlled ultrafast ππ*-πσ* dynamics in tryptophan-based peptides with tailored micro-environment

Ultrafast charge, energy and structural dynamics in molecules are driven by the topology of the multidimensional potential energy surfaces that determines the coordinated electronic and nuclear motion. These processes are also strongly influenced by the interaction with the molecular environment, making very challenging a general understanding of these dynamics on a microscopic level. Here we use electrospray and mass spectrometry technologies to produce isolated molecular ions with a controlled micro-environment. We measure ultrafast photo-induced ππ*-πσ* dynamics in tryptophan species in the presence of a single, charged adduct. A striking increase of the timescale by more than one order of magnitude is observed when changing the added adduct atom. A model is proposed to rationalize the results, based on the localized and delocalized effects of the adduct on the electronic structure of the molecule. These results offer perspectives to control ultrafast molecular processes by designing the micro-environment on the Angström length scale.

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.

Revealing the role of excited state proton transfer (ESPT) in excited state hydrogen transfer (ESHT): systematic study in phenol-(NH3) n clusters

Excited State Hydrogen Transfer (ESHT), proposed at the end of the 20th century by the corresponding authors, has been observed in many neutral or protonated molecules and become a new paradigm to understand excited state dynamics/photochemistry of aromatic molecules. For example, a significant number of photoinduced proton-transfer reactions from X–H bonds have been re-defined as ESHT, including those of phenol, indole, tryptophan, aromatic amino acid cations and so on. Photo-protection mechanisms of biomolecules, such as isolated nucleic acids of DNA, are also discussed in terms of ESHT. Therefore, a systematic and up-to-date description of ESHT mechanism is important for researchers in chemistry, biology and related fields. In this review, we will present a general model of ESHT which unifies the excited state proton transfer (ESPT) and the ESHT mechanisms and reveals the hidden role of ESPT in controlling the reaction rate of ESHT. For this purpose, we give an overview of experimental and theoretical work on the excited state dynamics of phenol–(NH₃)_n clusters and related molecular systems. The dynamics has a significant dependence on the number of solvent molecules in the molecular cluster. Three-color picosecond time-resolved IR/near IR spectroscopy has revealed that ESHT becomes an electron transfer followed by a proton transfer in highly solvated clusters. The systematic change from ESHT to decoupled electron/proton transfer according to the number of solvent molecules is rationalized by a general model of ESHT including the role of ESPT.

A general model of excited state hydrogen transfer (ESHT) which unifies ESHT and the excited state proton transfer (ESPT) is presented from experimental and theoretical works on phenol–(NH₃)_n. The hidden role of ESPT is revealed.

Effects of Ring Fluorination on the Ultraviolet Photodissociation Dynamics of Phenol

The dynamics of photoinduced O-H bond fission in five fluorinated phenols (2-fluorophenol, 3-fluorophenol, 2,6-difluorophenol, 3,4,5-trifluorophenol, and pentafluorophenol) have been investigated by H Rydberg atom photofragment translational spectroscopy following excitation at many wavelengths in the range 220 ≤ λ ≤ 275 nm. The presence of multiple fluorine substituents reduces the efficiency of O-H bond fission (by tunneling) from the first excited (1¹ππ*) electronic state, whereas all bar the perfluorinated species undergo O-H bond fission when excited at shorter wavelengths (to the 2¹ππ* state). As in bare phenol, O-H bond fission is deduced to occur by non-adiabatic coupling at conical intersections between the photoprepared "bright" ππ* states and the 1¹πσ* potential energy surface. In all cases, the fluorophenoxyl photoproducts are found to be formed in a range of vibrational levels, all of which include an odd number of quanta (typically one) in an out-of-plane (a″) vibrational mode; this product vibration is viewed as a legacy of the parent out-of-plane motions that promote non-adiabatic coupling to the dissociative 1¹πσ* potential. The radical products also show activity in in-plane vibrations involving coupled (both in- and out-of-phase) C-O and C-F wagging motions, which can be traced to the impulse between the recoiling O and H atoms and, in detail, are sensitive to the presence (or not) of an intramolecular F···H-O hydrogen bond. Upper limit values for the O-H bond dissociation energies are reported for all molecules studied apart from pentafluorophenol.

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.

Photolysis dynamics of m- and o-fluorophenol: Substitution effects on tunneling mechanism

The photolysis dynamics of m-fluorophenol (m-FPhOH) and o-fluorophenol (o-FPhOH) have been investigated with time-resolved velocity map imaging (TR-VMI) and time-resolved ion-yield (TR-IY) techniques. Following excitation to the origin of S₁ (ππ∗) state of m- and o-FPhOH, H atoms elimination mediated by tunneling through the potential barrier under the S₁ (ππ∗)/S₂ (πσ∗) conical intersection (CI) has been observed as a Gaussian feature signal centered at a total kinetic energy release (TKER) of ∼6000 cm^-1 for both molecules. The quantum tunneling mechanism has been identified as the main decay pathway of S₁ state for m-FPhOH, and the tunneling lifetime of 2.1 ns has been obtained from the TR-VMI measurements of H fragments. This tunneling mechanism is further confirmed by the studies on the selective O-H deuterated species, m-FPhOD, and consolidated by our theoretical calculations. However, the photolysis dynamics is quite different for the photoexcited o-FPhOH. The much lower yield of the H atoms originating from tunneling hinders the extraction of a reliable tunneling lifetime for o-FPhOH. Our theoretical calculations exhibit a broader and higher potential barrier exists beneath the S₁/S₂ CI of o-FPhOH, which increase the difficulty for tunneling. Furthermore, the special existence of intramolecular hydrogen bond in o-FPhOH is probably also the key factor that affects the tunneling rate, which would restrict the O-H stretch motion.

Photoinduced water oxidation in pyrimidine-water clusters: a combined experimental and theoretical study

The photocatalytic oxidation of water with molecular or polymeric N-heterocyclic chromophores is a topic of high current interest in the context of artificial photosynthesis, that is, the conversion of solar energy to clean fuels. Hydrogen-bonded clusters of N-heterocycles with water molecules in a molecular beam are simple model systems for which the basic mechanisms of photochemical water oxidation can be studied under well-defined conditions. In this work, we explored the photoinduced H-atom transfer reaction in pyrimidine-water clusters yielding pyrimidinyl and hydroxyl radicals with laser spectroscopy, mass spectrometry and trajectory-based ab initio molecular dynamics simulations. The oxidation of water by photoexcited pyrimidine is unequivocally confirmed by the detection of the pyrimidinyl radical. The dynamics simulations provide information on the time scales and branching ratios of the reaction. While relaxation to local minima of the S₁ potential-energy surface is the dominant reaction channel, the H-atom transfer reaction occurs on ultrafast time scales (faster than about 100 fs) with a branching ratio of a few percent.

Insights into the Mechanism of Nonadiabatic Photodissociation from Product Vibrational Distributions. The Remarkable Case of Phenol

The fate of a photoexcited molecule is often strongly influenced by electronic degeneracies, such as conical intersections, which break the Born-Oppenheimer separation of electronic and nuclear motion. Detailed information concerning internal energy redistribution in a nonadiabatic process can be extracted from the product state distribution of a photofragment in photodissociation. Here, we focus on the nonadiabatic photodissociation of phenol and discuss the internal excitation of the phenoxyl fragment using both symmetry analysis and wave packet dynamics. It is shown that unique and general selection rules exist, which can be attributed to the geometric phase in the adiabatic representation. Further, our results provide a reinterpretation of the experimental data, shedding light on the impact of conical intersections on the product state distribution.

Real-Time Observation of Fermi Resonances in the S1 State of Phenol

Fermi resonances in the first electronically excited (S₁) state of phenol have been observed in real time. Quantum beats associated with coherent superposition of Fermi resonant eigenstates are manifested as temporal oscillations of the ionization cross sections of which the amplitudes are strongly dependent on the total ionization energy. This indicates that coherently excited eigenstates are effectively decomposed into their zeroth-order states, providing the unique opportunity for the investigation of nonstationary state dynamics in real-time. Energy gaps (Δν̃) of eigenstates within the laser coherence width have been most precisely determined up to date, giving Δν̃ ∼ 3.302 ± 0.001 or 1.655 ± 0.001 cm^-1 for the 1¹/4¹10b¹ or 12²/8a¹ Fermi doublets, respectively. Dephasing rate suddenly increases as the S₁ internal energy becomes above ∼1500 cm^-1, revealing the important role of energy randomization dynamics during the H atom tunneling process of phenol in S₁.

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.

Multidimensional H Atom Tunneling Dynamics of Phenol: Interplay between Vibrations and Tunneling

Multidimensional facets of the hydrogen tunneling dynamics of phenol excited in S₁ (ππ*) have been unraveled to give particular S₁ vibronic states strongly coupled or actively decoupled to the O-H tunneling coordinate. Strong mode-dependent variation of the tunneling rate measured with picosecond lasers indicates that tunneling probability is extremely sensitive to low-frequency vibrational modes seemingly orthogonal to the O-H elongation coordinate unless the rate of energy randomization exceeds that of tunneling. The multidimensional nature of tunneling has also been manifested in efficient internal-to-translational energy transfers observed at S₁ vibronic modes strongly coupled to the tunneling coordinate, giving insights into otherwise the formidable multidimensional map of tunneling process. The nonadiabatic bifurcation dynamics in the later stage of the chemical reaction has been disentangled by analyzing picosecond time-resolved product state distributions, resolving a long controversial issue regarding the origin of high or low kinetic energy component of the product translational energy distributions.

Up to a Sign. The Insidious Effects of Energetically Inaccessible Conical Intersections on Unimolecular Reactions

It is now well established that conical intersections play an essential role in nonadiabatic radiationless decay where their double-cone topography causes them to act as efficient funnels channeling wave packets from the upper to the lower adiabatic state. Until recently, little attention was paid to the effect of conical intersections on dynamics on the lower state, particularly when the total energy involved is significantly below that of the conical intersection seam. This energetic deficiency is routinely used as a sufficient condition to exclude consideration of excited states in ground state dynamics. In this account, we show that, this energy criterion notwithstanding, energy inaccessible conical intersections can and do exert significant influence on lower state dynamics. The origin of this influence is the geometric phase, a signature property of conical intersections, which is the fact that the real-valued electronic wave function changes sign when transported along a loop containing a conical intersection, making the wave function double-valued. This geometric phase is permitted by an often neglected property of the real-valued adiabatic electronic wave function; namely, it is determined only up to an overall sign. Noting that in order to change sign a normalized, continuous function must go through zero, for loops of ever decreasing radii, demonstrating the need for an electronic degeneracy (intersection) to accompany the geometric phase. Since the total wave function must be single-valued a compensating geometry dependent phase needs to be included in the total electronic-nuclear wave function. This Account focuses on how this consequence of the geometric phase can modify nuclear dynamics energetically restricted to the lower state, including tunneling dynamics, in directly measurable ways, including significantly altering tunneling lifetimes, thus confounding the relation between measured lifetimes and barrier heights and widths, and/or completely changing product rotational distributions. Some progress has been made in understanding the origin of this effect. It has emerged that for a system where the lower adiabatic potential energy surface exhibits a topography comprised of two saddle points separated by a high energy conical intersection, the effect of the geometric phase can be quite significant. In this case topologically distinct paths through the two adiabatic saddle points may lead to interference. This was pointed out by Mead and Truhlar almost 50 years ago and denoted the Molecular Aharonov-Bohm effect. Still, the difficulty in anticipating a significant geometric phase effect in tunneling dynamics due to energetically inaccessible conical intersections leads to the attribute insidious that appears in the title of this Account. Since any theory is only as relevant as the prevalence of the systems it describes, we include in this Account examples of real systems where these effects can be observed. The accuracy of the reviewed calculations is high since we use fully quantum mechanical dynamics and construct the geometric phase using an accurate diabatic state fit of high quality ab initio data, energies, energy gradients, and interstate couplings. It remains for future work to establish the prevalence of this phenomenon and its deleterious effects on the conventional wisdom discussed in this work.

Solvent reorganization triggers photo-induced solvated electron generation in phenol

Charge-transfer states with large electron–hole separation, correlating to the formation of solvated electrons, are found below the maximum of the absorbing ππ* band of solvated phenol.

Excited-state dissociation dynamics of phenol studied by a new time-resolved technique

Phenol is an important model molecule for the theoretical and experimental investigation of dissociation in the multistate potential energy surfaces. Recent theoretical calculations [X. Xu et al., J. Am. Chem. Soc. 136, 16378 (2014)] suggest that the phenoxyl radical produced in both the X and A states from the O-H bond fission in phenol can contribute substantially to the slow component of photofragment translational energy distribution. However, current experimental techniques struggle to separate the contributions from different dissociation pathways. A new type of time-resolved pump-probe experiment is described that enables the selection of the products generated from a specific time window after molecules are excited by a pump laser pulse and can quantitatively characterize the translational energy distribution and branching ratio of each dissociation pathway. This method modifies conventional photofragment translational spectroscopy by reducing the acceptance angles of the detection region and changing the interaction region of the pump laser beam and the molecular beam along the molecular beam axis. The translational energy distributions and branching ratios of the phenoxyl radicals produced in the X, A, and B states from the photodissociation of phenol at 213 and 193 nm are reported. Unlike other techniques, this method has no interference from the undissociated hot molecules. It can ultimately become a standard pump-probe technique for the study of large molecule photodissociation in multistates.

Nuclear Ensemble Approach with Importance Sampling

We show that the importance sampling technique can effectively augment the range of problems where the nuclear ensemble approach can be applied. A sampling probability distribution function initially determines the collection of initial conditions for which calculations are performed, as usual. Then, results for a distinct target distribution are computed by introducing compensating importance sampling weights for each sampled point. This mapping between the two probability distributions can be performed whenever they are both explicitly constructed. Perhaps most notably, this procedure allows for the computation of temperature dependent observables. As a test case, we investigated the UV absorption spectra of phenol, which has been shown to have a marked temperature dependence. Application of the proposed technique to a range that covers 500 K provides results that converge to those obtained with conventional sampling. We further show that an overall improved rate of convergence is obtained when sampling is performed at intermediate temperatures. The comparison between calculated and the available measured cross sections is very satisfactory, as the main features of the spectra are correctly reproduced. As a second test case, one of Tully's classical models was revisited, and we show that the computation of dynamical observables also profits from the importance sampling technique. In summary, the strategy developed here can be employed to assess the role of temperature for any property calculated within the nuclear ensemble method, with the same computational cost as doing so for a single temperature.

From fundamental science to product: a bottom-up approach to sunscreen development

Despite the pivotal role of ultraviolet (UV) radiation in sustaining life on Earth, overexposure to this type of radiation can have catastrophic effects, such as skin cancer. Sunscreens, the most common form of artificial protection against such harmful effects, absorb UV radiation before it reaches vulnerable skin cells. Absorption of UV radiation prompts ultrafast molecular events in sunscreen molecules which, ideally, would allow for fast and safe dissipation of the excess energy. However, our knowledge of these mechanisms remains limited. In this article, we will review recent advances in the field of ultrafast photodynamics (light induced molecular processes occurring within femtoseconds, fs, 10-15 s to picoseconds, ps, 10-12 s) of sunscreens. We follow a bottom-up approach to common sunscreen active ingredients, analysing any emerging trends from the current literature on the subject. Moreover, we will identify the main questions that remain unanswered, pinpoint some of the main challenges and finally comment on the outlook of this exciting field of research.

A 'bottom up', ab initio computational approach to understanding fundamental photophysical processes in nitrogen containing heterocycles, DNA bases and base pairs

The availability of non-radiative decay mechanisms by which photoexcited molecules can revert to their ground electronic state, without experiencing potentially deleterious chemical transformation, is fundamental to molecular photostability. This Perspective Article combines results of new ab initio electronic structure calculations and prior experimental data in an effort to systematise trends in the non-radiative decay following UV excitation of selected families of heterocyclic molecules. We start with the prototypical uni- and bicyclic molecules phenol and indole, and explore the structural and photophysical consequences of incorporating progressively more nitrogen atoms within the respective ring structures en route to the DNA bases thymine, cytosine, adenine and guanine. For each of the latter, we identify low energy non-radiative decay pathways via conical intersections with the ground state potential energy surface accessed by out-of-plane ring deformations. This is followed by summary descriptions and illustrations of selected rival (electron driven H atom transfer) non-radiative excited state decay processes that demand consideration once the nucleobases are merely components in larger biomolecular systems like nucleosides, and both individual and stacked base-pairs.

Photoinduced water splitting in pyridine water clusters

Ab initio calculations predict that pyridine (Py) can act as a photo-catalyst to split water by the absorption of a UV photon following the reaction Py-H₂O + hν → PyH˙ + OH˙. To test this prediction, we performed two types of experiments: in the first, we characterize the electronic spectroscopy of the PyH˙ radical in the gas phase. In the second, we evidence the reaction through the UV excitation of molecular Py-(H₂O)_n clusters obtained in a supersonic expansion and monitoring the PyH˙ reaction product. The results show unambiguously that PyH˙ is produced, and thus that water is split using pyridine as a photo-catalyst.

Watching proton transfer in real time: Ultrafast photoionization-induced proton transfer in phenol-ammonia complex cation

In this paper, we give a full account of our previous work [C. C. Shen et al., J. Chem. Phys. 141, 171103 (2014)] on the study of an ultrafast photoionization-induced proton transfer (PT) reaction in the phenol-ammonia (PhOH-NH₃) complex using ultrafast time-resolved ion photofragmentation spectroscopy implemented by the photoionization-photofragmentation pump-probe detection scheme. Neutral PhOH-NH₃ complexes prepared in a free jet are photoionized by femtosecond 1 + 1 resonance-enhanced multiphoton ionization via the S₁ state. The evolving cations are then probed by delayed pulses that result in ion fragmentation, and the ionic dynamics is followed by measuring the parent-ion depletion as a function of the pump-probe delay time. By comparing with systems in which PT is not feasible and the steady-state ion photofragmentation spectra, we concluded that the observed temporal evolutions of the transient ion photofragmentation spectra are consistent with an intracomplex PT reaction after photoionization from the initial non-PT to the final PT structures. Our experiments revealed that PT in [PhOH-NH₃]⁺ cation proceeds in two distinct steps: an initial impulsive wave-packet motion in ∼70 fs followed by a slower relaxation of about 1 ps that stabilizes the system into the final PT configuration. These results indicate that for a barrierless PT system, even though the initial PT motions are impulsive and ultrafast, the time scale to complete the reaction can be much slower and is determined by the rate of energy dissipation into other modes.

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.

Sequential Proton-Coupled Electron Transfer Mediates Excited-State Deactivation of a Eumelanin Building Block

Skin photoprotection is commonly believed to rely on the photochemistry of 5,6-dihydroxyindole (DHI)- and 5,6-dihydroxyindole-2-carboxylic acid (DHICA)-based eumelanin building blocks. Attempts to elucidate the underlying excited-state relaxation mechanisms have been partly unsuccessful due to the marked instability to oxidation. We report a study of the excited-state deactivation of DHI using steady-state and time-resolved fluorescence accompanied by high-level quantum-chemistry calculations including solvent effects. Spectroscopic data show that deactivation of the lowest excited state of DHI in aqueous buffer proceeds on the 100 ps time scale and is 20 times faster than in methanol. Quantum-chemical calculations reveal that the excited-state decay mechanism is a sequential proton-coupled electron transfer, which involves the initial formation of a solvated electron from DHI, followed by the transfer of a proton to the solvent. This unexpected finding would prompt a revision of current notions about eumelanin photophysics and photobiology.