Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)πσ? excited state with the bound (1)ππ? excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)ππ?-(1)πσ? and (1)πσ?-S(0) conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S(0), (1)ππ?, (1)πσ?) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)ππ? and (1)πσ? states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)σ/(2)π electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (ν(16a)) and ring-distortion (ν(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor.     

Time-dependent quantum wave-packet description of the 1pi sigma* photochemistry of phenol

The photoinduced hydrogen elimination reaction in phenol via the conical intersections of the dissociative 1pi sigma* state with the 1pi pi* state and the electronic ground state has been investigated by time-dependent quantum wave-packet calculations. A model including three intersecting electronic potential-energy surfaces (S0, 1pi sigma*, and 1pi pi*) and two nuclear degrees of freedom (OH stretching and OH torsion) has been constructed on the basis of accurate ab initio multireference electronic-structure data. The electronic population transfer processes at the conical intersections, the branching ratio between the two dissociation channels, and their dependence on the initial vibrational levels have been investigated by photoexciting phenol from different vibrational levels of its ground electronic state. The nonadiabatic transitions between the excited states and the ground state occur on a time scale of a few tens of femtoseconds if the 1pi pi*-1pi sigma* conical intersection is directly accessible, which requires the excitation of at least one quantum of the OH stretching mode in the 1pi pi* state. It is shown that the node structure, which is imposed on the nuclear wave packet by the initial preparation as well as by the transition through the first conical intersection (1pi pi*-1pi sigma*), has a profound effect on the nonadiabatic dynamics at the second conical intersection (1pi sigma*-S0). These findings suggest that laser control of the photodissociation of phenol via IR mode-specific excitation of vibrational levels in the electronic ground state should be possible.     

Photochemistry of pyrrole: time-dependent quantum wave-packet description of the dynamics at the 1pi sigma*-S0 conical intersections

The photoinduced hydrogen-elimination reaction in pyrrole via the conical intersections of the two (1)pi sigma(*) excited states with the electronic ground states [(1)B(1)(pi sigma(*))-S(0) and (1)A(2)(pi sigma(*))-S(0)] have been investigated by time-dependent quantum wave-packet calculations. Model potential-energy surfaces of reduced dimensionality have been constructed on the basis of accurate multireference ab initio electronic-structure calculations. For the (1)B(1)-S(0) conical intersection, the model includes the NH stretching coordinate as the tuning mode and the hydrogen out-of-plane bending coordinate as the coupling mode. For the (1)A(2)-S(0) conical intersection, the NH stretching coordinate and the screwing coordinate of the ring hydrogens are taken into account. The latter is the dominant coupling mode of this conical intersection. The electronic population-transfer processes at the conical intersections, the branching ratio between the dissociation channels, and their dependence on the initial preparation of the system have been investigated for pyrrole and deuterated pyrrole. It is shown that the excitation of the NH stretching mode strongly enhances the reaction rate, while the excitation of the coupling mode influences the branching ratio of different dissociation channels. The results suggest that laser control of the photodissociation of pyrrole via mode-specific vibrational excitation should be possible. The calculations provide insight into the microscopic details of ultrafast internal-conversion processes in pyrrole via hydrogen-detachment processes, which are aborted at the (1)pi sigma(*)-S(0) conical intersections. These mechanisms are of relevance for the photostability of the building blocks of life (e.g., the DNA bases).     

Laser control of reactions of photoswitching functional molecules

Laser control schemes of reactions of photoswitching functional molecules are proposed based on the quantum mechanical wave-packet dynamics and the design of laser parameters. The appropriately designed quadratically chirped laser pulses can achieve nearly complete transitions of wave packet among electronic states. The laser parameters can be optimized by using the Zhu-Nakamura theory of nonadiabatic transition. This method is effective not only for the initial photoexcitation process but also for the pump and dump scheme in the middle of the overall photoswitching process. The effects of momentum of the wave packet crossing a conical intersection on the branching ratio of products have also been clarified. These control schemes mentioned above are successfully applied to the cyclohexadiene/hexatriene photoisomerization (ring-opening) process which is the reaction center of practical photoswitching molecules such as diarylethenes. The overall efficiency of the ring opening can be appreciably increased by using the appropriately designed laser pulses compared to that of the natural photoisomerization without any control schemes.     

Controlling product selection in the photodissociation of formaldehyde: direct quantum dynamics from the S1 barrier

Controlling the selectivity between H(2)+CO and H+HCO in the S(1)/S(0) nonadiabatic photodissociation of formaldehyde has been investigated using direct quantum dynamics. Simulations started from the S(1) transition state have suggested that a key feature for controlling the branching ratio of ground-state products is the size of the momentum given to the wavepacket along the transition vector. Our results show that letting the wavepacket fall down from the barrier to the conical intersection with no initial momentum leads to H(2)+CO, while extra momentum toward products favors the formation of H+HCO through the same conical intersection. Quantum dynamics results are interpreted in semiclassical terms with the aid of a Mulliken-like analysis of the final population distribution among both products and the reactant on each electronic state.     

A non-adiabatic quantum-classical dynamics study of the intramolecular excited state hydrogen transfer in ortho-nitrobenzaldehyde

Ab initio surface-hopping dynamics calculations have been performed to simulate the intramolecular excited state hydrogen transfer dynamics of ortho-nitrobenzaldehyde (o-NBA) in the gas phase from the electronic S(1) excited state. Upon UV excitation, the hydrogen is transferred from the aldehyde substituent to the nitro group, generating o-nitrosobenzoic acid through a ketene intermediate. The semiclassical propagations show that the deactivation from the S(1) is ultrafast, in agreement with the experimental measurements, which detect the ketene in less than 400 fs. The trajectories show that the deactivation mechanism involves two different conical intersections. The first one, a planar configuration with the hydrogen partially transferred, is responsible for the branching between the formation of a biradical intermediate and the regeneration of the starting material. The conversion of the biradical to the ketene corresponds to the passage through a second intersection region in which the ketene group is formed.     

Ultrafast dynamics of UV-excited imidazole

The ultrafast dynamics of UV-excited imidazole in the gas phase is investigated by theoretical nonadiabatic dynamics simulations and experimental time-resolved photoelectron spectroscopy. The results show that different electronic excited-state relaxation mechanisms occur, depending on the pump wavelength. When imidazole is excited at 239.6 nm, deactivation through the NH-dissociation conical intersection is observed on the sub-50 fs timescale. After 200.8 nm excitation, competition between NH-dissociation and NH-puckering conical intersections is observed. The NH-dissociation to NH-puckering branching ratio is predicted to be 21:4, and the total relaxation time is elongated by a factor of eight. A procedure for simulation of photoelectron spectra based on dynamics results is developed and employed to assign different features in the experimental spectra.     

Ab initio based surface-hopping dynamics study on ultrafast internal conversion in cyclopropanone

Cyclopropanone exhibits an intriguing phenomenon that the fluorescence from the S(1) state disappears below 365 nm. This is ascribed to the ultrafast S(1) → S(0) internal conversion process via conical intersection, which deprives opportunity of the fluorescence emission. In this work, we have used ab initio based surface hopping dynamics method to study vibrational-mode-dependent S(1) → S(0) internal conversion of cyclopropanone. A new conical intersection between the S(1) and S(0) states is determined by the state-averaged CASSCF/cc-pVDZ calculations, which is confirmed to play a critical role in the ultrafast S(1) → S(0) internal conversion by the nonadiabatic dynamics simulations. It is found that the internal conversion occurs more efficiently when the initial kinetic energies are distributed in the four vibrational modes related to the C═O group, especially in the C-O stretching and the O-C-C-C out-of-plane torsional modes. Meanwhile, the S(1) lifetime and the time scale of the S(1) → S(0) internal conversion are estimated by the ab initio based dynamics simulations, which is consistent with the ultrafast S(1) → S(0) internal conversion and provides further evidence that the ultrafast internal conversion is responsible for the fluorescence disappearance of cyclopropanone.     

Geometric phase effects in the coherent control of the branching ratio of photodissociation products of phenol

Optimal control simulation is used to examine the control mechanisms in the photodissociation of phenol within a two-dimensional, three-electronic-state model with two conical intersections. This model has two channels for H-atom elimination, which correspond to the (2)pi and (2)sigma states of the phenoxyl radical. The optimal pulse that enhances (2)sigma dissociation initially generates a wave packet on the S(1) potential-energy surface of phenol. This wave packet is bifurcated at the S(2)-S(1) conical intersection into two components with opposite phases because of the geometric phase effect. The destructive interference caused by the geometric phase effect reduces the population around the S(1)-S(0) conical intersection, which in turn suppresses nonadiabatic transitions and thus enhances dissociation to the (2)sigma limit. The optimal pulse that enhances S(0) dissociation, on the other hand, creates a wave packet on the S(2) potential-energy surface of phenol via an intensity borrowing mechanism, thus avoiding geometric phase effects at the S(2)-S(1) conical intersection. This wave packet hits the S(1)-S(0) conical intersection directly, resulting in preferred dissociation to the (2)pi limit. The optimal pulse that initially prepares the wave packet on the S(1) potential-energy surface (PES) has a higher carrier frequency than the pulse that prepares the wave packet on the S(2) PES. This counterintuitive effect is explained by the energy-level structure and the S(2)-S(1) vibronic coupling mechanism.     

Conical intersections and semiclassical trajectories: comparison to accurate quantum dynamics and analyses of the trajectories

Semiclassical trajectory methods are tested for electronically nonadiabatic systems with conical intersections. Five triatomic model systems are presented, and each system features two electronic states that intersect via a seam of conical intersections (CIs). Fully converged, full-dimensional quantum mechanical scattering calculations are carried out for all five systems at energies that allow for electronic de-excitation via the seam of CIs. Several semiclassical trajectory methods are tested against the accurate quantum mechanical results. For four of the five model systems, the diabatic representation is the preferred (most accurate) representation for semiclassical trajectories, as correctly predicted by the Calaveras County criterion. Four surface hopping methods are tested and have overall relative errors of 40%-60%. The semiclassical Ehrenfest method has an overall error of 66%, and the self-consistent decay of mixing (SCDM) and coherent switches with decay of mixing (CSDM) methods are the most accurate methods overall with relative errors of approximately 32%. Furthermore, the CSDM method is less representation dependent than both the SCDM and the surface hopping methods, making it the preferred semiclassical trajectory method. Finally, the behavior of semiclassical trajectories near conical intersections is discussed.     

Three-state conical intersections in cytosine and pyrimidinone bases

Three-state conical intersections have been located and characterized for cytosine and its analog 5-methyl-2-pyrimidinone using multireference configuration-interaction ab initio methods. The potential energy surfaces for each base contain three different three-state intersections: two different S(0)-S(1)-S(2) intersections (gs/pi pi(*)/n(N)pi(*) and gs/pi pi(*)/n(O)pi(*)) and an S(1)-S(2)-S(3) intersection (pi pi(*)/n(N)pi(*)/n(O)pi(*)). Two-state seam paths from these intersections are shown to be connected to previously reported two-state conical intersections. Nonadiabatic coupling terms have been calculated, and the effects of the proximal third state on these quantities are detailed. In particular, it is shown that when one of these loops incorporates more than one seam point, there is a profound and predictable effect on the phase of the nonadiabatic coupling terms, and as such provides a diagnostic for the presence and location of additional seams. In addition, it is shown that each of the three three-state conical intersections located on cytosine and 5-methyl-2-pyrimidinone is qualitatively similar between the two bases in terms of energies and character, implying that, like with the stationary points and two-state conical intersections previously reported for these two bases, there is an underlying pattern of energy surfaces for 2-pyrimidinone bases, in general, and this pattern also includes three-state conical intersections.     

Ab initio molecular dynamics and time-resolved photoelectron spectroscopy of electronically excited uracil and thymine

The reaction dynamics of excited electronic states in nucleic acid bases is a key process in DNA photodamage. Recent ultrafast spectroscopy experiments have shown multicomponent decays of excited uracil and thymine, tentatively assigned to nonadiabatic transitions involving multiple electronic states. Using both quantum chemistry and first principles quantum molecular dynamics methods we show that a true minimum on the bright S2 electronic state is responsible for the first step that occurs on a femtosecond time scale. Thus the observed femtosecond decay does not correspond to surface crossing as previously thought. We suggest that subsequent barrier crossing to the minimal energy S2/S1 conical intersection is responsible for the picosecond decay.     

Electronic and steric effects on the photo-induced C→E ring-opening of structurally modified furylfulgides

The ultrafast C→E ring-opening reactions of four selectively modified furylfulgides have been studied by means of ultrafast broadband transient absorption spectroscopy after femtosecond laser excitation at λ = 500 nm. A large difference in the dynamics was found in the case of benzannulation at the furyl moiety as an example for an electronic effect by extension of the conjugated π-electron system compared to furylfulgides carrying sterically different alkyl substituents at the central cyclohexadiene (CHD) ring. The measured very similar spectro-temporal absorption maps for the furylfulgides with a methyl or isopropyl group at the CHD ring or an intramolecular alkyl bridge from the CHD to the furyl moiety showed two distinctive excited-state absorptions with slightly different decay times. The first time constant (τ(1) = 0.39-0.57 ps) was assigned to the rapid departure of the excited wavepacket from the Franck-Condon region. The slightly longer second decay time of τ(2) = 0.66-0.92 ps, depending on the compound, was attributed to the electronic deactivation and ring-opening through a conical intersection to the S(0) state. In contrast, the benzannulation at the furyl moiety was found to lead to a bi-phasic excited-state decay with τ(2) = 4.7 ps and a much slower additional contribution of τ(3) = 17.4 ps, ≈25 times longer compared to the normal furylfulgides. The drastic change is attributed to a trapping of excited molecules in a local potential energy minimum en route to the conical intersection.     

Ab initio study of ultrafast photochemical reaction dynamics of phenol blue

Phenol blue (PB) is a primary skeletal structure part of indoaniline dyes and well-known as a solvatochromic dye. It has been recently observed by pump-probe (PP) transient absorption measurements that PB shows ultrafast ground state recovery within a few hundred femtoseconds after photoexcitation. In this work, the ultrafast photochemical reaction mechanism of PB has been investigated using direct ab initio (CASSCF) nonadiabatic molecular dynamics with the trajectory surface hopping (TSH) method. The swarm of trajectories starting from the S1 Franck-Condon (FC) point has mostly shown surface hops (nonadiabatic transitions) from the S1 state to the S0 state at 110-120 fs in the vicinity of an S1/S0 conical intersection and after decay to the S0 state bifurcated into two (Reverse and Forward) directions with almost the same branching ratio and reached the vicinity of the S0 minimum energy point at 200-300 fs, which is in good agreement with the fast time component of the ground state recovery in the PP measurements. After reaching the vicinity of the S0 minimum energy point, the trajectories showed a coherent vibration of bending motion between quinoneimine and aniline rings with a low frequency of 43 cm-1, which presumably corresponds to a coherently photoexcitation-induced vibrational mode with a low frequency recently observed by the PP measurements.     

Comment on "Optical conversion of conical intersection to avoided crossing" by Y. Arasaki and K. Takatsuka, Phys. Chem. Chem. Phys., 2010, 12, 1239

A recent paper in this journal proposed the conversion of conical intersections to avoided crossings by lowering the symmetry with an optical field. The article also claimed that the characters of nonadiabatic transitions caused by avoided crossings and conical intersections are qualitatively different. The present comment shows that this proposal and this claim result from an incorrect appreciation of the nature of conical intersections and avoided crossings. Conical intersections are moved, not removed, by almost all perturbations. Furthermore, there is no dichotomy between avoided crossing mechanisms and conical intersection mechanisms; as the parameters of the problem change and the typical locally avoided crossing involved in nonadiabatic dynamics becomes farther from the conical intersection, there is a gradual shift in the nature of the nonadiabatic transitions, with a continuum of possible behaviors, not just two.     

The non-adiabatic nanoreactor: towards the automated discovery of photochemistry

The ab initio nanoreactor has previously been introduced to automate reaction discovery for ground state chemistry. In this work, we present the nonadiabatic nanoreactor, an analogous framework for excited state reaction discovery. We automate the study of nonadiabatic decay mechanisms of molecules by probing the intersection seam between adiabatic electronic states with hyper-real metadynamics, sampling the branching plane for relevant conical intersections, and performing seam-constrained path searches. We illustrate the effectiveness of the nonadiabatic nanoreactor by applying it to benzene, a molecule with rich photochemistry and a wide array of photochemical products. Our study confirms the existence of several types of S₀/S₁ and S₁/S₂ conical intersections which mediate access to a variety of ground state stationary points. We elucidate the connections between conical intersection energy/topography and the resulting photoproduct distribution, which changes smoothly along seam space segments. The exploration is performed with minimal user input, and the protocol requires no previous knowledge of the photochemical behavior of a target molecule. We demonstrate that the nonadiabatic nanoreactor is a valuable tool for the automated exploration of photochemical reactions and their mechanisms.

The nonadiabatic nanoreactor is a tool for automated photochemical reaction discovery that extensively explores intersection seams and links conical intersections to photoproduct distributions.     

Nonadiabatic decay dynamics of 9H-guanine in aqueous solution

The nonadiabatic decay of the biologically relevant guanine tautomer (9H-guanine) in aqueous solution has been investigated by trajectory surface hopping simulations in a quantum mechanical-molecular mechanical (QM-MM) framework. The QM part (9H-guanine) was treated at the semiempirical OM2/MRCI level, while the MM part (water) was described by the TIP3P force field. The optimized geometries for the relevant minima and conical intersections are qualitatively similar for 9H-guanine in the gas phase and in aqueous solution, while there are notable solvent-induced shifts in the computed vertical excitation energies (up to about 0.4 eV). Overall, the results from the static OM2/MRCI-based calculations are in reasonable agreement with the available ab initio and experimental data. The dynamics simulations show ultrafast nonradiative decay for 9H-guanine in water that is even slightly faster than in the gas phase, with time constants of 20 fs and around 0.3 ps for the S(2)→ S(1) and S(1)→ S(0) internal conversions, respectively. They predict a change in the S(1)→ S(0) decay mechanism when going from the gas phase to aqueous solution: the major pathway for 9H-guanine in water involves a conical intersection with an out-of-plane distortion of the carbonyl oxygen atom, which does not play any significant role in the gas phase, where the decay mainly proceeds via two other conical intersections characterized by ring distortions and out-of-plane displacement of the amino group, respectively. Possible reasons for this change in the mechanism are analyzed.     

Nonadiabatic decay dynamics of a benzylidene malononitrile

The photoinduced nonadiabatic decay dynamics of 2-[4-(dimethylamino)benzylidene]malononitrile (DMN) in the gas phase is investigated at the semiempirical OM2/MRCI level using surface hopping simulations. A lifetime of 1.2 ps is predicted for the S(1) state, in accordance with experimental observation. The dominant reaction coordinate is found to be the twisting around the C7═C8 double bond accompanied by pronounced pyramidalization at the C8 atom. Motion along this coordinate leads to the lowest-energy conical intersection (CI(01α)). Several other S(0)/S(1) conical intersections have also been located by full optimization but play no role in the dynamics. The time-resolved fluorescence spectrum of DMN is simulated by computing emission energies and oscillator strengths along the trajectories. It compares well with the experimental spectrum. The use of different active spaces in the OM2/MRCI calculations yields similar results and thus demonstrates their internal consistency.     

Competitive decay at two- and three-state conical intersections in excited-state intramolecular proton transfer

We demonstrate the existence of a simultaneous degeneracy (not required by symmetry) of three electronic states in malonaldehyde. This is one of the first reports of such a triple degeneracy involving S0, S1, and S2 in a molecule with a closed-shell ground state. We further report on a two-state S2/S1 conical intersection which is higher in energy than the three-state intersection, but closer to the Franck-Condon point. First-principles quantum dynamics calculations of the photochemistry after excitation to S2 show that there is a competition between these intersections, with more than half of the population decaying to S1 through the higher energy S2/S1 intersection. Surprisingly, much of the population which makes it to the triple degeneracy point is not funneled directly to S0, but rather remains trapped on S1. We attribute this to the large dimensionality of the branching plane at a three-state intersection (the degeneracy is lifted along at least five distinct molecular displacements).