Probing highly efficient photoisomerization of a bridged azobenzene by a combination of CASPT2//CASSCF calculation with semiclassical dynamics simulation

Mechanism of phototriggered isomerization of azobenzene and its derivatives is of broad interest. In this paper, the S(0) and S(1) potential energy surfaces of the ethylene-bridged azobenzene (1) that was recently reported to have highly efficient photoisomerization were determined by ab initio electronic structure calculations at different levels and further investigated by a semiclassical dynamics simulation. Unlike azobenzene, the cis isomer of 1 was found to be more stable than the trans isomer, consistent with the experimental observation. The thermal isomerization between cis and trans isomers proceeds via an inversion mechanism with a high barrier. Interestingly, only one minimum-energy conical intersection was determined between the S(0) and S(1) states (CI) for both cis → trans and trans → cis photoisomerization processes and confirmed to act as the S(1) → S(0) decay funnel. The S(1) state lifetime is ～30 fs for the trans isomer, while that for the cis isomer is much longer, due to a redistribution of the initial excitation energies. The S(1) relaxation dynamics investigated here provides a good account for the higher efficiency observed experimentally for the trans → cis photoisomerization than the reverse process. Once the system decays to the S(0) state via CI, formation of the trans product occurs as the downhill motion on the S(0) surface, while formation of the cis isomer needs to overcome small barriers on the pathways of the azo-moiety isomerization and rotation of the phenyl ring. These features support the larger experimental quantum yield for the cis → trans photoisomerization than the trans → cis process.     

Nonadiabatic dynamics of a truncated indigo model

Indigo (1) is stable when exposed to ultraviolet light. We employ electronic structure calculations and nonadiabatic trajectory surface-hopping dynamics simulations to study the photoinduced processes and the photoprotection mechanism of an indigo model, bispyrroleindigo (2). Consistent with recent static ab initio calculations on 1 and 2 (Phys. Chem. Chem. Phys., 2011, 13, 1618), we find an efficient deactivation process that proceeds as follows. After vertical photoexcitation, the S(1)(ππ*) state undergoes an essentially barrierless intramolecular single proton transfer and relaxes to the minimum of an S(1) tautomer, which is structurally and energetically close to a nearby conical intersection that acts as a funnel to the S(0) state; after this internal conversion, a reverse single hydrogen transfer leads back to the equilibrium structure of the most stable S(0) tautomer. This deactivation process is completely dominant in our semiempirical OM2/MRCI nonadiabatic dynamics simulations. The other two mechanisms considered previously, namely excited-state intramolecular double proton transfer and trans-cis double bond isomerization, are not seen in any of the 325 trajectories of the present surface-hopping simulations. On the basis of the computed time-dependent populations of the S(1) state, we estimate an S(1) lifetime of about 700 fs for 2 in the gas phase.     

Molecular mechanisms of the photostability of indigo

The photophysics of indigo as well as of bispyrroleindigo, the basic chromophore of indigo, has been investigated with ab initio electronic-structure calculations. Vertical electronic excitation energies and excited-state potential-energy profiles have been calculated with the CASSCF, CASPT2 and CC2 methods. The calculations reveal that indigo and bispyrroleindigo undergo intramolecular single-proton transfer between adjacent N-H and C=O groups in the (1)ππ* excited state. The nearly barrierless proton transfer provides the pathway for a very efficient deactivation of the (1)ππ* state via a conical intersection with the ground state. While a low-lying S(1)-S(0) conical intersection exists also after double-proton transfer, the latter reaction path exhibits a much higher barrier. The reaction path for trans→cis photoisomerization via the twisting of the central C=C bond has been investigated for bispyrroleindigo. It has been found that the twisting of the central C=C bond is unlikely to play a role in the photochemistry of indigo, because of a large potential-energy barrier and a rather high energy of the S(1)-S(0) conical intersection of the twisted structure. These findings indicate that the exceptional photostability of indigo is the result of rapid internal conversion via intramolecular single-proton transfer, combined with the absence of a low-barrier reaction path for the generation of the cis isomer via trans→cis photoisomerization.     

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.     

Ab initio trajectory surface-hopping study on ultrafast deactivation process of thiophene

The ultrafast S(1)((1)ππ*) → S(0) deactivation process of thiophene in the gas phase has been simulated with the complete active space self-consistent field (CASSCF) based fewest switch surface hopping method. It was found that most of the calculated trajectories (～80%) decay to the ground state (S(0)) with an averaged time constant of 65 ± 5 fs. This is in good agreement with the experimental value of about 80 fs. Two conical intersections were determined to be responsible for the ultrafast S(1)((1)ππ*) → S(0) internal conversion process. After thiophene is excited to the S(1)((1)ππ*) state in the Franck-Condon region, it quickly relaxes to the minimum of the S(1)((1)ππ*) state, then overcomes a small barrier near the conical intersection (CI((1)ππ*/(1)πσ*)), and eventually arrives at the minimum of one C-S bond fission (S(1)((1)πσ*)). In the vicinity of this minimum, the conical intersection (CI((1)πσ*/S(0))) funnels the electron population to the ground state (S(0)), completing the ultrafast S(1)((1)ππ*) → S(0) internal conversion process. This decay mechanism matches well with previous experimental and theoretical studies.     

Vibronic spectra of jet-cooled 2-aminopurine·H2O clusters studied by UV resonant two-photon ionization spectroscopy and quantum chemical calculations

For understanding the major- and minor-groove hydration patterns of DNAs and RNAs, it is important to understand the local solvation of individual nucleobases at the molecular level. We have investigated the 2-aminopurine·H(2)O monohydrate by two-color resonant two-photon ionization and UV/UV hole-burning spectroscopies, which reveal two isomers, denoted A and B. The electronic spectral shift δν of the S(1) ← S(0) transition relative to bare 9H-2-aminopurine (9H-2AP) is small for isomer A (-70 cm(-1)), while that of isomer B is much larger (δν = -889 cm(-1)). B3LYP geometry optimizations with the TZVP basis set predict four cluster isomers, of which three are doubly H-bonded, with H(2)O acting as an acceptor to a N-H or -NH2 group and as a donor to either of the pyrimidine N sites. The "sugar-edge" isomer A is calculated to be the most stable form with binding energy D(e) = 56.4 kJ/mol. Isomers B and C are H-bonded between the -NH2 group and pyrimidine moieties and are 2.5 and 6.9 kJ/mol less stable, respectively. Time-dependent (TD) B3LYP/TZVP calculations predict the adiabatic energies of the lowest (1)ππ* states of A and B in excellent agreement with the observed 0(0)(0) bands; also, the relative intensities of the A and B origin bands agree well with the calculated S(0) state relative energies. This allows unequivocal identification of the isomers. The R2PI spectra of 9H-2AP and of isomer A exhibit intense low-frequency out-of-plane overtone and combination bands, which is interpreted as a coupling of the optically excited (1)ππ* state to the lower-lying (1)nπ* dark state. In contrast, these overtone and combination bands are much weaker for isomer B, implying that the (1)ππ* state of B is planar and decoupled from the (1)nπ* state. These observations agree with the calculations, which predict the (1)nπ* above the (1)ππ* state for isomer B but below the (1)ππ* for both 9H-2AP and isomer A.     

Isomerization of dewar benzenes involving electron donor-acceptor exciplexes

Valence photoisomerization of hexamethyl (Dewar benzene) (HMDB) is sensitized by aromatic singlet photosensitizers 1,4-dicyanobenzene, 1-cyanonaphthalene, 9-cyanoanthracene, and 9,10-dicyanoanthracene with a limiting quantum efficiency of 1.0 in cyclohexane solvent. Quenching of the fluorescence of the aromatic sensitizers leads to exciplex emission which is identical to that obtained by quenching with the isomer, hexamethylbenzene (HMB). The emission is identified as HMB exciplex emission on the basis of relative lifetime and dual quenching experiments. The relative yield of HMDB-derived (“adiabatic”) emission is 20–50% depending on the excitation energy of the HMB exciplex product. Neither biacetyl singlet or triplet nor 1-cyanonaphthalene triplet photosensitization is successful in bringing about isomerization of HMDB. Dimethyl 1,4,5,6-tetramethylbicy-clo[2.20]hexa-2,5-diene-2, 3-dicarboxylate undergoes valence isomerization on quenching electron donor fluorophores, with a quantum efficiency of 0.2. The aromatic valence isomer is not produced in an excited state in this case. Factors which govern the efficiency of adiabatic and diabatic isomerization of the Dewar benzenes are discussed, including sensitizer redox properties, configuration, and multiplicity, the excitation energy and binding characteristics of exciplexes, and the Dewar benzene substituent pattern.     

Photoisomerization dynamics of N-1-methyl-2-(tolylazo) imidazole and the effect of complexation with Cu(ii)

Azo-compounds containing an imidazole moiety have the potential to photoregulate biofunctions, such as gene-expression and enzymatic action. Photoinduced isomerization of the azo-backbone is the vital process for such applications, but the photoisomerization dynamics of azo-imidazole compounds has not been well explored. We investigated the photoisomerization dynamics of trans-N-1-methyl-2-(tolylazo) imidazole (trans-MTAI) using femtosecond transient absorption spectroscopy following photoexcitation to the S(2) state. Time evolution of the transient electronic spectra and the global analysis of the temporal profiles reveal a three state relaxation (S(2) → S(1) → S(0)) process in different kinds of solvents. The lifetime of the S(2) state is independent of the viscosity of solvent, whereas that of the S(1) state becomes longer with increasing solvent viscosity. This observation clearly indicates that the large amplitude motion that leads to the trans → cis isomerization occurs only in the S(1) state and relaxation of the S(2) state is not associated with the isomerization process. We have also investigated the excited state dynamics of [Cu(trans-MTAI)(2)]Cl(2) to examine the effect of complexation with the metal ion on the isomerization dynamics of trans-MTAI. It is observed that the photoinduced isomerization of the azo-backbone in trans-MTAI is completely inhibited upon complexation with Cu(ii).     

Trans-cis photoisomerization of azobenzene by n →π^＊ excitation： A semiclassical dynamics study

A realistic dynamics simulation study is reported for the trans-cis photoisomerization of azobenzene triggered by the n →π^＊ excitation and the results show that the formation ofcis isomer follows the rotational motion around the N=N bond. The simulation find that the CNN bond angle bending vibrations also play a significant role in the vibronic coupling between the HOMO and LUMO, which essentially leads a nonadiabatic transition of the molecule to the electronic ground state.     

Nonadiabatic molecular dynamics study of the cis-trans photoisomerization of azobenzene excited to the S1 state

Ab initio nonadiabatic dynamics simulations of cis-to-trans isomerization of azobenzene upon S(1) (n-π*) excitation are carried out employing the fewest-switches surface hopping method. Azobenzene photoisomerization occurs purely as a rotational motion of the central CNNC moiety. Two nonequivalent rotational pathways corresponding to clockwise or counterclockwise rotation are available. The course of the rotational motion is strongly dependent on the initial conditions. The internal conversion occurs via an S(0)/S(1) crossing seam located near the midpoint of both of these rotational pathways. Based on statistical analysis, it is shown that the occurrence of one or other pathway can be completely controlled by selecting adequate initial conditions.     

Time-dependent partitioning theory of the control of radiationless transitions in 24-mode pyrazine

We consider the control of internal conversion between the S(2)((1)B(2u)) excited electronic state of pyrazine and the S(1)((1)B(3u)) state. The study is performed both during and after the femtosecond excitation of the ground electronic state S(0)((1)A(g)) to form the S(2) state. The dynamics is examined using the newly developed "effective modes" technique which enables the full computation of quantum dynamics in multi-dimensional spaces. Using this technique, we also investigate the coherent control of population transfer from S(0) to the S(2) and S(1) electronic states. We find that the use of shaped laser pulses enables a significant delay of the internal conversion. For example, after 60 fs, the S(2) population amounts to ～60% of the initial S(0) population, and remains at ～20% after 100 fs, in contrast to the S(0) electronic state which is completely depopulated within 75 fs.     

Trans-cis photoisomerization of azobenzene by n→π* excitation:A semiclassical dynamics study

A realistic dynamics simulation study is reported for the trans-cis photoisomerization of azobenzene triggered by the n→π*excitation and the results show that the formation of cis isomer follows the rotational motion around the N=N bond.The simulation find that the CNN bond angle bending vibrations also play a significant role in the vibronic coupling between the HOMO and LUMO,which essentially leads a nonadiabatic transition of the molecule to the electronic ground state.     

Theoretical studies of the photochemical dynamics of acetylacetone: isomerzation, dissociation, and dehydration reactions

The potential energy surfaces of the C-O cleavage, rotational isomerization, keto-enolic tautomerization, and dehydration reactions of acetylacetone in the lowest triplet and ground states have been determined using the complete active space self-consistent field and density functional theory methods. The main photochemical mechanism obtained indicates that the acetylacetone molecule in the S(2)((1)pipi*) state can relax to the T(1)((3)pipi*) state via the S(2)-S(1) vibronic interaction and an S(1)/T(1)/T(2) intersection. The C-O fission pathway is the predominant dissociation process in the T(1)((3)pipi) state. Rotational isomerization reactions proceed difficultly in the ground state but very easily in the T(1)((3)pipi*) state. Keto-enolic tautomerization takes place with little probability for acetylacetone in the gas phase.     

电子效应调控丁烯二腈光分子开关的CASSCF和MS-CASPT2理论研究

采用多组态CASSCF方法和MS-CASPT2方法研究了丁烯二腈中性分子及阳离子和阴离子的顺-反异构化机理.结果表明,中性分子和离子态的光顺-反异构化反应经历不同的非绝热跃迁途径:中性丁烯二腈受光激发至S1态后,需克服一个不低于19.7 k J/mol的能垒才有机会到达基态和激发态间的圆锥交叉(S_1/S_0-CI),随后经非辐射跃迁回到基态,S_1/S_0-CI在结构上偏离C=C双键旋转路径,且能量较高,因此会降低旋转速度,阻碍旋转的单向性;丁烯二腈阳离子和阴离子自由基的D_0态和D_1态旋转势能面在90°处相交,优化的D_1/D_0-CI与D_1态中间体的结构和能量均相近,因此从D1态经由D_1/D_0-C_I无辐射跃迁到D_0态的过程无势垒,在此过程中C=C旋转方向性得到最大限度的保持.研究结果证实了电子诱导不仅能降低基态热旋转势垒,而且能够调控光旋转的非绝热跃迁机理.     

A fast photoswitch for minimally perturbed peptides: investigation of the trans-->cis photoisomerization of N-methylthioacetamide

Thio amino acids can be integrated into the backbone of peptides without significantly perturbing their structure. In this contribution we use ultrafast infrared and visible spectroscopy as well as state-of-the-art ab initio computations to investigate the photoisomerization of the trans form of N-methylthioacetamide (NMTAA) as a model conformational photoswitch. Following the S2 excitation of trans-NMTAA in water, the return of the molecule into the trans ground state and the formation of the cis isomer is observed on a dual time scale, with a fast component of 8-9 ps and a slow time constant of approximately 250 ps. On both time scales the probability of isomerization to the cis form is found to be 30-40%, independently of excitation wavelength. Ab initio CASPT2//CASSCF photochemical reaction path calculations indicate that, in vacuo, the trans-->cis isomerization event takes place on the S1 and/or T1 triplet potential energy surfaces and is controlled by very small energy barriers, in agreement with the experimentally observed picosecond time scale. Furthermore, the calculations identify one S2/S1 and four nearly isoenergetic S1/S0 conical intersection decay channels. In line with the observed isomerization probability, only one of the S1/S0 conical intersections yields the cis conformation upon S1-->S0 decay. A substantially equivalent excited-state relaxation results from four T1/S0 intersystem crossing points.     

Time-resolved IR spectroscopy of N-methylthioacetamide: trans-->cis isomerization upon n-pi and pi-pi excitation and cis-->trans photoreaction

Time-resolved infrared spectroscopy was used to study the photoisomerization of N-Methylthioacetamide (NMTAA) in D2O in both the cis-->trans and the trans-->cis direction upon selective excitation of the n-pi (S1) and pi-pi (S2) electronic transitions. While isomerization and the return to the ground state takes place on two distinct time scales (cis isomerization is 30-40%, independent of the electronic state excited, while the cis-->trans isomerization proceeds with a 60-70% quantum efficiency. These results support a mechanism by which isomerization takes place via one common intermediate state independent of electronic excitation energy and initial conformation.     

Does one‐photon photocyclization of trans‐diarylethylenes involve adiabatic trans‐to‐cis photoisomerization? potential energy surface calculations for 1‐styrylnaphthalene

Potential energy surface (PES) for 1‐styrylnaphthalene was calculated by PM3 method for the S₀ state and PM3‐CI(2x2) method with configuration interaction for the S₁ state. Scanning PES along both isomerization and cyclization reaction coordinates enabled to reveal the minimum energy path (MEP) with low barriers on the S₁ PES from E‐isomer to dihydrocyclophotoproduct (DHP). This is consistent with formation of the photocyclization product in one‐photon process during irradiation of E‐isomer. Additionally, the MEP was found to bypass the coordinate region of Z‐isomer, i.e. one‐photon E‐isomer‐to‐DHP photocyclization does not demand participation of the excited Z‐isomer. Therefore, adiabatic trans‐to‐cis isomerization is likely not an intermediate stage on the E‐isomer photocyclization pathway, and experimentally observed one‐photon formation of the DHP from the E‐isomer is likely not an evidence for adiabatic trans‐to‐cis photoisomerization, as it is usually assumed. According to the results obtained, two photochemical reactions of E‐isomer, photoisomerization to Z‐isomer and photocyclization to DHP, are not consecutive but parallel reactions with branching at perpendicular conformer on the S₁ PES. © 2012 Wiley Periodicals, Inc.     

Mechanism for the Nonadiabatic Photooxidation of Benzene to Phenol: Orientation‐Dependent Proton‐Coupled Electron Transfer

An efficient catalytic one‐step conversion of benzene to phenol was achieved recently by selective photooxidation under mild conditions with 2,3‐dichloro‐5,6‐dicyano‐p‐benzoquinone (DDQ) as the photocatalyst. Herein, high‐level electronic structure calculations in the gas phase and in acetonitrile solution are reported to explore the underlying mechanism. The initially populated ¹ππ* state of DDQ can relax efficiently through a nearby dark ¹nπ* doorway state to the ³ππ* state of DDQ, which is found to be the precursor state involved in the initial intermolecular electron transfer from benzene to DDQ. The subsequent triplet‐state reaction between DDQ radical anions, benzene radical cations, and water is computed to be facile. The formed DDQH and benzene‐OH radicals can undergo T₁→S₀ intersystem crossing and concomitant proton‐coupled electron transfer (PCET) to generate the products DDQH₂ and phenol. Two of the four considered nonadiabatic pathways involve an orientation‐dependent triplet PCET process, followed by intersystem crossing to the ground state (S₀). The other two first undergo a nonadiabatic T₁→S₀ transition to produce a zwitterionic S₀ complex, followed by a barrierless proton transfer. The present theoretical study identifies novel types of nonadiabatic PCET processes and provides detailed mechanistic insight into DDQ‐catalyzed photooxidation.     

Unraveling ultrafast dynamics in photoexcited aniline

A combination of ultrafast time-resolved velocity map imaging (TR-VMI) methods and complete active space self-consistent field (CASSCF) ab initio calculations are implemented to investigate the electronic excited-state dynamics in aniline (aminobenzene), with a perspective for modeling (1)πσ* mediated dynamics along the amino moiety in the purine derived DNA bases. This synergy between experiment and theory has enabled a comprehensive picture of the photochemical pathways/conical intersections (CIs), which govern the dynamics in aniline, to be established over a wide range of excitation wavelengths. TR-VMI studies following excitation to the lowest-lying (1)ππ* state (1(1)ππ*) with a broadband femtosecond laser pulse, centered at wavelengths longer than 250 nm (4.97 eV), do not generate any measurable signature for (1)πσ* driven N-H bond fission on the amino group. Between wavelengths of 250 and >240 nm (<5.17 eV), coupling from 1(1)ππ* onto the (1)πσ* state at a 1(1)ππ*/(1)πσ* CI facilitates ultrafast nonadiabatic N-H bond fission through a (1)πσ*/S(0) CI in <1 ps, a notion supported by CASSCF results. For excitation to the higher lying 2(1)ππ* state, calculations reveal a near barrierless pathway for CI coupling between the 2(1)ππ* and 1(1)ππ* states, enabling the excited-state population to evolve through a rapid sequential 2(1)ππ* → 1(1)ππ* → (1)πσ* → N-H fission mechanism, which we observe to take place in 155 ± 30 fs at 240 nm. We also postulate that an analogous cascade of CI couplings facilitates N-H bond scission along the (1)πσ* state in 170 ± 20 fs, following 200 nm (6.21 eV) excitation to the 3(1)ππ* surface. Particularly illuminating is the fact that a number of the CASSCF calculated CI geometries in aniline bear an exceptional resemblance with previously calculated CIs and potential energy profiles along the amino moiety in guanine, strongly suggesting that the results here may act as an excellent grounding for better understanding (1)πσ* driven dynamics in this ubiquitous genetic building block.     

Ab initio study on one-way photoisomerization of the maleic acid and fumaric acid anion radical system as a model system of their esters

Potential energy surfaces (PESs) of the maleic acid anion radical (MA-*: cis isomer)/fumaric acid anion radical (FA-*: trans isomer) system as a model system of their esters have been studied in detail using CASSCF method. The results suggest the following: The photoisomerization is initiated with the H-C-C-H dihedral angle distortion [hydrogen out of plain (HOOP) motion] on the D1 PES. The C-C-C-C dihedral angle distortion occurs on the D0 PES after the deactivation from D1 to D0. A large fraction of the net motion along the isomerization coordinate occurs on the D0 PES. The D0 state is responsible for the one-way nature of the photoisomerization.