Manipulating nonadiabatic conical intersection dynamics by optical cavities

Optical cavities hold great promise to manipulate and control the photochemistry of molecules. We demonstrate how molecular photochemical processes can be manipulated by strong light–matter coupling. For a molecule with an inherent conical intersection, optical cavities can induce significant changes in the nonadiabatic dynamics by either splitting the pristine conical intersections into two novel polaritonic conical intersections or by creating light-induced avoided crossings in the polaritonic surfaces. This is demonstrated by exact real-time quantum dynamics simulations of a three-state two-mode model of pyrazine strongly coupled to a single cavity photon mode. We further explore the effects of external environments through dissipative polaritonic dynamics computed using the hierarchical equation of motion method. We find that cavity-controlled photochemistry can be immune to external environments. We also demonstrate that the polariton-induced changes in the dynamics can be monitored by transient absorption spectroscopy.

Optical cavities hold great promise to manipulate and control the photochemistry of molecules.     

A novel conical intersection topography and its consequences: the 1, 2 2A conical intersection seam of the vinoxy radical

A region of the 1, 2 2A seam of accidental conical intersections in the vinoxy radical exhibits a novel topography which has important consequences for both upper-state to lower-state and lower-state to upper-state nonadiabatic transitions. The impact of this topography on these nonadiabatic transitions is described. We also considered the possibility that this conical intersection seam provides a dynamical bottleneck to the photodissociation of vinoxy to H+ketene by facilitating nonadiabatic recrossing. Our analysis of the conical topographies and the proximity of the conical intersections to the transition state for dissociation to H+ketene does not support nonadiabatic recrossing as an effective dynamical bottleneck blocking the H+ketene channel.     

Excited‐State Intramolecular Proton Transfer in a Blue Fluorescence Chromophore Induces Dual Emission

Compared with green fluorescence protein (GFP) chromophores, the recently synthesized blue fluorescence protein (BFP) chromophore variant presents intriguing photochemical properties, for example, dual fluorescence emission, enhanced fluorescence quantum yield, and ultra‐slow excited‐state intramolecular proton transfer (ESIPT; J. Phys. Chem. Lett., 2014 , 5, 92); however, its photochemical mechanism is still elusive. Herein we have employed the CASSCF and CASPT2 methods to study the mechanistic photochemistry of a truncated BFP chromophore variant in the S₀ and S₁ states. Based on the optimized minima, conical intersections, and minimum‐energy paths (ESIPT, photoisomerization, and deactivation), we have found that the system has two competitive S₁ relaxation pathways from the Franck–Condon point of the BFP chromophore variant. One is the ESIPT path to generate an S₁ tautomer that exhibits a large Stokes shift in experiments. The generated S₁ tautomer can further evolve toward the nearby S₁/S₀ conical intersection and then jumps down to the S₀ state. The other is the photoisomerization path along the rotation of the central double bond. Along this path, the S₁ system runs into an S₁/S₀ conical intersection region and eventually hops to the S₀ state. The two energetically allowed S₁ excited‐state deactivation pathways are responsible for the in‐part loss of fluorescence quantum yield. The considerable S₁ ESIPT barrier and the sizable barriers that separate the S₁ tautomers from the S₁/S₀ conical intersections make these two tautomers establish a kinetic equilibrium in the S₁ state, which thus results in dual fluorescence emission.     

Curl equations as substratum for the derivation of the electronic nonadiabatic coupling terms

The idea to derive the nonadiabatic coupling terms by solving the Curl equations (Avery, J.; Baer, M.; Billing, G. D. Mol Phys 2002, 100, 1011) is extended to a three‐state system where the first and second states form one conical intersection, i.e., τ₁₂ and the second and the third states form another conical intersection, i.e., τ₂₃. Whereas the two‐state Curl equations form a set of linear differential equations, the extension to a three‐state system not only increases the number of equations but also leads to nonlinear terms. In the present study, we developed a perturbative scheme, which guarantees convergence if the overlap between the two interacting conical intersections is not too strong. Among other things, we also revealed that the nonadiabatic coupling term between the first and third states, i.e., τ₁₃ (such interactions do not originate from conical intersection) is formed due to the interaction between τ₁₂ and τ₂₃. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

A computational study of the mechanisms of the photoisomerization reactions of monocyclic and bicyclic olefins

The mechanisms of the photochemical isomerization reactions were investigated theoretically using a model system of cyclohexene (1), cycloheptene (2), norbornene (3), and two bicyclic olefins (4 and 5) using the CASSCF (six-electron/six-orbital active space) and MP2-CAS methods with the 6-311(d,p) basis set. The structures of the conical intersections, which play a decisive role in such photoisomerizations, were obtained. The intermediates and transition structures of the ground state were also calculated to assist in providing a qualitative explanation of the reaction pathways. Two photoreaction pathways were examined in the present work. The first can produce a photoproduct with an extra ring. The other can yield a photoproduct with a smaller ring with an external double bond. Both pathways involve cyclic carbene intermediates. Also, our model investigations suggest that both reaction pathways follow a similar photochemical pattern as follows: reactant → Franck-Condon region → conical intersection → cyclic carbene intermediate → transition state → photoproduct. Moreover, these two reaction pathways can compete with each other since the energetics of their conical intersection points are quite similar. Our present theoretical results agree with the available experimental observations.     

Photochemical isomerization reactions of cyanopyrroles: a theoretical study

The mechanisms of the photochemical isomerization reactions were investigated theoretically using a model system of 2-cyanopyrrole and 2-cyano-5-methylpyrrole with the CASSCF (eight-electron/seven-orbital active space) and MP2-CAS methods and the 6-311(d,p) basis set. The structures of the conical intersections, which play a decisive role in such phototranspositions, were obtained. The intermediates and transition structures of the ground state were also calculated to assist in providing a qualitative explanation of the reaction pathways. Our model investigations suggest that the preferred reaction route for the cyanopyrroles is as follows: reactant --> Franck-Condon region --> conical intersection --> photoproduct. In particular, the conical intersection mechanism found in this work gives a better explanation than the previously proposed internal cyclization-isomerization mechanism and supports the experimental observations. In addition, we suggest a simple p-pi orbital topology model, which can be used as a guide tool to predict the location at which conical intersections are likely to occur, as well as the conformations of the phototransposition products of various heterocycles.     

飞秒时间分辨的光电子影像技术研究2-氯吡啶超快内转换动力学

利用时间分辨的飞秒光电子影像技术结合时间分辨的质谱技术, 研究了2-氯吡啶分子激发态的超快过程. 实时观察到了2-氯吡啶分子第二激发态(S₂)向第一激发态(S₁)高振动能级的的超快内转换过程,该内转换的时间常数为(162±5)fs. 实验结果表明, 通过S₂/S₀的锥形交叉衰减到基态的衰减通道也是退布居的重要通道, 其时间尺度为(5.5±0.3) ps.     

Internal conversion of the anionic GFP chromophore: in and out of the I-twisted S1/S0 conical intersection seam

The functional diversity of the green fluorescent protein (GFP) family is intimately connected to the interplay between competing photo-induced transformations of the chromophore motif, anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI⁻). Its ability to undergo Z/E-isomerization is of particular importance for super-resolution microscopy and emerging opportunities in optogenetics. Yet, key dynamical features of the underlying internal conversion process in the native HBDI⁻ chromophore remain largely elusive. We investigate the intrinsic excited-state behavior of isolated HBDI⁻ to resolve competing decay pathways and map out the factors governing efficiency and the stereochemical outcome of photoisomerization. Based on non-adiabatic dynamics simulations, we demonstrate that non-selective progress along the two bridge-torsional (i.e., phenolate, P, or imidazolinone, I) pathways accounts for the three decay constants reported experimentally, leading to competing ultrafast relaxation primarily along the I-twisted pathway and S₁ trapping along the P-torsion. The majority of the population (∼70%) is transferred to S₀ in the vicinity of two approximately enantiomeric minima on the I-twisted intersection seam (MECI-Is). Despite their sloped, reactant-biased topographies (suggesting low photoproduct yields), we find that decay through these intersections leads to products with a surprisingly high quantum yield of ∼30%. This demonstrates that E-isomer generation results at least in part from direct isomerization on the excited state. A photoisomerization committor analysis reveals a difference in intrinsic photoreactivity of the two MECI-Is and that the observed photoisomerization is the combined result of two effects: early, non-statistical dynamics around the less reactive intersection followed by later, near-statistical behavior around the more reactive MECI-I. Our work offers new insight into internal conversion of HBDI⁻ that both establishes the intrinsic properties of the chromophore and enlightens principles for the design of chromophore derivatives and protein variants with improved photoswitching properties.

The Z–E photoisomerization quantum yield of the HBDI⁻ chromophore is a result of early, non-statistical dynamics around a less reactive I-twisted intersection and later, statistical behavior around the more reactive, near-enantiomeric counterpart.     

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.     

Exploring molecular complexity: conical intersections and NH3 photodissociation

The role of conical intersections in the photodissociation of the A 1A2" state of NH3 is investigated using extended atomic basis sets and a configuration state function expansion of approximately 8.5 million terms. A previously unknown portion of the 1 1A-2 1A seam of conical intersections with only C(s) symmetry is located. This portion of the seam is readily accessible from the equilibrium geometry of the A 1A2" state. These conical intersections are expected to play a role in the competition between adiabatic and nonadiabatic pathways for NH3(A 1A2") photodissociation.     

Photoinduced multi-mode quantum dynamics of pyrrole at the – conical intersections

The photoinduced dynamics of pyrrole at the – S₀ and – S₀ conical intersections has been investigated by multi-mode time-dependent quantum wave-packet calculations. Diabatic potential-energy surfaces have been constructed for both conical intersection using accurate multi-reference ab initio electronic-structure calculations. In addition to the NH stretching coordinate, the three (four) symmetry-allowed coupling modes of A₂ (B₁) symmetry have been considered for the – S₀ (– S₀) conical intersections. Wave-packet dynamics calculations have been performed for three-dimensional models, taking account of the two dominant coupling modes of each conical intersection. The electronic population-transfer processes at the conical intersections, the branching ratio for the dissociation to the ground and excited states of the pyrrolyl radical, and their dependence on the initial preparation of the system have been investigated. It is shown that the excitation of the NH stretching mode strongly enhances the photodissociation rate, while the excitation of the strongest coupling mode has a pronounced effect on the branching ratio of the photodissociation process. Although the inclusion of the second (weaker) coupling mode has little effect on the electronic population dynamics, it leads to interesting changes of the nodal pattern of the wave packet at the conical intersections. The calculations provide insight into the effect of the multiple coupling modes on the process of direct photodissociation through a conical intersection.     

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.     

CASCSF study on the photochemical transposition reactions of pyrazines

Several reaction pathways for the photochemical transformations of methyl-substituted pyrazine in its first excited state 1(pi --> pi*) have been determined using the CASSCF (six-orbital/six-electron active space) and MP2-CAS methods with the 6-311G(d) basis set. Our model investigations suggest that conical intersections play a crucial role in the photoisomerization of pyrazines. Moreover, the present theoretical findings indicate that all of the photoisomerizations of pyrazines adopt the same reaction path as follows: pyrazine --> Franck-Condon region --> conical intersection --> pyrimidine. That is, although an excited-state pyrazine molecule can initiate a phototransposition process easily, this process can be completed on the ground-state potential energy surface after passage through a conical intersection where a fast, radiationless decay is possible. The existence of these nonadiabatic reaction pathways is consistent with the available experimental observations of the photochemistry and photophysics of pyrazine and its methyl derivatives. In the present work, we propose a simple p-pi orbital topology model, which can be used as a diagnostic tool to predict the location of the conical intersections, as well as the geometries of the phototransposition products of various heterocycles.     

Photoisomerization of Arylazopyrazole Photoswitches: Stereospecific Excited‐State Relaxation

Electronic structure calculations and nonadiabatic dynamics simulations (more than 2000 trajectories) are used to explore the Z–E photoisomerization mechanism and excited‐state decay dynamics of two arylazopyrazole photoswitches. Two chiral S₁/S₀ conical intersections with associated enantiomeric S₁ relaxation paths that are barrierless and efficient (timescale of ca. 50 fs) were found. For the parent arylazopyrazole (Z8) both paths contribute evenly to the S₁ excited‐state decay, whereas for the dimethyl derivative (Z11) each of the two chiral cis minima decays almost exclusively through one specific enantiomeric S₁ relaxation path. To our knowledge, the Z11 arylazopyrazole is thus the first example for nearly stereospecific unidirectional excited‐state relaxation.     

Ultrafast Dynamics Through Conical Intersections in 2,6-dimethylpyridine Studied with Time-resolved Photoelectron Imaging

The ultrafast dynamics through conical intersections in 2,6-dimethylpyridine has been stud-ied by femtosecond time-resolved photoelectron imaging coupled with time-resolved mass spectroscopy. Upon absorption of 266 nm pump laser, 2,6-dimethylpyridine is excited to the S₂ state with a ππ^* character from S₀state. The time evolution of the parent ion sig-nals consists of two exponential decays. One is a fast component on a timescale of 635 fs and the other is a slow component with a timescale of 4.37 ps. Time-dependent photo-electron angular distributions and energy-resolved photoelectron spectroscopy are extracted from time-resolved photoelectron imaging and provide the evolutive information of S₂ state. In brief, the ultrafast component is a population transfer from S₂ to S₁ through the S₂/S₁ conical intersections, the slow component is attributed to simultaneous IC from the S₂ state and the higher vibrational levels of S₁ state to S₀ state, which involves the coupling of S₂/S₀ and S₁/S₀ conical intersections. Additionally, the observed ultrafast S₂→S₁ transition occurs only with an 18% branching ratio.     

Nonradiative Excited-State Decay via Conical Intersection in Graphene Nanostructures

Chemical groups are known to tune the luminescent efficiencies of graphene-related nanomaterials, but some species, including the epoxide group (−COC−), are suspected to act as emission-quenching sites. Herein, by performing nonadiabatic excited-state dynamics simulations, we reveal a fast (within 300 fs) nonradiative excited-state decay of a graphene epoxide nanostructure from the lowest excited singlet (S₁) state to the ground (S₀) state via a conical intersection (CI), at which the energy difference between the S₁ and S₀ states is approximately zero. This CI is induced after breaking one C−O bond at the −COC− moiety during excited-state structural relaxation. This study ascertains the role of epoxide groups in inducing the nonradiative recombination of the excited electron-hole, providing important insights into the CI-promoted nonradiative de-excitations and the luminescence tuning of relevant materials. In addition, it shows the feasibility of utilizing nonadiabatic excited-state dynamics simulations to investigate the photophysical processes of the excited states of graphene nanomaterials.     

Solvent dielectric delimited nitro–nitrito photorearrangement in a perylenediimide derivative

The discovery of vibrant excited-state dynamics and distinctive photochemistry has established nitrated polycyclic aromatic hydrocarbons as an exhilarating class of organic compounds. Herein, we report the atypical photorearrangement of nitro-perylenediimide (NO₂-PDI) to nitrito-perylenediimide (ONO-PDI), triggered by visible-light excitation and giving rise to linkage isomers in the polar aprotic solvent acetonitrile. ONO-PDI has been isolated and unambiguously characterized using standard spectroscopic, spectrometric, and elemental composition techniques. Although nitritoaromatic compounds are conventionally considered to be crucial intermediates in the photodissociation of nitroaromatics, experimental evidence for this has not been observed heretofore. Ultrafast transient absorption spectroscopy combined with computational investigations revealed the prominence of a conformationally relaxed singlet excited-state (S^CR₁) of NO₂-PDI in the photoisomerization pathway. Theoretical transition state (TS) analysis indicated the presence of a six-membered cyclic TS, which is pivotal in connecting the S^CR₁ state to the photoproduct state. This article addresses prevailing knowledge gaps in the field of organic linkage isomers and provides a comprehensive understanding of the unprecedented photoisomerization mechanism operating in the case of NO₂-PDI.

The unprecedented photorearrangement of nitro-perylenediimide (NO₂-PDI) to nitrito-perylenediimide (ONO-PDI) is shown to occur through a cyclic six-membered transition state triggered by visible-light excitation.     

Conical intersections in solution: formulation, algorithm, and implementation with combined quantum mechanics/molecular mechanics method

The significance of conical intersections in photophysics, photochemistry, and photodissociation of polyatomic molecules in gas phase has been demonstrated by numerous experimental and theoretical studies. Optimization of conical intersections of small- and medium-size molecules in gas phase has currently become a routine optimization process, as it has been implemented in many electronic structure packages. However, optimization of conical intersections of small- and medium-size molecules in solution or macromolecules remains inefficient, even poorly defined, due to large number of degrees of freedom and costly evaluations of gradient difference and nonadiabatic coupling vectors. In this work, based on the sequential quantum mechanics and molecular mechanics (QM/MM) and QM/MM-minimum free energy path methods, we have designed two conical intersection optimization methods for small- and medium-size molecules in solution or macromolecules. The first one is sequential QM conical intersection optimization and MM minimization for potential energy surfaces; the second one is sequential QM conical intersection optimization and MM sampling for potential of mean force surfaces, i.e., free energy surfaces. In such methods, the region where electronic structures change remarkably is placed into the QM subsystem, while the rest of the system is placed into the MM subsystem; thus, dimensionalities of gradient difference and nonadiabatic coupling vectors are decreased due to the relatively small QM subsystem. Furthermore, in comparison with the concurrent optimization scheme, sequential QM conical intersection optimization and MM minimization or sampling reduce the number of evaluations of gradient difference and nonadiabatic coupling vectors because these vectors need to be calculated only when the QM subsystem moves, independent of the MM minimization or sampling. Taken together, costly evaluations of gradient difference and nonadiabatic coupling vectors in solution or macromolecules can be reduced significantly. Test optimizations of conical intersections of cyclopropanone and acetaldehyde in aqueous solution have been carried out successfully.