A theoretical characterization of the photoisomerization channels of 1,2-cyclononadienes on both singlet and triplet potential-energy surfaces

The ground-, (1)(pipi*)-, and (3)(pipi*)-state potential-energy surfaces of 1,2-cyclononadiene and isomeric C(9)H(14) species, as well as 1-methyl-1,2-cyclononadiene and isomeric C(10)H(16) species were all mapped using CASSCF and the 6-31G(d) basis set. Theoretical results were found to be in good agreement with the available experimental observations for both 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene isomerization reactions under singlet and triplet direct or sensitized irradiation. Extremely efficient decay occurs from the first singlet excited state to the ground state through at least three different conical intersections (surface crossings). The first of these crossing points is accessed by a one-bond ring closure. From this conical intersection point (CI-A or CI-C), some possible subsequent ground-state reaction paths have been identified: 1) intramolecular C--H bond insertion to form the bicyclic photoproduct and 2) intramolecular C--H bond insertion to form tricyclic photoproducts. An excited state [1,3]-sigmatropic shift leads to the second conical intersection (CI-B or CI-E), which can give a three-bond cyclononyne species. Besides these, in the singlet photochemical reactions of 1-methyl-1,2-cyclononadiene, excited-state, one allenic C--H bond insertion leads to a third conical intersection (CI-D). Possible ground-state reaction pathways from this structure lead to the formation of a diene photoproduct or to transannular insertion photoproducts. Moreover, in the case of triplet 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene photoisomerization reactions, both chemical reactions will adopt a 1,3-biradical (T(1)/S(0)-1, T(1)/S(0)-2, and T(1)/S(0)-3), which may undergo intersystem crossings leading to the formation of tricyclic or bicyclic photoproducts. The results obtained allow a number of predictions to be made.     

The [2?+?2?+?2] mechanisms of trimerization of three ethynes and monosilaethylenes

The mechanisms of [2?+?2?+?2] reactions of three ethynes and monosilaethylenes to form benzene and 1,3,5-trisilacyclohexane were studied by ab initio MO methods. The reaction mechanisms were analyzed by configuration interaction/localized molecular orbital/CASSCF calculations. Although the [2?+?2?+?2] reaction of ethyne is typically ??homologous?? concerted, that of monosilaethylene is polarized (ionic-cyclic) one-step reaction. In addition, the aromaticity along the intrinsic reaction coordinate pathway was studied using the index of deviation from aromaticity. Although the transition state of trimerization of ethyne does not have an aromatic nature for the ??- and ??-bonds formation system, the crossing point of the ??-bond formation and ??-bond breaking shows an aromatic nature.     

Advances in the metallotropic [1,3]-shift of alkynyl carbenoids

Metal complexes of alkynyl carbenes undergo a [1,3]-bond shift known as a metallotropic shift. Since the discovery of the metallotropic [1,3]-shift of Rh-carbenoid, many more alkynyl carbene complexes with Ti, Cr, Mn, Mo, Ru, W, Re, Pt, and Au have been discovered to show metallotropic shift behavior. This article briefly summarizes the [1,3]-bond shift of alkynyl carbenes and their metal complexes.     

Computational study on the origin of the stereoselectivity for the photochemical denitrogenation of diazabicycloheptene

The origin of the inversion stereoselectivity of housane formation via photochemical nitrogen extrusion of diazabicycloheptene (DBH) has been investigated using reaction path computations and multireference second-order perturbation theory within a CASPT2//CASSCF scheme. We show that the primary photoproduct of the reaction is an exo-axial conformer of the diazenyl diradical ((1) DZ) which displays a cyclopenta-1,3-diyl moiety with a Cs-like structure. (1) DZ is selectively generated via decay at a linear-axial conical intersection, and it is located in a shallow region of the ground state potential energy surface that provides access to five different reaction pathways. Reaction path analysis (including probing with classical trajectories) indicates that production of inverted housane can only occur via impulsive population of an axial-to-equatorial pathway, and it is thus inconsistent with thermal equilibration of the primary (1) DZ conformer. Similarly, according to the same analysis, the decrease of inversion stereoselectivity and even the retention (stereochemical memory effect) observed for suitably substituted DBHs are explained by dynamics effects where the axial-to-equatorial impulsive motion is restrained by the inertia and/or steric hindrance of the substituents. These results shade light on the poorly understood mechanisms that allow a photochemical reaction, in which a large amount of energy is deposited in the reactant by photon absorption, to show a high degree of stereoselectivity.     

Pseudo-bimolecular [2+2] cycloaddition studied by time-resolved photoelectron spectroscopy

The first study of pseudo‐bimolecular cycloaddition reaction dynamics in the gas phase is presented. We used femtosecond time‐resolved photoelectron spectroscopy (TRPES) to study the [2+2] photocycloaddition in the model system pseudo‐gem‐divinyl[2.2]paracyclophane. From X‐ray crystal diffraction measurements we found that the ground‐state molecule can exist in two conformers; a reactive one in which the vinyl groups are immediately situated for [2+2] cycloaddition and a nonreactive conformer in which they point in opposite directions. From the measured S₁ lifetimes we assigned a clear relation between the conformation and the excited‐state reactivity; the reactive conformer has a lifetime of 13 ps, populating the ground state through a conical intersection leading to [2+2] cycloaddition, whereas the nonreactive conformer has a lifetime of 400 ps. Ab initio calculations were performed to locate the relevant conical intersection (CI) and calculate an excited‐state [2+2] cycloaddition reaction path. The interpretation of the results is supported by experimental results on the similar but nonreactive pseudo‐para‐divinyl[2.2]paracyclophane, which has a lifetime of more than 500 ps in the S₁ state.     

Computational study of the mechanisms of the photoisomerization reactions of bicycloalkene

The mechanisms of the photochemical isomerization reactions were investigated theoretically using a model system, bicyclo[4,1,0]hept-3-ene (1), with the CASSCF (six-electron/six-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 bicyclo[4,1,0]hept-3-ene is as follows: reactant → Franck-Condon region → conical intersection → intermediate → transition state → photoproduct. Two reaction paths, which lead to final photoproducts, have been identified: (path I) ring expansion and (path II) ring closure. The former is more favorable than the latter. Also, our theoretical findings strongly indicate that there is a substantial interaction between the cyclopropane moiety and the isolated carbon-carbon double bond in the excited state of 1.     

Optimal control of ultrafast cis-trans photoisomerization of retinal in rhodopsin via a conical intersection

Optimal control simulation is applied to the cis-trans photoisomerization of retinal in rhodopsin within a two-dimensional, two-electronic-state model with a conical intersection [S. Hahn and G. Stock, J. Phys. Chem. B 104, 1146 (2000)]. For this case study, we investigate coherent control mechanisms, in which laser pulses work cooperatively with a conical intersection that acts as a "wave-packet cannon." Optimally designed pulses largely consist of shaping subpulses that prepare a wave packet, which is localized along a reaction coordinate and has little energy in the coupling mode, through multiple electronic transitions. This shaping process is shown to be essential for achieving a high target yield although the envelopes of the calculated pulses depend on the local topography of the potential-energy surfaces around the conical intersection and the choice of target. The control mechanisms are analyzed by considering the motion of reduced wave packets in a nuclear configuration space as well as by snapshots of probability current-density maps.     

Quenching of tryptophan (1)(pi,pi*) fluorescence induced by intramolecular hydrogen abstraction via an aborted decarboxylation mechanism

CASSCF computations show that the hydrogen-transfer-induced fluorescence quenching of the (1)(pi,pi*) excited state of zwitterionic tryptophan occurs in three steps: (1) formation of an intramolecular excited-state complex, (2) hydrogen transfer from the amino acid side chain to the indole chromophore, and (3) radiationless decay through a conical intersection, where the reaction path bifurcates to a photodecarboxylation and a phototautomerization route. We present a general model for fluorescence quenching by hydrogen donors, where the radiationless decay occurs at a conical intersection (real state crossing). At the intersection, the reaction responsible for the quenching is aborted, because the reaction path bifurcates and can proceed forward to the products or backward to the reactants. The position of the intersection along the quenching coordinate depends on the nature of the states and, in turn, affects the formation of photoproducts during the quenching. For a (1)(n,pi*) model system reported earlier (Sinicropi, A.; Pogni, R.; Basosi, R.; Robb, M. A.; Gramlich, G.; Nau, W. M.; Olivucci, M. Angew. Chem., Int. Ed. 2001, 40, 4185-4189), the ground and the excited state of the chromophore are hydrogen acceptors, and the excited-state hydrogen transfer is nonadiabatic and leads directly to the intersection point. There, the hydrogen transfer is aborted, and the reaction can return to the reactant pair or proceed further to the hydrogen-transfer products. In the tryptophan case, the ground state is not a hydrogen acceptor, and the excited-state hydrogen transfer is an adiabatic, sequential proton and electron transfer. The decay to the ground state occurs along a second reaction coordinate associated with decarboxylation of the amino acid side chain and the corresponding aborted conical intersection. The results show that, for (1)(pi,pi*) states, the hydrogen transfer alone is not sufficient to induce the quenching, and explain why fluorescence quenching induced by hydrogen donors is less general for (1)(pi,pi*) than for (1)(n,pi*) states.     

"H/Vinyl" and "Alkyl/Vinyl" conical intersections leading to carbene formation from the excited states of cyclohexene and norbornene

The photochemical [1,2]-shifts leading to carbene intermediates from cyclohexene and norbornene have been studied using CASSCF calculations and a 6-31G* basis set. In each case, a S(1)/S(0) conical intersection hyperline was identified that extends from the region of the reactant excited state to the carbene product. It is traditionally thought that the Rydberg R(pi,3s) state is responsible for carbene formation on photolysis of cyclic alkenes, but these new results indicate an efficient mechanism for carbene formation following excitation to the (1)(pipi*) state. This pathway is essentially barrierless and involves internal conversion to the ground-state surface via conical intersections between the excited zwitterionic valence state and the ground-state surface, similar to those responsible for cis-trans isomerization in ethene and other acyclic alkenes. These results are in excellent agreement with recent experimental data obtained using femtosecond spectroscopy.     

Cyclooctatetraene computational photo- and thermal chemistry: a reactivity model for conjugated hydrocarbons

We use ab initio CASSCF and CASPT2 computations to construct the composite multistate relaxation path relevant to cycloocta-1,3,5,7-tetraene singlet photochemistry. The results show that an efficient population of the dark excited state (S(1)) takes place after ultrafast decay from the spectroscopic excited state (S(2)). A planar D(8)(h)-symmetric minimum represents the collecting point on S(1). Nonadiabatic transitions to S(0) appear to be controlled by two different tetraradical-type conical intersections, which are directly accessible from the S(1) minimum following specific excited-state reaction paths. The higher-energy conical intersection belongs to the same type of intersections previously documented in linear and cyclic conjugated hydrocarbons and features a triangular -(CH)(3)- kink. This point mediates both cis --> trans photoisomerization and cyclopropanation reactions. The lowest energy conical intersection has a boat-shaped structure. This intersection accounts for production of semibullvalene or for double-bond shifting. The mapping of both photochemical and thermal reaction paths (including also Cope rearrangements, valence isomerizations, ring inversions, and double-bond shifting) has allowed us to draw a comprehensive reactivity scheme for cyclooctatetraene, which rationalizes the experimental observations and documents the complex network of photochemical and thermal reaction path interconnections. The factors controlling the selection and accessibility of a number of conjugated hydrocarbon prototype conical intersections and ground-state relaxation channels are discussed.     

A unified perspective on the hydrogen atom transfer and proton-coupled electron transfer mechanisms in terms of topographic features of the ground and excited potential energy surfaces as exemplified by the reaction between phenol and radicals

The relation between the hydrogen atom transfer (HAT) and proton-coupled electron transfer (PCET) mechanisms is discussed and is illustrated by multiconfigurational electronic structure calculations on the ArOH + R(*) --> ArO(*) + RH reactions. The key topographic features of the Born-Oppenheimer potential energy surfaces that determine the predominant reaction mechanism are the conical intersection seam of the two lowest states and reaction saddle points located on the shoulders of this seam. The saddle point corresponds to a crossing of two interacting valence bond states corresponding to the reactant and product bonding patterns, and the conical intersection corresponds to the noninteracting intersection of the same two diabatic states. The locations of mechanistically relevant conical intersection structures and relevant saddle point structures are presented for the reactions between phenol and the N- and O-centered radicals, (*)NH2 and (*)OOCH3. Points on the conical intersection of the ground doublet D0 and first excited doublet D1 states are found to be in close geometric and energetic proximity to the reaction saddle points. In such systems, either the HAT mechanism or both the HAT mechanism and the proton-coupled electron transfer (PCET) mechanism can take place, depending on the relative energetic accessibility of the reaction saddle points and the D0/D1 conical intersection seams. The discussion shows how the two mechanisms are related and how they blend into each other along intermediate reaction paths. The recognition that the saddle point governing the HAT mechanism is on the shoulder of the conical intersection governing the PCET mechanism is used to provide a unified view of the competition between the two mechanisms (and the blending of the two mechanisms) in terms of the prominent and connected features of the potential energy surface, namely the saddle point and the conical intersection. The character of the dual mechanism may be understood in terms of the dominant valence bond configurations of the intersecting states, which are zero-order approximations to the diabatic states.     

Theoretical study of the formation of naphthalene from the radical/π-bond addition between single-ring aromatic hydrocarbons

The experimental investigations performed in the 1960s on the o-benzyne + benzene reaction as well as the more recent studies on reactions involving π-electrons highlight the importance of π-bonding for different combustion processes related to PAH's and soot formation. In the present investigation radical/π-bond addition reactions between single-ring aromatic compounds have been proposed and computationally investigated as possible pathways for the formation of two-ring fused compounds, such as naphthalene, which serve as precursors to soot formation. The computationally generated optimized structures for the stationary points were obtained with uB3LYP/6-311+G(d,p) calculations, while the energies of the optimized complexes were refined using the uCCSD(T) method and the cc-pVDZ basis set. The computations have addressed the relevance of a number of radical/π-bond addition reactions including the singlet benzene + o-benzyne reaction, which leads to formation of naphthalene and acetylene through fragmentation of the benzobicyclo[2,2,2]octatriene intermediate. For this reaction, the high-pressure limit rate constants for the individual elementary reactions involved in the overall process were evaluated using transition state theory analysis. Other radical/π-bond addition reactions studied were between benzene and triplet o-benzyne, between benzene and phenyl radical, and between phenyl radicals, for all of which potential energy surfaces were produced. On the basis of the results of these reaction studies, it was found necessary to propose and subsequently confirm additional, alternative pathways for the formation of the types of PAH compounds found in combustion systems. The potential energy surface for one reaction in particular, the phenyl + phenyl addition, is shown to contain a low-energy channel leading to formation of naphthalene that is energetically comparable to the other examined conventional pathways leading to formation of biphenyl compounds. This channel is the first evidence of a reaction which involves an aromatic radical adding to the nonradical π-bond site of another aromatic radical which leads directly to a fused ring structure.     

Gold-catalyzed intermolecular reactions of propiolic acids with alkenes: [4 + 2] annulation and enyne cross metathesis

A gold-catalyzed intermolecular reaction of propiolic acids with alkenes led to a [4 + 2] annulation or enyne cross metathesis. The [4 + 2] annulation proceeds with net cis-addition with respect to alkenes and provides an expedient route to α,β-unsaturated δ-lactones, for which preliminary asymmetric reactions were also demonstrated. For 1,2-disubstituted alkenes, unprecedented enyne cross metathesis occurred to give 1,3-dienes in a completely stereospecific fashion. DFT calculations and experiments indicated that the cyclobutene derivatives are not viable intermediates and that the steric interactions during concerted σ-bond rearrangements are responsible for the observed unique stereospecificity.     

Thermal [3,3]-rearrangement of 1,1-disubstituted allyl carboxylates: lone pair participation and the geminal bond participation

The diastereoselectivity of the [3,3]-rearrangement of 1,1-disubsstituted allyl carboxylates was studied. In this heteroatom-containing system, the transition state has a boat-like transition structure (TS) because of the participation of the lone pairs and the secondary orbital interaction. Although the TS for the [1,3]-rearrangement has a far higher barrier, it does not proceed in the usual antarafacial manner due to the cyclic orbital interaction among two lone pairs of the carboxylate and the allylic lumo. In conjunction with the geminal bond participation, delocalization to the σ-bond at the Z-position shows a bonding character in the transition state of the [3,3]-rearrangement. Therefore, we predicted that a more electron-withdrawing σ-bond prefers the Z-position in the product. We designed the 1,1-disubstituted substrates with trimethylsilyl and pentyl groups, and found that the trimethylsilyl group prefers the Z-position despite its steric bulkiness. We confirmed our prediction by experimentation. This Z-selectivity was improved when a trimethylgermyl group was used.     

Photoreactivity of furocoumarins and DNA in PUVA therapy: formation of psoralen-thymine adducts

The mechanism of the [2+2] cycloaddition photoreaction of psoralen and a DNA nucleobase, thymine, cornerstone of the furocoumarin-based PUVA (psoralen+UVA radiation) phototherapy, has been studied by the quantum-chemical multiconfigurational CASPT2 method. Triplet- and singlet-mediated mono- and diadduct formations have been determined to take place via singlet-triplet crossings and conical intersections, correlated with the initially promoted triplet or singlet states in different possible reactive orientations. Pyrone-side monoadducts are suggested to be formed in the triplet manifold of the system, and to be less prone to yield diadducts because of the properties of the monoadduct lowest triplet state and the minor accessibility of its excited singlet states. Furan-side monoadducts are better produced in the singlet manifold after reaching a conical intersection with the ground state of the system. From there, the absorption of a second photon would in this case trigger the formation of the diadduct. The proposed mechanisms enable rationalizing the phototherapeutic behavior of several furocoumarins.     

Photophysics of (1-butyl-4-(1H-inden-1-ylidene)-1,4-dihydropyridine (BIDP): an experimental test for conical intersections

Fluorescence experiments on (1-butyl-4-(1H-inden-1-ylidene)-1,4-dihydropyridine (BIDP) are reported in liquid and glassy solutions. The data indicate a fast decay in the fluid nonpolar, nonprotic solutions (decay times approximately 10(-12) s) and rapid but considerably slower decay in polar ones. In frozen solutions (polar and nonpolar), the fluorescence quantum yield is much higher (near 0.5 and around 0.1 in polar and nonpolar glasses, respectively). The rapid nonradiative transitions in fluid solutions are assigned to internal conversion in both solvent classes, as intersystem crossing is much slower and no net reaction is observed. These results are in agreement with predictions made for the closely related (in terms of electronic structure) but simpler molecule cyclopentadienyl-1,4-dihydropyridine (CPDHP) for which an S1/S0 conical intersection was recently proposed [Int. J. Quant. Chem. 2005, 102, 961]. The crossing of the two lowest singlet states is calculated to vanish in polar solvents such as methyl cyanide, leading to longer lifetime of S1 of CPDHP. As BIDP has a very similar electronic structure, the model predicts a corresponding change in this larger molecule. The strong fluorescence observed in the glassy environments is rationalized by the hindering of the internal torsion required to reach the geometry of the conical intersection.     

The mechanism of photochemical 1,3-silyl migration of allylsilane

The photochemical reaction mechanisms of model compounds for 4-tert-butyl-1-(4-phenylphenyl)-1-(1,1-dimethylallyl)silacyclohexane are investigated using a complete active space comprised of six electrons in six orbitals with the standard 6-31G(d) basis set. It is concluded that the stereochemistry in the photochemical 1,3-silyl migrations of allylsilanes has a retention preference, in accord with the Woodward-Hoffmann rules. The calculated conical intersection (CI) structure suggests a dissociation path to radicals in addition to a 1,3-shift path. The bulkiness and rigidness of a silacyclohexane moiety does not affect the stereochemistry, but a slightly elongated Si-C bond length in the CI structure would promote the dissociation path.     

Nonradiative decay mechanisms of the biologically relevant tautomer of guanine

We present the excited-state potential energy profiles of the biologically relevant 9H-keto-amino tautomer of guanine with respect to the radiationless decay via the out-of-plane deformation of the six-membered ring as well as the dissociation of NH bonds. The CASPT2//CASSCF method is employed for the reaction-path calculations. The reaction path for the out-of-plane deformation in the (1)pi pi* state leads in a barrierless way to a conical intersection with the electronic ground state. For the NH dissociation via the (1)pi sigma* state, the 9H-keto-amino tautomer is shown to have lower energy barriers than the 7H tautomers which we have studied recently. These two radiationless decay mechanisms explain the unexpected missing of the biologically relevant form in the resonant two-photon ionization spectrum of guanine in a supersonic jet. It is suggested that these ultrafast deactivation processes provide the biologically relevant tautomer of guanine with a high degree of photostability.     

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.     

Photochemical reactivity of 2-vinylbiphenyl and 2-vinyl-1,3-terphenyl: the balance between nonadiabatic and adiabatic photocyclization

A mechanism for the photochemical conversion of 2-vinyl-1,3-terphenyl to 8,9a-dihydrophenanthrene (Lewis, F. D.; Zuo, X.; Gevorgyan, V.; Rubin, M. J. Am. Chem. Soc. 2002, 124, 13664-13665) is presented in this study, based on ab initio restricted active space self-consistent field calculations and a molecular mechanics-valence bond dynamics simulation of a model system: the syn isomer of 2-vinylbiphenyl. An extended crossing seam between the ground and first excited electronic states was found to be largely responsible for the efficient photocyclization of the photochemically active syn isomer. This mechanism is nonadiabatic in nature, with an excited-state reaction pathway approaching the crossing region during the initial stage of cyclization. Dynamics simulation shows that this seam is easily accessible by vibrational motion in the branching space, once a small barrier is passed on the S1 excited-state potential energy surface. Ultrafast radiationless decay to the ground state then follows, and the cyclization is completed on this surface. A second possible mechanism was identified, which involves complete adiabatic cyclization on the S1 surface, with decay to the ground state (at a different conical intersection) only taking place once the product is formed. Thus, there is a competition between these two mechanisms-nonadiabatic and adiabatic-governed by the dynamics of the system. A large quantum yield is predicted for the photocyclization of the syn isomer of 2-vinylbiphenyl and 2-vinyl-1,3-terphenyl, in agreement with experimental observations.