1
|
Guo X, Li G, Shi Z, Wang L. Surface Hopping with Reliable Wave Function by Introducing Auxiliary Wave Packets to Trajectory Branching. J Phys Chem Lett 2024:3345-3353. [PMID: 38498301 DOI: 10.1021/acs.jpclett.4c00437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
It is well-known that the widely utilized fewest switches surface hopping method suffers from the severe overcoherence problem, and thus adiabatic populations calculated by wave functions are generally inferior to those based on active states. More importantly, to achieve a complete description of nonadiabatic dynamics, the density matrix is essential. In this paper, we present an auxiliary branching corrected surface hopping (A-BCSH) method that introduces auxiliary wave packets (WPs) on the adiabatic potential energy surfaces for trajectory branching. Both rapid and gradual separation of WP components on different surfaces are characterized, and thus the correct decoherence time along each trajectory is captured. As demonstrated in the three standard Tully models, A-BCSH exhibits excellent internal consistency. Namely, close adiabatic populations are obtained based on both wave functions and active states. In particular, A-BCSH successfully obtains a reliable time-dependent spatial distribution of the density matrix, which relies only on electronic wave functions. Due to its high performance, our A-BCSH method provides a new and highly promising perspective on further development of more consistent surface hopping with reliable wave function.
Collapse
Affiliation(s)
- Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Guijie Li
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
2
|
Salvadori G, Mazzeo P, Accomasso D, Cupellini L, Mennucci B. Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response. J Mol Biol 2024; 436:168358. [PMID: 37944793 DOI: 10.1016/j.jmb.2023.168358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/28/2023] [Accepted: 11/02/2023] [Indexed: 11/12/2023]
Abstract
In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.
Collapse
Affiliation(s)
- Giacomo Salvadori
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Patrizia Mazzeo
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Davide Accomasso
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| |
Collapse
|
3
|
Liu Z, Hu H, Sun X. Multistate Reaction Coordinate Model for Charge and Energy Transfer Dynamics in the Condensed Phase. J Chem Theory Comput 2023; 19:7151-7170. [PMID: 37815937 PMCID: PMC10601487 DOI: 10.1021/acs.jctc.3c00770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Indexed: 10/12/2023]
Abstract
Constructing multistate model Hamiltonians from all-atom electronic structure calculations and molecular dynamics simulations is crucial for understanding charge and energy transfer dynamics in complex condensed phases. The most popular two-level system model is the spin-boson Hamiltonian, where the nuclear degrees of freedom are represented as shifted normal modes. Recently, we proposed the general multistate nontrivial extension of the spin-boson model, i.e., the multistate harmonic (MSH) model, which is constructed by extending the spatial dimensions of each nuclear mode so as to satisfy the all-atom reorganization energy restrictions for all pairs of electronic states. In this work, we propose the multistate reaction coordinate (MRC) model with a primary reaction coordinate and secondary bath modes as in the Caldeira-Leggett form but in extended spatial dimensions. The MRC model is proven to be equivalent to the MSH model and offers an intuitive physical picture of the nuclear-electronic feedback in nonadiabatic processes such as the inherent trajectory of the reaction coordinate. The reaction coordinate is represented in extended dimensions, carrying the entire reorganization energies and bilinearly coupled to the secondary bath modes. We demonstrate the MRC model construction for photoinduced charge transfer in an organic photovoltaic caroteniod-porphyrin-C60 molecular triad dissolved in tetrahydrofuran as well as excitation energy transfer in a photosynthetic light-harvesting Fenna-Matthews-Olson complex. The MRC model provides an effective and robust platform for investigating quantum dissipative dynamics in complex condensed-phase systems since it allows a consistent description of realistic spectral density, state-dependent system-bath couplings, and heterogeneous environments due to static disorder in reorganization energies.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
| | - Haorui Hu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
- Shanghai
Frontiers Science Center of Artificial Intelligence and Deep Learning, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| |
Collapse
|
4
|
Shao C, Shi Z, Xu J, Wang L. Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics. J Phys Chem Lett 2023; 14:7680-7689. [PMID: 37606199 DOI: 10.1021/acs.jpclett.3c02019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
5
|
Zhu HY, Li QS. Theoretical Understanding on the Facilitated Photoisomerization of a Carbonyl Supported Borane System. Chemphyschem 2023:e202300435. [PMID: 37646234 DOI: 10.1002/cphc.202300435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 08/29/2023] [Accepted: 08/30/2023] [Indexed: 09/01/2023]
Abstract
Boron compound BOMes2 containing an internal B-O bond undergoes highly efficient photoisomerization, followed by sequential structural transformations, resulting in a rare eight-membered B, O-heterocycle (S. Wang, et al. Org. Lett. 2019, 21, 5285-5289). In this work, the detailed reaction mechanisms of such a unique carbonyl-supported tetracoordinate boron system in the first excited singlet (S1 ) state and the ground (S0 ) state were investigated by using the complete active space self-consistent field and its second-order perturbation (MS-CASPT2//CASSCF) method combined with time-dependent density functional theory (TD-DFT). Moreover, an imine-substituted tetracoordinated organic boron system (BNMes2 ) was selected for comparative study to explore the intrinsic reasons for the difference in reactivity between the two types of compounds. Steric factor was found to influence the photoisomerization activity of BNMes2 and BOMes2 . These results rationalize the experimental observations and can provide helpful insights into understanding the excited-state dynamics of heteroatom-doped tetracoordinate organoboron compounds, which facilitates the rational design of boron-based materials with superior photoresponsive performances.
Collapse
Affiliation(s)
- Hong-Yang Zhu
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, 5 South Zhongguancun Street, Beijing, P. R. China
| | - Quan-Song Li
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, 5 South Zhongguancun Street, Beijing, P. R. China
| |
Collapse
|
6
|
Zhao X, Li J, Luo J, Liu J. Significant Acceleration of E-Z Photoisomerization induced by Molecular Planarity Breaking. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2023]
|
7
|
Chang XP, Zhao G, Zhang TS, Xie BB. Quantum mechanics/molecular mechanics studies on mechanistic photophysics of cytosine aza-analogues: 2,4-diamino-1,3,5-triazine and 2-amino-1,3,5-triazine in aqueous solution. Phys Chem Chem Phys 2023; 25:7669-7680. [PMID: 36857660 DOI: 10.1039/d2cp05639a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
The excited-state properties and photophysics of cytosine aza-analogues, i.e., 2,4-diamino-1,3,5-triazine (2,4-DT) and 2-amino-1,3,5-triazine (2-AT) in solution have been systematically explored using the QM(MS-CASPT2//CASSCF)/MM approach. The excited-state nonradiative relaxation mechanisms for the initially photoexcited S1(ππ*) state decay back to the S0 state are proposed in terms of the present computed minima, surface crossings (conical intersections and singlet-triplet crossings), and excited-state decay paths in the S1, S2, T1, T2, and S0 states. Upon photoexcitation to the bright S1(ππ*) state, 2,4-DT quickly relaxes to its S1 minimum and then overcomes a small energy barrier of 5.1 kcal mol-1 to approach a S1/S0 conical intersection, where the S1 system hops to the S0 state through S1 → S0 internal conversion (IC). In addition, at the S1 minimum, the system could partially undergo intersystem crossing (ISC) to the T1 state, followed by further ISC to the S0 state via the T1/S0 crossing point. In the T1 state, an energy barrier of 7.9 kcal mol-1 will trap 2,4-DT for a while. In parallel, for 2-AT, the system first relaxes to the S1 minimum and then S1 → S0 IC or S1 → T1 → S0 ISCs take place to the S0 state by surmounting a large barrier of 15.3 kcal mol-1 or 11.9 kcal mol-1, respectively, which heavily suppress electronic transition to the S0 state. Different from 2,4-DT, upon photoexcitation in the Franck-Condon region, 2-AT can quickly evolve in an essentially barrierless manner to nearby S2/S1 conical intersection, where the S2 and T1 states can be populated. Once it hops to the S2 state, the system will overcome a relatively small barrier (6.6 kcal mol-1vs. 15.3 kcal mol-1) through IC to the S0 state. Similarly, an energy barrier of 11.9 kcal mol-1 heavily suppresses the T1 state transformation to the S0 state. The present work manifests that the amination/deamination of the triazine rings can affect some degree of different vertical and adiabatic excitation energies and nonradiative decay pathways in solution. It not only rationalizes excited-state decay dynamics of 2,4-DT and 2-AT in aqueous solution but could also provide insights into the understanding of the photophysics of aza-nucleobases.
Collapse
Affiliation(s)
- Xue-Ping Chang
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, China.
| | - Geng Zhao
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, China.
| | - Teng-Shuo Zhang
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou 310014, China
| | - Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou 311231, China
| |
Collapse
|
8
|
Chang XP, Yu L, Zhang TS, Cui G. Quantum mechanics/molecular mechanics studies on the mechanistic photophysics of sunscreen oxybenzone in methanol solution. Phys Chem Chem Phys 2022; 24:13293-13304. [PMID: 35607908 DOI: 10.1039/d2cp01263d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Herein, we have employed the QM(CASPT2//CASSCF)/MM method to explore the photophysical and photochemical mechanism of oxybenzone (OB) in methanol solution. Based on the optimized minima, conical intersections and crossing points, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decay paths in the 1ππ*, 1nπ*, 3ππ*, 3nπ*, and S0 states, we have identified several feasible excited-state relaxation pathways for the initially populated S2(1ππ*) state to decay to the initial enol isomer' S0 state. The major one is the singlet-mediated and stretch-torsion coupled ESIPT pathway, in which the system first undergoes an essentially barrierless 1ππ* ESIPT process to generate the 1ππ* keto species, and finally realizes its ground state recovery through the subsequent carbonyl stretch-torsion facilitating S1 → S0 internal conversion (IC) and the reverse ground-state intramolecular proton transfer (GSIPT) process. The minor ones are related to intersystem crossing (ISC) processes. At the S2(1ππ*) minimum, an S2(1ππ*)/S1(1nπ*)/T2(3nπ*) three-state intersection region helps the S2 system branch into the T1 state through a S2 → S1 → T1 or S2 → T2 → T1 process. Once it has reached the T1 state, the system may relax to the S0 state via direct ISC or via subsequent nearly barrierless 3ππ* ESIPT to yield the T1 keto tautomer and ISC. The resultant S0 keto species significantly undergoes reverse GSIPT and only a small fraction yields the trans-keto form that relaxes back more slowly. However, due to small spin-orbit couplings at T1/S0 crossing points, the ISC to S0 state occurs very slowly. The present work rationalizes not only the ultrafast excited-state decay dynamics of OB but also its phosphorescence emission at low temperature.
Collapse
Affiliation(s)
- Xue-Ping Chang
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, China.
| | - Li Yu
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, China.
| | - Teng-Shuo Zhang
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou 310014, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
9
|
Chang XP, Zhang TS, Cui G. Theoretical Studies on the Excited-State Decay Mechanism of Homomenthyl Salicylate in a Gas Phase and an Acetonitrile Solution. J Phys Chem A 2021; 126:16-28. [PMID: 34963284 DOI: 10.1021/acs.jpca.1c07108] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Here, we employ the CASPT2//CASSCF and QM(CASPT2//CASSCF)/MM approaches to explore the photochemical mechanism of homomenthyl salicylate (HMS) in vacuum and an acetonitrile solution. The results show that in both cases, the excited-state relaxation mainly involves a spectroscopically "bright" S1(1ππ*) state and the lower-lying T1 and T2 states. In the major relaxation pathway, the photoexcited S1 keto system first undergoes an essentially barrierless excited-state intramolecular proton transfer (ESIPT) to generate the S1 enol minimum, near which a favorable S1/S0 conical intersection decays the system to the S0 state followed by a reverse ground-state intramolecular proton transfer (GSIPT) to repopulate the initial S0 keto species. In the minor one, an S1/T2/T1 three-state intersection in the keto region makes the T1 state populated via direct and T2-mediated intersystem crossing (ISC) processes. In the T1 state, an ESIPT occurs, which is followed by ISC near a T1/S0 crossing point in the enol region to the S0 state and finally back to the S0 keto species. In addition, a T1/S0 crossing point near the T1 keto minimum can also help the system decay to the S0 keto species. However, small spin-orbit couplings between T1 and S0 at these T1/S0 crossing points make ISC to the S0 state very slow and make the system trapped in the T1 state for a while. The present work rationalizes not only the ultrafast excited-state decay dynamics of HMS but also its low quantum yield of phosphorescence at 77 K.
Collapse
Affiliation(s)
- Xue-Ping Chang
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, P. R. China
| | - Teng-Shuo Zhang
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou 310014, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
10
|
Chang XP, Zhang TS, Fang YG, Cui G. Quantum Mechanics/Molecular Mechanics Studies on the Photophysical Mechanism of Methyl Salicylate. J Phys Chem A 2021; 125:1880-1891. [PMID: 33645980 DOI: 10.1021/acs.jpca.0c10589] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Methyl salicylate (MS) as a subunit of larger salicylates found in commercial sunscreens has been shown to exhibit keto-enol tautomerization and dual fluorescence emission via excited-state intramolecular proton transfer (ESIPT) after the absorption of ultraviolet (UV) radiation. However, its excited-state relaxation mechanism is unclear. Herein, we have employed the quantum mechanics(CASPT2//CASSCF)/molecular mechanics method to explore the ESIPT and excited-state relaxation mechanism of MS in the lowest three electronic states, that is, S0, S1, and T1 states, in a methanol solution. Based on the optimized geometric and electronic structures, conical intersections and crossing points, and minimum-energy paths combined with the computed linearly interpolated Cartesian coordinate paths, the photophysical mechanism of MS has been proposed. The S1 state is a spectroscopically bright 1ππ* state in the Franck-Condon region. From the initially populated S1 state, there exist three nonradiative relaxation paths to repopulate the S0 state. In the first one, the S1 system (i.e., ketoB form) first undergoes an ESIPT path to generate an S1 tautomer (i.e., enol form) that exhibits a large Stokes shift in experiments. The generated S1 enol tautomer further evolves toward the nearby S1/S0 conical intersection and then hops to the S0 state, followed by the backward ground-state intramolecular proton transfer (GSIPT) to the initial ketoB form S0 state. In the second one, the S1 system first hops through the S1 → T1 intersystem crossing (ISC) to the T1 state, which then further decays to the S0 state via T1 → S0 ISC at the T1/S0 crossing point. In the third path, the T1 system that stems from the S1 → T1 ISC process via the S1/T1 crossing point first takes place a T1 ESIPT to generate a T1 enol tautomer, which can further decay to the S0 state via T1-to-S0 ISC. Finally, the GSIPT occurs to back the system to the initial ketoB form S0 state. Our present work could contribute to understanding the photophysics of MS and its derivatives.
Collapse
Affiliation(s)
- Xue-Ping Chang
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, P. R. China
| | - Teng-Shuo Zhang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Ye-Guang Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
11
|
Hu Y, Xu C, Ye L, Gu FL, Zhu C. Nonadiabatic molecular dynamics simulation for the ultrafast photoisomerization of dMe-OMe-NAIP based on TDDFT on-the-fly potential energy surfaces. Phys Chem Chem Phys 2021; 23:5236-5243. [PMID: 33629668 DOI: 10.1039/d0cp06104b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Global switching on-the-fly trajectory surface hopping molecular dynamics simulation was performed on the accurate TD-B3LYP/6-31G* potential energy surfaces for E-to-Z and Z-to-E photoisomerization of dMe-OMe-NAIP up to S1(ππ*) excitation. The present TD(DFT) simulation provides accurate calculation for conical intersections between the first-excited and ground states. Thus, simulated quantum yield and lifetime of 0.23 and 620 fs (0.15 and 600 fs) for E-to-Z (Z-to-E) isomerization are in good (relatively good) agreement with experimental observation of 0.25 and 480 fs (0.24 and 430 fs), respectively. Simulated results reveal that photoisomerization pathways are initially uphill to conical intersection zones on the S1 potential energy surface and then downhill to product zones. Three types of representative conical intersections are found for determining photoisomerization mechanisms: one is the rotation type responsible for reactive isomerization and the other two are close to E and Z configurations, respectively, only for nonreactive isomerization. The present conclusions can be held in general for similar large NAIP systems of photoinduced isomerization based on E and Z configurations.
Collapse
Affiliation(s)
- Ying Hu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry & Environment of South China Normal University, Guangzhou 51006, P. R. China.
| | - Chao Xu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry & Environment of South China Normal University, Guangzhou 51006, P. R. China.
| | - Linfeng Ye
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry & Environment of South China Normal University, Guangzhou 51006, P. R. China. and Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, P. R. China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry & Environment of South China Normal University, Guangzhou 51006, P. R. China.
| | - Chaoyuan Zhu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry & Environment of South China Normal University, Guangzhou 51006, P. R. China. and Department of Applied Chemistry and Institute of Molecular Science, National Chiao-Tung University, Hsinchu 30010, Taiwan. and Department of Applied Chemistry and Center for Emergent Functional Matter Science, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan
| |
Collapse
|
12
|
Abstract
When nonadiabatic dynamics are described on the basis of trajectories, severe trajectory branching occurs when the nuclear wave packets on some potential energy surfaces are reflected while those on the remaining surfaces are not. As a result, the traditional Ehrenfest mean field (EMF) approximation breaks down. In this study, two versions of the branching corrected mean field (BCMF) method are proposed. Namely, when trajectory branching is identified, BCMF stochastically selects either the reflected or the nonreflected group to build the new mean field trajectory or splits the mean field trajectory into two new trajectories with the corresponding weights. As benchmarked in six standard model systems and an extensive model base with two hundred diverse scattering models, BCMF significantly improves the accuracy while retaining the high efficiency of the traditional EMF. In fact, BCMF closely reproduces the exact quantum dynamics in all investigated systems, thus highlighting the essential role of branching correction in nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
- Jiabo Xu
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
13
|
Xiao BY, Xu JB, Wang LJ. New energy-based decoherence correction approaches for trajectory surface hopping. CHINESE J CHEM PHYS 2020. [DOI: 10.1063/1674-0068/cjcp2006098] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Bing-yang Xiao
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jia-bo Xu
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Lin-jun Wang
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|