Herein we have employed high-level multi-reference CASSCF and MS-CASPT2 electronic structure methods to systematically study the photochemical mechanism of intramolecularly hydrogen-bonded 2-（2＇-hydroxyphenyl）-4-me...Herein we have employed high-level multi-reference CASSCF and MS-CASPT2 electronic structure methods to systematically study the photochemical mechanism of intramolecularly hydrogen-bonded 2-（2＇-hydroxyphenyl）-4-methyloxazole. At the CASSCF level, we have optimized minima, conical intersections, minimum-energy reaction paths relevant to the excited-state intramolecular proton transfer （ESIPT）, rotation, photoisomerization, and the excited-state deactivation pathways. The energies of all structures and paths are refined by the MS-CASPT2 method. On the basis of the present results, we found that the ESIPT process in a conformer with the OH... N hydrogen bond is essentially barrierless process; whereas, the ESIPT process is inhibited in the other conformer with the OH... O hydrogen bond. The central single-bond rotation of the S1 enol species is energetically unfavorable due to a large barrier. In addition, the excited-state deactivation of the S1 keto species, as a result of the ultrafast ESIPT, is very efficient because of the existence of two easily-approached keto S1/S0 conical intersections. In stark contrast to the S1 keto species, the decay of the S1 enol species is almostly blocked. The present theoretical study contributes valuable knowledge to the understanding of photochemistry of similar intramolecularly hydrogen-bonded molecular and biological systems.展开更多
In this work,we employ electronic structure calculations and nonadiabatic dynamics simulations based on many-body Green function and BetheSalpeter equation(GW/BSE)methods to study excited-state properties of a zinc ph...In this work,we employ electronic structure calculations and nonadiabatic dynamics simulations based on many-body Green function and BetheSalpeter equation(GW/BSE)methods to study excited-state properties of a zinc phthalocyanine-fullerene(ZnPcC_(60))dyad with 6-6 and 5-6 configurations.In the former,the initially populated locally excited(LE)state of ZnPc is the lowest S1 state and thus,its subsequent charge separation is relatively slow.In contrast,in the latter,the S1 state is the LE state of C_(60)while the LE state of ZnPc is much higher in energy.There also exist several charge-transfer(CT)states between the LE states of ZnPc and C_(60).Thus,one can see apparent charge separation dynamics during excited-state relaxation dynamics from the LE state of ZnPc to that of C_(60).These points are verified in dynamics simulations.In the first 200 fs,there is a rapid excitation energy transfer from ZnPc to C_(60),followed by an ultrafast charge separation to form a CT intermediate state.This process is mainly driven by hole transfer from C_(60)to ZnPc.The present work demonstrates that different bonding patterns(i.e.5-6 and 6-6)of the C−N linker can be used to tune excited-state properties and thereto optoelectronic properties of covalently bonded ZnPc-C_(60)dyads.Methodologically,it is proven that combined GW/BSE nonadiabatic dynamics method is a practical and reliable tool for exploring photoinduced dynamics of nonperiodic dyads,organometallic molecules,quantum dots,nanoclusters,etc.展开更多
Understanding the excited state dynamics of donor-acceptor(D-A)complexes is of fundamental importance both experimentally and theoretically.Herein,we have first explored the photoinduced dynamics of a recently synthes...Understanding the excited state dynamics of donor-acceptor(D-A)complexes is of fundamental importance both experimentally and theoretically.Herein,we have first explored the photoinduced dynamics of a recently synthesized paddle-wheel BODIPY-hexaoxatriphenylene(BODIPY is the abbreviation for BF_(2)-chelated dipyrromethenes)conjugates D-A complexes with the combination of both electronic structure calculations and nonadiabatic dynamics simulations.On the basis of computational results,we concluded that the BODIPY-hexaoxatriphenylene(BH)conjugates will be promoted to the local excited(LE)states of the BODIPY fragments upon excitation,which is followed by the ultrafast exciton transfer from LE state to charge transfer(CT).Instead of the photoinduced electron transfer process proposed in previous experimental work,such a exciton transfer process is accompanied with the photoinduced hole transfer from BODIPY to hexaoxatriphenylene.Additionally,solvent effects are found to play an important role in the photoinduced dynamics.Specifically,the hole transfer dynamics is accelerated by the acetonitrile solvent,which can be ascribed to significant influences of the solvents on the charge transfer states,i.e.the energy gaps between LE and CT excitons are reduced greatly and the non-adiabatic couplings are increased in the meantime.Our present work not only provides valuable insights into the underlying photoinduced mechanism of BH,but also can be helpful for the future design of novel donor-acceptor conjugates with better optoelectronic performance.展开更多
As the major and abundant type of glucosinolates(GL)in plants,sinigrin has potential functions in promoting health and insect defense.The final step in the biosynthesis of sinigrin core structure is highly representat...As the major and abundant type of glucosinolates(GL)in plants,sinigrin has potential functions in promoting health and insect defense.The final step in the biosynthesis of sinigrin core structure is highly representative in GL compounds,which corresponds to the process from 3-methylthiopropyl ds-GL to 3-methylthiopropyl GL catalyzed by sulfotransferase(SOT).However,due to the lack of the crystallographic structure of SOT complexed with the 3-methylthiopropyl GL,little is known about this sulfonation process.Fortunately,the crystal structure of SOT 18 from Arabidopsis thaliana(At SOT18)containing the substance(sinigrin)similar to 3-methylthiopropyl GL has been determined.To understand the enzymatic mechanism,we employed molecular dynamics(MD)simulation and quantum mechanics combined with molecular mechanics(QM/MM)methods to study the conversion from ds-sinigrin to sinigrin catalyzed by AtSOT18.The calculated results demonstrate that the reaction occurs through a concerted dissociative mechanism.Moreover,Lys93,Thr96,Thr97,Tyr130,His155,and two enzyme peptide chains(Pro92-Lys93 and Gln95-Thr96-Thr97)play a role in positioning the substrates and promoting the catalytic reaction by stabilizing the transition state geometry.Particularly,His155 acts as a catalytic base while Lys93 acts as a catalytic acid in the reaction process.The presently proposed concerted dissociative mechanism explains the role of At SOT18 in sinigrin biosynthesis,and could be instructive for the study of GL biosynthesis catalyzed by other SOTs.展开更多
Herein we have employed the MS-CASPT2//CASSCF method to study the S1 excited-state intramolecular proton transfers （ESIPTs） of recently synthesized ortho-hydroxyl GFP core chromophores, i.e. OHIM, CHBDI, and MHBID, ...Herein we have employed the MS-CASPT2//CASSCF method to study the S1 excited-state intramolecular proton transfers （ESIPTs） of recently synthesized ortho-hydroxyl GFP core chromophores, i.e. OHIM, CHBDI, and MHBID, and their excited-state relaxation pathways. We have found that in OHIM and CHBDI, the ESIPT process is associated with small barriers of 3.4 and 4.2 kcal/mol; while, in MHBDI, it becomes essentially barrierless. Moreover, we have found two main S1 excited-state radiationless channels. In the first one, the enol S1 species decays to the So state via the enol S1/S0 conical intersection after overcoming considerable barriers of 7.0 and 7.7 kcal/mol in OHIM and CHBDI （however, in MHBDI, it is nearly barrierless）. In the second one, the keto S] species is first generated through the ESIPT event; then, it is de-excited into the So state in the vicinity of the keto S1/S0 conical intersection. These energetically allowed excited-state decay channels rationalize ex- perimentally observed ultralow fluorescence quantum yields. The insights gained from the present work may help to guide the design of new ortho-hydroxyl GFP core chromophores with improved fluorescence emission and brightness.展开更多
In(relativistic)electronic structure methods,the quaternion matrix eigenvalue problem and the linear response(Bethe-Salpeter)eigenvalue problem for excitation energies are two frequently encoun-tered structured eigenv...In(relativistic)electronic structure methods,the quaternion matrix eigenvalue problem and the linear response(Bethe-Salpeter)eigenvalue problem for excitation energies are two frequently encoun-tered structured eigenvalue problems.While the former problem was thoroughly studied,the later problem in its most general form,namely,the complex case without assuming the positive definiteness of the electronic Hessian,was not fully understood.In view of their very similar mathematical structures,we examined these two problems from a unified point of view.We showed that the identification of Lie group structures for their eigenvectors provides a framework to design diagonalization algorithms as well as numerical optimizations techniques on the corresponding manifolds.By using the same reduction algorithm for the quaternion matrix eigenvalue problem,we provided a necessary and sufficient condition to characterize the different scenarios,where the eigenvalues of the original linear response eigenvalue problem are real,purely imaginary,or complex.The result can be viewed as a natural generalization of the well-known condition for the real matrix case.展开更多
Excited-state double proton transfer(ESDPT)is a controversial issue which has long been plagued with theoretical and experimental communities.Herein,we took 1,8-dihydroxy-2-naphthaldehyde(DHNA)as a prototype and used ...Excited-state double proton transfer(ESDPT)is a controversial issue which has long been plagued with theoretical and experimental communities.Herein,we took 1,8-dihydroxy-2-naphthaldehyde(DHNA)as a prototype and used combined complete active space selfconsistent field(CASSCF)and multi-state complete active-space second-order perturbation(MS-CASPT2)methods to investigate ESDPT and excited-state deactivation pathways of DHNA.Three different tautomer minima of S1-ENOL,S1-KETO-1,and S1-KETO-2 and two crucial conical intersections of S1 S0-KETO-1 and S1 S0-KETO-2 in and between the S0 and S1 states were obtained.S1-KETO-1 and S1-KETO-2 should take responsibility for experimentally observing dual-emission bands.In addition,two-dimensional potential energy surfaces(2 D-PESs)and linear interpolated internal coordinate paths connecting relevant structures were calculated at the MS-CASPT2//CASSCF level and confirmed a stepwise ESDPT mechanism.Specifically,the first proton transfer from S1-ENOL to S1-KETO-1 is barrierless,whereas the second one from S1-KETO-1 to S1-KETO-2 demands a barrier of ca.6.0 kcal/mol.The linear interpolated internal coordinate path connecting S1-KETO-1(S1-KETO-2)and S_(1) S0-KETO-1(S1 S0-KETO-2)is uphill with a barrier of ca.12.0 kcal/mol,which will trap DHNA in the S_(1) state while therefore enabling dual-emission bands.On the other hand,the S1/S0 conical intersections would also prompt the S_(1) system to decay to the S_(0) state,which could be to certain extent suppressed by locking the rotation of the C5-C8-C9-O10 dihedral angle.These mechanistic insights are not only helpful for understanding ESDPT but also useful for designing novel molecular materials with excellent photoluminescent performances.展开更多
文摘Herein we have employed high-level multi-reference CASSCF and MS-CASPT2 electronic structure methods to systematically study the photochemical mechanism of intramolecularly hydrogen-bonded 2-（2＇-hydroxyphenyl）-4-methyloxazole. At the CASSCF level, we have optimized minima, conical intersections, minimum-energy reaction paths relevant to the excited-state intramolecular proton transfer （ESIPT）, rotation, photoisomerization, and the excited-state deactivation pathways. The energies of all structures and paths are refined by the MS-CASPT2 method. On the basis of the present results, we found that the ESIPT process in a conformer with the OH... N hydrogen bond is essentially barrierless process; whereas, the ESIPT process is inhibited in the other conformer with the OH... O hydrogen bond. The central single-bond rotation of the S1 enol species is energetically unfavorable due to a large barrier. In addition, the excited-state deactivation of the S1 keto species, as a result of the ultrafast ESIPT, is very efficient because of the existence of two easily-approached keto S1/S0 conical intersections. In stark contrast to the S1 keto species, the decay of the S1 enol species is almostly blocked. The present theoretical study contributes valuable knowledge to the understanding of photochemistry of similar intramolecularly hydrogen-bonded molecular and biological systems.
基金support from the National Natural Science Foundation of China(No.21688102,No.21590801,and No.21520102005)support from Sichuan Science and Technology Program Grant(2020YJ0161)。
文摘In this work,we employ electronic structure calculations and nonadiabatic dynamics simulations based on many-body Green function and BetheSalpeter equation(GW/BSE)methods to study excited-state properties of a zinc phthalocyanine-fullerene(ZnPcC_(60))dyad with 6-6 and 5-6 configurations.In the former,the initially populated locally excited(LE)state of ZnPc is the lowest S1 state and thus,its subsequent charge separation is relatively slow.In contrast,in the latter,the S1 state is the LE state of C_(60)while the LE state of ZnPc is much higher in energy.There also exist several charge-transfer(CT)states between the LE states of ZnPc and C_(60).Thus,one can see apparent charge separation dynamics during excited-state relaxation dynamics from the LE state of ZnPc to that of C_(60).These points are verified in dynamics simulations.In the first 200 fs,there is a rapid excitation energy transfer from ZnPc to C_(60),followed by an ultrafast charge separation to form a CT intermediate state.This process is mainly driven by hole transfer from C_(60)to ZnPc.The present work demonstrates that different bonding patterns(i.e.5-6 and 6-6)of the C−N linker can be used to tune excited-state properties and thereto optoelectronic properties of covalently bonded ZnPc-C_(60)dyads.Methodologically,it is proven that combined GW/BSE nonadiabatic dynamics method is a practical and reliable tool for exploring photoinduced dynamics of nonperiodic dyads,organometallic molecules,quantum dots,nanoclusters,etc.
基金supported by the National Natural Science Foundation of China(No.22003043 for Xiang-Yang Liu)the National Natural Science Foundation of China(No.21688102,No.21590801,and No.21520102005 for Ganglong Cui)+1 种基金Sichuan Science and Technology Program(No.2020YJ0161 for Xiang-Yang Liu)the High Performance Computing Center of Sichuan Normal University。
文摘Understanding the excited state dynamics of donor-acceptor(D-A)complexes is of fundamental importance both experimentally and theoretically.Herein,we have first explored the photoinduced dynamics of a recently synthesized paddle-wheel BODIPY-hexaoxatriphenylene(BODIPY is the abbreviation for BF_(2)-chelated dipyrromethenes)conjugates D-A complexes with the combination of both electronic structure calculations and nonadiabatic dynamics simulations.On the basis of computational results,we concluded that the BODIPY-hexaoxatriphenylene(BH)conjugates will be promoted to the local excited(LE)states of the BODIPY fragments upon excitation,which is followed by the ultrafast exciton transfer from LE state to charge transfer(CT).Instead of the photoinduced electron transfer process proposed in previous experimental work,such a exciton transfer process is accompanied with the photoinduced hole transfer from BODIPY to hexaoxatriphenylene.Additionally,solvent effects are found to play an important role in the photoinduced dynamics.Specifically,the hole transfer dynamics is accelerated by the acetonitrile solvent,which can be ascribed to significant influences of the solvents on the charge transfer states,i.e.the energy gaps between LE and CT excitons are reduced greatly and the non-adiabatic couplings are increased in the meantime.Our present work not only provides valuable insights into the underlying photoinduced mechanism of BH,but also can be helpful for the future design of novel donor-acceptor conjugates with better optoelectronic performance.
基金supported by the National Natural Science Foundation of China(No.21973005)。
文摘As the major and abundant type of glucosinolates(GL)in plants,sinigrin has potential functions in promoting health and insect defense.The final step in the biosynthesis of sinigrin core structure is highly representative in GL compounds,which corresponds to the process from 3-methylthiopropyl ds-GL to 3-methylthiopropyl GL catalyzed by sulfotransferase(SOT).However,due to the lack of the crystallographic structure of SOT complexed with the 3-methylthiopropyl GL,little is known about this sulfonation process.Fortunately,the crystal structure of SOT 18 from Arabidopsis thaliana(At SOT18)containing the substance(sinigrin)similar to 3-methylthiopropyl GL has been determined.To understand the enzymatic mechanism,we employed molecular dynamics(MD)simulation and quantum mechanics combined with molecular mechanics(QM/MM)methods to study the conversion from ds-sinigrin to sinigrin catalyzed by AtSOT18.The calculated results demonstrate that the reaction occurs through a concerted dissociative mechanism.Moreover,Lys93,Thr96,Thr97,Tyr130,His155,and two enzyme peptide chains(Pro92-Lys93 and Gln95-Thr96-Thr97)play a role in positioning the substrates and promoting the catalytic reaction by stabilizing the transition state geometry.Particularly,His155 acts as a catalytic base while Lys93 acts as a catalytic acid in the reaction process.The presently proposed concerted dissociative mechanism explains the role of At SOT18 in sinigrin biosynthesis,and could be instructive for the study of GL biosynthesis catalyzed by other SOTs.
文摘Herein we have employed the MS-CASPT2//CASSCF method to study the S1 excited-state intramolecular proton transfers （ESIPTs） of recently synthesized ortho-hydroxyl GFP core chromophores, i.e. OHIM, CHBDI, and MHBID, and their excited-state relaxation pathways. We have found that in OHIM and CHBDI, the ESIPT process is associated with small barriers of 3.4 and 4.2 kcal/mol; while, in MHBDI, it becomes essentially barrierless. Moreover, we have found two main S1 excited-state radiationless channels. In the first one, the enol S1 species decays to the So state via the enol S1/S0 conical intersection after overcoming considerable barriers of 7.0 and 7.7 kcal/mol in OHIM and CHBDI （however, in MHBDI, it is nearly barrierless）. In the second one, the keto S] species is first generated through the ESIPT event; then, it is de-excited into the So state in the vicinity of the keto S1/S0 conical intersection. These energetically allowed excited-state decay channels rationalize ex- perimentally observed ultralow fluorescence quantum yields. The insights gained from the present work may help to guide the design of new ortho-hydroxyl GFP core chromophores with improved fluorescence emission and brightness.
基金supported by the National Natural Science Foundation of China (No.21973003)the Beijing Normal University Startup Package
文摘In(relativistic)electronic structure methods,the quaternion matrix eigenvalue problem and the linear response(Bethe-Salpeter)eigenvalue problem for excitation energies are two frequently encoun-tered structured eigenvalue problems.While the former problem was thoroughly studied,the later problem in its most general form,namely,the complex case without assuming the positive definiteness of the electronic Hessian,was not fully understood.In view of their very similar mathematical structures,we examined these two problems from a unified point of view.We showed that the identification of Lie group structures for their eigenvectors provides a framework to design diagonalization algorithms as well as numerical optimizations techniques on the corresponding manifolds.By using the same reduction algorithm for the quaternion matrix eigenvalue problem,we provided a necessary and sufficient condition to characterize the different scenarios,where the eigenvalues of the original linear response eigenvalue problem are real,purely imaginary,or complex.The result can be viewed as a natural generalization of the well-known condition for the real matrix case.
基金supported by the National Key Research and Development Program of China for BinBin Xie(No.2019YFA0709400)the National Natural Science Foundation of China for Bin-Bin Xie(No.21903068)+1 种基金Xiang-Yang Liu(No.22003043)Natural Science Foundation of Zhejiang Province for Bin-Bin Xie(No.LQ19B030007)。
文摘Excited-state double proton transfer(ESDPT)is a controversial issue which has long been plagued with theoretical and experimental communities.Herein,we took 1,8-dihydroxy-2-naphthaldehyde(DHNA)as a prototype and used combined complete active space selfconsistent field(CASSCF)and multi-state complete active-space second-order perturbation(MS-CASPT2)methods to investigate ESDPT and excited-state deactivation pathways of DHNA.Three different tautomer minima of S1-ENOL,S1-KETO-1,and S1-KETO-2 and two crucial conical intersections of S1 S0-KETO-1 and S1 S0-KETO-2 in and between the S0 and S1 states were obtained.S1-KETO-1 and S1-KETO-2 should take responsibility for experimentally observing dual-emission bands.In addition,two-dimensional potential energy surfaces(2 D-PESs)and linear interpolated internal coordinate paths connecting relevant structures were calculated at the MS-CASPT2//CASSCF level and confirmed a stepwise ESDPT mechanism.Specifically,the first proton transfer from S1-ENOL to S1-KETO-1 is barrierless,whereas the second one from S1-KETO-1 to S1-KETO-2 demands a barrier of ca.6.0 kcal/mol.The linear interpolated internal coordinate path connecting S1-KETO-1(S1-KETO-2)and S_(1) S0-KETO-1(S1 S0-KETO-2)is uphill with a barrier of ca.12.0 kcal/mol,which will trap DHNA in the S_(1) state while therefore enabling dual-emission bands.On the other hand,the S1/S0 conical intersections would also prompt the S_(1) system to decay to the S_(0) state,which could be to certain extent suppressed by locking the rotation of the C5-C8-C9-O10 dihedral angle.These mechanistic insights are not only helpful for understanding ESDPT but also useful for designing novel molecular materials with excellent photoluminescent performances.
