Mechanistic Insights into the Excited-State Intramolecular Proton Transfer (ESIPT) Process of 2-(2-Aminophenyl)naphthalene.

The 2-(2-aminophenyl)naphthalene molecule attracted much attention due to excited-state intramolecular proton transfer (ESIPT) from an amino NH2 group to a carbon atom of an adjacent aromatic ring. The ESIPT mechanisms of 2-(2-aminophenyl)naphthalene are still unclear. Herein, the decay pathways of this molecule in vacuum were investigated by combining static electronic structure calculations and nonadiabatic dynamics simulations. The calculations indicated the existence of two stable structures (S0-1 and S0-2) in the S0 and S1 states. For the S0-1 isomer, upon excitation to the Franck-Condon point, the system relaxed to the S1 minimum quickly, and then there exist four decay pathways (two ESIPT ones and two decay channels with C atom pyramidalization). In the ESIPT decay pathways, the system encounters the S1S0-PT-1 or S1S0-PT-2 conical intersection, which funnels the system rapidly to the S0 state. In the other two pathways, the system de-excited from the S1 to the S0 state via the S1S0-1 or S1S0-2 conical intersection. For the S0-2 structure, the decay pathways were similar to those of S0-1. The dynamics simulations showed that 75 and 69% of trajectories experienced the two ESIPT conical intersections for the S0-1 and S0-2 structures, respectively. Our simulations showed that the lifetime of the S1 state of S0-1 (S0-2) is estimated to be 358 (400) fs. Notably, we not only found the detailed reaction mechanism of the system but also found that the different ground-state configurations of this system have little effect on the reaction mechanism in vacuum.

Photochemical Mechanisms of Hydroxyquinoline Benzimidazole: Insights from Electronic Structure Calculations and Nonadiabatic Dynamics Simulations.

Excited-state intramolecular double proton transfer (ESIDPT) has received much attention because of its widespread existence in the life reactions of living organisms, and materials with this property are significant for their special luminescent properties. In this work, the complete active space self-consistent field (CASSCF) and OM2/multireference configuration interaction (OM2/MRCI) methods have been employed to study the static electronic structure calculations of the photochemistry and the possibility of ESIDPT process of hydroxyquinoline benzimidazole (HQB) molecule, along with the nonadiabatic dynamics simulations. The computational results show that the HQB molecule is relaxed to the S1-ENOL minimum after being excited to the Franck-Condon point in the S1 state. Subsequently, during the nonadiabatic deactivation process, the OH···N proton transfer and the twisting of benzimidazole occur before arriving at the single proton transfer conical intersection S1S0-KETO. Finally, the system can either return to the initial ground-state structure S0-ENOL or to the single proton transfer ground-state structure S0-KETO, both of which have almost the same probability. The dynamics simulations also show that no double proton transfer occurs. The excited-state lifetime of HQB is fitted to 1.1 ps, and only 64% of the dynamic trajectories return to the ground state within the 2.0 ps simulation time. We hope the detailed reaction mechanism of the HQB molecule will provide new insights into similar systems.

Structures and Energetics of E2H3+ (E = As, Sb, and Bi) Cations.

E2H2 (E = As, Sb, Bi) structures involving multiple bonds have attracted much attention recently. The E2H3+ cations (protonated E2H2) are predicted to be viable with substantial proton affinities (>180 kcal/mol). Herein, the bonding characters and energetics of a number of E2H3+ isomers are explored through CCSD(T) and DFT methods. For the As2H3+ system, the CCSD(T)/cc-pVQZ-PP method predicts that the vinylidene-like structure lies lowest in energy, with the trans and cis isomers higher by 6.7 and 9.3 kcal/mol, respectively. However, for Sb2H3+ and Bi2H3+ systems, the trans isomer is the global minimum, while the energies of the cis and vinylidene-like structures are higher, respectively, by 2.0 and 2.4 kcal/mol for Sb2H3+ and 1.6 and 15.0 kcal/mol for Bi2H3+. Thus, the vinyledene-like structure is the lowest energy for the arsenic system but only a transition state of the bismuth system. With permanent dipole moments, all minima may be observable in microwave experiments. Besides, we have also obtained transition states and planar-cis structures with higher energies. The current results should provide new insights into the various isomers and provide a number of predictions for future experiments.

Excited-state relaxation mechanisms of 2,2'-(1-phenyl-1H-1,2,4-triazole-3,5-diyl)diphenol: single- or double-proton transfer?

Triazole compounds are important organic systems with excellent electronic properties, which have diagnostic potential in the fields of organic electronics and organic photovoltaics. The important photophysical nature of these systems is the transformation between the enol and keto forms after excited-state proton transfer. In this study, the IR vibrational spectrum, ESIPT mechanism, and excited-state decay dynamics of 2,2'-(1-phenyl-1H-1,2,4-triazole-3,5-diyl)diphenol (ExPh) were explored using electronic structure calculations and non-adiabatic dynamics simulations. Two S1/S0 conical intersections with distinct proton transfer (ESIPT-I and ESIPT-II) involved were obtained. The associated two-dimensional S1 minimum-energy potential energy surface indicated that the dynamical roles of these two S1/S0 conical intersections in the S1 excited-state decay were quite different. The ESIPT-I reaction was more favorable to occur than the ESIPT-II process. Our dynamics simulations supported this hypothesis with the whole trajectories decaying to the ground state via the S1S0-1 conical intersection, which involved the ESIPT-I process. The ESIPT-Involved efficient deactivation pathway could be partially responsible for the decrease in fluorescence emission. These results and ESIPT mechanisms are helpful for understanding the decay pathways of similar systems.

Non-adiabatic dynamics simulations of the S1 excited-state relaxation of diacetyl phenylenediamine.

The small molecule built around the benzene ring, diacetyl phenylenediamine (DAPA), has attracted much attention due to its synthesis accessibility, large Stokes shift, etc. However, its meta structure m-DAPA does not fluoresce. In a previous investigation, it was found that such a property is due to the fact that it undergoes an energy-reasonable double proton transfer conical intersection during the deactivation of the S1 excited-state, then returns to the ground state by a nonradiative relaxation process eventually. However, our static electronic structure calculations and non-adiabatic dynamics analysis results indicate that only one reasonable non-adiabatic deactivation channel exists: after being excited to the S1 state, m-DAPA undergoes an ultrafast and barrierless ESIPT process and reaches the single-proton-transfer conical intersection. Subsequently, the system either returns to the keto-form S0 state minimum with proton reversion or returns to the single-proton-transfer S0 minimum after undergoing a slight twist of the acetyl group. The dynamics results show that the S1 excited-state lifetime of m-DAPA is 139 fs. In other words, we propose an efficient single-proton-transfer non-adiabatic deactivation channel of m-DAPA that is different from previous work, which can provide important mechanistic information of similar fluorescent materials.

Combined QM (MS-CASPT2)/MM studies on photocyclization and photoisomerization of a fulgide derivative in toluene solution.

Photocyclization and photoisomerization of fulgides have been extensively studied experimentally and computationally due to their significant potential applications for example as photoswitches in memory devices. However, the reported excited-state decay mechanisms of fulgides do not include the effects of solvation explicitly to date. Herein, calculations using the high-level MS-CASPT2//CASSCF method were conducted to explore the photoinduced excited-state decay processes of the Eα conformer of a fulgide derivative in toluene with solvent effects treated by implicit PCM and explicit QM/MM models, respectively. Several minima and conical intersections were optimized successfully in and between the S0 and S1 states; then, two nonadiabatic excited-state decay channels that could efficiently drive the system to the ground state were proposed based on the excited-state ring-closure and isomerization paths. In addition, we also found that in the ring-closure path, the potential energy surface is essentially barrierless before approaching the conical intersection, while it needs to overcome a small energy barrier along the E → Z photoisomerization path for the nonadiabatic S1 → S0 internal conversion process. The present computational results could provide useful mechanistic insights into the photoinduced cyclization and isomerization reactions of fulgide and its derivatives.

Excited-state photochemistry dynamics of 2-(1-naphthyl) phenol: electronic structure calculations and non-adiabatic dynamics simulations.

The excited-state proton transfer processes and the formation mechanism of quinone methide of (1-naphthyl)phenol were investigated by combining static electronic structure calculations and non-adiabatic dynamics simulations in vacuum. The results indicated the existence of two minimum energy structures (S0-ENOL-1 and S0-ENOL-2) in the ground and excited states, which correspond to two ESIPT pathways. Upon excitation of S0-ENOL-1 to the bright S1 state, the system relaxes to the S1 minimum quickly in the enol region, for which two decay pathways have been described. The first is a barrierless ESIPT-1 process that generates keto species. Afterwards, the system encounters a keto conical intersection, which funnels the system to the ground state. The generated keto species, in the S0 state, either regenerated the starting material via ground-state proton transfer or yielded the keto product at the end of the simulations. In the other pathway, the system de-excites from the S1 state to the S0 state via one enol-type conical intersection. The dynamics simulations showed that 58.8% of trajectories experience keto-type conical intersection and the rest undergo enol-type conical intersection. Besides the ESIPT-1 process, a new-type ESIPT (ESIPT-2), which was not observed experimentally, was found with the irradiation of S0-ENOL-2. The ESIPT-2 process occurs after overcoming a small barrier (0.9 kcal mol-1) and yields a distinct quinone methide. Our simulation results also showed that the S1 lifetime of S0-ENOL-1 (S0-ENOL-2) would be 437 (617) fs in the gas phase. These results provide detailed and important mechanistic insights into the systems in which ESPT to carbon atoms occurs.

Excited-State Deactivation Mechanism of 3,5-bis(2-hydroxyphenyl)-1H-1,2,4-triazole: Electronic Structure Calculations and Nonadiabatic Dynamics Simulations.

3,5-bis(2-Hydroxyphenyl)-1H-1,2,4-triazole (bis-HPTA) has attracted wide attention due to the important application in the detection of microorganisms and insecticidal activity. However, the mechanisms of excited-state intramolecular proton transfer (ESIPT) process and decay pathways are still a matter of debate. In this work, we have comprehensively investigated the photodynamics of bis-HPTA by executing combined electronic structure calculations and nonadiabatic surface-hopping dynamics simulations. Based on the computed electronic structure and dynamics information, we propose two nonadiabatic deactivation channels that efficiently populate the ground state from the Franck-Condon region. In the first one, after being excited to the bright S1 state, bis-HPTA molecule undergoes an ultrafast and barrierless ESIPT-1 process. Then, the system encounters with an energetically accessible S1/S0 conical intersection (CI), which funnels the system to the ground state speedily. Afterward, the keto species either arrives at the keto product or return to its enol species via a ground-state proton transfer in the S0 state. In the other excited-state decay channel, the S1 system hops to the ground state through a different CI, which involves the ESIPT-2 process. In our dynamics simulations, about 79.6% of the trajectories decay to the S0 state via the first CI, while the remaining ones employ the second conical intersection. The results of dynamics simulations also demonstrated that the lifetime of the S1 state is estimated as 315 fs. The present work will give elaborating mechanistic information of similar compounds in various environments.

Comprehensive analysis of the differential expression profile of microRNAs in missed abortion.

To screen the key circulating microRNAs (miRNAs) involved in missed abortion (MA) and explore their role in MA process. We examined the miRNA profile from the serum of three MA patients and three early pregnancy induced abortion patients (controls) by next-generation sequencing. We analyzed the target genes of the differentially expressed (DE) miRNAs to analyze the function and pathway enrichment using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, respectively. We validated five candidate miRNAs by real time-qPCR. Integrated miRNA-mRNA-pathway network analysis was performed to show the interaction network of the candidate miRNAs and their target genes of interest with the involved pathways. It was observed that 227 miRNAs were differently expressed between the MA group and the early pregnancy control group, with 58 miRNAs downregulated and 169 miRNAs upregulated in the MA group. Real-time qPCR results revealed that expression of the five candidate miRNAs, namely hsa-miR-22-3p, hsa-miR-145-3p, hsa-miR-107, hsa-miR-361-3p, and hsa-miR-378c, was consistent with the miRNA data obtained by sequencing. Integrated miRNA-mRNA-pathway network analysis illustrated that target genes of the candidate miRNAs were mainly involved in the PI3K-Akt signaling pathway, HIF-1 signaling pathway, and VEGF signaling pathway, which would have potential significance in pregnancy and MA. We are the first to reveal the DE miRNAs involved in MA and illustrate their functional interaction network. These results might provide potential circulating biomarkers and new therapeutic targets for MA.

Photochemical mechanism of 1,5-benzodiazepin-2-one: electronic structure calculations and nonadiabatic surface-hopping dynamics simulations.

Due to the significant applications in bioimaging, sensing, optoelectronics etc., photoluminescent materials have attracted more and more attention in recent years. 1,5-Benzodiazepin-2-one and its derivatives have been used as fluorogenic probes for the detection of biothiols. However, their photochemical and photophysical properties have remained ambiguous until now. In this work, we have adopted combined static electronic structure calculations and nonadiabatic surface-hopping dynamics simulations to study the photochemical mechanism of 1,5-benzodiazepin-2-one. Firstly, we optimized minima and conical intersections in S0 and S1 states; then, we proposed three nonadiabatic decay pathways that efficiently populate the ground state from the Franck-Condon region based on computed electronic structure information and dynamics simulations. In the first pathway, upon photoexcitation to the S1 state, the system proceeds with an ultrafast excited-state intramolecular proton transfer (ESIPT) process. Then, the molecule tends to rotate around the C-C bond until it encounters keto conical intersections, from which the system can easily decay to the ground state. The other two pathways involve the enol channels in which the S1 system hops to the ground state via two enol S1/S0 conical intersections, respectively. These three energetically allowed S1 excited-state deactivation pathways are responsible for the decrease of fluorescence quantum yield. The present work will provide detailed mechanistic information of similar systems.

Electronic structure calculations and nonadiabatic dynamics simulations of excited-state relaxation of Pigment Yellow 101.

Pigment Yellow 101 (PY101) is widely used as a typical pigment due to its excellent excited-state properties. However, the origin of its photostability is still elusive. In this work, we have systematically investigated the photodynamics of PY101 by performing combined electronic structure calculations and trajectory-based nonadiabatic dynamics simulations. On the basis of the results, we have found that upon photoexcitation to the S1 state, PY101 undergoes an essentially barrierless excited-state intramolecular single proton transfer generating an S1 keto species. In the keto region, there is an energetically accessible S1/S0 conical intersection that funnels the system to the S0 state quickly. In the S0 state, the keto species either goes back to its trans-enol species through a ground-state reverse hydrogen transfer or arrives at the cis-keto region. In addition, we have found an additional excited-state decay channel for the S1 enol species, which is directly linked to an S1/S0 conical intersection located in the enol region. This mechanism has also been confirmed by our dynamics simulations, in which about 54% of the trajectories decay to the S0 state via the enol S1/S0 conical intersection; while the remaining ones employ the keto S1/S0 conical intersection. The gained mechanistic information helps us understand the photostability of the PY101 chromophore and its variants with the same molecular scaffold.

A theoretical study of ruthenium complexes with 2,2'-biimidazole-like ligands: structural, optical and emissive properties.

The structural and optical properties of five ruthenium complexes, recently synthesized for their photooxidative and photophysical properties, have been studied by means of density functional theory (DFT) and time-dependent DFT (TD-DFT). The structures of [Ru(bpy)2(BiimH2)](2+) (bpy = 2,2'-bipyridine; BiimH2 = 2,2'-biimidazole) 1, [Ru(bpy)2(TMBiimH2)](2+) (TM BiimH2 = 4,5,4',5'-tetramethyl-2,2'-biimidazole) 5, [Ru(bpy)2(L1H2)](2+) (L1H2 = 4,5-dimethyl-2(N,N-diacetyl)(carboximidamide-1H-imidazole)) 6, [Ru(bpy)2(L2H2)](2+) (L2H2 = N(1),N(1),N(2),N(2)-tetrakis(acetyl)ethanediimidamide) 7 and [Ru(phen)2(TMBiimH2)](2+) (phen = 1,10'-phenanthroline) 8 have been fully optimized in the electronic ground state as well as in the lowest triplet T1 excited state. The theoretical absorption spectra of the five complexes that compare rather well with the experimental spectra have been analyzed on the basis of TD-DFT calculations without and with spin-orbit coupling (SOC). The deprotonated form [Ru(bpy)2(L2H)](+)7d contributes mostly to the experimental absorption spectrum of complex 7. The spectra of all molecules are characterized by the presence of low-lying metal-to-ligand charge transfer (MLCT) excited states between 500 and 400 nm, ligand-centered (LC) excited states on the biimidazole-like ligands between 350 and 300 nm and on the bpy ligands between 300 and 250 nm. The theoretical emission wavelengths deduced from the lowest triplet T1 properties calculated at 661 nm (1), 690 nm (5) and 660 nm (8) reproduce the experimental emission spectra of these molecules characterized by a maximum at 638 nm (1), 646 nm (5) and 652 nm (8). In contrast the low theoretical emission wavelengths (>1000 nm) obtained for complexes 6, 7 and 7d favorable to non-radiative decays explain the low intensity of the experimental emission spectra of these two complexes. The SOC is of little effect in this class of molecules where metal-centered (MC) excited states do not perturb the lowest part of the absorption spectra leading to negligible splitting of low-lying triplet states.

Theoretical study on photooxidation mechanism of ruthenium complex [Ru(II)-(bpy)2 (TMBiimH2 )](2+) with molecular oxygen.

Photoinduced reactions of ruthenium complexes with molecular oxygen have attracted a lot of experimental attention; however, the reaction mechanism remains elusive. In this work, we have used the density functional theory method to scrutinize the visible-light induced photooxidation mechanism of the ruthenium complex [Ru(II)-(bpy)2 (TMBiimH2 )](2+) (bpy: 2, 2-bipyridine and TMBiimH2 : 4, 5, 4, 5-tetramethyl-2, 2-biimidazole) initiated by the attack of molecular oxygen. The present computational results not only explain very well recent experiments, also provide new mechanistic insights. We found that: (1) the triplet energy transfer process between the triplet molecular oxygen and the metal-ligand charge transfer triplet state of the ruthenium complex, which leads to singlet molecular oxygen, is thermodynamically favorable; (2) the singlet oxygen addition process to the S0 ruthenium complex is facile in energy; (3) the chemical transformation from endoperoxide to epidioxetane intermediates can be either two- or one-step reaction (the latter is energetically favored). These findings contribute important mechanistic information to photooxidation reactions of ruthenium complexes with molecular oxygen. © 2016 Wiley Periodicals, Inc.

How Photoisomerization Drives Peptide Folding and Unfolding: Insights from QM/MM and MM Dynamics Simulations.

Photoswitchable azobenzene cross-linkers can control the folding and unfolding of peptides by photoisomerization and can thus regulate peptide affinities and enzyme activities. Using quantum mechanics/molecular mechanics (QM/MM) methods and classical MM force fields, we report the first molecular dynamics simulations of the photoinduced folding and unfolding processes in the azobenzene cross-linked FK-11 peptide. We find that the interactions between the peptide and the azobenzene cross-linker are crucial for controlling the evolution of the secondary structure of the peptide and responsible for accelerating the folding and unfolding events. They also modify the photoisomerization mechanism of the azobenzene cross-linker compared with the situation in vacuo or in solution.

Excited-State Proton-Transfer-Induced Trapping Enhances the Fluorescence Emission of a Locked GFP Chromophore.

The chemical locking of the central single bond in core chromophores of green fluorescent proteins (GFPs) influences their excited-state behavior in a distinct manner. Experimentally, it significantly enhances the fluorescence quantum yield of GFP chromophores with an ortho-hydroxyl group, while it has almost no effect on the photophysics of GFP chromophores with a para-hydroxyl group. To unravel the underlying physical reasons for this different behavior, we report static electronic structure calculations and nonadiabatic dynamics simulations on excited-state intramolecular proton transfer, cis-trans isomerization, and excited-state deactivation in a locked ortho-substituted GFP model chromophore (o-LHBI). On the basis of our previous and present results, we find that the S1 keto species is responsible for the fluorescence emission of the unlocked o-HBI and the locked o-LHBI species. Chemical locking does not change the parts of the S1 and S0 potential energy surfaces relevant to enol-keto tautomerization; hence, in both chromophores, there is an ultrafast excited-state intramolecular proton transfer that takes only 35 fs on average. However, the locking effectively hinders the S1 keto species from approaching the keto S1/S0 conical intersections so that most of trajectories are trapped in the S1 keto region for the entire 2 ps simulation time. Therefore, the fluorescence quantum yield of o-LHBI is enhanced compared with that of unlocked o-HBI, in which the S1 excited-state decay is efficient and ultrafast. In the case of the para-substituted GFP model chromophores p-HBI and p-LHBI, chemical locking hardly affects their efficient excited-state deactivation via cis-trans isomerization; thus, the fluorescence quantum yields in these chromophores remain very low. The insights gained from the present work may help to guide the design of new GFP chromophores with improved fluorescence emission and brightness.

Photophysics of Auramine-O: electronic structure calculations and nonadiabatic dynamics simulations.

Diphenylmethane dyes are very useful photoinduced molecular rotors; however, their photophysical mechanisms are still elusive until now. In this work, we adopted combined static electronic structure calculations (MS-CASPT2//CASSCF) and trajectory-based surface-hopping dynamics simulations (OM2/MRCI) to study the S1 excited-state relaxation mechanism of a representative diphenylmethane dye Auramine-O. On the basis of the optimized S1 minima and the computed emission bands, we have for the first time assigned experimentally proposed three transient states (i.e. S1-LE, S1-I1 or S1-I2, and S1-II). Mechanistically, upon irradiation to the S1 state, the system first relaxes to the locally excited S1 minimum (S1-LE). Starting from this point, there exist two kinds of relaxation paths to S1-II. In the sequential path, the system first evolves into S1-I1 or S1-I2 and then runs into S1-II; in the concerted one, the system, bypassing S1-I1 and S1-I2, directly runs into S1-II. In addition, the system can decay to the S0 state in the vicinity of three S1/S0 conical intersections i.e. S1S0-I1, S1S0-I2, and S1S0-II. In the S1 dynamic simulations, 54% trajectories decay to the S0 state via S1S0-II; the remaining trajectories are de-excited to the S0 state via S1S0-I1 (11%) and S1S0-I2 (35%). Our present theoretical investigation does not support the experimentally proposed S1 excited-state hypothesis that the intramolecular rotation of the two dimethyl groups around the C-N bond is responsible for the rapid decay of the emission band at about 500 nm; instead, it should be heavily interrelated with the rotation of the two dimethylanilino groups. Finally, this work provides important mechanistic insights into similar diphenylmethane dyes.

Photodissociation dynamics of CH3C(O)SH in argon matrix: A QM/MM nonadiabatic dynamics simulation.

In this work, we have first employed the combined quantum mechanics/molecular mechanics (QM/MM) method to study the photodissociation mechanism of thioacetic acid CH3C(O)SH in the S1, T1, and S0 states in argon matrix. CH3C(O)SH is treated quantum mechanically using the complete active space self-consistent field and complete active space second-order perturbation theory methods; argon matrix is described classically using Lennard-Jones potentials. We find that the C-S bond fission is predominant due to its small barriers of ca. 3.0 and 1.0 kcal/mol in the S1 and T1 states. It completely suppresses the nearby C-C bond fission. After the bond fission, the S1 radical pair of CH3CO and SH can decay to the S0 and T1 states via internal conversion and intersystem crossing, respectively. In the S0 state, the radical pair can either recombine to form CH3C(O)SH or proceed to form molecular products of CH2CO and H2S. We have further employed our recently developed QM/MM generalized trajectory-based surface-hopping method to simulate the photodissociation dynamics of CH3C(O)SH. In 1 ps dynamics simulation, 56% trajectories stay at the Franck-Condon region; the S1 C-S bond fission takes place in the remaining 44% trajectories. Among all nonadiabatic transitions, the S1 → S0 internal conversion is major (55%) but the S1 → T1 intersystem crossing is still comparable and cannot be ignored, which accounts for 28%. Finally, we have found a radical channel generating the molecular products of CH2CO and H2S, which is complementary to the concerted molecular channel. The present work sets the stage for simulating photodissociation dynamics of similar thio-carbonyl systems in matrix.

Surface-hopping dynamics simulations of malachite green: a triphenylmethane dye.

Malachite green is a typical triphenylmethane dye widely used in fundamental and industrial research; however, its excited-state relaxation dynamics remains elusive. In this work we simulate its photodynamics from the S2 and S1 states using the fewest-switches surface-hopping scheme. In the S2 photodynamics, the system first relaxes to the S2 minimum, which immediately hops to the S1 state via an S2/S1 conical intersection. In the S1 state, 90% trajectories evolve into a structurally symmetric S1 minimum; the remaining ones proceed toward two propeller-like S1 minima. Two kinds of S1 minima then decay to the S0 state via the S1/S0 conical intersections. The S1 photodynamics is overall similar to the S1 excited-state dynamics as a result of the ultrafast S2 → S1 internal conversion in the S2 photodynamics, but the weights of the trajectories that decay to the S0 state via three different S1/S0 conical intersections are variational. Moreover, the S2 relaxation dynamics mainly happens in a concerted synchronous rotation of three phenyl rings. In comparison, in the S1 relaxation dynamics, the rotations of two aminophenyl rings can proceed in the same and opposite directions. In certain trajectories, only the rotation of an aminophenyl ring is active. On the basis of the results, the S2 and S1 excited-state lifetimes of malachite green in vacuo are calculated to be 424 fs and 1.2 ps, respectively. The present work provides important mechanistic insights for similar triphenylmethane dyes.

Excited-state ring-opening mechanism of cyclic ketones: a MS-CASPT2//CASSCF study.

We have employed complete active space self-consistent field (CASSCF) and its second-order perturbation (MS-CASPT2) methods to study the S1 and T1 excited-state ring-opening mechanisms and S1 excited-state deactivation channels of cyclopropanone, cyclobutanone, cyclopentanone, and cyclohexanone. On the basis of optimized minima, transition states, conical intersections, refined energies, and relaxed two-dimensional S1 and T1 potential energy surfaces, we find that, with the ring-strain decrease from cyclopropanone to cyclohexanone, (1) the ring-opening S1 and T1 barrier increases from 0.0 and 0.0 to 19.7 and 10.4 kcal/mol, respectively; (2) the electronic state responsible for the dominant ring-opening reaction varies from the S1 state of cyclopropanone to the T1 state of cyclopentanone and cyclohexanone; and (3) the S1 ring opening gradually becomes inefficient even blocked in cyclopentanone and cyclohexanone. This work shows that these dissimilar excited-state dynamics could originate from different ring strain of small cyclic ketones.

Excited-state intramolecular proton transfer to carbon atoms: nonadiabatic surface-hopping dynamics simulations.

Excited-state intramolecular proton transfer (ESIPT) between two highly electronegative atoms, for example, oxygen and nitrogen, has been intensely studied experimentally and computationally, whereas there has been much less theoretical work on ESIPT to other atoms such as carbon. We have employed CASSCF, MS-CASPT2, RI-ADC(2), OM2/MRCI, DFT, and TDDFT methods to study the mechanistic photochemistry of 2-phenylphenol, for which such an ESIPT has been observed experimentally. According to static electronic structure calculations, irradiation of 2-phenylphenol populates the bright S1 state, which has a rather flat potential in the Franck-Condon region (with a shallow enol minimum at the CASSCF level) and may undergo an essentially barrierless ESIPT to the more stable S1 keto species. There are two S1/S0 conical intersections that mediate relaxation to the ground state, one in the enol region and one in the keto region, with the latter one substantially lower in energy. After S1 → S0 internal conversion, the transient keto species can return back to the S0 enol structure via reverse ground-state hydrogen transfer in a facile tautomerization. This mechanistic scenario is verified by OM2/MRCI-based fewest-switches surface-hopping simulations that provide detailed dynamic information. In these trajectories, ESIPT is complete within 118 fs; the corresponding S1 excited-state lifetime is computed to be 373 fs in vacuum. Most of the trajectories decay to the ground state via the S1/S0 conical intersection in the keto region (67%), and the remaining ones via the enol region (33%). The combination of static electronic structure computations and nonadiabatic dynamics simulations is expected to be generally useful for understanding the mechanistic photophysics and photochemistry of molecules with intramolecular hydrogen bonds.