1
|
Zhu Y, Peng J, Xu C, Lan Z. Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation. J Phys Chem Lett 2024:9601-9619. [PMID: 39270134 DOI: 10.1021/acs.jpclett.4c01751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2024]
Abstract
The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.
Collapse
Affiliation(s)
- Yifei Zhu
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
2
|
Hokyo H, Ito K, Kalia RK, Kapadia R, Nakano A, Shimamura K, Shimojo F, Vashishta P. Photoinduced Negative Differential Resistance at a Graphene/Silicon Interface: A Nonadiabatic Quantum Molecular Dynamics Study. J Phys Chem Lett 2024; 15:9226-9232. [PMID: 39225635 DOI: 10.1021/acs.jpclett.4c02272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
The oscillatory retinal neuron (ORN) is a promising technology for achieving in-sensor cognitive image computing without external power. While its operation is based on photoinduced negative differential resistance (NDR) at a graphene/silicon interface to directly convert the incident optical signal into voltage oscillations, the optoelectronic mechanism of the NDR remains elusive. Here, nonadiabatic quantum molecular dynamics simulations show that the interplay of band alignment and charge transfer rates of photoexcited carriers at varying applied voltages gives rise to NDR at a graphene/silicon interface under illumination. Such intrinsic NDR at an interface, along with extrinsic circuit-level factors, could enable the much needed rational design of desired image computing functionality of ORN devices in the era of ubiquitous AI on edge devices.
Collapse
Affiliation(s)
- Hinata Hokyo
- Collaboratory for Advanced Computing and Simulation, University of Southern California, Los Angeles, California 90089-0242, United States
- Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan
| | - Kai Ito
- Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan
| | - Rajiv K Kalia
- Collaboratory for Advanced Computing and Simulation, University of Southern California, Los Angeles, California 90089-0242, United States
| | - Rehan Kapadia
- Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California 90089-2560, United States
| | - Aiichiro Nakano
- Collaboratory for Advanced Computing and Simulation, University of Southern California, Los Angeles, California 90089-0242, United States
| | - Kohei Shimamura
- Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan
| | - Fuyuki Shimojo
- Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan
| | - Priya Vashishta
- Collaboratory for Advanced Computing and Simulation, University of Southern California, Los Angeles, California 90089-0242, United States
| |
Collapse
|
3
|
Chen W, Xie H, Chen Z, Qu Z, Ren H, Li X. Influence of Singlet Oxygen in H + O 2 Collision Reaction. J Phys Chem A 2024; 128:7226-7234. [PMID: 39140992 DOI: 10.1021/acs.jpca.4c03317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2024]
Abstract
The collision reaction of H + O2 = OH + O is a pivotal step in combustion. To investigate the influence of singlet oxygen on this reaction, we computed potential energy surfaces (PESs) for all six lowest states using high-level ab initio methods and coupled them with embedded atom neural network (EANN) fitting. By integrating quasi-classical trajectory (QCT) with trajectory surface hopping (TSH) based on the fitted PESs, we simulated the dynamics of both ground- and excited-states to derive the reaction rate constants for the forward and reverse processes. The results reveal that the forward reaction facilitates radical generation, promoting combustion reactions. Furthermore, calculations of reverse reaction rate constants indicate that all electronic states ultimately yield ground-state oxygen, leading to radical deactivation and exerting an inhibitory effect on combustion processes.
Collapse
Affiliation(s)
- Wenlan Chen
- School of Chemical Engineering, Sichuan University, Chengdu 610065, China
| | - Haohan Xie
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Zhizhou Chen
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Zexing Qu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Haisheng Ren
- School of Chemical Engineering, Sichuan University, Chengdu 610065, China
- Engineering Research Center of Combustion and Cooling for Aerospace Power, Ministry of Education, Sichuan University, Chengdu 610065, China
| | - Xiangyuan Li
- School of Chemical Engineering, Sichuan University, Chengdu 610065, China
- Engineering Research Center of Combustion and Cooling for Aerospace Power, Ministry of Education, Sichuan University, Chengdu 610065, China
| |
Collapse
|
4
|
Myers CA, Miyazaki K, Trepl T, Isborn CM, Ananth N. GPU-accelerated on-the-fly nonadiabatic semiclassical dynamics. J Chem Phys 2024; 161:084114. [PMID: 39193942 DOI: 10.1063/5.0223628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2024] [Accepted: 08/11/2024] [Indexed: 08/29/2024] Open
Abstract
GPU-accelerated on-the-fly nonadiabatic dynamics is enabled by interfacing the linearized semiclassical dynamics approach with the TeraChem electronic structure program. We describe the computational workflow of the "PySCES" code interface, a Python code for semiclassical dynamics with on-the-fly electronic structure, including parallelization over multiple GPU nodes. We showcase the abilities of this code and present timings for two benchmark systems: fulvene solvated in acetonitrile and a charge transfer system in which a photoexcited zinc-phthalocyanine donor transfers charge to a fullerene acceptor through multiple electronic states on an ultrafast timescale. Our implementation paves the way for an efficient semiclassical approach to model the nonadiabatic excited state dynamics of complex molecules, materials, and condensed phase systems.
Collapse
Affiliation(s)
- Christopher A Myers
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Ken Miyazaki
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - Thomas Trepl
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
5
|
Lacroix T, Le Dé B, Riva A, Dunnett AJ, Chin AW. MPSDynamics.jl: Tensor network simulations for finite-temperature (non-Markovian) open quantum system dynamics. J Chem Phys 2024; 161:084116. [PMID: 39206827 DOI: 10.1063/5.0223107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2024] [Accepted: 08/08/2024] [Indexed: 09/04/2024] Open
Abstract
The MPSDynamics.jl package provides an easy-to-use interface for performing open quantum systems simulations at zero and finite temperatures. The package has been developed with the aim of studying non-Markovian open system dynamics using the state-of-the-art numerically exact Thermalized-Time Evolving Density operator with Orthonormal Polynomials Algorithm based on environment chain mapping. The simulations rely on a tensor network representation of the quantum states as matrix product states (MPS) and tree tensor network states. Written in the Julia programming language, MPSDynamics.jl is a versatile open-source package providing a choice of several variants of the Time-Dependent Variational Principle method for time evolution (including novel bond-adaptive one-site algorithms). The package also provides strong support for the measurement of single and multi-site observables, as well as the storing and logging of data, which makes it a useful tool for the study of many-body physics. It currently handles long-range interactions, time-dependent Hamiltonians, multiple environments, bosonic and fermionic environments, and joint system-environment observables.
Collapse
Affiliation(s)
- Thibaut Lacroix
- Institut für Theoretische Physik und IQST, Universität Ulm, Albert-Einstein-Allee 11, D-89081 Ulm, Germany
| | - Brieuc Le Dé
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 Place Jussieu, 75005 Paris, France
| | - Angela Riva
- LPENS, Département de Physique, École Normale Supérieure, Centre Automatique et Systèmes (CAS), MINES ParisTech, Université PSL, Sorbonne Université, CNRS, Inria, 75005 Paris, France
| | - Angus J Dunnett
- Multiverse Computing, 7 rue de la Croix Martre, 91120 Palaiseau, France
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 Place Jussieu, 75005 Paris, France
| |
Collapse
|
6
|
Curchod BFE, Orr-Ewing AJ. Perspective on Theoretical and Experimental Advances in Atmospheric Photochemistry. J Phys Chem A 2024; 128:6613-6635. [PMID: 39021090 PMCID: PMC11331530 DOI: 10.1021/acs.jpca.4c03481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2024] [Revised: 07/02/2024] [Accepted: 07/03/2024] [Indexed: 07/20/2024]
Abstract
Research that explores the chemistry of Earth's atmosphere is central to the current understanding of global challenges such as climate change, stratospheric ozone depletion, and poor air quality in urban areas. This research is a synergistic combination of three established domains: earth observation, for example, using satellites, and in situ field measurements; computer modeling of the atmosphere and its chemistry; and laboratory measurements of the properties and reactivity of gas-phase molecules and aerosol particles. The complexity of the interconnected chemical and photochemical reactions which determine the composition of the atmosphere challenges the capacity of laboratory studies to provide the spectroscopic, photochemical, and kinetic data required for computer models. Here, we consider whether predictions from computational chemistry using modern electronic structure theory and nonadiabatic dynamics simulations are becoming sufficiently accurate to supplement quantitative laboratory data for wavelength-dependent absorption cross-sections, photochemical quantum yields, and reaction rate coefficients. Drawing on presentations and discussions from the CECAM workshop on Theoretical and Experimental Advances in Atmospheric Photochemistry held in March 2024, we describe key concepts in the theory of photochemistry, survey the state-of-the-art in computational photochemistry methods, and compare their capabilities with modern experimental laboratory techniques. From such considerations, we offer a perspective on the scope of computational (photo)chemistry methods based on rigorous electronic structure theory to become a fourth core domain of research in atmospheric chemistry.
Collapse
|
7
|
Liu D, Wang B, Vasenko AS, Prezhdo OV. Decoherence ensures convergence of non-adiabatic molecular dynamics with number of states. J Chem Phys 2024; 161:064104. [PMID: 39120030 DOI: 10.1063/5.0222557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2024] [Accepted: 07/24/2024] [Indexed: 08/10/2024] Open
Abstract
Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems.
Collapse
Affiliation(s)
| | - Bipeng Wang
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Andrey S Vasenko
- HSE University, 101000 Moscow, Russia
- Donostia International Physics Center (DIPC), 20018 San Sebastián-Donostia, Euskadi, Spain
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
8
|
Salvadori G, Mennucci B. Analogies and Differences in the Photoactivation Mechanism of Bathy and Canonical Bacteriophytochromes Revealed by Multiscale Modeling. J Phys Chem Lett 2024; 15:8078-8084. [PMID: 39087732 PMCID: PMC11376688 DOI: 10.1021/acs.jpclett.4c01823] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2024]
Abstract
Bacteriophytochromes are light-sensing biological machines that switch between two photoreversible states, Pr and Pfr. Their relative stability is opposite in canonical and bathy bacteriophytochromes, but in both cases the switch between them is triggered by the photoisomerization of an embedded bilin chromophore. We applied an integrated multiscale strategy of excited-state QM/MM nonadiabatic dynamics and (QM/)MM molecular dynamics simulations with enhanced sampling techniques to the Agrobacterium fabrum bathy phytochrome and compared the results with those obtained for the canonical phytochrome Deinococcus radiodurans. Contrary to what recently suggested, we found that photoactivation in both phytochromes is triggered by the same hula-twist motion of the bilin chromophore. However, only in the bathy phytochrome, the bilin reaches the final rotated structure already in the first intermediate. This allows a reorientation of the binding pocket in a microsecond time scale, which can propagate through the entire protein causing the spine to tilt.
Collapse
Affiliation(s)
- Giacomo Salvadori
- Institute for Computational Biomedicine (INM-9/IAS-5), Forschungszentrum Jülich, 52428 Jülich, Germany
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, 56124 Pisa, Italy
| |
Collapse
|
9
|
Pandey P, Durga Prasad M. Time-Dependent Multireference Coupled-Cluster Method (TDMRCCM) for the Bath-Dynamics in System-Bath Approach to Nonadiabatic Dynamics. J Phys Chem A 2024; 128:5990-5998. [PMID: 39012785 DOI: 10.1021/acs.jpca.4c01432] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2024]
Abstract
The time-dependent multireference coupled-cluster method (TDMRCCM) fits well in the scheme of the system-bath separation used to study the nonadiabatic dynamics. In TDMRCCM, a projection operator is defined as one that projects the full Hilbert space onto the space spanned by the collection of system degrees of freedom, called the model space. The inverse of this projection operator is a wave operator that acts on the model space and takes its projection back to the full wave function. This wave operator is defined as an exponential of the excitation terms and, hence, can be expanded into a Taylor series, which we have truncated in this work at the second-order of excitations. The attraction of using TDMRCCM for describing the bath dynamics is due to the exponential nature of the ansatz used in the method, which makes it possible for the higher-order excitations to be absorbed by the lower-order terms, even upon truncating the series. This improves the accuracy of the numerical calculations using TDMRCCM as an approximation for the bath-mode dynamics in the system-bath framework for nonadiabatic dynamics. We present the theoretical details of TDMRCCM and the numerical results for implementing this method to study the dynamics in the butatriene cation.
Collapse
Affiliation(s)
- Prachi Pandey
- School of Chemistry, University of Hyderabad, Hyderabad, Telangana 500046, India
| | - M Durga Prasad
- School of Chemistry, University of Hyderabad, Hyderabad, Telangana 500046, India
| |
Collapse
|
10
|
Meng G, Gardner J, Hertl N, Dou W, Maurer RJ, Jiang B. First-Principles Nonadiabatic Dynamics of Molecules at Metal Surfaces with Vibrationally Coupled Electron Transfer. PHYSICAL REVIEW LETTERS 2024; 133:036203. [PMID: 39094165 DOI: 10.1103/physrevlett.133.036203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Accepted: 06/12/2024] [Indexed: 08/04/2024]
Abstract
Accurate description of nonadiabatic dynamics of molecules at metal surfaces involving electron transfer has been a long-standing challenge for theory. Here, we tackle this problem by first constructing high-dimensional neural network diabatic potentials including state crossings determined by constrained density functional theory, then applying mixed quantum-classical surface hopping simulations to evolve coupled electron-nuclear motion. Our approach accurately describes the nonadiabatic effects in CO scattering from Au(111) without empirical parameters and yields results agreeing well with experiments under various conditions for this benchmark system. We find that both adiabatic and nonadiabatic energy loss channels have important contributions to the vibrational relaxation of highly vibrationally excited CO(v_{i}=17), whereas relaxation of low vibrationally excited states of CO(v_{i}=2) is weak and dominated by nonadiabatic energy loss. The presented approach paves the way for accurate first-principles simulations of electron transfer mediated nonadiabatic dynamics at metal surfaces.
Collapse
|
11
|
Yang H, Dong H, Martens CC, Zheng Y. Nonadiabatic Coupling-Induced Quantum Coherence in Two-Dimensional Materials. J Phys Chem Lett 2024; 15:6363-6369. [PMID: 38857307 PMCID: PMC11194825 DOI: 10.1021/acs.jpclett.4c01140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Revised: 05/20/2024] [Accepted: 06/05/2024] [Indexed: 06/12/2024]
Abstract
Two-dimensional materials provide a rich platform demonstrating quantum effects, and the process of electron-hole recombination occurring in them has significant applications in the fields of the photocatalytic and optoelectronic community. Here, we present nonadiabatic coupling-induced quantum coherence and quantum beats in Al-doped blue phosphorene. The work improves our understanding and utilization of nonadiabatic coupling in low-dimensional materials from a new perspective. In addition, our investigations provide meaningful guidance for manipulating quantum coherence in low-dimensional materials and promoting their novel optoelectronic properties.
Collapse
Affiliation(s)
- Huan Yang
- School
of Physics, Shandong University, Jinan 250100, China
| | - Hao Dong
- School
of Physics, Shandong University, Jinan 250100, China
| | - Craig C. Martens
- Department
of Chemistry, University of California,
Irvine, Irvine, California 92697-2025, United States
| | - Yujun Zheng
- School
of Physics, Shandong University, Jinan 250100, China
| |
Collapse
|
12
|
Gu K, Lin S. Advances in the Dynamics of Adsorbate Diffusion on Metal Surfaces: Focus on Hydrogen and Oxygen. Chemphyschem 2024; 25:e202400083. [PMID: 38511509 DOI: 10.1002/cphc.202400083] [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: 01/28/2024] [Revised: 02/26/2024] [Accepted: 03/21/2024] [Indexed: 03/22/2024]
Abstract
Adsorbates on metal surfaces are typically formed from the dissociative chemisorption of molecules occurring at gas-solid interfaces. These adsorbed species exhibit unique diffusion behaviors on metal surfaces, which are influenced by their translational energy. They play crucial roles in various fields, including heterogeneous catalysis and corrosion. This review examines recent theoretical advancements in understanding the diffusion dynamics of adsorbates on metal surfaces, with a specific emphasis on hydrogen and oxygen atoms. The diffusion processes of adsorbates on metal surfaces involve two energy transfer mechanisms: surface phonons and electron-hole pair excitations. This review also surveys new theoretical methods, including the characterization of the electron-hole pair excitation within electronic friction models, the acceleration of quantum chemistry calculations through machine learning, and the treatment of atomic nuclear motion from both quantum mechanical and classical perspectives. Furthermore, this review offers valuable insights into how energy transfer, nuclear quantum effects, supercell sizes, and the topography of potential energy surfaces impact the diffusion behavior of hydrogen and oxygen species on metal surfaces. Lastly, some preliminary research proposals are presented.
Collapse
Affiliation(s)
- Kaixuan Gu
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou, 350002, China
| | - Sen Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou, 350002, China
| |
Collapse
|
13
|
Chen BJ, Pradhan E, Nooijen M, Zeng T. Adiabat-to-Diabat Angle in Seam Space: Renner-Teller-Type and Pseudo-Jahn-Teller-Type Problems. Chemphyschem 2024; 25:e202400130. [PMID: 38427966 DOI: 10.1002/cphc.202400130] [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: 02/05/2024] [Revised: 02/29/2024] [Accepted: 03/01/2024] [Indexed: 03/03/2024]
Abstract
In this study, we examine the adiabat-to-diabat (ATD) angles for trajectories in 2-dimensional vibrational subspace of the seam space of two degenerate states. In circulating around the tangential touching degeneracy center, the ATD angle is changed by2 π ${2\pi }$ or 0, similar to the Renner-Teller problem and the pseudo-Jahn-Teller problem, respectively. These ATD angle profiles may be indistinguishable from those of circulating multiple conical intersections or a pseudo-Jahn-Teller center. Methods to discern those seemingly indistinguishable cases are proposed. A sharp zigzag variation of the ATD angle is seen as a feature for trajectories that graze a pseudo-Jahn-Teller-type tangential touching center, in contrast to the monotonic steep variation for grazing a conical intersection or a Renner-Teller-type tangential touching center.
Collapse
Affiliation(s)
- Benny Jun Chen
- Department of Chemistry, University of Waterloo, Waterloo, ON, N2 L3G1, Canada
| | - Ekadashi Pradhan
- Department of Chemistry, York University, Toronto, Ontario, M3 J1P3, Canada
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, ON, N2 L3G1, Canada
| | - Tao Zeng
- Department of Chemistry, York University, Toronto, Ontario, M3 J1P3, Canada
| |
Collapse
|
14
|
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
|
15
|
Wu W, Scholes GD. Foundations of Quantum Information for Physical Chemistry. J Phys Chem Lett 2024; 15:4056-4069. [PMID: 38587240 DOI: 10.1021/acs.jpclett.4c00180] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
Quantum information, a field in which great advances have been made in the past decades, now presents opportunities for advanced chemistry. One roadblock to progress, especially for experimental chemical science, is that new concepts and technical definitions need to be learned. In this paper, we review some basic, but sometimes misunderstood, concepts of quantum information based on the mathematical formulation of quantum mechanics that will be useful for chemists interested in discovering ways that chemistry can contribute to the quantum information field. We cover topics including qubits and their density matrix formalism, quantum measurement as a quantum operation, information theory, and entanglement. We focus on the difference between the concepts in the quantum context and the classic context. We also discuss the relation and distinction among entanglement, correlation, and coherence. We aim to clarify the rigorous definition of these concepts and then indicate some examples in physical chemistry.
Collapse
Affiliation(s)
- Weijun Wu
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Gregory D Scholes
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
16
|
Janoš J, Figueira Nunes JP, Hollas D, Slavíček P, Curchod BFE. Predicting the photodynamics of cyclobutanone triggered by a laser pulse at 200 nm and its MeV-UED signals-A trajectory surface hopping and XMS-CASPT2 perspective. J Chem Phys 2024; 160:144305. [PMID: 38591685 DOI: 10.1063/5.0203105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 03/14/2024] [Indexed: 04/10/2024] Open
Abstract
This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born-Oppenheimer molecular dynamics with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from Born-Oppenheimer molecular dynamics coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S2 can follow two pathways for its nonradiative decay: (i) a ring-opening in S2 and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S1, where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function-ΔPDF(r)] for each (interpolated) pathway as well as a time-resolved signal [ΔPDF(r,t) and ΔI/I(s,t)] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals [both ΔPDF(r,t) and ΔI/I(s,t)] and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals.
Collapse
Affiliation(s)
- Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | | | - Daniel Hollas
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
17
|
Xu J, Shi Z, Wang L. Consistent Construction of the Density Matrix from Surface Hopping Trajectories. J Chem Theory Comput 2024; 20:2349-2361. [PMID: 38490993 DOI: 10.1021/acs.jctc.4c00137] [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/2024]
Abstract
Proper construction of the density matrix based on surface hopping trajectories remains a difficult problem. Due to the well-known overcoherence in traditional surface hopping simulations, the electronic wave function cannot be used directly. In this work, we propose a consistent density matrix construction method, which takes the advantage of occupation of active states to rescale the coherence calculated by wave functions and ensures the intrinsic consistency of the density matrix. This new trajectory analysis method can be used for both Tully's fewest switches surface hopping (FSSH) and our recently proposed branching corrected surface hopping (BCSH). As benchmarked in both one- and two-dimensional standard scattering models, the new approach combined with BCSH trajectories achieves highly accurate time-dependent spatial distributions of adiabatic populations and coherence compared to exact quantum results.
Collapse
Affiliation(s)
- Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
18
|
Guo X, Li G, Shi Z, Wang L. Surface Hopping with Reliable Wave Function by Introducing Auxiliary Wave Packets to Trajectory Branching. J Phys Chem Lett 2024:3345-3353. [PMID: 38498301 DOI: 10.1021/acs.jpclett.4c00437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
It is well-known that the widely utilized fewest switches surface hopping method suffers from the severe overcoherence problem, and thus adiabatic populations calculated by wave functions are generally inferior to those based on active states. More importantly, to achieve a complete description of nonadiabatic dynamics, the density matrix is essential. In this paper, we present an auxiliary branching corrected surface hopping (A-BCSH) method that introduces auxiliary wave packets (WPs) on the adiabatic potential energy surfaces for trajectory branching. Both rapid and gradual separation of WP components on different surfaces are characterized, and thus the correct decoherence time along each trajectory is captured. As demonstrated in the three standard Tully models, A-BCSH exhibits excellent internal consistency. Namely, close adiabatic populations are obtained based on both wave functions and active states. In particular, A-BCSH successfully obtains a reliable time-dependent spatial distribution of the density matrix, which relies only on electronic wave functions. Due to its high performance, our A-BCSH method provides a new and highly promising perspective on further development of more consistent surface hopping with reliable wave function.
Collapse
Affiliation(s)
- Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Guijie Li
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
19
|
Bao S, Raymond N, Nooijen M. Time dependent vibrational electronic coupled cluster (VECC) theory for non-adiabatic nuclear dynamics. J Chem Phys 2024; 160:094105. [PMID: 38426527 DOI: 10.1063/5.0190034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2023] [Accepted: 01/28/2024] [Indexed: 03/02/2024] Open
Abstract
A time-dependent vibrational electronic coupled-cluster (VECC) approach is proposed to simulate photo-electron/UV-VIS absorption spectra as well as time-dependent properties for non-adiabatic vibronic models, going beyond the Born-Oppenheimer approximation. A detailed derivation of the equations of motion and a motivation for the ansatz are presented. The VECC method employs second-quantized bosonic construction operators and a mixed linear and exponential ansatz to form a compact representation of the time-dependent wave-function. Importantly, the method does not require a basis set, has only a few user-defined inputs, and has a classical (polynomial) scaling with respect to the number of degrees of freedom (of the vibronic model), resulting in a favorable computational cost. In benchmark applications to small models and molecules, the VECC method provides accurate results compared to multi-configurational time-dependent Hartree calculations when predicting short-time dynamical properties (i.e., photo-electron/UV-VIS absorption spectra) for non-adiabatic vibronic models. To illustrate the capabilities, the VECC method is also successfully applied to a large vibronic model for hexahelicene with 14 electronic states and 63 normal modes, developed in the group by Aranda and Santoro [J. Chem. Theory Comput. 17, 1691, (2021)].
Collapse
Affiliation(s)
- Songhao Bao
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Neil Raymond
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
20
|
Makri N. Kink Sum for Long-Memory Small Matrix Path Integral Dynamics. J Phys Chem B 2024. [PMID: 38437738 DOI: 10.1021/acs.jpcb.3c08282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
The small matrix decomposition of the real-time path integral (SMatPI) allows for numerically exact and efficient propagation of the reduced density matrix (RDM) for system-bath Hamiltonians. Its high efficiency lies in the small size of the SMatPI matrices employed in the iterative algorithm, whose size is equal to that of the full RDM. By avoiding the storage and multiplication of large tensors, the SMatPI algorithm is applicable in multistate systems under a variety of conditions. The main computational effort is the evaluation of path sums within the entangled memory length to construct the SMatPI matrices. A number of methods are available for this task, each with its own favorable parameter regime, but calculations with strong system-bath coupling and long memory at low temperatures remain out of reach. The present paper evaluates the path sums by binning the paths (in forward time only) based on their amplitudes, which depend on the number and type of kinks they contain. The algorithm is very efficient, leading to a dramatic acceleration of path sums and significantly extending the accessible memory length in the most challenging regimes.
Collapse
Affiliation(s)
- Nancy Makri
- Departments of Chemistry and Physics, University of Illinois, 505 S. Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
21
|
Salvadori G, Mazzeo P, Accomasso D, Cupellini L, Mennucci B. Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response. J Mol Biol 2024; 436:168358. [PMID: 37944793 DOI: 10.1016/j.jmb.2023.168358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/28/2023] [Accepted: 11/02/2023] [Indexed: 11/12/2023]
Abstract
In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.
Collapse
Affiliation(s)
- Giacomo Salvadori
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Patrizia Mazzeo
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Davide Accomasso
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| |
Collapse
|
22
|
Cao Y, Halls MD, Friesner RA. Highly efficient implementation of analytic nonadiabatic derivative couplings within the pseudospectral method. J Chem Phys 2024; 160:084106. [PMID: 38385510 DOI: 10.1063/5.0188277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 01/31/2024] [Indexed: 02/23/2024] Open
Abstract
A pseudospectral implementation of nonadiabatic derivative couplings in the Tamm-Dancoff approximation is reported, and the accuracy and efficiency of the pseudospectral nonadiabatic derivative couplings are studied. Our results demonstrate that the pseudospectral method provides mean absolute errors of 0.2%-1.9%, while providing a significant speedup. Benchmark calculations on fullerenes (Cn, n up to 100) using B3LYP achieved 10- to 15-fold, 8- to 17-fold, and 43- to 75-fold speedups for 6-31G**, 6-31++G**, and cc-pVTZ basis sets, respectively, when compared to the conventional spectral method.
Collapse
Affiliation(s)
- Yixiang Cao
- Schrödinger Inc., 1540 Broadway, 24th Floor, New York, New York 10036, USA
| | - Mathew D Halls
- Schrödinger Inc., 9868 Scranton, Suite 3200, San Diego, California 92121, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
23
|
Ramos P, Friedman H, Li BY, Garcia C, Sletten E, Caram JR, Jang SJ. Nonadiabatic Derivative Couplings through Multiple Franck-Condon Modes Dictate the Energy Gap Law for Near and Short-Wave Infrared Dye Molecules. J Phys Chem Lett 2024; 15:1802-1810. [PMID: 38329913 DOI: 10.1021/acs.jpclett.3c02629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Near infrared (NIR, 700-1000 nm) and short-wave infrared (SWIR, 1000-2000 nm) dye molecules exhibit significant nonradiative decay rates from the first singlet excited state to the ground state. While these trends can be empirically explained by a simple energy gap law, detailed mechanisms of nearly universal behavior have remained unsettled for many cases. Theoretical and experimental results for two representative NIR/SWIR dye molecules reported here clarify the key mechanism for the observed energy gap law behavior. It is shown that the first derivative nonadiabatic coupling terms serve as major coupling pathways for nonadiabatic decay processes from the first excited singlet state to the ground state for these NIR and SWIR dye molecules and that vibrational modes other than the highest frequency modes also make significant contributions to the rate. This assessment is corroborated by further theoretical comparison with possible alternative mechanisms of intersystem crossing to triplet states and also by comparison with experimental data for deuterated molecules.
Collapse
Affiliation(s)
- Pablo Ramos
- Department of Chemistry and Biochemistry, Queens College, City University of New York, 65-30 Kissena Boulevard, New York, New York 11367, United States
| | - Hannah Friedman
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Barry Y Li
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Cesar Garcia
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Ellen Sletten
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Justin R Caram
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Seogjoo J Jang
- Department of Chemistry and Biochemistry, Queens College, City University of New York, 65-30 Kissena Boulevard, New York, New York 11367, United States
- Chemistry and Physics PhD programs, Graduate Center, City University of New York, New York, New York 10016, United States
| |
Collapse
|
24
|
Zhou JG, Shu Y, Wang Y, Leszczynski J, Prezhdo O. Dissociation Time, Quantum Yield, and Dynamic Reaction Pathways in the Thermolysis of trans-3,4-Dimethyl-1,2-dioxetane. J Phys Chem Lett 2024; 15:1846-1855. [PMID: 38334951 PMCID: PMC10895692 DOI: 10.1021/acs.jpclett.3c03578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 01/26/2024] [Accepted: 02/06/2024] [Indexed: 02/10/2024]
Abstract
The thermolysis of trans-3,4-dimethyl-1,2-dioxetane is studied by trajectory surface hopping. The significant difference between long and short dissociation times is rationalized by frustrated dissociations and the time spent in triplet states. If the C-C bond breaks through an excited state channel, then the trajectory passes over a ridge of the potential energy surface of that state. The calculated triplet quantum yields match the experimental results. The dissociation half-times and quantum yields follow the same ascending order as per the product states, justifying the conjecture that the longer dissociation time leads to a higher quantum yield, proposed in the context of the methylation effect. The populations of the molecular Coulomb Hamiltonian and diagonal states reach equilibrium, but the triplet populations with different Sz components fluctuate indefinitely. Certain initial velocities, leading the trajectories to given product states, can be identified as the most characteristic features for sorting trajectories according to their product states.
Collapse
Affiliation(s)
- Jian-Ge Zhou
- Interdisciplinary
Nanotoxicity Center, Department of Chemistry, Physics and Atmospheric
Sciences, Jackson State University, Jackson, Mississippi 39217, United States
| | - Yinan Shu
- Department
of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Yuchen Wang
- Department
of Chemistry and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Jerzy Leszczynski
- Interdisciplinary
Nanotoxicity Center, Department of Chemistry, Physics and Atmospheric
Sciences, Jackson State University, Jackson, Mississippi 39217, United States
| | - Oleg Prezhdo
- Department
of Chemistry and Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
25
|
Li X, Lubbers N, Tretiak S, Barros K, Zhang Y. Machine Learning Framework for Modeling Exciton Polaritons in Molecular Materials. J Chem Theory Comput 2024; 20:891-901. [PMID: 38168674 DOI: 10.1021/acs.jctc.3c01068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
A light-matter hybrid quasiparticle, called a polariton, is formed when molecules are strongly coupled to an optical cavity. Recent experiments have shown that polariton chemistry can manipulate chemical reactions. Polariton chemistry is a collective phenomenon, and its effects increase with the number of molecules in a cavity. However, simulating an ensemble of molecules in the excited state coupled to a cavity mode is theoretically and computationally challenging. Recent advances in machine learning (ML) techniques have shown promising capabilities in modeling ground-state chemical systems. This work presents a general protocol to predict excited-state properties, such as energies, transition dipoles, and nonadiabatic coupling vectors with the hierarchically interacting particle neural network. ML predictions are then applied to compute the potential energy surfaces and electronic spectra of a prototype azomethane molecule in the collective coupling scenario. These computational tools provide a much-needed framework to model and understand many molecules' emerging excited-state polariton chemistry.
Collapse
Affiliation(s)
- Xinyang Li
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nicholas Lubbers
- Information Sciences, Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Kipton Barros
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
26
|
Shimray SA, Ningthoujam A, Khaidem DKS, Chipem FAS. Theoretical studies on the photo protective mechanism of curcuminoids. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 305:123449. [PMID: 37774584 DOI: 10.1016/j.saa.2023.123449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Revised: 09/16/2023] [Accepted: 09/20/2023] [Indexed: 10/01/2023]
Abstract
In this work, the deactivation pathways of curcuminoids after photoexcitation was studied by employing density functional theory to explore their UVA radiation screening capacity. A comprehensive computational characterization of the excited-state processes of curcumin, demethoxycurcumin, and bis-demethoxycurcumin was done. The molecules exist in diketo and enol forms which are in equilibrium and interconvertible through keto-enol tautomerism. The enolic forms of each of the studied molecules have eight geometric cis-trans isomers as a result of torsion rotation about three different carbon-carbon double bonds across the aliphatic chain. For each geometric isomer, sixteen possible rotamers are found to exist due to rotation about five different carbon-carbon single bond rotations, also across the skeleton of the aliphatic chain. Upon photoexcitation, the studied molecules follow three main pathways of radiationless decay: (a) rotamerism and interconversion between rotamers of comparable energies which are in equilibrium, (b) interconversion between the cis-trans geometrical isomers where an efficient vibrational relaxation path is formed at ∼90° during torsion rotation about carbon-carbon double bond, and (c) excited state intramolecular proton transfer in a single O-H stretching vibration through a cyclic intramolecular hydrogen bonded ring formed at the centre of the molecule giving back the original structure. The absorption and emission spectra of the molecules were also simulated where the theoretically obtained absorption and emission maxima are close to the reported experimental values.
Collapse
Affiliation(s)
- Sophy A Shimray
- Department of Chemistry, Manipur University, Canchipur 795 003, India
| | - Amar Ningthoujam
- Department of Chemistry, Manipur University, Canchipur 795 003, India
| | | | | |
Collapse
|
27
|
Yang A, Luo J, Xie Z, Chen Q, Xie Q. Ab Initio Quantum Dynamics Simulation of the Impact of Graphene on the Carrier Lifetime of the ZnV 2O 6 Photocatalyst. J Phys Chem Lett 2024; 15:23-33. [PMID: 38127901 DOI: 10.1021/acs.jpclett.3c02387] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
We used a nonadiabatic molecular dynamics simulation to determine the carrier dynamics of a graphene/ZnV2O6 heterostructure in the search for an effective photocatalyst material. The C-2p orbital promotes the wave function overlap, guiding electrons to move between graphene and ZnV2O6, successfully achieving good mixing with the valence and conduction bands in ZnV2O6 materials, which is conducive to supporting carrier migration. The overlap between graphene/ZnV2O6 electrons and hole wave functions is less than that of ZnV2O6, and there is small absolute nonadiabatic coupling. The charge separation caused by graphene increases the carrier lifetime and prevents nonradiative electron-hole recombination. This study reveals the microscopic mechanism of extending the carrier lifetime of ZnV2O6 by introducing graphene, providing useful insights for regulating the electronic structure, promoting electron transfer and ultrafast electron and hole transfer. This strategy provides design considerations for advanced photocatalytic materials.
Collapse
Affiliation(s)
- Anqi Yang
- Institute of New Type Optoelectronic Materials and Technology, College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, China
- Special and Key Laboratory of Guizhou Provincial Higher Education for Green Energy-Saving Materials, Guiyang 550025, China
| | - Jiaolian Luo
- Special and Key Laboratory of Guizhou Provincial Higher Education for Green Energy-Saving Materials, Guiyang 550025, China
- School of materials science and engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Zhenyu Xie
- Special and Key Laboratory of Guizhou Provincial Higher Education for Green Energy-Saving Materials, Guiyang 550025, China
- Architectural Engineering College, Guizhou Minzu University, Guiyang 550025, China
| | - Qian Chen
- Institute of New Type Optoelectronic Materials and Technology, College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, China
- State Key Laboratory of Public Big Data, Guizhou University, Guiyang 550025, China
| | - Quan Xie
- Institute of New Type Optoelectronic Materials and Technology, College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, China
| |
Collapse
|
28
|
Hu D, Ying W, Huo P. Resonance Enhancement of Vibrational Polariton Chemistry Obtained from the Mixed Quantum-Classical Dynamics Simulations. J Phys Chem Lett 2023; 14:11208-11216. [PMID: 38055902 PMCID: PMC10726371 DOI: 10.1021/acs.jpclett.3c02985] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 11/21/2023] [Accepted: 12/01/2023] [Indexed: 12/08/2023]
Abstract
We applied a variety of mixed quantum-classical (MQC) approaches to simulate the VSC-influenced reaction rate constant. All of these MQC simulations treat the key vibrational levels associated with the reaction coordinate in the quantum subsystem (as quantum states), whereas all other degrees of freedom (DOFs) are treated inside the classical subsystem. We find that, as long as we have the quantum state descriptions for the vibrational DOFs, one can correctly describe the VSC resonance condition when the cavity frequency matches the bond vibrational frequency. This correct resonance behavior can be obtained regardless of the detailed MQC methods that one uses. The results suggest that the MQC approaches can generate semiquantitative agreement with the exact results for rate constant changes when changing the cavity frequency, the light-matter coupling strength, or the cavity lifetime. The finding of this work suggests that one can use computationally economic MQC approaches to explore the collective coupling scenario when many molecules are collectively coupled to many cavity modes in the future.
Collapse
Affiliation(s)
- Deping Hu
- Center
for Advanced Materials Research, Beijing
Normal University, Zhuhai 519087, China
| | - Wenxiang Ying
- Department
of Chemistry, University of Rochester, Rochester, New York 14627, United States
| | - Pengfei Huo
- Department
of Chemistry, University of Rochester, Rochester, New York 14627, United States
- Institute
of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
29
|
Amati G, Mannouch JR, Richardson JO. Detailed balance in mixed quantum-classical mapping approaches. J Chem Phys 2023; 159:214114. [PMID: 38054513 DOI: 10.1063/5.0176291] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Accepted: 11/07/2023] [Indexed: 12/07/2023] Open
Abstract
The violation of detailed balance poses a serious problem for the majority of current quasiclassical methods for simulating nonadiabatic dynamics. In order to analyze the severity of the problem, we predict the long-time limits of the electronic populations according to various quasiclassical mapping approaches by applying arguments from classical ergodic theory. Our analysis confirms that regions of the mapping space that correspond to negative populations, which most mapping approaches introduce in order to go beyond the Ehrenfest approximation, pose the most serious issue for reproducing the correct thermalization behavior. This is because inverted potentials, which arise from negative electronic populations entering the nuclear force, can result in trajectories unphysically accelerating off to infinity. The recently developed mapping approach to surface hopping (MASH) provides a simple way of avoiding inverted potentials while retaining an accurate description of the dynamics. We prove that MASH, unlike any other quasiclassical approach, is guaranteed to describe the exact thermalization behavior of all quantum-classical systems, confirming it as one of the most promising methods for simulating nonadiabatic dynamics in real condensed-phase systems.
Collapse
Affiliation(s)
- Graziano Amati
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jonathan R Mannouch
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| |
Collapse
|
30
|
Huang DM, Green AT, Martens CC. A first principles derivation of energy-conserving momentum jumps in surface hopping simulations. J Chem Phys 2023; 159:214108. [PMID: 38047505 DOI: 10.1063/5.0178534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2023] [Accepted: 11/06/2023] [Indexed: 12/05/2023] Open
Abstract
The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]-along with its many later variations-forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH.
Collapse
Affiliation(s)
| | - Austin T Green
- University of California, Irvine, California 92697-2025, USA
| | - Craig C Martens
- University of California, Irvine, California 92697-2025, USA
| |
Collapse
|
31
|
Weight BM, Li X, Zhang Y. Theory and modeling of light-matter interactions in chemistry: current and future. Phys Chem Chem Phys 2023; 25:31554-31577. [PMID: 37842818 DOI: 10.1039/d3cp01415k] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2023]
Abstract
Light-matter interaction not only plays an instrumental role in characterizing materials' properties via various spectroscopic techniques but also provides a general strategy to manipulate material properties via the design of novel nanostructures. This perspective summarizes recent theoretical advances in modeling light-matter interactions in chemistry, mainly focusing on plasmon and polariton chemistry. The former utilizes the highly localized photon, plasmonic hot electrons, and local heat to drive chemical reactions. In contrast, polariton chemistry modifies the potential energy curvatures of bare electronic systems, and hence their chemistry, via forming light-matter hybrid states, so-called polaritons. The perspective starts with the basic background of light-matter interactions, molecular quantum electrodynamics theory, and the challenges of modeling light-matter interactions in chemistry. Then, the recent advances in modeling plasmon and polariton chemistry are described, and future directions toward multiscale simulations of light-matter interaction-mediated chemistry are discussed.
Collapse
Affiliation(s)
- Braden M Weight
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
- Department of Physics and Astronomy, University of Rochester, Rochester, NY, 14627, USA
| | - Xinyang Li
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| |
Collapse
|
32
|
Adrion DM, Karunaratne WV, Lopez SA. Multiconfigurational photodynamics simulations reveal the mechanism of photodecarbonylations of cyclopropenones in explicit aqueous environments. Chem Sci 2023; 14:13205-13218. [PMID: 38023495 PMCID: PMC10664470 DOI: 10.1039/d3sc03805j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 11/02/2023] [Indexed: 12/01/2023] Open
Abstract
Gas-evolving photochemical reactions use light and mild conditions to access strained organic compounds irreversibly. Cyclopropenones are a class of light-responsive molecules used in bioorthogonal photoclick reactions; their excited-state decarbonylation reaction mechanisms are misunderstood due to their ultrafast (<100 femtosecond) lifetimes. We have combined multiconfigurational quantum mechanical (QM) calculations and non-adiabatic molecular dynamics (NAMD) simulations to uncover the excited-state mechanism of cyclopropenone and a photoprotected cyclooctyne-(COT)-precursor in gaseous and explicit aqueous environments. We explore the role of H-bonding with fully quantum mechanical explicitly solvated NAMD simulations for the decarbonylation reaction. The cyclopropenones pass through asynchronous conical intersections and have dynamically concerted photodecarbonylation mechanisms. The COT-precursor has a higher quantum yield of 55% than cyclopropenone (28%) because these trajectories prefer to break a σCC bond to avoid the strained trans-cyclooctene geometries. Our solvated simulations show an increased quantum yield (58%) for the systems studied here.
Collapse
Affiliation(s)
- Daniel M Adrion
- Department of Chemistry and Chemical Biology, Northeastern University Boston Massachusetts 02115 USA
| | - Waruni V Karunaratne
- Department of Chemistry and Chemical Biology, Northeastern University Boston Massachusetts 02115 USA
| | - Steven A Lopez
- Department of Chemistry and Chemical Biology, Northeastern University Boston Massachusetts 02115 USA
| |
Collapse
|
33
|
Gu K, Lin S. Sustained Hydrogen Spillover on Pt/Cu(111) Single-Atom Alloy: Dynamic Insights into Gas-Induced Chemical Processes. Angew Chem Int Ed Engl 2023; 62:e202312796. [PMID: 37830406 DOI: 10.1002/anie.202312796] [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: 08/30/2023] [Revised: 10/01/2023] [Accepted: 10/13/2023] [Indexed: 10/14/2023]
Abstract
Hydrogen spillover, involving the surface migration of dissociated hydrogen atoms from active metal sites to the relatively inert catalyst support, plays a crucial role in hydrogen-involved catalytic processes. However, a comprehensive understanding of how H atoms are driven to spill over from active sites onto the catalyst support is still lacking. Here, we examine the atomic-scale perspective of the H spillover process on a Pt/Cu(111) single atom alloy surface using machine-learning accelerated molecular dynamics calculations based on density functional theory. Our results show that when an impinging H2 dissociates at an active Pt site, the Pt atom undergoes deactivation due to the dissociated hydrogen atoms that attach to it. Interestingly, collisions between H2 and sticking H atoms facilitate H spillover onto the host Cu, leading to the reactivation of the Pt atom and the realization of a continuous H spillover process. This work underscores the importance of the interaction between gas molecules and adsorbates as a driving force in elucidating chemical processes under a gaseous atmosphere, which has so far been underappreciated in thermodynamic studies.
Collapse
Affiliation(s)
- Kaixuan Gu
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou, 350002, China
| | - Sen Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou, 350002, China
| |
Collapse
|
34
|
Talbot JJ, Arias-Martinez JE, Cotton SJ, Head-Gordon M. Fantastical excited state optimized structures and where to find them. J Chem Phys 2023; 159:171102. [PMID: 37916588 DOI: 10.1063/5.0172015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2023] [Accepted: 10/16/2023] [Indexed: 11/03/2023] Open
Abstract
The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.
Collapse
Affiliation(s)
- Justin J Talbot
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
| | - Juan E Arias-Martinez
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Stephen J Cotton
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
35
|
Mutsuji A, Saita K, Maeda S. An energy decomposition and extrapolation scheme for evaluating electron transfer rate constants: a case study on electron self-exchange reactions of transition metal complexes. RSC Adv 2023; 13:32097-32103. [PMID: 37920761 PMCID: PMC10619204 DOI: 10.1039/d3ra05784d] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Accepted: 10/23/2023] [Indexed: 11/04/2023] Open
Abstract
A simple approach to the analysis of electron transfer (ET) reactions based on energy decomposition and extrapolation schemes is proposed. The present energy decomposition and extrapolation-based electron localization (EDEEL) method represents the diabatic energies for the initial and final states using the adiabatic energies of the donor and acceptor species and their complex. A scheme for the efficient estimation of ET rate constants is also proposed. EDEEL is semi-quantitative by directly evaluating the seam-of-crossing region of two diabatic potentials. In a numerical test, EDEEL successfully provided ET rate constants for electron self-exchange reactions of thirteen transition metal complexes with reasonable accuracy. In addition, its energy decomposition and extrapolation schemes provide all the energy values required for activation-strain model (ASM) analysis. The ASM analysis using EDEEL provided rational interpretations of the variation of the ET rate constants as a function of the transition metal complexes. These results suggest that EDEEL is useful for efficiently evaluating ET rate constants and obtaining a rational understanding of their magnitudes.
Collapse
Affiliation(s)
- Akihiro Mutsuji
- Graduate School of Chemical Sciences and Engineering, Hokkaido University Sapporo Hokkaido 060-8628 Japan
| | - Kenichiro Saita
- Department of Chemistry, Graduate School of Science, Hokkaido University Sapporo Hokkaido 060-0810 Japan
| | - Satoshi Maeda
- Department of Chemistry, Graduate School of Science, Hokkaido University Sapporo Hokkaido 060-0810 Japan
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University Sapporo Hokkaido 001-0021 Japan
- ERATO Maeda Artificial Intelligence for Chemical Reaction Design and Discovery Project, Hokkaido University Sapporo Hokkaido 060-0810 Japan
- Research and Services Division of Materials Data and Integrated System (MaDIS), National Institute for Materials Science (NIMS) Tsukuba Ibaraki 305-0044 Japan
| |
Collapse
|
36
|
Zhu Y, Prezhdo OV, Long R, Fang WH. Twist Angle-Dependent Intervalley Charge Carrier Transfer and Recombination in Bilayer WS 2. J Am Chem Soc 2023; 145:22826-22835. [PMID: 37796526 DOI: 10.1021/jacs.3c09170] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
A twist angle at a van der Waals junction provides a handle to tune its optoelectronic properties for a variety of applications, and a comprehensive understanding of how the twist modulates electronic structure, interlayer coupling, and carrier dynamics is needed. We employ time-dependent density functional theory and nonadiabatic molecular dynamics to elucidate angle-dependent intervalley carrier transfer and recombination in bilayer WS2. Repulsion between S atoms in twisted configurations weakens interlayer coupling, increases the interlayer distance, and softens layer breathing modes. Twisting has a minor influence on K valleys while it lowers Γ valleys and raises Q valleys because their wave functions are delocalized between layers. Consequently, the reduced energy gaps between the K and Γ valleys accelerate the hole transfer in the twisted structures. Intervalley electron transfer proceeds nearly an order of magnitude faster than hole transfer. The more localized wave functions at K than Q values and larger bandgaps result in smaller nonadiabatic couplings for intervalley recombination, making it 3-4 times slower in twisted than high-symmetry structures. B2g breathing, E2g in-plane, and A1g out-of-plane modes are most active during intervalley carrier transfer and recombination. The faster intervalley transfer and extended carrier lifetimes in twisted junctions are favorable for optoelectronic device performance.
Collapse
Affiliation(s)
- Yonghao Zhu
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P.R. China
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P.R. China
| | - Wei-Hai Fang
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P.R. China
| |
Collapse
|
37
|
Zhao X, Merritt ICD, Lei R, Shu Y, Jacquemin D, Zhang L, Xu X, Vacher M, Truhlar DG. Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation. J Chem Theory Comput 2023; 19:6577-6588. [PMID: 37772732 DOI: 10.1021/acs.jctc.3c00813] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/30/2023]
Abstract
Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.
Collapse
Affiliation(s)
- Xiaorui Zhao
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, P. R. China
| | | | - Ruiqing Lei
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
| | - Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
- Institut Universitaire de France (IUF), Paris 75005, France
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Xuefei Xu
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, P. R. China
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
38
|
Park W, Komarov K, Lee S, Choi CH. Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory. J Phys Chem Lett 2023; 14:8896-8908. [PMID: 37767969 PMCID: PMC10561896 DOI: 10.1021/acs.jpclett.3c02296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Accepted: 09/21/2023] [Indexed: 09/29/2023]
Abstract
The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.
Collapse
Affiliation(s)
- Woojin Park
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Konstantin Komarov
- Center
for Quantum Dynamics, Pohang University
of Science and Technology, Pohang 37673, South Korea
| | - Seunghoon Lee
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Cheol Ho Choi
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
39
|
Kim J, Woo KC, Kang M, Kim SK. Dynamic Role of the Intramolecular Hydrogen Bonding in the S 1 State Relaxation Dynamics Revealed by the Direct Measurement of the Mode-Dependent Internal Conversion Rate of 2-Chlorophenol and 2-Chlorothiophenol. J Phys Chem Lett 2023; 14:8428-8436. [PMID: 37712655 DOI: 10.1021/acs.jpclett.3c02208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/16/2023]
Abstract
The dynamic role of the intramolecular hydrogen bond in the S1 relaxation of cis-2-chlorophenol (2-CP) or cis-2-chlorothiophenol (2-CTP) has been investigated in a state-specific manner. Whereas ultrafast internal conversion is dominant for 2-CP, the H-tunneling competes with internal conversion for 2-CTP even at the S1 origin. The S0-S1 internal conversion rate of 2-CTP could be directly measured from the S1 lifetimes of 2-CTP-d1 (Cl-C6H4-SD) as the D-tunneling is kinetically blocked, allowing distinct estimations of tunneling and internal conversion rates with increasing the energy. The internal conversion rate of 2-CTP increases by two times at the out-of-plane torsional mode excitation, suggesting that the internal conversion is facilitated at the nonplanar geometry. It then sharply increases at ∼600 cm-1, indicating that the S1/S0 conical intersection is readily accessible at the extended C-Cl bond length. The strength of the intramolecular hydrogen bond should be responsible for the distinct dynamic behaviors of 2-CP and 2-CTP.
Collapse
Affiliation(s)
- Junggil Kim
- Department of Chemistry, KAIST, Daejeon 34141, Republic of Korea
| | - Kyung Chul Woo
- Department of Chemistry, KAIST, Daejeon 34141, Republic of Korea
| | - Minseok Kang
- Department of Chemistry, KAIST, Daejeon 34141, Republic of Korea
| | - Sang Kyu Kim
- Department of Chemistry, KAIST, Daejeon 34141, Republic of Korea
| |
Collapse
|
40
|
Shao C, Shi Z, Xu J, Wang L. Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics. J Phys Chem Lett 2023; 14:7680-7689. [PMID: 37606199 DOI: 10.1021/acs.jpclett.3c02019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
41
|
Mandal A, Taylor MA, Weight BM, Koessler ER, Li X, Huo P. Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics. Chem Rev 2023; 123:9786-9879. [PMID: 37552606 PMCID: PMC10450711 DOI: 10.1021/acs.chemrev.2c00855] [Citation(s) in RCA: 31] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Indexed: 08/10/2023]
Abstract
When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.
Collapse
Affiliation(s)
- Arkajit Mandal
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Department
of Chemistry, Columbia University, New York, New York 10027, United States
| | - Michael A.D. Taylor
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| | - Braden M. Weight
- Department
of Physics and Astronomy, University of
Rochester, Rochester, New York 14627, United
States
| | - Eric R. Koessler
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Xinyang Li
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Pengfei Huo
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
42
|
Gardner J, Habershon S, Maurer RJ. Assessing Mixed Quantum-Classical Molecular Dynamics Methods for Nonadiabatic Dynamics of Molecules on Metal Surfaces. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:15257-15270. [PMID: 37583439 PMCID: PMC10424245 DOI: 10.1021/acs.jpcc.3c03591] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/27/2023] [Revised: 06/30/2023] [Indexed: 08/17/2023]
Abstract
Mixed quantum-classical (MQC) methods for simulating the dynamics of molecules at metal surfaces have the potential to accurately and efficiently provide mechanistic insight into reactive processes. Here, we introduce simple two-dimensional models for the scattering of diatomic molecules at metal surfaces based on recently published electronic structure data. We apply several MQC methods to investigate their ability to capture how nonadiabatic effects influence molecule-metal energy transfer during the scattering process. Specifically, we compare molecular dynamics with electronic friction, Ehrenfest dynamics, independent electron surface hopping, and the broadened classical master equation approach. In the case of independent electron surface hopping, we implement a simple decoherence correction approach and assess its impact on vibrationally inelastic scattering. Our results show that simple, low-dimensional models can be used to qualitatively capture experimentally observed vibrational energy transfer and provide insight into the relative performance of different MQC schemes. We observe that all approaches predict similar kinetic energy dependence but return different vibrational energy distributions. Finally, by varying the molecule-metal coupling, we can assess the coupling regime in which some MQC methods become unsuitable.
Collapse
Affiliation(s)
- James Gardner
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Reinhard J. Maurer
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
- Department
of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
43
|
Lassmann Y, Curchod BFE. Probing the sensitivity of ab initio multiple spawning to its parameters. Theor Chem Acc 2023; 142:66. [PMID: 37520272 PMCID: PMC10382418 DOI: 10.1007/s00214-023-03004-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 07/13/2023] [Indexed: 08/01/2023]
Abstract
Full multiple spawning (FMS) offers a strategy to simulate the nonadiabatic dynamics of molecular systems by describing their nuclear wavefunctions by a linear combination of coupled trajectory basis functions (TBFs). Applying a series of controlled approximations to the full multiple spawning (FMS) equations leads to the ab initio multiple spawning (AIMS), which is compatible with an on-the-fly propagation of the TBFs and an accurate description of nonadiabatic processes. The AIMS strategy and its numerical implementations, however, rely on a series of user-defined parameters. Herein, we investigate the influence of these parameters on the electronic-state population of two molecular systems- trans-azomethane and a two-dimensional model of the butatriene cation. This work highlights the stability of AIMS with respect to most of its parameters, underlines the specific parameters that require particular attention from the user of the method, and offers prescriptions for an informed selection of their value. Supplementary Information The online version contains supplementary material available at 10.1007/s00214-023-03004-w.
Collapse
Affiliation(s)
- Yorick Lassmann
- Centre for Computational Chemistry, School of Chemistry, Cantock’s Close, University of Bristol, Bristol, BS8 1TS UK
| | - Basile F. E. Curchod
- Centre for Computational Chemistry, School of Chemistry, Cantock’s Close, University of Bristol, Bristol, BS8 1TS UK
| |
Collapse
|
44
|
Akher FB, Shu Y, Varga Z, Truhlar DG. Semiclassical Multistate Dynamics for Six Coupled 5A' States of O + O 2. J Chem Theory Comput 2023. [PMID: 37441750 DOI: 10.1021/acs.jctc.3c00517] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2023]
Abstract
Dynamics simulations of high-energy O2-O collisions play an important role in simulating thermal energy content and heat flux in flows around hypersonic vehicles. To carry out such dynamics simulations efficiently requires accurate global potential energy surfaces and (in most algorithms) state couplings for many energetically accessible electronic states. The ability to treat collisions involving many coupled electronic states has been a challenge for decades. Very recently, a new diabatization method, the parametrically managed diabatization by deep neural network (PM-DDNN), has been developed. The PM-DDNN method uses a deep neural network architecture with an activation function parametrically dependent on input data to discover and fit the diabatic potential energy matrix (DPEM) as a function of geometry, and the adiabatic potential energy surfaces are obtained by diagonalization of a small matrix with analytic matrix elements. Here, we applied the PM-DDNN method to the six lowest-energy potential energy surfaces in the 5A' manifold of O3 to perform simultaneous diabatization and fitting; the data are obtained by extended multistate complete-active-space second-order perturbation theory. We then used the adiabatic surfaces for dynamics calculations with three methods: coherent switching with decay of mixing (CSDM), curvature-driven CSDM (κCSDM), and electronically curvature-driven CSDM (eκCSDM). The κCSDM calculations require only adiabatic potential energies and gradients. The three dynamical methods are in good agreement. We then calculated electronically nonadiabatic, electronically inelastic, and dissociative cross sections for seven initial collision energies, five initial vibrational levels, and four initial rotational levels. Trends in the electronically inelastic cross sections as functions of the initial collision energy and vibrational level were rationalized in terms of the coordinate ranges where the gaps between the second and third potential energy surfaces are small.
Collapse
Affiliation(s)
- Farideh Badichi Akher
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Yinan Shu
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- 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
|
45
|
Lu D, Galvão BRL, Varandas AJC, Guo H. Quantum and semiclassical studies of nonadiabatic electronic transitions between N( 4S) and N( 2D) by collisions with N 2. Phys Chem Chem Phys 2023; 25:15656-15665. [PMID: 37278325 DOI: 10.1039/d3cp01429k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The dynamics and kinetics of spin-forbidden transitions between N(2D) and N(4S) via collisions with N2 molecules are investigated using a quantum wave packet (WP) method and the semi-classical coherent switches with decay of mixing (CSDM) method. These electronic transition processes are competing with exchange reaction channels on both the doublet and quartet potential energy surfaces. The WP and CSDM quenching rate coefficients are found in reasonable agreement with each other, and both reproduce the previous theoretical results. For the excitation process, the agreement between the two approaches is dependent on the treatment of the zero-point energy (ZPE) in the product, because the high endoergicity of this process leads to severe violation of the vibrational ZPE. The Gaussian-binning (GB) method is found to improve the agreement with the quantum result. The excitation rate coefficients are found to be two orders of magnitude smaller than that of the adiabatic exchange reaction, underscoring the inefficient intersystem crossing due to the weak spin-orbit coupling between the two spin manifolds of the N3 system.
Collapse
Affiliation(s)
- Dandan Lu
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, 87131, New Mexico, USA.
| | - Breno R L Galvão
- Centro Federal de Educação Tecnológica de Minas Gerais, CEFET-MG, Av. Amazonas 5253, (30421-169), Belo Horizonte, Minas Gerais, Brazil
| | - Antonio J C Varandas
- Departamento de Física, Universidade Federal do Espírito Santo, Vitória, Brazil
- Coimbra Chemistry Centre and Chemistry Department, University of Coimbra, Coimbra, Portugal
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, 87131, New Mexico, USA.
| |
Collapse
|
46
|
Abstract
The topological properties of an object, associated with an integer called the topological invariant, are global features that cannot change continuously but only through abrupt variations, hence granting them intrinsic robustness. Engineered metamaterials (MMs) can be tailored to support highly nontrivial topological properties of their band structure, relative to their electronic, electromagnetic, acoustic and mechanical response, representing one of the major breakthroughs in physics over the past decade. Here, we review the foundations and the latest advances of topological photonic and phononic MMs, whose nontrivial wave interactions have become of great interest to a broad range of science disciplines, such as classical and quantum chemistry. We first introduce the basic concepts, including the notion of topological charge and geometric phase. We then discuss the topology of natural electronic materials, before reviewing their photonic/phononic topological MM analogues, including 2D topological MMs with and without time-reversal symmetry, Floquet topological insulators, 3D, higher-order, non-Hermitian and nonlinear topological MMs. We also discuss the topological aspects of scattering anomalies, chemical reactions and polaritons. This work aims at connecting the recent advances of topological concepts throughout a broad range of scientific areas and it highlights opportunities offered by topological MMs for the chemistry community and beyond.
Collapse
Affiliation(s)
- Xiang Ni
- Photonics Initiative, Advanced Science Research Center, City University of New York, New York, New York 10031, United States
- School of Physics and Electronics, Central South University, Changsha, Hunan 410083, China
| | - Simon Yves
- Photonics Initiative, Advanced Science Research Center, City University of New York, New York, New York 10031, United States
| | - Alex Krasnok
- Department of Electrical and Computer Engineering, Florida International University, Miami, Florida 33174, USA
| | - Andrea Alù
- Photonics Initiative, Advanced Science Research Center, City University of New York, New York, New York 10031, United States
- Department of Electrical Engineering, City College, The City University of New York, 160 Convent Avenue, New York, New York 10031, United States
- Physics Program, The Graduate Center, The City University of New York, 365 Fifth Avenue, New York, New York 10016, United States
| |
Collapse
|
47
|
Chen WK, Wang SR, Liu XY, Fang WH, Cui G. Nonadiabatic Derivative Couplings Calculated Using Information of Potential Energy Surfaces without Wavefunctions: Ab Initio and Machine Learning Implementations. Molecules 2023; 28:molecules28104222. [PMID: 37241962 DOI: 10.3390/molecules28104222] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Revised: 05/16/2023] [Accepted: 05/18/2023] [Indexed: 05/28/2023] Open
Abstract
In this work, we implemented an approximate algorithm for calculating nonadiabatic coupling matrix elements (NACMEs) of a polyatomic system with ab initio methods and machine learning (ML) models. Utilizing this algorithm, one can calculate NACMEs using only the information of potential energy surfaces (PESs), i.e., energies, and gradients as well as Hessian matrix elements. We used a realistic system, namely CH2NH, to compare NACMEs calculated by this approximate PES-based algorithm and the accurate wavefunction-based algorithm. Our results show that this approximate PES-based algorithm can give very accurate results comparable to the wavefunction-based algorithm except at energetically degenerate points, i.e., conical intersections. We also tested a machine learning (ML)-trained model with this approximate PES-based algorithm, which also supplied similarly accurate NACMEs but more efficiently. The advantage of this PES-based algorithm is its significant potential to combine with electronic structure methods that do not implement wavefunction-based algorithms, low-scaling energy-based fragment methods, etc., and in particular efficient ML models, to compute NACMEs. The present work could encourage further research on nonadiabatic processes of large systems simulated by ab initio nonadiabatic dynamics simulation methods in which NACMEs are always required.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Hebei Key Laboratory of Inorganic Nano-Materials, College of Chemistry and Materials Science, Hebei Normal University, Shijiazhuang 050024, China
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Sheng-Rui Wang
- 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
- Hefei National Laboratory, Hefei 230088, 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
| |
Collapse
|
48
|
Troß J, Carter-Fenk K, Cole-Filipiak NC, Schrader P, Word M, McCaslin LM, Head-Gordon M, Ramasesha K. Excited-State Dynamics during Primary C-I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy. J Phys Chem A 2023; 127:4103-4114. [PMID: 37103479 DOI: 10.1021/acs.jpca.3c01414] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2023]
Abstract
In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to nπC═O* states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed nπC═O* and nσC-I* character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin-orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C-I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N4,5 edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C-I homolysis. By examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C-I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental-theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin-orbit coupled system.
Collapse
Affiliation(s)
- Jan Troß
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Neil C Cole-Filipiak
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Paul Schrader
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Mi'Kayla Word
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Laura M McCaslin
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Krupa Ramasesha
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| |
Collapse
|
49
|
Mi X, Zhang M, Zhang L, Wu C, Zhou T, Xu H, Xie C, Li Z, Liu Y. Geometric Phase Effect in Attosecond Stimulated X-ray Raman Spectroscopy. J Phys Chem A 2023; 127:3608-3613. [PMID: 37053512 DOI: 10.1021/acs.jpca.3c00594] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Conical intersections (CIs) are diabolical points in the potential energy surfaces generally caused by point-wise degeneracy of different electronic states, and give rise to the geometric phases (GPs) of molecular wave functions. Here we theoretically propose and demonstrate that the transient redistribution of ultrafast electronic coherence in attosecond Raman signal (TRUECARS) spectroscopy is capable of detecting the GP effect in excited state molecules by applying two probe pulses including an attosecond and a femtosecond X-ray pulse. The mechanism is based on a set of symmetry selection rules in the presence of nontrivial GPs. The model of this work can be realized for probing the geometric phase effect in the excited state dynamics of complex molecules with appropriate symmetries, using attosecond light sources such as free-electron X-ray lasers.
Collapse
Affiliation(s)
- Xiaoyu Mi
- State Key Laboratory for Mesoscopic Physics and Collaborative Innovation Center of Quantum Matter, School of Physics, Peking University, Beijing 100871, China
| | - Ming Zhang
- State Key Laboratory for Mesoscopic Physics and Collaborative Innovation Center of Quantum Matter, School of Physics, Peking University, Beijing 100871, China
| | - Linfeng Zhang
- State Key Laboratory for Mesoscopic Physics and Collaborative Innovation Center of Quantum Matter, School of Physics, Peking University, Beijing 100871, China
| | - Chengyin Wu
- State Key Laboratory for Mesoscopic Physics and Collaborative Innovation Center of Quantum Matter, School of Physics, Peking University, Beijing 100871, China
- Peking University Yangtze Delta Institute of Optoelectronics, Nantong, Jiangsu 226000 China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Tianyu Zhou
- State Key Laboratory for Mesoscopic Physics and Collaborative Innovation Center of Quantum Matter, School of Physics, Peking University, Beijing 100871, China
| | - Haitan Xu
- School of Materials Science and Intelligent Engineering, Nanjing University, Suzhou, Jiangsu 215163, China
- School of Physical Sciences, University of Science and Technology of China, Hefei 230026, China
| | - Changjian Xie
- Institute of Modern Physics, Shanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an 710127, China
| | - Zheng Li
- State Key Laboratory for Mesoscopic Physics and Collaborative Innovation Center of Quantum Matter, School of Physics, Peking University, Beijing 100871, China
- Peking University Yangtze Delta Institute of Optoelectronics, Nantong, Jiangsu 226000 China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Yunquan Liu
- State Key Laboratory for Mesoscopic Physics and Collaborative Innovation Center of Quantum Matter, School of Physics, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
50
|
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
|