GW/BSE Nonadiabatic Dynamics Simulations on Excited-State Relaxation Processes of Zinc Phthalocyanine-Fullerene Dyads:Roles of Bridging Chemical Bonds

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.

Solvent Effects on Excited-State Relaxation Dynamics of Paddle-Wheel BODIPY-Hexaoxatriphenylene Conjugates:Insights from Non-adiabatic Dynamics Simulations

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.

QM/MM Study on Enzymatic Mechanism in Sinigrin Biosynthesis

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.

Structured Eigenvalue Problems in Electronic Structure Methods from a Unified Perspective

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 of 1,8-Dihydroxy-2-Naphthaldehyde:a MS-CASPT2//CASSCF Study

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.

多原子体系非绝热动力学的近似理论方法

非绝热动力学普遍存在于光物理和光化学过程中。描述非绝热跃迁需要处理电子-原子核间的相互耦合运动。由于计算量随体系尺度增大剧烈增长,准确的量子动力学计算目前只适用于描述小分子体系。为了研究多原子分子体系的非绝热过程,近年来发展了一些基于量子-经典动力学近似方法。本文将对典型的这类方法包括经典Ehrenfest方法、面跳跃方法、基于Wigner表示的混合量子-经典方法进行简要的介绍,并讨论如何将量子-经典动力学方法与电子结构从头算手段相结合,模拟非绝热过程。重点阐明各种方法的基本思想和优缺点,并对该领域的发展进行展望。