1
|
Huang KY, Li GY, Liang X, Li K, Li L, Cui G, Liu XY. "On-the-Fly" Nonadiabatic Dynamics Simulation on the Ultrafast Photoisomerization of a Molecular Photoswitch Iminothioindoxyl: An RMS-CASPT2 Investigation. J Phys Chem A 2024; 128:7145-7157. [PMID: 39145596 DOI: 10.1021/acs.jpca.4c03685] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/16/2024]
Abstract
Iminothioindoxyl (ITI) is a new class of photoswitch that exhibits many excellent properties including well-separated absorption bands in the visible region for both conformers, ultrafast Z to E photoisomerization as well as the millisecond reisomerization at room temperature for the E isomer, and switchable ability in both solids and various solvents. However, the underlying ultrafast photoisomerization mechanism at the atomic level remains unclear. In this work, we have employed a combination of high-level RMS-CASPT2-based static electronic structure calculations and nonadiabatic dynamics simulations to investigate the ultrafast photoisomerization dynamics of ITI. Based on the minimum-energy structures, minimum-energy conical intersections, linear interpolation internal coordinate paths, and nonadiabatic dynamics simulations, the overall photoisomerization scenario of ITI upon excitation is established. Upon excitation around 416 nm, the molecule will be excited to the S2 state considering its close energy to the experimentally measured absorption maximum and larger oscillator strength, from which ultrafast decay of S2 to S1 state can take place efficiently with a time constant of 62 fs. However, the photoisomerization is not likely to complete in the S2 state since the dihedral associated with the Z to E isomerization changes little during the relaxation. Upon relaxing to the S1 state, the molecule will decay to the S0 state ultrafast with a time constant of 232 fs. In contrast, the decay of the S1 state is important for the isomerization considering that the dihedral related to the isomerization of the hopping structures is close to 90°. Therefore, the S1/S0 intersection region should be important for the isomerization of ITI. Arriving at the S0 state, the molecule can either go back to the original Z reactant or isomerize to the E products. At the end of the 500 fs simulation time, the E configuration accounts for nearly 37% of the final structures. Moreover, the photoisomerization mechanism is different from the isomerization mechanism in the ground state; i.e., instead of the inversion mechanism in the ground state, the photoisomerization prefers the rotation mechanism. Our results not only agree well with previous experimental studies but also provide some novel insights that could be helpful for future improvements in the performance of the ITI photoswitches.
Collapse
Affiliation(s)
- Kai-Yue Huang
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Gao-Yi Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Xiaoqin Liang
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Kai Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Laicai Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| |
Collapse
|
2
|
Zanchet A, Roncero O, Karabulut E, Solem N, Romanzin C, Thissen R, Alcaraz C. The role of intersystem crossing in the reactive collision of S+(4S) with H2. J Chem Phys 2024; 161:044302. [PMID: 39037135 DOI: 10.1063/5.0214447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2024] [Accepted: 07/05/2024] [Indexed: 07/23/2024] Open
Abstract
We report a study on the reactive collision of S+(4S) with H2, HD, and D2 combining guided ion beam experiments and quantum-mechanical calculations. It is found that the reactive cross sections reflect the existence of two different mechanisms, one being spin-forbidden. Using different models, we demonstrate that the spin-forbidden pathway follows a complex mechanism involving three electronic states instead of two as previously thought. The good agreement between theory and experiment validates the methodology employed and allows us to fully understand the reaction mechanism. This study also provides new fundamental insights into the intersystem crossing process.
Collapse
Affiliation(s)
- Alexandre Zanchet
- Instituto de Física Fundamental, CSIC, Serrano 123, 28006 Madrid, Spain
| | - Octavio Roncero
- Instituto de Física Fundamental, CSIC, Serrano 123, 28006 Madrid, Spain
| | - Ezman Karabulut
- Vocational School of Health Services, Bitlis Eren University, 13000 Bitlis, Turkey
| | - Nicolas Solem
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France and Synchrotron SOLEIL, L'Orme des Merisiers, 91192 Saint Aubin, Gif-sur-Yvette, France
| | - Claire Romanzin
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France and Synchrotron SOLEIL, L'Orme des Merisiers, 91192 Saint Aubin, Gif-sur-Yvette, France
| | - Roland Thissen
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France and Synchrotron SOLEIL, L'Orme des Merisiers, 91192 Saint Aubin, Gif-sur-Yvette, France
| | - Christian Alcaraz
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France and Synchrotron SOLEIL, L'Orme des Merisiers, 91192 Saint Aubin, Gif-sur-Yvette, France
| |
Collapse
|
3
|
Huang H, Peng J, Zhang Y, Gu FL, Lan Z, Xu C. The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications. J Chem Phys 2024; 160:234101. [PMID: 38884395 DOI: 10.1063/5.0215036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/15/2024] [Indexed: 06/18/2024] Open
Abstract
Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.
Collapse
Affiliation(s)
- Haiyi Huang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yulin Zhang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
4
|
Shu Y, Truhlar DG. Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients. J Chem Theory Comput 2024; 20:4396-4426. [PMID: 38819014 DOI: 10.1021/acs.jctc.4c00424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
5
|
He X, Cheng X, Wu B, Liu J. Nonadiabatic Field with Triangle Window Functions on Quantum Phase Space. J Phys Chem Lett 2024; 15:5452-5466. [PMID: 38747729 PMCID: PMC11129318 DOI: 10.1021/acs.jpclett.4c00793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2024] [Revised: 04/29/2024] [Accepted: 05/01/2024] [Indexed: 05/24/2024]
Abstract
Recent progress on the constraint coordinate-momentum phase space (CPS) formulation of finite-state quantum systems has revealed that the triangle window function approach is an isomorphic representation of the exact population-population correlation function of the two-state system. We use the triangle window (TW) function and the CPS mapping kernel element to formulate a novel useful representation of discrete electronic degrees of freedom (DOFs). When it is employed with nonadiabatic field (NaF) dynamics, a new variant of the NaF approach (i.e., NaF-TW) is proposed. The NaF-TW expression of the population of any adiabatic state is always positive semidefinite. Extensive benchmark tests of model systems in both the condensed phase and gas phase demonstrate that the NaF-TW approach is able to faithfully capture the dynamical interplay between electronic and nuclear DOFs in a broad region, including where the states remain coupled all the time, as well as where the bifurcation characteristic of nuclear motion is important.
Collapse
Affiliation(s)
- Xin He
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Xiangsong Cheng
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Baihua Wu
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
6
|
Liu XY, Wang SR, Fang WH, Cui G. Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation. J Chem Theory Comput 2024; 20:3426-3439. [PMID: 38656202 DOI: 10.1021/acs.jctc.4c00068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Sheng-Rui Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
7
|
Zhang Y, Jia MR, Liu XY, Fang WH, Cui G. Photoinduced Dynamics of a Single-Walled Carbon Nanotube with a sp 3 Defect: The Importance of Excitonic Effects. J Phys Chem A 2024; 128:3311-3320. [PMID: 38654690 DOI: 10.1021/acs.jpca.4c00803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Herein, we employed linear-response time-dependent functional theory nonadiabatic dynamic simulations to explore the photoinduced exciton dynamics of a chiral single-walled carbon nanotube CNT(6,5) covalently doped with a 4-nitrobenzyl group (CNT65-NO2). The results indicate that the introduction of a sp3 defect leads to the splitting of the degenerate VBM/VBM-1 and CBM/CBM+1 states. Both the VBM upshift and the CBM downshift are responsible for the experimentally observed redshifted E11* trapping state. The simulations reveal that the photoinduced exciton relaxation dynamics completes within 500 fs, which is consistent with the experimental work. On the other hand, we also conducted the nonadiabatic carrier (electron and hole) dynamic simulations, which completely ignore the excitonic effects. The comparison demonstrates that excitonic effects are indispensable. Deep analyses show that such effects induce several dark states, which play an important role in regulating the photoinduced dynamics of CNT65-NO2. The present work demonstrates the importance of including excitonic effects in simulating photoinduced processes of carbon nanotubes. In addition, it not only rationalizes previous experiments but also provides valuable insights that will help in the future rational design of novel covalently doped carbon nanotubes with superior photoluminescent properties.
Collapse
Affiliation(s)
- Yang Zhang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Meng-Ru Jia
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
8
|
Jaiswal VK, Montorsi F, Aleotti F, Segatta F, Keefer D, Mukamel S, Nenov A, Conti I, Garavelli M. Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level. J Chem Phys 2024; 160:164316. [PMID: 38686819 DOI: 10.1063/5.0203624] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 04/10/2024] [Indexed: 05/02/2024] Open
Abstract
We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.
Collapse
Affiliation(s)
- V K Jaiswal
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - F Montorsi
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - F Aleotti
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - F Segatta
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Daniel Keefer
- Molecular Spectroscopy Department, Max Planck Institute for Polymer Research, Mainz 55128, Germany
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California Irvine, Irvine, California 92697, USA
| | - A Nenov
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - I Conti
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - M Garavelli
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| |
Collapse
|
9
|
Wada S, Tsutsumi T, Saita K, Taketsugu T. Ab initio molecular dynamics study of intersystem crossing dynamics for MH 2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces. J Comput Chem 2024; 45:552-562. [PMID: 38009451 DOI: 10.1002/jcc.27271] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 10/18/2023] [Accepted: 11/14/2023] [Indexed: 11/28/2023]
Abstract
Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.
Collapse
Affiliation(s)
- Satoi Wada
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo, Japan
| | - Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
| | - Kenichiro Saita
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Japan
| |
Collapse
|
10
|
Bian X, Subotnik JE. Spin-Dependent Stereochemistry: A Nonadiabatic Quantum Dynamics Case Study of S + H 2 → SH + H Reaction. J Phys Chem Lett 2024:3434-3440. [PMID: 38507592 DOI: 10.1021/acs.jpclett.3c03344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/22/2024]
Abstract
We study the spin-dependent stereodynamics of the S + H2 → SH + H reaction by using full-dimensional quantum dynamics calculations with zero total nuclear angular momentum along the triplet 3A″ states and singlet 1A' states. We find that the interplay between the electronic spin direction and the molecular geometry has a measurable influence on the singlet-triplet intersystem crossing reaction probabilities. Our results show that for some incident scattering angles in the body-fixed frame, the relative difference in intersystem crossing reaction probabilities (as determined between spin up and spin down initial states) can be as large as 15%. Our findings are an ab initio demonstration of spin-dependent nonadiabatic dynamics, which we hope will shine light as far as understanding the chiral-induced spin selectivity effect.
Collapse
Affiliation(s)
- Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
11
|
Chang XP, Wang JL, Peng LY, Cen XJ, Yin BW, Xie BB. Mechanistic photophysics of tellurium-substituted cytosine: Electronic structure calculations and nonadiabatic dynamics simulations. Photochem Photobiol 2024; 100:339-354. [PMID: 37435854 DOI: 10.1111/php.13835] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Revised: 06/24/2023] [Accepted: 06/26/2023] [Indexed: 07/13/2023]
Abstract
Previously, the MS-CASPT2 method was performed to study the static and qualitative photophysics of tellurium-substituted cytosine (TeC). To get quantitative information, we used our recently developed QTMF-FSSH dynamics method to simulate the excited-state decay of TeC. The CASSCF method was adopted to reduce the calculation costs, which was confirmed to provide reliable structures and energies as those of MS-CASPT2. A detailed structural analysis showed that only 5% trajectories will hop to the lower triplet or singlet state via the twisted (S2 /S1 /T2 )T intersection, while 67% trajectories will choose the planar intersections of (S2 /S1 /T3 /T2 /T1 )P and (S2 /S1 /T2 /T1 )P but subsequently become twisted in other electronic states. By contrast, ~28% trajectories will maintain in a plane throughout dynamics. Electronic population revealed that the S2 population will ultrafast transfer to the lower triplet or singlet state. Later, the TeC system will populate in the spin-mixed electronic states composed of S1 , T1 and T2 . At the end of 300 fs, most trajectories (~74%) will decay to the ground state and only 17.4% will survive in the triplet states. Our dynamics simulation verified that tellurium substitution will enhance the intersystem crossings, but the very short triplet lifetime (ca. 125 fs) will make TeC a less effective photosensitizer.
Collapse
Affiliation(s)
- Xue-Ping Chang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou, China
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang, China
| | - Jie-Lei Wang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou, China
| | - Ling-Ya Peng
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, China
| | - Xu-Jiang Cen
- Ningbo Zhongtian Engineering Co., Ltd., Ningbo, China
| | - Bo-Wen Yin
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou, China
| | - Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou, China
| |
Collapse
|
12
|
Wu B, He X, Liu J. Nonadiabatic Field on Quantum Phase Space: A Century after Ehrenfest. J Phys Chem Lett 2024; 15:644-658. [PMID: 38205956 DOI: 10.1021/acs.jpclett.3c03385] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Nonadiabatic transition dynamics lies at the core of many electron/hole transfer, photoactivated, and vacuum field-coupled processes. About a century after Ehrenfest proposed "Phasenraum" and the Ehrenfest theorem, we report a conceptually novel trajectory-based nonadiabatic dynamics approach, nonadiabatic field (NAF), based on a generalized exact coordinate-momentum phase space formulation of quantum mechanics. It does not employ the conventional Born-Oppenheimer or Ehrenfest trajectory in the nonadiabatic coupling region. Instead, in NAF the equations of motion of the independent trajectory involve a nonadiabatic nuclear force term in addition to an adiabatic nuclear force term of a single electronic state. A few benchmark tests for gas phase and condensed phase systems indicate that NAF offers a practical tool to capture the correct correlation of electronic and nuclear dynamics for processes where the states remain coupled all the time as well as for the asymptotic region where the coupling of electronic states vanishes.
Collapse
Affiliation(s)
- Baihua Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Xin He
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
13
|
Doležel J, Poryvai A, Slanina T, Filgas J, Slavíček P. Spin-Vibronic Coupling Controls the Intersystem Crossing of Iodine-Substituted BODIPY Triplet Chromophores. Chemistry 2024; 30:e202303154. [PMID: 37905588 DOI: 10.1002/chem.202303154] [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: 09/27/2023] [Revised: 10/30/2023] [Accepted: 10/30/2023] [Indexed: 11/02/2023]
Abstract
4,4-Difluoro-4-borata-3a-azonia-4a-aza-s-indacene (BODIPY) dyes are extensively used in various applications of their triplet states, ranging from photoredox catalysis, through triplet sensitization to photodynamic therapy. However, the rational design of BODIPY triplet chromophores by ab initio modelling is limited by their strong interactions of spin, electronic and vibrational dynamics. In particular, spin-vibronic coupling is often overlooked when estimating intersystem crossing (ISC) rates. In this study, a combined experimental and theoretical approach using spin-vibronic coupling to correctly describe ISC in BODIPY dyes was developed. For this purpose, seven π-extended BODIPY derivatives with iodine atoms in different positions were examined. It was found that the heavy-atom effect of iodine atoms is site specific, causing high triplet yields in only some positions. This site-specific ISC was explained by El-Sayed rules, so both the contribution and character of the molecular orbitals involved in the excitation must be considered when predicting the ISC rates. Overall, the rational design of BODIPY triplet chromophores requires using (i) the high-quality electronic structure theory, including both static and dynamical correlations; and (ii) the two-component wave function Hamiltonian, and rationalizing; and (iii) ISC based on the character of the molecular orbitals of heavy atoms involved in the excitation, expanding El-Sayed rules beyond their traditional applications.
Collapse
Affiliation(s)
- Jiří Doležel
- Institute of Organic Chemistry and Biochemistry of the Czech Academy, Flemingovo nám. 542/2, Prague 6, 160 00, Czech Republic
| | - Anna Poryvai
- Institute of Organic Chemistry and Biochemistry of the Czech Academy, Flemingovo nám. 542/2, Prague 6, 160 00, Czech Republic
| | - Tomáš Slanina
- Institute of Organic Chemistry and Biochemistry of the Czech Academy, Flemingovo nám. 542/2, Prague 6, 160 00, Czech Republic
| | - Josef Filgas
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28, Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28, Prague 6, Czech Republic
| |
Collapse
|
14
|
Wu Y, Rawlinson J, Littlejohn RG, Subotnik JE. Linear and angular momentum conservation in surface hopping methods. J Chem Phys 2024; 160:024119. [PMID: 38205852 DOI: 10.1063/5.0179599] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2023] [Accepted: 11/29/2023] [Indexed: 01/12/2024] Open
Abstract
We demonstrate that, for systems with spin-orbit coupling and an odd number of electrons, the standard fewest switches surface hopping algorithm does not conserve the total linear or angular momentum. This lack of conservation arises not so much from the hopping direction (which is easily adjusted) but more generally from propagating adiabatic dynamics along surfaces that are not time reversible. We show that one solution to this problem is to run along eigenvalues of phase-space electronic Hamiltonians H(R, P) (i.e., electronic Hamiltonians that depend on both nuclear position and momentum) with an electronic-nuclear coupling Γ · P [see Eq. (25)], and we delineate the conditions that must be satisfied by the operator Γ. The present results should be extremely useful as far as developing new semiclassical approaches that can treat systems where the nuclear, electronic orbital, and electronic spin degrees of freedom altogether are all coupled together, hopefully including systems displaying the chiral-induced spin selectivity effect.
Collapse
Affiliation(s)
- Yanze Wu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Jonathan Rawlinson
- Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
15
|
Wenzel M, Mitric R. Prediction of fluorescence quantum yields using the extended thawed Gaussian approximation. J Chem Phys 2023; 159:234113. [PMID: 38108487 DOI: 10.1063/5.0178106] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 11/26/2023] [Indexed: 12/19/2023] Open
Abstract
Spontaneous emission and internal conversion rates are calculated within harmonic approximations and compared to the results obtained within the semi-classical extended thawed Gaussian approximation (ETGA). This is the first application of the ETGA in the calculation of internal conversion and emission rates for real molecular systems, namely, formaldehyde, fluorobenzene, azulene, and a dicyano-squaraine dye. The viability of the models as black-box tools for prediction of spontaneous emission and internal conversion rates is assessed. All calculations were done using a consistent protocol in order to investigate how different methods perform without previous experimental knowledge using density functional theory (DFT) and time-dependent DFT (TD-DFT) with B3LYP, PBE0, ωB97XD, and CAM-B3LYP functionals. Contrasting the results with experimental data shows that there are further improvements required before theoretical predictions of emission and internal conversion rates can be used as reliable indicators for the photo-luminescence properties of molecules. We find that the ETGA performs rather similar to the vertical harmonical model. Including anharmonicities in the calculation of internal conversion rates has a moderate effect on the quantitative results in the studied systems. The emission rates are fairly stable with respect to computational parameters, but the internal conversion rate reveals itself to be highly dependent on the choice of the spectral line shape function, particularly the width of the Lorentzian function, associated with homogeneous broadening.
Collapse
Affiliation(s)
- Michael Wenzel
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| | - Roland Mitric
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| |
Collapse
|
16
|
Liu XY, Chen WK, Fang WH, Cui G. Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications. J Chem Theory Comput 2023. [PMID: 37984502 DOI: 10.1021/acs.jctc.3c00960] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2023]
Abstract
Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, P. R. China
| | - Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Hefei National Laboratory, Hefei 230088, 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
- Hefei National Laboratory, Hefei 230088, P. R. China
| |
Collapse
|
17
|
Liu S, Liu SS, Tang XM, Liu XY, Yang JJ, Cui G, Li L. Solvent effects on the photoinduced charge separation dynamics of directly linked zinc phthalocyanine-perylenediimide dyads: a nonadiabatic dynamics simulation with an optimally tuned screened range-separated hybrid functional. Phys Chem Chem Phys 2023; 25:28452-28464. [PMID: 37846460 DOI: 10.1039/d3cp03517d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
Herein, we have employed a combination of the optimally tuned screened range-separated hybrid (OT-SRSH) functional, the polarizable continuum model (PCM), and nonadiabatic dynamics (NAMD) simulations to investigate the photoinduced dynamics of directly linked donor-acceptor dyads formed using zinc phthalocyanine (ZnPc) and perylenediimide (PDI), in which ZnPc is the donor while PDI is the acceptor. Our simulations aim to analyze the behavior of these dyads upon local excitation of the ZnPc moiety in the gas phase and in benzonitrile. Our findings indicate that the presence of a solvent can significantly influence the excited state dynamics of ZnPc-PDI dyads. Specifically, the polar solvent benzonitrile effectively lowers the vertical excitation energies of the charge transfer (CT) state from ZnPc to PDI. As a result, the energetic order of the locally excited (LE) states of ZnPc and the CT states is reversed compared to the gas phase. Consequently, the photoinduced electron transfer (PET) dynamics from ZnPc to PDI, which is absent in the gas phase, takes place in benzonitrile with a time constant of 10.4 ps. Importantly, our present work not only qualitatively agrees with experimental results but also provides in-depth insights into the underlying mechanisms responsible for the photoinduced dynamics of ZnPc-PDI. Moreover, this study emphasizes the importance of appropriately considering solvent effects in NAMD simulation of organic donor-acceptor systems, taking into account the distinct excited state dynamics observed in the gas phase and benzonitrile. Furthermore, the combination of the OT-SRSH functional, the PCM solvent model, and nonadiabatic dynamics simulations shows promise as a strategy for investigating the complex excited state dynamics of organic donor-acceptor systems in solvents. These findings will be valuable for the future design of novel organic donor-acceptor structures with improved performance.
Collapse
Affiliation(s)
- Shuai Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China.
| | - Sha-Sha Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China.
| | - Xiao-Mei Tang
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China.
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China.
| | - Jia-Jia Yang
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China.
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Hefei National Laboratory, Hefei 230088, China
| | - Laicai Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China.
| |
Collapse
|
18
|
Polonius S, Zhuravel O, Bachmair B, Mai S. LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution. J Chem Theory Comput 2023; 19:7171-7186. [PMID: 37788824 PMCID: PMC10601485 DOI: 10.1021/acs.jctc.3c00805] [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/25/2023] [Indexed: 10/05/2023]
Abstract
We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - Oleksandra Zhuravel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
19
|
Wang Y, Gu Z, Cao L, Zhang B, Zhang S. Channel Competition and Control of Relaxation Pathways in S 1 State of Acrolein: Role of Conical Intersection and Surface Crossing. J Phys Chem A 2023; 127:8595-8601. [PMID: 37801298 DOI: 10.1021/acs.jpca.3c05376] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
Channel competition and further photochemical control of relaxation pathways in excited molecules are of primary importance in photochemistry and related areas. Acrolein, as the simplest and most typical α,β-enone, is suitable to provide a model for understanding the photochemistry and photophysics of α,β-enones. Here, the ultrafast dynamics in acrolein following S1(nπ*) excitation has been studied by time-resolved photoelectron imaging (TRPEI) and mass spectroscopy. The competition between intersystem crossing (ISC) and internal conversion (IC) is investigated. The key factor influencing the decay pathways and the relative contributions are revealed to be the position of the excitation relative to the energy of the S1/S0 conical intersection (CI), which is obtained to be 3.65-3.76 eV experimentally. If the excitation is above the CI, IC is superior to ISC and most excited molecules go back to the ground. Otherwise, ISC will dominate the relaxation and lead the triplet products formation. These results show the potential of affecting the dynamics and governing the fate of excited molecules by adjusting the excitation conditions from the point of view of chemical control.
Collapse
Affiliation(s)
- Yanmei Wang
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, P. R. China
| | - Zhenfei Gu
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, P. R. China
| | - Ling Cao
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, P. R. China
| | - Bing Zhang
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, P. R. China
| | - Song Zhang
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, P. R. China
| |
Collapse
|
20
|
Danilov D, Jenkins AJ, Bearpark MJ, Worth GA, Robb MA. Coherent Mixing of Singlet and Triplet States in Acrolein and Ketene: A Computational Strategy for Simulating the Electron-Nuclear Dynamics of Intersystem Crossing. J Phys Chem Lett 2023:6127-6134. [PMID: 37364275 DOI: 10.1021/acs.jpclett.3c01187] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023]
Abstract
We present a theoretical study of intersystem crossing (ISC) in acrolein and ketene with the Ehrenfest method that can describe a superposition of singlet and triplet states. Our simulations illustrate a new mechanistic effect of ISC, namely, that a superposition of singlets and triplets yields nonadiabatic dynamics characteristic of that superposition rather than the constituent state potential energy surfaces. This effect is particularly significant in ketene, where mixing of singlet and triplet states along the approach to a singlet/singlet conical intersection occurs, with the spin-orbit coupling (SOC) remaining small throughout. In both cases, the effects require many recrossings of the singlet/triplet state crossing seam, consistent with the textbook treatment of ISC.
Collapse
Affiliation(s)
- Don Danilov
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 82 Wood Lane, W12 0BZ London, United Kingdom
| | - Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J Bearpark
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 82 Wood Lane, W12 0BZ London, United Kingdom
| | - Graham A Worth
- Department of Chemistry, University College London, 20 Gordon St., WC1H 0AJ London, United Kingdom
| | - Michael A Robb
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 82 Wood Lane, W12 0BZ London, United Kingdom
| |
Collapse
|
21
|
Chakraborty P, Couto RC, List NH. Deciphering Methylation Effects on S 2( ππ*) Internal Conversion in the Simplest Linear α,β-Unsaturated Carbonyl. J Phys Chem A 2023. [PMID: 37331016 DOI: 10.1021/acs.jpca.3c02582] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2023]
Abstract
Chemical substituents can influence photodynamics by altering the location of critical points and the topography of the potential energy surfaces (electronic effect) and by selectively modifying the inertia of specific nuclear modes (inertial effects). Using nonadiabatic dynamics simulations, we investigate the impact of methylation on S2(ππ*) internal conversion in acrolein, the simplest linear α,β-unsaturated carbonyl. Consistent with time constants reported in a previous time-resolved photoelectron spectroscopy study, S2 → S1 deactivation occurs on an ultrafast time scale (∼50 fs). However, our simulations do not corroborate the sequential decay model used to fit the experiment. Instead, upon reaching the S1 state, the wavepacket bifurcates: a portion undergoes ballistic S1 → S0 deactivation (∼90 fs) mediated by fast bond-length alternation motion, while the remaining decays on the picosecond time scale. Our analysis reveals that methyl substitution, generally assumed to mainly exert inertial influence, is also manifested in important electronic effects due to its weak electron-donating ability. While methylation at the β C atom gives rise to effects principally of an inertial nature, such as retarding the twisting motion of the terminal -CHCH3 group and increasing its coupling with pyramidalization, methylation at the α or carbonyl C atom modifies the potential energy surfaces in a way that also contributes to altering the late S1-decay behavior. Specifically, our results suggest that the observed slowing of the picosecond component upon α-methylation is a consequence of a tighter surface and reduced amplitude along the central pyramidalization, effectively restricting the access to the S1/S0-intersection seam. Our work offers new insight into the S2(ππ*) internal conversion mechanisms in acrolein and its methylated derivatives and highlights site-selective methylation as a tuning knob to manipulate photochemical reactions.
Collapse
Affiliation(s)
- Pratip Chakraborty
- Department of Chemistry, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
| | - Rafael C Couto
- Department of Chemistry, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
| | - Nanna H List
- Department of Chemistry, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
| |
Collapse
|
22
|
Zheng Z, Shi Y, Zhou JJ, Prezhdo OV, Zheng Q, Zhao J. Ab initio real-time quantum dynamics of charge carriers in momentum space. NATURE COMPUTATIONAL SCIENCE 2023; 3:532-541. [PMID: 38177418 DOI: 10.1038/s43588-023-00456-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Accepted: 04/21/2023] [Indexed: 01/06/2024]
Abstract
Application of the non-adiabatic molecular dynamics (NAMD) approach is limited to studying carrier dynamics in the momentum space, as a supercell is required to sample the phonon excitation and electron-phonon (e-ph) interaction at different momenta in a molecular dynamics simulation. Here we develop an ab initio approach for the real-time charge carrier quantum dynamics in the momentum space (NAMD_k) by directly introducing e-ph coupling into the Hamiltonian based on the harmonic approximation. The NAMD_k approach maintains the zero-point energy and includes memory effects of carrier dynamics. The application of NAMD_k to the hot carrier dynamics in graphene reveals the phonon-specific relaxation mechanism. An energy threshold of 0.2 eV-defined by two optical phonon modes-separates the hot electron relaxation into fast and slow regions with lifetimes of pico- and nanoseconds, respectively. The NAMD_k approach provides an effective tool to understand real-time carrier dynamics in the momentum space for different materials.
Collapse
Affiliation(s)
- Zhenfa Zheng
- Department of Physics, ICQD/Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, China
| | - Yongliang Shi
- Department of Physics, ICQD/Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, China.
- Center for Spintonics and Quantum Systerms, State Key Laboratory for Mechanical Behavior of Materials, School of Materials Science and Engineering, Xi'an Jiaotong University, Xi'an, China.
- State Key Laboratory of Surface Physics and Department of Physics, Fudan University, Shanghai, China.
| | - Jin-Jian Zhou
- School of Physics, Beijing Institute of Technology, Beijing, China
| | - Oleg V Prezhdo
- Departments of Chemistry, Physics, and Astronomy, University of Southern California, Los Angeles, CA, USA
| | - Qijing Zheng
- Department of Physics, ICQD/Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, China.
| | - Jin Zhao
- Department of Physics, ICQD/Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, China.
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA, USA.
- Hefei National Laboratory, University of Science and Technology of China, Hefei, China.
| |
Collapse
|
23
|
Shu Y, Zhang L, Wu D, Chen X, Sun S, Truhlar DG. New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States. J Chem Theory Comput 2023; 19:2419-2429. [PMID: 37079755 DOI: 10.1021/acs.jctc.2c01173] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/22/2023]
Abstract
It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Dihua Wu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Xiye Chen
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
24
|
Dergachev ID, Dergachev VD, Rooein M, Mirzanejad A, Varganov SA. Predicting Kinetics and Dynamics of Spin-Dependent Processes. Acc Chem Res 2023; 56:856-866. [PMID: 36926853 DOI: 10.1021/acs.accounts.2c00843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
Abstract
ConspectusPredicting mechanisms and rates of nonadiabatic spin-dependent processes including photoinduced intersystem crossings, thermally activated spin-forbidden reactions, and spin crossovers in metal centers is a very active field of research. These processes play critical roles in transition-metal-based and metalloenzymatic catalysis, molecular magnets, light-harvesting materials, organic light-emitting diodes, photosensitizers for photodynamic therapy, and many other applications. Therefore, accurate modeling of spin-dependent processes in complex systems and on different time scales is important for many problems in chemistry, biochemistry, and materials sciences.Nonadiabatic statistical theory (NAST) and nonadiabatic molecular dynamics (NAMD) are two complementary approaches to modeling the kinetics and dynamics of spin-dependent processes. NAST predicts the probabilities and rate constants of nonradiative transitions between electronic states with different spin multiplicities using molecular properties at only few critical points on the potential energy surfaces (PESs), including the reactant minimum and the minimum energy crossing point (MECP) between two spin states. This makes it possible to obtain molecular properties for NAST calculations using accurate but often computationally expensive electronic structure methods, which is critical for predicting the rate constants of spin-dependent processes. Alternatively, NAST can be used to study spin-dependent processes in very large complex molecular systems using less computationally expensive electronic structure methods. The nuclear quantum effects, such as zero-point vibrational energy, tunneling, and interference between reaction paths can be easily incorporated. However, the statistical and local nature of NAST makes it more suitable for large systems and slow kinetics. In contrast, NAMD explores entire PESs of interacting electronic states, making it ideal for modeling fast barrierless spin-dependent processes. Because the knowledge of large portions of PESs is often needed, the simulations require a very large number of electronic structure calculations, which limits the NAMD applicability to relatively small molecular systems and ultrafast kinetics.In this Account, we discuss our contribution to the development of the NAST and NAMD approaches for predicting the rates and mechanism of spin-dependent processes. First, we briefly describe our NAST and NAMD implementations. The NAST implementation is an extension of the transition state theory to the processes involving two crossing potential energy surfaces of different spin multiplicities. The NAMD approach includes the trajectory surface hopping (TSH) and ab initio multiple spawning (AIMS) methods. Second, we discuss several applications of NAST and NAMD to model spin-dependent processes in different systems. The NAST applicability to large complex systems is demonstrated by the studies of the spin-forbidden isomerization of the active sites of metal-sulfur proteins. Our implementation of the MECP search algorithm within the fully ab initio fragment molecular orbital method allows applying NAST to systems with thousands of atoms, such as the solvated protein rubredoxin. Applications of NAMD to ultrafast spin-dependent processes are represented by the generalized AIMS simulations utilizing the fast GPU-based TeraChem electronic structure program to gain insight into the complex photoexcited state relaxation in 2-cyclopentenone.
Collapse
Affiliation(s)
- Ilya D Dergachev
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Vsevolod D Dergachev
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Mitra Rooein
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Amir Mirzanejad
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
25
|
Xie RF, Zhang JB, Wu Y, Li L, Liu XY, Cui G. Non-negligible roles of charge transfer excitons in ultrafast excitation energy transfer dynamics of a double-walled carbon nanotube. J Chem Phys 2023; 158:054108. [PMID: 36754819 DOI: 10.1063/5.0134353] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Herein, we employed a developed linear response time dependent density functional theory-based nonadiabatic dynamics simulation method that explicitly takes into account the excitonic effects to investigate photoinduced excitation energy transfer dynamics of a double-walled carbon nanotube (CNT) model with different excitation energies. The E11 excitation of the outer CNT will generate a local excitation (LE) |out*〉 exciton due to its low energy, which does not induce any charge separation. In contrast, the E11 excitation of the inner CNT can generate four kinds of excitons with the LE exciton |in*〉 dominates. In the 500-fs dynamics simulation, the LE exciton |in*〉 and charge transfer (CT) excitons |out-in+〉 and |out+in-〉 are all gradually converted to the |out*〉 exciton, corresponding to a photoinduced excitation energy transfer, which is consistent with experimental studies. Finally, when the excitation energy is close to the E22 state of the outer CNT (∼1.05 eV), a mixed population of different excitons, with the |out*〉 exciton dominated, is generated. Then, photoinduced energy transfer from the outer to inner CNTs occurs in the first 50 fs, which is followed by an inner to outer excitation energy transfer that is completed in 400 fs. The present work not only sheds important light on the mechanistic details of wavelength-dependent excitation energy transfer of a double-walled CNT model but also demonstrates the roles and importance of CT excitons in photoinduced excitation energy transfer. It also emphasized that explicitly including the excitonic effects in electronic structure calculations and nonadiabatic dynamics simulations is significant for correct understanding/rational design of optoelectronic properties of periodically extended systems.
Collapse
Affiliation(s)
- Rui-Fang Xie
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Jing-Bin Zhang
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Yang Wu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Laicai Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
26
|
Wenzel M, Mitric R. Internal conversion rates from the extended thawed Gaussian approximation: Theory and validation. J Chem Phys 2023; 158:034105. [PMID: 36681643 DOI: 10.1063/5.0130340] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
The theoretical prediction of the rates of nonradiative processes in molecules is fundamental in assessing their emissive properties. In this context, global harmonic models have been widely used to simulate vibronic spectra as well as internal conversion rates and to predict photoluminescence quantum yields. However, these simplified models suffer from the limitations that are inherent to the harmonic approximation and can have a severe effect on the calculated internal conversion rates. Therefore, the development of more accurate semiclassical methods is highly desirable. Here, we introduce a procedure for the calculation of nonradiative rates in the framework of the time-dependent semi-classical Extended Thawed Gaussian Approximation (ETGA). We systematically investigate the performance of the ETGA method by comparing it to the adiabatic and vertical harmonic methods, which belong to the class of widely used global harmonic models. Its performance is tested in potentials that cannot be treated adequately by global harmonic models, beginning with Morse potentials of varying anharmonicity followed by a double well potential. The calculated radiative and nonradiative internal conversion rates are compared to reference values based on exact quantum dynamics. We find that the ETGA has the capability to predict internal conversion rates in anharmonic systems with an appreciable energy gap, whereas the global harmonic models prove to be insufficient.
Collapse
Affiliation(s)
- Michael Wenzel
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| | - Roland Mitric
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| |
Collapse
|
27
|
Suda K, Yokogawa D. Spin-Orbit Coupling Calculation Combined with the Reference Interaction Site Model Self-Consistent Field Explicitly Including Constrained Spatial Electron Density Distribution. J Chem Theory Comput 2022; 18:6043-6051. [PMID: 36069633 DOI: 10.1021/acs.jctc.2c00580] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Studying the radiative and non-radiative decay processes of molecules in a solution is an important issue in the design of organic and functional molecules. Theoretical approaches have great potential for revealing this decay process through computation of various parameters, such as the energy surfaces at the excited state and spin-orbit coupling (SOC). The development of quantum chemical programs has enabled the calculation of SOC values to become popular for the gas phase. However, SOC calculations in solution have some difficulties that need to be overcome. In the present study, the authors combined the SOC calculations with the reference interaction site model self-consistent field explicitly including constrained spatial electron density distribution. To validate the reliability of our method, the decay process of dimethylaminobenzonitrile in cyclohexane and acetonitrile was studied. By computing the SOC values in both solution systems, the authors were able to investigate the decay process at the atomistic level. Furthermore, a natural transition orbital analysis and the measurement of the decomposed SOC values were found to provide a clear understanding of intersystem crossing.
Collapse
Affiliation(s)
- Kayo Suda
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
28
|
Cofer-Shabica DV, Menger MFSJ, Ou Q, Shao Y, Subotnik JE, Faraji S. INAQS, a Generic Interface for Nonadiabatic QM/MM Dynamics: Design, Implementation, and Validation for GROMACS/Q-CHEM simulations. J Chem Theory Comput 2022; 18:4601-4614. [PMID: 35901266 DOI: 10.1021/acs.jctc.2c00204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The accurate description of large molecular systems in complex environments remains an ongoing challenge for the field of computational chemistry. This problem is even more pronounced for photoinduced processes, as multiple excited electronic states and their corresponding nonadiabatic couplings must be taken into account. Multiscale approaches such as hybrid quantum mechanics/molecular mechanics (QM/MM) offer a balanced compromise between accuracy and computational burden. Here, we introduce an open-source software package (INAQS) for nonadiabatic QM/MM simulations that bridges the sampling capabilities of the GROMACS MD package and the excited-state infrastructure of the Q-CHEM electronic structure software. The interface is simple and can be adapted easily to other MD codes. The code supports a variety of different trajectory-based molecular dynamics, ranging from Born-Oppenheimer to surface hopping dynamics. To illustrate the power of this combination, we simulate electronic absorption spectra, free-energy surfaces along a reaction coordinate, and the excited-state dynamics of 1,3-cyclohexadiene in solution.
Collapse
Affiliation(s)
- D Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, 231 S. 34 Street, Cret Wing 141D, Philadelphia, Pennsylvania 19104-6243, United States
| | - Maximilian F S J Menger
- Zernike Institute for Advanced Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Qi Ou
- AI for Science Institute, Beijing 100080, China
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, 231 S. 34 Street, Cret Wing 141D, Philadelphia, Pennsylvania 19104-6243, United States
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| |
Collapse
|
29
|
Abstract
Intersystem crossing (ISC), a vital component of the electronic and nuclear transitions that compose photophysics, has been successfully simulated in light elements and transition metal complexes. Derived from the Z-dependent spin-orbit coupling (SOC), ISC is expected to be of greater importance in heavier elements, but few attempts have been made at the simulation of ISC in lanthanides or actinides. In this work, we explore several of the challenges that will need to be overcome in order to treat ISC in late-row elements, including the loss of spin as a good quantum number, the need to include SOC variationally via two- or four-component electronic structure, and the high density of states present in late-row complexes. Density functional theory (DFT) calculations are used to illustrate several of these effects, while a model Hamiltonian is used to illustrate the importance of momentum rescaling in surface hopping simulations of strongly coupled states.
Collapse
Affiliation(s)
- Andrew J S Valentine
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
30
|
Xie BB, Jia PK, Wang KX, Chen WK, Liu XY, Cui G. Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems. J Phys Chem A 2022; 126:1789-1804. [PMID: 35266391 DOI: 10.1021/acs.jpca.1c10195] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.
Collapse
Affiliation(s)
- Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Pei-Ke Jia
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Ke-Xin Wang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, Sichuan, 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
|
31
|
Chen WK, Cui G, Liu XY. Solvent effects on excited-state relaxation dynamics of paddle-wheel BODIPY-Hexaoxatriphenylene conjugates: Insights from non-adiabatic dynamics simulations. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2110214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
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 BF2-chelated dipyrromethenes) conjugates D-A complexes with the combination of both electronic structure calculations and non-adiabatic 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.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| |
Collapse
|
32
|
Mao D, Chen XR, Li DH, Liu XY, Cui G, Li L. Ultrafast charge transfer in a nonfullerene all-small-molecule organic solar cell: a nonadiabatic dynamics simulation with optimally tuned range-separated functional. Phys Chem Chem Phys 2022; 24:27173-27183. [DOI: 10.1039/d2cp03822f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The combination of nonadiabatic dynamics simulation and optimally tuned range-separated functional might be a powerful tool for elucidating the ultrafast charge transfer in nonfullerene all-small-molecule organic solar cells.
Collapse
Affiliation(s)
- Dan Mao
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, 610068, China
| | - Xin-Rui Chen
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, 610068, China
| | - Dong-Heng Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, 610068, China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, 610068, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, 100875, China
| | - Laicai Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, 610068, China
| |
Collapse
|
33
|
Computational Characterization of Nanosystems. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
34
|
Zhang S, Zeng YP, Wan XJ, Xu DH, Liu XY, Cui G, Li L. Ultrafast Exciton Delocalization and Localization Dynamics of a Perylene Bisimide Quadruple π-Stack: A Nonadiabatic Dynamics Simulation. Phys Chem Chem Phys 2022; 24:7293-7302. [DOI: 10.1039/d2cp00018k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Unraveling the photogenerated exciton dynamics of πstacked molecular aggregates is of great importance for both fundamental studies and industrial applications. Among various πstacked molecular aggregates, perylene tetracarboxylic acid bisimides (PBI)...
Collapse
|
35
|
Xie BB, Tang XF, Liu XY, Chang XP, Cui G. Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations. Phys Chem Chem Phys 2021; 23:27124-27149. [PMID: 34849517 DOI: 10.1039/d1cp03994f] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Photophysics and photochemistry are basic subjects in the study of light-matter interactions and are ubiquitous in diverse fields such as biology, energy, materials, and environment. A full understanding of mechanistic photophysics and photochemistry underpins many recent advances and applications. This contribution first provides a short discussion on the theoretical calculation methods we have used in relevant studies, then we introduce our latest progress on the mechanistic photophysics and photochemistry of two classes of molecular systems, namely unnatural bases and sunscreens. For unnatural bases, we disclose the intrinsic driving forces for the ultrafast population to reactive triplet states, impacts of the position and degree of chalcogen substitutions, and the effects of complex environments. For sunscreen molecules, we reveal the photoprotection mechanisms that dissipate excess photon energy to the surroundings by ultrafast internal conversion to the ground state. Finally, relevant theoretical challenges and outlooks are discussed.
Collapse
Affiliation(s)
- Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China.
| | - Xiu-Fang Tang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China.
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, Sichuan, 610068, China
| | - Xue-Ping Chang
- College of Chemistry and Chemical Engineering, Xinyang Normal University, Xinyang 464000, 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
|
36
|
Chen WK, Sun XW, Fang Q, Liu XY, Cui GL. GW/BSE nonadiabatic dynamics simulations on excited-state relaxation processes of zinc phthalocyanine-fullerene dyads: Roles of bridging chemical bonds. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2109162] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Wen-kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Xin-wei Sun
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Qiu Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Xiang-yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Gang-long Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
37
|
Wang Y, Yarkony DR. Conical intersection seams in spin-orbit coupled systems with an even number of electrons: A numerical study based on neural network fit surfaces. J Chem Phys 2021; 155:174115. [PMID: 34742185 DOI: 10.1063/5.0067660] [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/14/2022] Open
Abstract
In this work, we consider the existence and topography of seams of conical intersections (CIs) for two key singlet-triplet systems, including a uniformly scaled spin-orbit interaction. The basic one triplet and one singlet state system denoted as (S0,T1) and the two singlets and one triplet system denoted as (S0,S1,T1) are treated. Essential to this analysis are realistic electronic structure data taken from a recently reported neural network fit for the 1,21A and 13A states of NH3, including Hsf (spin-free) and Hso (spin-orbit) surfaces derived from high quality ab initio wavefunctions. Three types of seams for the (S0,S1,T1) system are reported, which depend on the choice of the electronic Hamiltonian, He. The nonrelativistic CI seam [He = Hsf, (S0,S1)], the energy minimized nonrelativistic singlet-triplet intersection seam [He = Hsf, (S0,T1)], and the fully relativistic seam in the spin-diabatic representation (He = Htot = Hsf + Hso) are reported as functions of R(N-H). The derivative couplings are computed using He = Htot and Hsf from the fit data. The line integral of the derivative coupling is employed to juxtapose the geometric phase in the relativistic, He = Htot, and nonrelativistic, He = Hsf, cases. It is found for the (S0,T1) system that there is no CI in the spin-adiabatic representation, while for the (S0,S1,T1) system, CI can only be formed for two pairs of spin-adiabatic electronic states. The geometric phase effect thus needs to be handled with care when it comes to spin-nonconserving dynamics simulations.
Collapse
Affiliation(s)
- Yuchen Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
38
|
Mukherjee S, Varganov SA. Intersystem crossing and internal conversion dynamics with GAIMS-TeraChem: Excited state relaxation in 2-cyclopentenone. J Chem Phys 2021; 155:174107. [PMID: 34742200 DOI: 10.1063/5.0068040] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Excited states relaxation in complex molecules often involves two types of nonradiative transitions, internal conversion (IC) and intersystem crossing (ISC). In the situations when the timescales of IC and ISC are comparable, an interplay between these two types of transitions can lead to complex nonadiabatic dynamics on multiple electronic states of different characters and spin multiplicities. We demonstrate that the generalized ab initio multiple spawning (GAIMS) method interfaced with the fast graphics processing unit-based TeraChem electronic structure code can be used to model such nonadiabatic dynamics involving both the IC and ISC transitions in molecules of moderate size. We carried out 1500 fs GAIMS simulations leading to the creation of up to 2500 trajectory basis functions to study the excited states relaxation in 2-cyclopentenone. After a vertical excitation from the ground state to the bright S2 state, the molecule quickly relaxes to the S1 state via conical intersection. The following relaxation proceeds along two competing pathways: one involves IC to the ground state, and the other is dominated by ISC to the low-lying triplet states. The time constants describing the population transfer between the six lowest singlet and triplet states predicted by the GAIMS dynamics are in good agreement with the characteristic times of IC and ISC obtained from the analysis of the time-resolved photoelectron spectrum.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA
| |
Collapse
|
39
|
Liu J, Lan Z, Yang J. An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics. Phys Chem Chem Phys 2021; 23:22313-22323. [PMID: 34591049 DOI: 10.1039/d1cp03477d] [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
We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.
Collapse
Affiliation(s)
- Jie Liu
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| |
Collapse
|
40
|
Ren JJ, Wang YH, Li WT, Jiang T, Shuai ZG. Time-dependent density matrix renormalization group coupled with n-mode representation potentials for the excited state radiationless decay rate: Formalism and application to azulene. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2108138] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Jia-jun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Yuan-heng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Wei-tang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Zhi-gang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| |
Collapse
|
41
|
Avagliano D, Bonfanti M, Garavelli M, González L. QM/MM Nonadiabatic Dynamics: the SHARC/COBRAMM Approach. J Chem Theory Comput 2021; 17:4639-4647. [PMID: 34114454 DOI: 10.1021/acs.jctc.1c00318] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the SHARC/COBRAMM approach to enable easy and efficient excited-state dynamics simulations at different levels of electronic structure theory in the presence of complex environments using a quantum mechanics/molecular mechanics (QM/MM) setup. SHARC is a trajectory surface-hoping method that can incorporate the simultaneous effects of nonadiabatic and spin-orbit couplings in the excited-state dynamics of molecular systems. COBRAMM allows ground- and excited-state QM/MM calculations using a subtractive scheme, with electrostatic embedding and a hydrogen link-atom approach. The combination of both free and open-source program packages provides a modular and extensive framework to model nonadiabatic processes after light irradiation from the atomistic scale to the nano-scale. As an example, the relaxation of acrolein from S1 to T1 in solution is provided.
Collapse
Affiliation(s)
- Davide Avagliano
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
| | - Matteo Bonfanti
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Viale Del Risorgimento, 4, I-40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Viale Del Risorgimento, 4, I-40136 Bologna, Italy
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria.,Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
| |
Collapse
|
42
|
Guan Y, Xie C, Guo H, Yarkony DR. Enabling a Unified Description of Both Internal Conversion and Intersystem Crossing in Formaldehyde: A Global Coupled Quasi-Diabatic Hamiltonian for Its S 0, S 1, and T 1 States. J Chem Theory Comput 2021; 17:4157-4168. [PMID: 34132545 DOI: 10.1021/acs.jctc.1c00370] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
In our recent work, a diabatic Hamiltonian that couples the S0 and S1 states of formaldehyde was constructed using a robust fitting-and-diabatizing procedure with artificial neural networks, which is capable of representing adiabatic energies, energy gradients, and derivative couplings over a wide range of geometries including seams of conical intersection. In this work, based on the diabatization of S0 and S1, the spin-orbit couplings between singlet states (S0, S1) and triplet state T1 are also determined in the same diabatic representation. The diabatized spin-orbit couplings are then fit with a symmetrized neural-network functional form. The ab initio spin-orbit couplings are well reproduced in large configuration space. Together with the neural-network-based potential energy surface for T1, the full quasi-diabatic Hamiltonian for the S0, S1, and T1 states is completed, enabling a unified description of both internal conversion and intersystem crossing in formaldehyde. The vibrational levels on the three adiabatic states are found to be in good agreement with known experimental band origins.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, People's Republic of China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
43
|
Wang Y, Ren J, Shuai Z. Evaluating the anharmonicity contributions to the molecular excited state internal conversion rates with finite temperature TD-DMRG. J Chem Phys 2021; 154:214109. [PMID: 34240969 DOI: 10.1063/5.0052804] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we propose a new method to calculate molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory. This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts toward the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts toward the repulsive side. In both cases, a higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from S1 to S0 on the ab initio anharmonic PES approximated by the one-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on the anharmonic PES are 30%-40% higher than the rates under HA.
Collapse
Affiliation(s)
- Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
44
|
Zhou HJ, Xu DH, Yang QH, Liu XY, Cui G, Li L. Rational design of monolayer transition metal dichalcogenide@fullerene van der Waals photovoltaic heterojunctions with time-domain density functional theory simulations. Dalton Trans 2021; 50:6725-6734. [PMID: 33912883 DOI: 10.1039/d1dt00291k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
van der Waals heterojunctions formed by transition metal dichalcogenides (TMDs) and fullerenes are promising candidates for novel photovoltaic devices due to the excellent optoelectronic properties of both TMDs and fullerenes. However, relevant experimental and theoretical investigations remain scarce to the best of our knowledge. Herein, we have first employed static density functional theory (DFT) calculations in combination with time-domain density functional theory (TDDFT) based nonadiabatic dynamics simulations to rationally evaluate the photovoltaic performances of four TMD@fullerene heterostructures, i.e. WSe2@C60, WSe2@C70, MoTe2@C60 and MoTe2@C70, respectively. Our simulation results indicate that the C70-based heterostructures overall have better photoinduced electron transfer efficiencies than their C60-based counterparts, among which the performance of the WSe2@C70 heterostructure is the best and the electron transfer from WSe2 to C70 almost accomplishes within 1 ps. In addition, the large build-in potential of about 0.75 eV of WSe2@C70 is beneficial for the charge separation processes. Our present work not only selects the van der Waals TMD@fullerene heterojunctions that might have excellent photovoltaic properties, but also paves the way for the rational design of novel heterojunctions with better optoelectronic performances with DFT and TDDFT simulations in the future.
Collapse
Affiliation(s)
- Hong-Jun Zhou
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China.
| | | | | | | | | | | |
Collapse
|
45
|
"On-The-Fly" Non-Adiabatic Dynamics Simulations on Photoinduced Ring-Closing Reaction of a Nucleoside-Based Diarylethene Photoswitch. Molecules 2021; 26:molecules26092724. [PMID: 34066431 PMCID: PMC8125013 DOI: 10.3390/molecules26092724] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2021] [Revised: 04/29/2021] [Accepted: 05/03/2021] [Indexed: 11/17/2022] Open
Abstract
Nucleoside-based diarylethenes are emerging as an especial class of photochromic compounds that have potential applications in regulating biological systems using noninvasive light with high spatio-temporal resolution. However, relevant microscopic photochromic mechanisms at atomic level of these novel diarylethenes remain to be explored. Herein, we have employed static electronic structure calculations (MS-CASPT2//M06-2X, MS-CASPT2//SA-CASSCF) in combination with non-adiabatic dynamics simulations to explore the related photoinduced ring-closing reaction of a typical nucleoside-based diarylethene photoswitch, namely, PS-IV. Upon excitation with UV light, the open form PS-IV can be excited to a spectroscopically bright S1 state. After that, the molecule relaxes to the conical intersection region within 150 fs according to the barrierless relaxed scan of the C1–C6 bond, which is followed by an immediate deactivation to the ground state. The conical intersection structure is very similar to the ground state transition state structure which connects the open and closed forms of PS-IV, and therefore plays a crucial role in the photochromism of PS-IV. Besides, after analyzing the hopping structures, we conclude that the ring closing reaction cannot complete in the S1 state alone since all the C1–C6 distances of the hopping structures are larger than 2.00 Å. Once hopping to the ground state, the molecules either return to the original open form of PS-IV or produce the closed form of PS-IV within 100 fs, and the ring closing quantum yield is estimated to be 56%. Our present work not only elucidates the ultrafast photoinduced pericyclic reaction of the nucleoside-based diarylethene PS-IV, but can also be helpful for the future design of novel nucleoside-based diarylethenes with better performance.
Collapse
|
46
|
Abstract
This review article focuses on the understanding of intersystem crossing (ISC) in molecules. It addresses readers who are interested in the phenomenon of intercombination transitions between states of different electron spin multiplicities but are not familiar with relativistic quantum chemistry. Among the spin-dependent interaction terms that enable a crossover between states of different electron spin multiplicities, spin-orbit coupling (SOC) is by far the most important. If SOC is small or vanishes by symmetry, ISC can proceed by electronic spin-spin coupling (SSC) or hyperfine interaction (HFI). Although this review discusses SSC- and HFI-based ISC, the emphasis is on SOC-based ISC. In addition to laying the theoretical foundations for the understanding of ISC, the review elaborates on the qualitative rules for estimating transition probabilities. Research on the mechanisms of ISC has experienced a major revival in recent years owing to its importance in organic light-emitting diodes (OLEDs). Exemplified by challenging case studies, chemical substitution and solvent environment effects are discussed with the aim of helping the reader to understand and thereby get a handle on the factors that steer the efficiency of ISC.
Collapse
Affiliation(s)
- Christel M Marian
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University, Düsseldorf 40204, Germany;
| |
Collapse
|
47
|
Abstract
In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Institut de Chimie Radicalaire, CNRS 7273, Aix-Marseille University, 13013 Marseille, France;
| | - Dmitry A Fedorov
- Oak Ridge Associated Universities, Oak Ridge, Tennessee 37830, USA;
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, Nevada 89557-0216, USA;
| |
Collapse
|
48
|
Leng C, You S, Si Y, Qin HM, Liu J, Huang WQ, Li K. Unraveling the Mechanism of Near-Infrared Thermally Activated Delayed Fluorescence of TPA-Based Molecules: Effect of Hydrogen Bond Steric Hindrance. J Phys Chem A 2021; 125:2905-2912. [PMID: 33822612 DOI: 10.1021/acs.jpca.1c00739] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A recently synthesized novel molecule (named CAT-1) exhibits intriguing near-infrared (NIR) thermally activated delayed fluorescence (TADF) close to 1000 nm wavelength; however, the mechanism behind these intrinsic properties is not fully understood. Herein, we unravel that the fluorescence emission spectrum with a broad wavelength range (770-950 nm) of CAT-1 is primarily induced by hydrogen bond steric hindrance based on density functional theory and Marcus theory. It is found that the hydrogen bond steric hindrance plays a critical role in inhibiting the twist of the configuration of different excited states, which leads to the minor driving force for fast electron trapping between the excited states, as well as small internal reorganization energy caused by less changed geometric configuration. Furthermore, such steric hindrance will cause a more distorted plane, resulting in a less favorable electron delocalization. A faster reverse intersystem crossing (RISC) rate is then obtained due to the nearly unchanged conformation between excited states caused by steric hindrance, although the spin-orbit coupling is small. Consequently, the NIR TADF with a longer wavelength can be emitted in CAT-1. This work shows that the hydrogen bond steric hindrance can fine-tune the electronic interactions of the donor and acceptor units to control the TADF.
Collapse
Affiliation(s)
- Can Leng
- Science and Technology on Parallel and Distributed Processing Laboratory, National University of Defense Technology, Changsha 410073, China.,Laboratory of Software Engineering for Complex Systems, National University of Defense Technology, Changsha 410073, China.,National Supercomputer Center in Changsha, Changsha 410082, China
| | - Sheng You
- National Supercomputer Center in Changsha, Changsha 410082, China
| | - Yubing Si
- College of Chemistry, Zhengzhou University, Zhengzhou 450001, China
| | - Hai-Mei Qin
- College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen 361005, China
| | - Jie Liu
- Science and Technology on Parallel and Distributed Processing Laboratory, National University of Defense Technology, Changsha 410073, China.,Laboratory of Software Engineering for Complex Systems, National University of Defense Technology, Changsha 410073, China
| | - Wei-Qing Huang
- Department of Applied Physics, School of Physics and Electronics, Hunan University, Changsha 410082, China
| | - Keqin Li
- College of Computer Science and Electronic Engineering, Hunan University, Changsha 410082, China.,Department of Computer Science, State University of New York, New Paltz, New York 12561, United States
| |
Collapse
|
49
|
Bian X, Wu Y, Teh HH, Zhou Z, Chen HT, Subotnik JE. Modeling nonadiabatic dynamics with degenerate electronic states, intersystem crossing, and spin separation: A key goal for chemical physics. J Chem Phys 2021; 154:110901. [DOI: 10.1063/5.0039371] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yanze Wu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Hung-Hsuan Teh
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Zeyu Zhou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Hsing-Ta Chen
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
50
|
Li ZW, Yang JJ, Liu XY, Fang WH, Wang H, Cui G. Chemical Bonding as a New Avenue for Controlling Excited-State Properties and Excitation Energy-Transfer Processes in Zinc Phthalocyanine-Fullerene Dyads. Chemistry 2021; 27:4159-4167. [PMID: 33372312 DOI: 10.1002/chem.202004850] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Indexed: 11/08/2022]
Abstract
Whether chemical bonding can regulate the excited-state and optoelectronic properties of donor-acceptor dyads has been largely elusive. In this work, we used electronic structure and nonadiabatic dynamics methods to explore the excited-state properties of covalently bonded zinc phthalocyanine (ZnPc)-fullerene (C60 ) dyads with a 6-6 (or 5-6) bonding configuration in which ZnPc is bonded to two carbon atoms shared by the two hexagonal rings (or a pentagonal and a hexagonal ring) in C60 . In both cases, the locally excited (LE) states on ZnPc are spectroscopically bright. However, their different chemical bonding differentiates the electronic interactions between ZnPc and C60 . In the 5-6 bonding configuration, the LE states on ZnPc are much higher in energy than the LE states on C60 . Thus, the excitation energy transfer from ZnPc to C60 is thermodynamically favorable. On the other hand, in the 6-6 bonding configuration, such a process is inhibited because the LE states on ZnPc are the lowest ones. More detailed mechanisms are elucidated from nonadiabatic dynamics simulations. In the 6-6 bonding configuration, no excitation energy transfer was observed. In contrast, in the 5-6 bonding configuration, several LE and charge-transfer (CT) excitons were shown to participate in the energy-transfer process. Further analysis reveals that the photoinduced energy transfer is mediated by a CT exciton, such that electron- and hole-transfer processes take place in a concerted but asynchronous manner in the excitation energy transfer. It is also found that high-level electronic structure methods including exciton effects are indispensable to accurately describe photoinduced energy- and electron-transfer processes. Furthermore, this work opens up new avenues for regulating the excited-state properties of molecular donor-acceptor dyads by means of chemical bonding.
Collapse
Affiliation(s)
- Zi-Wen Li
- College of Chemistry, Beijing Normal University, Beijing, 100875, P. R. China
| | - Jia-Jia Yang
- College of Chemistry, Beijing Normal University, Beijing, 100875, P. R. China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu, 610068, P. R. China
| | - Wei-Hai Fang
- College of Chemistry, Beijing Normal University, Beijing, 100875, P. R. China
| | - Haobin Wang
- Department of Chemistry, University of Colorado Denver, Denver, Colorado, 80217-3364, USA
| | - Ganglong Cui
- College of Chemistry, Beijing Normal University, Beijing, 100875, P. R. China
| |
Collapse
|