Ganglong Cui

Find an error

Name: 崔刚龙; Cui, GangLong
Organization: Beijing Normal University , China
Department: College of Chemistry
Title: (PhD)
Co-reporter:Qian Wang, Xiao-Ying Xie, Juan Han, and Ganglong Cui
The Journal of Physical Chemistry B November 22, 2017 Volume 121(Issue 46) pp:10467-10467
Publication Date(Web):October 30, 2017
DOI:10.1021/acs.jpcb.7b09046
Semisynthetic alphabet can potentially increase the genetic information stored in DNA through the formation of unusual base pairs such as d5SICS:dNaM. However, recent experiments show that near-visible-light irradiation on the d5SICS and dNaM chromophores could lead to genetic mutations and damages. Until now, their photophysical mechanisms remain elusive. Herein, we have employed MS-CASPT2//CASSCF and QM(MS-CASPT2//CASSCF)/MM methods to explore the spectroscopic properties and excited-state relaxation mechanisms of d5SICS, dNaM, and d5SICS:dNaM in DNA. We have found that (1) the S2 state of d5SICS, the S1 state of dNaM, and the S2 state of d5SICS:dNaM are initially populated upon near-visible-light irradiation and (2) for d5SICS and d5SICS:dNaM, there are several parallel relaxation pathways to populate the lowest triplet state, but for dNaM, a main relaxation pathway is uncovered. Moreover, we have found that the excited-state relaxation mechanism of d5SICS:dNaM in DNA is similar to that of the isolated d5SICS chromophore. These mechanistic insights contribute to the understanding of photophysics and photochemistry of unusual base pairs and to the design of better semisynthetic genetic alphabet.
Co-reporter:Pin Xiao, Qian Wang, Wei-Hai Fang, and Ganglong Cui
The Journal of Physical Chemistry A June 8, 2017 Volume 121(Issue 22) pp:4253-4253
Publication Date(Web):May 17, 2017
DOI:10.1021/acs.jpca.7b03123
Photoinduced chemical reactions of organic compounds at the marine boundary layer have recently attracted significant experimental attention because this kind of photoreactions has been proposed to have substantial impact on local new particle formation and their photoproducts could be a source of secondary organic aerosols. In this work, we have employed first-principles density functional theory method combined with cluster models to systematically explore photochemical reaction pathways of nonanoic acids (NAs) to form volatile saturated and unsaturated C9 and C8 aldehydes at air–water interfaces. On the basis of the results, we have found that the formation of C9 aldehydes is not initiated by intermolecular Norrish type II reaction between two NAs but by intramolecular T1 C—O bond fission of NA generating acyl and hydroxyl radicals. Subsequently, saturated C9 aldehydes are formed through hydrogenation reaction of acyl radical by another intact NA. Following two dehydrogenation reactions, unsaturated C9 aldehydes are generated. In parallel, the pathway to C8 aldehydes is initiated by T1 C—C bond fission of NA, which generates octyl and carboxyl radicals; then, an octanol is formed through recombination reaction of octyl with hydroxyl radical. In the following, two dehydrogenation reactions result into an enol intermediate from which saturated C8 aldehydes are produced via NA-assisted intermolecular hydrogen transfer. Finally, two dehydrogenation reactions generate unsaturated C8 aldehydes. In these reactions, water and NA molecules are found to play important roles. They significantly reduce relevant reaction barriers. Our work has also explored oxygenation reactions of NA with molecular oxygen and radical–radical dimerization reactions.
Co-reporter:Xiao-Ying Xie, Qian Wang, Wei-Hai Fang, and Ganglong Cui
The Journal of Physical Chemistry C August 3, 2017 Volume 121(Issue 30) pp:16373-16373
Publication Date(Web):July 11, 2017
DOI:10.1021/acs.jpcc.7b04811
Nitric oxide (NO) is an important air pollutant. Its chemical conversion to ammonia (NH3) and water (H2O) molecules has recently attracted a lot of experimental attention. In this work, we have employed a periodic density functional theory method combined with a slab model to study the catalytic reaction of NO adsorbed on a hydroxylated rutile TiO2(110) surface. We have obtained two favorable NO adsorption structures: in the first one, the terminal N atom is bonded with a Ti5c surface atom (NadO); in the second one, both the N and O atoms are bonded with two nearby Ti5c surface atoms (NadOad). Interestingly, NadOad becomes more stable than NadO with the increasing coverage of hydroxyl groups, i.e., more than three hydroxyl groups in our slab model, which demonstrates that hydroxyls can seriously influence surface electronic structures and, thus, surface catalysis. Mechanistically, we have found that the N–O bond should be weakened prior to its dissociation. In the NadO adsorption structure, this weakening is achieved through a hydrogen atom transfer to the N atom of the NO molecule; in the NadOad adsorption structure, this N–O bond is already activated upon adsorption on the surface. After the N–O bond is broken, a series of hydrogen atom transfers to either the N or O atom take place, which eventually produces the final products. Our present computational results provide important mechanistic insights into NO removal from TiO2 surfaces.
Co-reporter:Bin-Bin Xie;Qian Wang;Wei-Wei Guo
Physical Chemistry Chemical Physics 2017 vol. 19(Issue 11) pp:7689-7698
Publication Date(Web):2017/03/15
DOI:10.1039/C7CP00478H
The photophysics of thiothymines has been extensively studied computationally in the past few years due to their significant potential as photosensitizers in photodynamic therapy. However, the corresponding computational studies of the photophysical mechanism of 2,4-dithiothymine are scarce. Herein we have employed the CASPT2//CASSCF and QM(CASPT2//CASSCF)/MM methods to systematically explore the excited-state decay mechanism of 2,4-dithiothymine in isolated, microsolvated, and aqueous surroundings. First, we have optimized minima and conical intersections in and between the lowest six excited singlet and triplet states i.e., , , , , and ; then, based on computed excited-state decay paths and spin–orbit couplings, we have proposed several nonadiabatic pathways that efficiently populate the lowest triplet state to explain the experimentally observed ultrahigh triplet-state quantum yield. Moreover, we have found that the excited-state decay mechanism in microsolvated and aqueous environments is more complicated than that in the gas phase. The solute–solvent interaction has significant effects on the excited-state potential energy surfaces of 2,4-dithiothymine and eventually on its excited-state decay mechanism. Finally, the present computational efforts contribute important mechanistic knowledge to the understanding of the photophysics of thiothymine-based photosensitizers.
Co-reporter:Pin Xiao;Dan Wu;Wei-Hai Fang
Catalysis Science & Technology (2011-Present) 2017 vol. 7(Issue 13) pp:2763-2771
Publication Date(Web):2017/07/03
DOI:10.1039/C7CY00785J
A novel light-triggered hydrogen evolution reaction from formic acid mediated by an Ir(III) photocatalyst has been experimentally reported recently. However, its reaction mechanism remains elusive. Herein, we have employed the density functional theory (DFT) method to explore this photocatalytic reaction in detail. On the basis of the results, we have proposed a possible photocatalytic reaction mechanism. In the formation of the metal hydride [Cp*Ir(bpy)(H)]+ (5), formic acid acts as a bridge assisting proton shuttling. Upon irradiation, two nonadiabatic excited-state decay pathways quickly populate the lowest triplet T1 state of the metal hydride from its initially populated excited singlet S1 state. In the T1 state, water and formic acid facilitate excited-state hydride/proton transfers from the Ir center to Cp* and bpy ligands producing several energetically lower triplet-state isomers demonstrating that the triplet-state metal hydride 5* could not be the only precursor for the photocatalysis. Adiabatic H2 evolution in the T1 state is energetically unfavorable. These T1 isomers hop, through radiationless T1 → S0 intersystem crossings via T1/S0 crossing points, to the S0 state in which H2 evolution takes place. In these reactions, solvents acting as assistants and catalysts reduce reaction barriers, thereby accelerating H2 release and enhancing the overall photocatalytic performance. Our current work provides significant mechanistic insights into light-induced hydrogen-evolution reactions of iridium-containing photocatalysts.
Co-reporter:Xue-Ping Chang;Pin Xiao;Juan Han;Wei-Hai Fang
Physical Chemistry Chemical Physics 2017 vol. 19(Issue 21) pp:13524-13533
Publication Date(Web):2017/05/31
DOI:10.1039/C7CP01511A
In contrast to photophysics of thio-substituted nucleobases, their photoinduced cross-linking reactions with canonical nucleobases remain scarcely investigated computationally. In this work, we have adopted combined CASPT2/PCM//CASSCF and B3LYP-D3/PCM electronic structure methods to study this kind of photochemical reaction of 5-fluoro-4-thiouridine (truncated 5-fluoro-1-methyl-4-thiouracil used in calculations) and 1-methylthymine (referred to as thymine for clarity hereinafter). On the basis of CASPT2/PCM computed results, we have proposed two efficient excited-state relaxation pathways to populate the lowest T1 state of the complex of 5-fluoro-1-methyl-4-thiouracil and thymine from its initially populated S2(1ππ*) state. In the first one, the S2 system first hops to the S1 state via an S2/S1 conical intersection, followed by a direct S1 → T1 intersystem crossing process enhanced by large S1/T1 spin–orbit coupling. In the second path, the resultant S1 system first jumps to the T2 state, from which an efficient T2 → T1 internal conversion occurs. The T1 cross-linking reaction is overall divided into two phases. The first phase is a stepwise and nonadiabatic photocyclization reaction, which starts from the T1 complex and ends up with an S0 thietane intermediate. The second phase is a thermal reaction. The system first rearranges its four- and six-membered rings to form three new rings; then, an S0 fluorine atom transfer occurs, followed by the formation of photoproducts. Finally, the present work paves the way for studying light-induced cross-linking reactions of thionucleobases with canonical bases in DNA and RNA.
Co-reporter:Heyang Lin;Xueping Chang;Dongpeng Yan;Wei-Hai Fang
Chemical Science (2010-Present) 2017 vol. 8(Issue 3) pp:2086-2090
Publication Date(Web):2017/02/28
DOI:10.1039/C6SC04354B
The formation of two-component molecular cocrystals can lead to the tunable excited state intramolecular proton transfer (ESIPT) process and emission, as first confirmed by both experimental and computational studies.
Co-reporter:Xiang-Yang Liu;Wei-Hai Fang
Theoretical Chemistry Accounts 2017 Volume 136( Issue 1) pp:
Publication Date(Web):2017 January
DOI:10.1007/s00214-016-2029-z
Two-state conical intersection optimization methods at both QM and QM/MM levels have been extensively implemented in many commercial and noncommercial packages in the past decade. In contrast, three-state conical intersection optimization methods are less concerned, in particular the QM/MM-based ones. In this work, we have developed a penalty function-based three-state conical intersection optimization approach in the framework of the QM/MM method. We first present the fundamental formulation of this approach, and its algorithm and implementation in our package; then, we have carried out several pilot applications on molecular systems in vacuo and aqueous solution to demonstrate the efficiency of our implemented method. Our current developments enable efficient determination of three-state conical intersection structures of molecules in solution and biological systems, which is at the heart of understanding the photophysical and photochemical mechanisms of large systems at the atomistic level.
Co-reporter:Ya-Ting Wang, Yuan-Jun Gao, Qian Wang, and Ganglong Cui
The Journal of Physical Chemistry A 2017 Volume 121(Issue 4) pp:
Publication Date(Web):January 4, 2017
DOI:10.1021/acs.jpca.6b11682
Intramolecularly bridged diarylethenes exhibit improved photocyclization quantum yields because the anti-syn isomerization that originally suppresses photocyclization in classical diarylethenes is blocked. Experimentally, three possible channels have been proposed to interpret experimental observation, but many details of photochromic mechanism remain ambiguous. In this work we have employed a series of electronic structure methods (OM2/MRCI, DFT, TDDFT, RI-CC2, DFT/MRCI, and CASPT2) to comprehensively study excited state properties, photocyclization, and photoreversion dynamics of 1,2-dicyano[2,2]metacyclophan-1-ene. On the basis of optimized stationary points and minimum-energy conical intersections, we have refined experimentally proposed photochromic mechanism. Only an S1/S0 minimum-energy conical intersection is located; thus, we can exclude the third channel experimentally proposed. In addition, we find that both photocyclization and photoreversion processes use the same S1/S0 conical intersection to decay the S1 system to the S0 state, so we can unify the remaining two channels into one. These new insights are verified by our OM2/MRCI nonadiabatic dynamics simulations. The S1 excited-state lifetimes of photocyclization and photoreversion are estimated to be 349 and 453 fs, respectively, which are close to experimentally measured values: 240 ± 60 and 250 fs in acetonitrile solution. The present study not only interprets experimental observations and refines previously proposed mechanism but also provides new physical insights that are valuable for future experiments.
Co-reporter:Bin-Bin Xie, Xiang-Yang Liu, Qiu Fang, Wei-Hai Fang, and Ganglong Cui
The Journal of Physical Chemistry Letters 2017 Volume 8(Issue 5) pp:
Publication Date(Web):February 14, 2017
DOI:10.1021/acs.jpclett.6b02897
We have employed combined electronic structure calculations and nonadiabatic dynamics simulations to study the S1 radiationless deactivation mechanism of pyrazole. In terms of MS-CASPT2 computed results, we propose that the 1πσ* state-driven nonadiabatic N–N dissociation is a major relaxation path; the ring-puckering deformation path as well as the 1πσ* state-driven N–H dissociation are less favorable. This excited-state decay mechanism is supported by MS-CASPT2 nonadiabatic dynamics simulations. The present study demonstrates that pyrazole has a different excited-state radiationless deactivation mechanism compared with its structural isomer imidazole, in which the 1πσ* state-driven nonadiabatic N–H dissociation plays a more important role. However, such a channel is suppressed in pyrazole; instead, the 1πσ* state-driven nonadiabatic N–N dissociation is dominant.
Co-reporter:Xiang-Yang Liu, Xue-Ping Chang, Shu-Hua Xia, Ganglong Cui, and Walter Thiel
Journal of Chemical Theory and Computation 2016 Volume 12(Issue 2) pp:753-764
Publication Date(Web):January 8, 2016
DOI:10.1021/acs.jctc.5b00894
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.
Co-reporter:Bin-Bin Xie, Shu-Hua Xia, Xue-Ping Chang and Ganglong Cui  
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 1) pp:403-413
Publication Date(Web):06 Nov 2015
DOI:10.1039/C5CP05312A
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.
Co-reporter:Dan Wu;Wei-Wei Guo;Xiang-Yang Liu ;Dr. Ganglong Cui
ChemPhysChem 2016 Volume 17( Issue 15) pp:2340-2347
Publication Date(Web):
DOI:10.1002/cphc.201600386

Abstract

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 S0 and S1 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 S1 relaxation pathways from the Franck–Condon point of the BFP chromophore variant. One is the ESIPT path to generate an S1 tautomer that exhibits a large Stokes shift in experiments. The generated S1 tautomer can further evolve toward the nearby S1/S0 conical intersection and then jumps down to the S0 state. The other is the photoisomerization path along the rotation of the central double bond. Along this path, the S1 system runs into an S1/S0 conical intersection region and eventually hops to the S0 state. The two energetically allowed S1 excited-state deactivation pathways are responsible for the in-part loss of fluorescence quantum yield. The considerable S1 ESIPT barrier and the sizable barriers that separate the S1 tautomers from the S1/S0 conical intersections make these two tautomers establish a kinetic equilibrium in the S1 state, which thus results in dual fluorescence emission.

Co-reporter:Xue-Ping Chang, Xiao-Ying Xie, Shi-Yun Lin, and Ganglong Cui
The Journal of Physical Chemistry A 2016 Volume 120(Issue 31) pp:6129-6136
Publication Date(Web):July 15, 2016
DOI:10.1021/acs.jpca.6b02669
Compared with isoalloxazine, the core chromophore of biologically important flavins, alloxazine exhibits much lower fluorescence quantum yield and larger intersystem-crossing quantum yield. However, its efficient radiationless relaxation pathways are still elusive. In this work, we have used the QM(MS-CASPT2//CASSCF)/MM method to explore the mechanistic photophysics of alloxazine chromophore in aqueous solution. On the basis of the optimized minima, conical intersections, and crossing points in the lowest 1ππ*, 1nπ*, 3ππ*, and 3nπ* states, we have proposed three energetically possible nonadiabatic relaxation pathways populating the lowest 3ππ* triplet state from the initially populated excited 1ππ* singlet state. The first is the direct 1ππ*→ 3ππ* intersystem crossing via the 1ππ*/3ππ* crossing point. The second is an indirect 1ππ* → 3ππ* intersystem crossing relayed by the dark 1nπ* singlet state. In this route, the 1ππ* system first decays to the 1nπ* state via the 1ππ*/1nπ* conical intersection, followed by an 1nπ*→ 3ππ* intersystem crossing at the 1nπ*/3ππ* crossing point to arrive at the final 3ππ* state. The third is similar to the second one; but its intersystem crossing is relayed by the 3nπ* triplet state. The 1ππ* system first decays to the 3nπ* state via the 1ππ*/3nπ* crossing point; the generated 3nπ* state is then de-excited to the 3ππ* state through the 3nπ*→ 3ππ* internal conversion at the 3nπ*/3ππ* conical intersection. According to the classical El-Sayed rule, we suggest the second and third paths play a much more important role than the first one in the formation of the lowest 3ππ* state.
Co-reporter:Xiao-Ying Xie, Chun-Xiang Li, Qiu Fang, and Ganglong Cui
The Journal of Physical Chemistry A 2016 Volume 120(Issue 30) pp:6014-6022
Publication Date(Web):July 11, 2016
DOI:10.1021/acs.jpca.6b05899
Herein we computationally studied the excited-state properties and decay dynamics of methyl-4-hydroxycinnamate (OMpCA) in the lowest three electronic states, that is, 1ππ*, 1nπ*, and S0 using combined MS-CASPT2 and CASSCF electronic structure methods. We found that one-water hydration can significantly stabilize and destabilize the vertical excitation energies of the spectroscopically bright 1ππ* and dark 1nπ* excited singlet states, respectively; in contrast, it has a much smaller effect on the 1ππ* and 1nπ* adiabatic excitation energies. Mechanistically, we located two 1ππ* excited-state relaxation channels. One is the internal conversion to the dark 1nπ* state, and the other is the 1ππ* photoisomerization that eventually leads the system to a 1ππ*/S0 conical intersection region, near which the radiationless internal conversion to the S0 state occurs. These two 1ππ* relaxation pathways play distinct roles in OMpCA and its two one-water complexes (OMpCA-W1 and OMpCA-W2). In OMpCA, the predominant 1ππ* decay route is the state-switching to the dark 1nπ* state, while in one-water complexes, the importance of the 1ππ* photoisomerization is significantly enhanced because the internal conversion to the 1nπ* state is heavily suppressed due to the one-water hydration.
Co-reporter:Dr. Shu-Hua Xia;Dr. Ganglong Cui;Dr. Wei-Hai Fang;Dr. Walter Thiel
Angewandte Chemie 2016 Volume 128( Issue 6) pp:2107-2112
Publication Date(Web):
DOI:10.1002/ange.201509622

Abstract

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.

Co-reporter:Dr. Shu-Hua Xia;Dr. Ganglong Cui;Dr. Wei-Hai Fang;Dr. Walter Thiel
Angewandte Chemie International Edition 2016 Volume 55( Issue 6) pp:2067-2072
Publication Date(Web):
DOI:10.1002/anie.201509622

Abstract

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.

Co-reporter:Shu-Hua Xia, Bin-Bin Xie, Qiu Fang, Ganglong Cui and Walter Thiel  
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 15) pp:9687-9697
Publication Date(Web):12 Feb 2015
DOI:10.1039/C5CP00101C
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.
Co-reporter:Dr. Xue-Ping Chang; Dr. Ganglong Cui; Dr. Wei-Hai Fang; Dr. Walter Thiel
ChemPhysChem 2015 Volume 16( Issue 5) pp:933-937
Publication Date(Web):
DOI:10.1002/cphc.201402897

Abstract

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 1ππ* state of DDQ can relax efficiently through a nearby dark 1nπ* doorway state to the 3ππ* 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 T1S0 intersystem crossing and concomitant proton-coupled electron transfer (PCET) to generate the products DDQH2 and phenol. Two of the four considered nonadiabatic pathways involve an orientation-dependent triplet PCET process, followed by intersystem crossing to the ground state (S0). The other two first undergo a nonadiabatic T1S0 transition to produce a zwitterionic S0 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.

Co-reporter:Pei-Jie Guan; Ganglong Cui;Dr. Qiu Fang
ChemPhysChem 2015 Volume 16( Issue 4) pp:805-811
Publication Date(Web):
DOI:10.1002/cphc.201402743

Abstract

We employed the complete active space self-consistent field (CASSCF) and its multistate second-order perturbation (MS-CASPT2) methods to explore the photochemical mechanism of 2-hydroxyazobenzene, the molecular scaffold of Sudan I and Orange II dyes. It was found that the excited-state intramolecular proton transfer (ESIPT) along the bright diabatic 1ππ* state is barrierless and ultrafast. Along this diabatic 1ππ* relaxation path, the system can jump to the dark 1nπ* state via the 1ππ*/1nπ* crossing point. However, ESIPT in this dark state is largely inhibited owing to a sizeable barrier. We also found two deactivation channels that decay 1ππ* keto and 1nπ* enol species to the ground state via two energetically accessible S1/S0 conical intersections. Finally, we encountered an interesting phenomenon in the excited-state hydrogen-bonding strength: it is reinforced in the 1ππ* state, whereas it is reduced in the 1nπ* state. The present work sets the stage for understanding the photophysics and photochemistry of Sudan I–IV, Orange II, Ponceau 2R, Ponceau 4R, and azo violet.

Co-reporter:Xue-Ping Chang, Chun-Xiang Li, Bin-Bin Xie, and Ganglong Cui
The Journal of Physical Chemistry A 2015 Volume 119(Issue 47) pp:11488-11497
Publication Date(Web):October 29, 2015
DOI:10.1021/acs.jpca.5b08434
p-Methoxy methylcinnamate (p-MMC) shares the same molecular skeleton with octyl methoxycinnamate sunscreen. It is recently found that adding one water to p-MMC can significantly enhance the photoprotection efficiency. However, the physical origin is elusive. Herein we have employed multireference complete active space self-consistent field (CASSCF) and multistate complete active-space second-order perturbation (MS-CASPT2) methods to scrutinize the photophysical and photochemical mechanism of p-MMC and its one-water complex p-MMC–W. Specifically, we optimize the stationary-point structures on the 1ππ*, 1nπ*, and S0 potential energy surfaces to locate the 1ππ*/S0 and 1ππ*/1nπ* conical intersections and to map 1ππ* and 1nπ* excited-state relaxation paths. On the basis of the results, we find that, for the trans p-MMC, the major 1ππ* deactivation path is decaying to the dark 1nπ* state via the in-plane 1ππ*/1nπ* crossing point, which only need overcome a small barrier of 2.5 kcal/mol; the minor one is decaying to the S0 state via the 1ππ*/S0 conical intersection induced by out-of-plane photoisomerization. For the cis p-MMC, these two decay paths are comparable 1ππ* deactivation paths: one is decaying to the dark 1nπ* state via the 1ππ*/1nπ* crossing point, and the second is decaying to the ground state via the 1ππ*/S0 conical intersection. One-water hydration stabilizes the 1ππ* state and meanwhile destabilizes the 1nπ* state. As a consequence, the 1ππ* deactivation path to the dark 1nπ* state is heavily inhibited. The related barriers are increased to 5.8 and 3.3 kcal/mol for the trans and cis p-MMC–W, respectively. In comparison, the barriers associated with the photoisomerization-induced 1ππ* decay paths are reduced to 2.5 and 1.3 kcal/mol for the trans and cis p-MMC–W. Therefore, the 1ππ* decay paths to the S0 state are dominant relaxation channels when adding one water molecule. Finally, the present work contributes a lot of knowledge to understanding the photoprotection mechanism of methylcinnamate derivatives.
Co-reporter:Bin-Bin Xie, Shu-Hua Xia, Li-Hong Liu, and Ganglong Cui
The Journal of Physical Chemistry A 2015 Volume 119(Issue 22) pp:5607-5617
Publication Date(Web):May 6, 2015
DOI:10.1021/acs.jpca.5b02549
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.
Co-reporter:Shu-Hua Xia, Xiang-Yang Liu, Qiu Fang, and Ganglong Cui
The Journal of Physical Chemistry A 2015 Volume 119(Issue 15) pp:3569-3576
Publication Date(Web):March 25, 2015
DOI:10.1021/acs.jpca.5b00302
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.
Co-reporter:Ganglong Cui, Pei-Jie Guan, and Wei-Hai Fang
The Journal of Physical Chemistry A 2014 Volume 118(Issue 26) pp:4732-4739
Publication Date(Web):June 6, 2014
DOI:10.1021/jp503865y
Intramolecularly hydrogen-bonded aromatic azo compound 1-cyclopropyldiazo-2-naphthol (CPDNO) exhibits complicated excited-state behaviors, e.g., wavelength-dependent photoinduced proton transfer and photoproducts. Its photochemistry differs from that of common aromatic azo compounds in which cis–trans photoisomerization is dominant. To rationalize the intriguing photochemistry of CPDNO at the atomic level, we have in this work employed the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods to explore the S0, S1, and S2 potential-energy profiles relevant to enol–keto proton transfer and isomerization reactions. It is found that the proton transfer along the bright diabatic 1ππ* potential-energy profile is almost barrierless, quickly forming the fluorescent 1ππ* keto minimum. In this process, the dark 1nπ* state is populated via a 1ππ*/1nπ* crossing point, but the proton transfer on this dark state is suppressed heavily as a result of a large barrier. In addition, two deactivation paths that decay the S1 enol and keto minima to the S0 state, respectively, were uncovered. For the former, it is exoenergetic and thereby thermodynamically favorable; for the latter, it is a little endothermic (ca. 5 kcal/mol). Both are energetically allowable concerning the available total energy. Finally, on the basis of the present results, the experimentally observed wavelength-dependent photoproducts were explained very well.
Co-reporter:Ganglong Cui and Walter Thiel
The Journal of Physical Chemistry Letters 2014 Volume 5(Issue 15) pp:2682-2687
Publication Date(Web):July 22, 2014
DOI:10.1021/jz501159j
Motivated by its potential use as a photosensitizer in photodynamic therapy, we report the first ab initio quantum mechanics/molecular mechanics (QM/MM) study of 4-thiothymidine in aqueous solution. The core chromophore 4-thiothymine was described using the multiconfigurational CASSCF and CASPT2 QM methods, while the ribose and the solvent water molecules were treated at the MM level (CHARMM and TIP3P, respectively). The minima of the five lowest electronic states (S0, S1, S2, T1, and T2) and six minimum-energy intersections were fully optimized at the QM(CASSCF)/MM level, and their energies were further refined by single-point QM(CASPT2)/MM and CASPT2 calculations. The relevant spin–orbit couplings were also computed. We find that (1) there are three efficient photophysical pathways that account for the experimentally observed ultrafast formation of the lowest triplet state with a quantum yield of nearly unity, (2) the striking qualitative differences in the photophysical behavior of 4-thiothymine and thymine originate from the different electronic structure of their S1 states, and (3) environmental effects play an important role. The present QM/MM calculations provide mechanistic insight that may guide the design of improved photosensitizers for photodynamic therapy.Keywords: 4-thiothymidine; intersystem crossing; photodynamic therapy; photosensitizers; quantum mechanics/molecular mechanics; spin−orbit coupling;
Co-reporter:Bin-Bin Xie, Shu-Hua Xia, Xue-Ping Chang and Ganglong Cui
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 1) pp:NaN413-413
Publication Date(Web):2015/11/06
DOI:10.1039/C5CP05312A
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.
Co-reporter:Pin Xiao, Dan Wu, Wei-Hai Fang and Ganglong Cui
Catalysis Science & Technology (2011-Present) 2017 - vol. 7(Issue 13) pp:NaN2771-2771
Publication Date(Web):2017/05/17
DOI:10.1039/C7CY00785J
A novel light-triggered hydrogen evolution reaction from formic acid mediated by an Ir(III) photocatalyst has been experimentally reported recently. However, its reaction mechanism remains elusive. Herein, we have employed the density functional theory (DFT) method to explore this photocatalytic reaction in detail. On the basis of the results, we have proposed a possible photocatalytic reaction mechanism. In the formation of the metal hydride [Cp*Ir(bpy)(H)]+ (5), formic acid acts as a bridge assisting proton shuttling. Upon irradiation, two nonadiabatic excited-state decay pathways quickly populate the lowest triplet T1 state of the metal hydride from its initially populated excited singlet S1 state. In the T1 state, water and formic acid facilitate excited-state hydride/proton transfers from the Ir center to Cp* and bpy ligands producing several energetically lower triplet-state isomers demonstrating that the triplet-state metal hydride 5* could not be the only precursor for the photocatalysis. Adiabatic H2 evolution in the T1 state is energetically unfavorable. These T1 isomers hop, through radiationless T1 → S0 intersystem crossings via T1/S0 crossing points, to the S0 state in which H2 evolution takes place. In these reactions, solvents acting as assistants and catalysts reduce reaction barriers, thereby accelerating H2 release and enhancing the overall photocatalytic performance. Our current work provides significant mechanistic insights into light-induced hydrogen-evolution reactions of iridium-containing photocatalysts.
Co-reporter:Shu-Hua Xia, Bin-Bin Xie, Qiu Fang, Ganglong Cui and Walter Thiel
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 15) pp:NaN9697-9697
Publication Date(Web):2015/02/12
DOI:10.1039/C5CP00101C
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.
Co-reporter:Xue-Ping Chang, Yiying Zheng, Ganglong Cui, Wei-Hai Fang and Walter Thiel
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 35) pp:NaN24721-24721
Publication Date(Web):2016/08/16
DOI:10.1039/C6CP04919B
We report a density functional study on the mechanism of the [2+2] photocyclization of atropisomeric maleimides. Experimentally, the reaction is known to proceed through the triplet state. We have located all relevant S0 and T1 minima and transition states, as well as the T1/S0 crossing points, and mapped eight stepwise photocyclization pathways for four different conformers in the T1 state that lead to distinct regioisomers. In the preferred four pathways (one for each conformer) the initially formed C–C bond involves the terminal carbon atom of the alkene moiety. This regioselectivity originates from electrostatic preferences (arising from the charge distribution in the polarized CC double bonds) and from the different thermodynamic stability of the resulting triplet diradical intermediates (caused by electron donation effects that stabilize the radical centers). The formation of the second C–C bond is blocked in the T1 state by prohibitively high barriers and thus occurs after intersystem crossing to the ground state. Furthermore, we rationalize substitution effects on enantioselectivity and diastereoselectivity and identify their origin.
Co-reporter:Bin-Bin Xie, Qian Wang, Wei-Wei Guo and Ganglong Cui
Physical Chemistry Chemical Physics 2017 - vol. 19(Issue 11) pp:NaN7698-7698
Publication Date(Web):2017/02/16
DOI:10.1039/C7CP00478H
The photophysics of thiothymines has been extensively studied computationally in the past few years due to their significant potential as photosensitizers in photodynamic therapy. However, the corresponding computational studies of the photophysical mechanism of 2,4-dithiothymine are scarce. Herein we have employed the CASPT2//CASSCF and QM(CASPT2//CASSCF)/MM methods to systematically explore the excited-state decay mechanism of 2,4-dithiothymine in isolated, microsolvated, and aqueous surroundings. First, we have optimized minima and conical intersections in and between the lowest six excited singlet and triplet states i.e., , , , , and ; then, based on computed excited-state decay paths and spin–orbit couplings, we have proposed several nonadiabatic pathways that efficiently populate the lowest triplet state to explain the experimentally observed ultrahigh triplet-state quantum yield. Moreover, we have found that the excited-state decay mechanism in microsolvated and aqueous environments is more complicated than that in the gas phase. The solute–solvent interaction has significant effects on the excited-state potential energy surfaces of 2,4-dithiothymine and eventually on its excited-state decay mechanism. Finally, the present computational efforts contribute important mechanistic knowledge to the understanding of the photophysics of thiothymine-based photosensitizers.
Co-reporter:Xue-Ping Chang, Pin Xiao, Juan Han, Wei-Hai Fang and Ganglong Cui
Physical Chemistry Chemical Physics 2017 - vol. 19(Issue 21) pp:NaN13533-13533
Publication Date(Web):2017/04/26
DOI:10.1039/C7CP01511A
In contrast to photophysics of thio-substituted nucleobases, their photoinduced cross-linking reactions with canonical nucleobases remain scarcely investigated computationally. In this work, we have adopted combined CASPT2/PCM//CASSCF and B3LYP-D3/PCM electronic structure methods to study this kind of photochemical reaction of 5-fluoro-4-thiouridine (truncated 5-fluoro-1-methyl-4-thiouracil used in calculations) and 1-methylthymine (referred to as thymine for clarity hereinafter). On the basis of CASPT2/PCM computed results, we have proposed two efficient excited-state relaxation pathways to populate the lowest T1 state of the complex of 5-fluoro-1-methyl-4-thiouracil and thymine from its initially populated S2(1ππ*) state. In the first one, the S2 system first hops to the S1 state via an S2/S1 conical intersection, followed by a direct S1 → T1 intersystem crossing process enhanced by large S1/T1 spin–orbit coupling. In the second path, the resultant S1 system first jumps to the T2 state, from which an efficient T2 → T1 internal conversion occurs. The T1 cross-linking reaction is overall divided into two phases. The first phase is a stepwise and nonadiabatic photocyclization reaction, which starts from the T1 complex and ends up with an S0 thietane intermediate. The second phase is a thermal reaction. The system first rearranges its four- and six-membered rings to form three new rings; then, an S0 fluorine atom transfer occurs, followed by the formation of photoproducts. Finally, the present work paves the way for studying light-induced cross-linking reactions of thionucleobases with canonical bases in DNA and RNA.
Co-reporter:Heyang Lin, Xueping Chang, Dongpeng Yan, Wei-Hai Fang and Ganglong Cui
Chemical Science (2010-Present) 2017 - vol. 8(Issue 3) pp:NaN2090-2090
Publication Date(Web):2016/11/14
DOI:10.1039/C6SC04354B
The formation of two-component molecular cocrystals can lead to the tunable excited state intramolecular proton transfer (ESIPT) process and emission, as first confirmed by both experimental and computational studies.
1,4-Butanediyl, 1-oxo-
Phenol, 2-(4-pyridinyl)-
1,3-Propanediyl, 1-oxo-
TRICYCLO[9.3.1.14,8]HEXADECA-1(15),2,4,6,8(16),11,13-HEPTAENE-2,3-DICARBONITRILE, 6,13-BIS(1,1-DIMETHYLETHYL)-15,16-DIMETHYL-
2-NAPHTHALENOL, 1-(CYCLOPROPYLAZO)-
Cyclopropanone
4H-Imidazol-4-one,3,5-dihydro-5-[(4-hydroxyphenyl)methylene]-2,3-dimethyl-
3H-Diazirin-3-one
(E)-p-Methoxycinnamic acid methyl ester
ETHANONE