1
|
Han D, Akimov AV. Nonadiabatic Dynamics with Exact Factorization: Implementation and Assessment. J Chem Theory Comput 2024; 20:5022-5042. [PMID: 38837952 DOI: 10.1021/acs.jctc.4c00343] [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/2024]
Abstract
In this work, we report our implementation of several independent-trajectory mixed-quantum-classical (ITMQC) nonadiabatic dynamics methods based on exact factorization (XF) in the Libra package for nonadiabatic and excited-state dynamics. Namely, the exact factorization surface hopping (SHXF), mixed quantum-classical dynamics (MQCXF), and mean-field (MFXF) are introduced. Performance of these methods is compared to that of several traditional surface hopping schemes, such as the fewest-switches surface hopping (FSSH), branching-corrected surface hopping (BCSH), and the simplified decay of mixing (SDM), as well as conventional Ehrenfest (mean-field, MF) method. Based on a comprehensive set of 1D model Hamiltonians, we find the ranking SHXF ≈ MQCXF > BCSH > SDM > FSSH ≫ MF, with the BCSH sometimes outperforming the XF methods in terms of describing coherences. Although the MFXF method can yield reasonable populations and coherences for some cases, it does not conserve the total energy and is therefore not recommended. We also find that the branching correction for auxiliary trajectories is important for the XF methods to yield accurate populations and coherences. However, the branching correction can worsen the quality of the energy conservation in the MQCXF. Finally, we find that using the time-dependent Gaussian width approximation used in the XF methods for computing decoherence correction can improve the quality of energy conservation in the MQCXF dynamics. The parameter-free scheme of Subotnik for computing the Gaussian widths is found to deliver the best performance in situations where such widths are not known a priori.
Collapse
Affiliation(s)
- Daeho Han
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
2
|
Ghosh A, Pramanik A, Pal S, Sarkar P. Emergence of Z-Scheme Photocatalysis for Total Water Splitting: An Improvised Route to High Efficiency. J Phys Chem Lett 2024:6841-6851. [PMID: 38917061 DOI: 10.1021/acs.jpclett.4c01162] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/27/2024]
Abstract
Photocatalytic water splitting to spontaneously produce H2 and O2 is a long-standing goal in solar energy conversion, presenting a significant challenge without using sacrificial electron donors or external biases. Inspired by natural photosynthesis, the design of artificial Z-scheme photocatalytic systems is at the forefront of this field. These systems achieve higher redox potential by separating photogenerated electrons and holes through a fast interlayer recombination process between valence and conduction band edges. Z-scheme photocatalysis involves using two different semiconductors with distinct bandgap energies. Here, we explore potential systems based on two-dimensional (2D) heterostructures composed of carbon, nitrogen, or similar main group elements. The advantages and disadvantages of these systems are discussed, with a focus on enhancing their efficiency through strategic design. Special emphasis is placed on the dynamics of excited charge carrier transfer and recombination processes, which are crucial for developing efficient photocatalytic systems for overall water splitting.
Collapse
Affiliation(s)
- Atish Ghosh
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| | - Anup Pramanik
- Department of Chemistry, Sidho-Kanho-Birsha University, Purulia 723104, India
| | - Sougata Pal
- Department of Chemistry, University of Gour Banga, Malda 732103, India
| | - Pranab Sarkar
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| |
Collapse
|
3
|
Huang H, Peng J, Zhang Y, Gu FL, Lan Z, Xu C. The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications. J Chem Phys 2024; 160:234101. [PMID: 38884395 DOI: 10.1063/5.0215036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/15/2024] [Indexed: 06/18/2024] Open
Abstract
Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.
Collapse
Affiliation(s)
- Haiyi Huang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yulin Zhang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
4
|
Zhang Q, Shao X, Li W, Mi W, Pavanello M, Akimov AV. Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:385901. [PMID: 38866023 DOI: 10.1088/1361-648x/ad577d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Accepted: 06/12/2024] [Indexed: 06/14/2024]
Abstract
In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.
Collapse
Affiliation(s)
- Qingxin Zhang
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| | - Xuecheng Shao
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China
| | - Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, People's Republic of China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People's Republic of China
| | - Michele Pavanello
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| |
Collapse
|
5
|
Sanz García J, Maskri R, Mitrushchenkov A, Joubert-Doriol L. Optimizing Conical Intersections without Explicit Use of Non-Adiabatic Couplings. J Chem Theory Comput 2024. [PMID: 38888629 DOI: 10.1021/acs.jctc.4c00326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
We present two alternative methods for optimizing minimum energy conical intersection (MECI) molecular geometries without knowledge of the derivative coupling (DC). These methods are based on the utilization of Lagrange multipliers: (i) one method uses an approximate calculation of the DC, while the other (ii) do not require the DC. Both methods use the fact that information on the DC is contained in the Hessian of the squared energy difference. Tests done on a set of small molecular systems, in comparison with other methods, show the ability of the proposed methods to optimize MECIs. Finally, we apply the methods to the furimamide molecule, to optimize and characterize its S1/S2 MECI, and to optimizing the S0/S1 MECI of the silver trimer.
Collapse
Affiliation(s)
- Juan Sanz García
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Rosa Maskri
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Alexander Mitrushchenkov
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Loïc Joubert-Doriol
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| |
Collapse
|
6
|
Polonius S, Lehrner D, González L, Mai S. Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations. J Chem Theory Comput 2024; 20:4738-4750. [PMID: 38768386 PMCID: PMC11171268 DOI: 10.1021/acs.jctc.4c00169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 04/22/2024] [Accepted: 05/07/2024] [Indexed: 05/22/2024]
Abstract
Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - David Lehrner
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
7
|
Shu Y, Truhlar DG. Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients. J Chem Theory Comput 2024; 20:4396-4426. [PMID: 38819014 DOI: 10.1021/acs.jctc.4c00424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
8
|
Chen J, Wang Y, Dou W. Floquet nonadiabatic mixed quantum-classical dynamics in periodically driven solid systems. J Chem Phys 2024; 160:214101. [PMID: 38828807 DOI: 10.1063/5.0204158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2024] [Accepted: 05/16/2024] [Indexed: 06/05/2024] Open
Abstract
In this paper, we introduce the Floquet mean-field dynamics and Floquet surface hopping approaches to study the nonadiabatic dynamics in periodically driven solid systems. We demonstrate that these two approaches can be formulated in both real and reciprocal spaces. Using the two approaches, we are able to simulate the interaction between electronic carriers and phonons under periodic drivings, such as strong light-matter interactions. Employing the Holstein and Peierls models, we show that strong light-matter interactions can effectively modulate the dynamics of electronic population and mobility. Notably, our study demonstrates the feasibility and effectiveness of modeling low-momentum carriers' interactions with phonons using a truncated reciprocal space basis, an approach impractical in real space frameworks. Moreover, we reveal that even with a significant truncation, carrier populations derived from surface hopping maintain greater accuracy compared to those obtained via mean-field dynamics. These results underscore the potential of our proposed methods in advancing the understanding of carrier-phonon interactions in various periodically driven materials.
Collapse
Affiliation(s)
- Jingqi Chen
- Fudan University, 220 Handan Road, Shanghai 200433, China
- Department of Chemistry, School of Science, Westlake University, Hangzhou 310024, Zhejiang, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou 310024, Zhejiang, China
| | - Yu Wang
- Department of Chemistry, School of Science, Westlake University, Hangzhou 310024, Zhejiang, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou 310024, Zhejiang, China
| | - Wenjie Dou
- Department of Chemistry, School of Science, Westlake University, Hangzhou 310024, Zhejiang, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou 310024, Zhejiang, China
- Department of Physics, School of Science, Westlake University, Hangzhou 310024, Zhejiang, China
| |
Collapse
|
9
|
Fu CF, Zheng Q, Li X, Yang J. Vertical Dipole Dominates Charge Carrier Lifetime in Monolayer Janus MoSSe. NANO LETTERS 2024; 24:6425-6432. [PMID: 38747348 DOI: 10.1021/acs.nanolett.4c01577] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
Two-dimensional semiconductor materials with vertical dipoles are promising photocatalysts as vertical dipoles not only promote the electron-hole separation but also enhance the carrier redox ability. However, the influence of vertical dipoles on carrier recombination in such materials, especially the competing relationship between vertical dipoles and band gaps, is not yet clear. Herein, first-principles calculations and nonadiabatic molecular dynamics simulations were combined to clarify the influence of band gap and vertical dipole on the carrier lifetime in Janus MoSSe monolayer. By comparing with the results of MoS2 and MoSe2 as well as exploring the carrier lifetime of MoSSe under strain regulation, it has been demonstrated that the vertical dipole, rather than the band gap, is the dominant factor affecting the carrier lifetime. Strikingly, a linear relationship between the carrier lifetime and vertical dipole is revealed. These findings have important implications for the design of high-performance photocatalysts and optoelectronic devices.
Collapse
Affiliation(s)
- Cen-Feng Fu
- Department of Chemical Physics and Hefei National Research Center for Physical Science at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Qijing Zheng
- Department of Physics, and ICQD/Hefei National Research Center for Physical Science at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| | - Xingxing Li
- Department of Chemical Physics and Hefei National Research Center for Physical Science at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| | - Jinlong Yang
- Department of Chemical Physics and Hefei National Research Center for Physical Science at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| |
Collapse
|
10
|
Sindhu A, Jain A. Electronic energy transfer in molecular wire: Coherences in the presence of anharmonicity. J Chem Phys 2024; 160:204117. [PMID: 38814012 DOI: 10.1063/5.0196606] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 05/06/2024] [Indexed: 05/31/2024] Open
Abstract
Electronic energy transfer in molecular wires is usually theoretically investigated with a harmonic bath to model the environment. The present study is a continuation of our previous work [A. Sindhu and A. Jain, Chem. Phys. Chem. 23, e2022003 (2022)] on studying the dynamics of molecular wires using surface hopping simulations. We extend our study to a 7-site model Hamiltonian and investigate the effects of an anharmonic bath on coherent energy transfer in molecular wires. We show that oscillatory and coherent population dynamics remain intact even in the presence of the anharmonic bath and further highlight the multiple channels available for energy flow in molecular wires.
Collapse
Affiliation(s)
- Aarti Sindhu
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Amber Jain
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| |
Collapse
|
11
|
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
|
12
|
Chen CG, Amadei A, D'Abramo M. Modeling the temperature dependence of the fluorescence properties of Indole in aqueous solution. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 313:124096. [PMID: 38442616 DOI: 10.1016/j.saa.2024.124096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 02/05/2024] [Accepted: 02/27/2024] [Indexed: 03/07/2024]
Abstract
In a recent paper, we proposed a scheme to describe the relaxation mechanism of the excited Indole in aqueous solution, involving the fluctuations among the diabatic electronic states 1Lb, 1La and 1πσ∗. Such a theoretical and computational model reproduced accurately the available experimental data at room temperature. Following these results, in the present work, we model the complex temperature dependence of the fluorescence properties of Indole in aqueous solution, with results further validating the proposed relaxation scheme. This scheme is able to explain the temperature effects on the fluorescence behavior indicating the water fluctuations as the main cause of (i) the stabilization of the dark state (1πσ∗) and (ii) the increase in temperature of the kinetics of the irreversible transition towards such a state.
Collapse
Affiliation(s)
- Cheng Giuseppe Chen
- Department of Chemistry, Sapienza University of Rome, Piazzale Aldo Moro, 5, Rome, 00185, Italy
| | - Andrea Amadei
- Department of Technological and Chemical Sciences, Tor Vergata University of Rome, Via della Ricerca Scientifica, 1, Rome, 00133, Italy.
| | - Marco D'Abramo
- Department of Chemistry, Sapienza University of Rome, Piazzale Aldo Moro, 5, Rome, 00185, Italy.
| |
Collapse
|
13
|
Kang M, Nuomin H, Chowdhury SN, Yuly JL, Sun K, Whitlow J, Valdiviezo J, Zhang Z, Zhang P, Beratan DN, Brown KR. Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics. Nat Rev Chem 2024; 8:340-358. [PMID: 38641733 DOI: 10.1038/s41570-024-00595-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/07/2024] [Indexed: 04/21/2024]
Abstract
Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.
Collapse
Affiliation(s)
- Mingyu Kang
- Duke Quantum Center, Duke University, Durham, NC, USA.
- Department of Physics, Duke University, Durham, NC, USA.
| | - Hanggai Nuomin
- Department of Chemistry, Duke University, Durham, NC, USA
| | | | - Jonathon L Yuly
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, USA
| | - Ke Sun
- Duke Quantum Center, Duke University, Durham, NC, USA
- Department of Physics, Duke University, Durham, NC, USA
| | - Jacob Whitlow
- Duke Quantum Center, Duke University, Durham, NC, USA
- Department of Electrical and Computer Engineering, Duke University, Durham, NC, USA
| | - Jesús Valdiviezo
- Kenneth S. Pitzer Theory Center, University of California, Berkeley, CA, USA
- Department of Chemistry, University of California, Berkeley, CA, USA
- Departamento de Ciencias, Sección Química, Pontificia Universidad Católica del Perú, Lima, Peru
| | - Zhendian Zhang
- Department of Chemistry, Duke University, Durham, NC, USA
| | - Peng Zhang
- Department of Chemistry, Duke University, Durham, NC, USA
| | - David N Beratan
- Department of Physics, Duke University, Durham, NC, USA.
- Department of Chemistry, Duke University, Durham, NC, USA.
- Department of Biochemistry, Duke University, Durham, NC, USA.
| | - Kenneth R Brown
- Duke Quantum Center, Duke University, Durham, NC, USA.
- Department of Physics, Duke University, Durham, NC, USA.
- Department of Chemistry, Duke University, Durham, NC, USA.
- Department of Electrical and Computer Engineering, Duke University, Durham, NC, USA.
| |
Collapse
|
14
|
Ghosh A, Das P, Kumar S, Sarkar P. Hot carrier relaxation dynamics of an aza-covalent organic framework during photoexcitation: An insight from ab initio quantum dynamics. J Chem Phys 2024; 160:164707. [PMID: 38647311 DOI: 10.1063/5.0200834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2024] [Accepted: 04/01/2024] [Indexed: 04/25/2024] Open
Abstract
In order to develop an efficient metal-free solar energy harvester, we herein performed the electronic structure calculation, followed by the hot carrier relaxation dynamics of two dimensional (2D) aza-covalent organic framework by time domain density functional calculations in conjunction with non-adiabatic molecular dynamics (NAMD) simulation. The electronic structure calculation shows that the aza-covalent organic framework (COF) is a direct bandgap semiconductor with acute charge separation and effective optical absorption in the UV-visible region. Our study of non-adiabatic molecular dynamics simulation predicts the sufficiently prolonged electron-hole recombination process (6.8 nanoseconds) and the comparatively faster electron (22.48 ps) and hole relaxation (0.51 ps) dynamics in this two-dimensional aza-COF. According to our theoretical analysis, strong electron-phonon coupling is responsible for the rapid charge relaxation, whereas the electron-hole recombination process is slowed down by relatively weak electron-phonon coupling, relatively lower non-adiabatic coupling, and quick decoherence time. We do hope that our results of NAMD simulation on exciton relaxation dynamics will be helpful for designing photovoltaic devices based on this two dimensional aza-COF.
Collapse
Affiliation(s)
- Atish Ghosh
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| | - Priya Das
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| | - Subhash Kumar
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| | - Pranab Sarkar
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| |
Collapse
|
15
|
Feng Z, Guo W, Kong WY, Chen D, Wang S, Tantillo DJ. Analogies between photochemical reactions and ground-state post-transition-state bifurcations shed light on dynamical origins of selectivity. Nat Chem 2024; 16:615-623. [PMID: 38216753 DOI: 10.1038/s41557-023-01410-y] [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/03/2023] [Accepted: 11/27/2023] [Indexed: 01/14/2024]
Abstract
Revealing the origins of kinetic selectivity is one of the premier tasks of applied theoretical organic chemistry, and for many reactions, doing so involves comparing competing transition states. For some reactions, however, a single transition state leads directly to multiple products, in which case non-statistical dynamic effects influence selectivity control. The selectivity of photochemical reactions-where crossing between excited-state and ground-state surfaces occurs near ground-state transition structures that interconvert competing products-also should be controlled by the momentum of the reacting molecules as they return to the ground state in addition to the shape of the potential energy surfaces involved. Now, using machine-learning-assisted non-adiabatic molecular dynamics and multiconfiguration pair-density functional theory, these factors are examined for a classic photochemical reaction-the deazetization of 2,3-diazabicyclo[2.2.2]oct-2-ene-for which we demonstrate that momentum dominates the selectivity for hexadiene versus [2.2.2] bicyclohexane products.
Collapse
Affiliation(s)
- Zhitao Feng
- Department of Chemistry, University of California, Davis, Davis, CA, USA
| | - Wentao Guo
- Department of Chemistry, University of California, Davis, Davis, CA, USA
| | - Wang-Yeuk Kong
- Department of Chemistry, University of California, Davis, Davis, CA, USA
| | - Dongjie Chen
- Department of Electrical and Computer Engineering, University of California, Davis, Davis, CA, USA
| | - Shunyang Wang
- Department of Chemistry, University of California, Davis, Davis, CA, USA
| | - Dean J Tantillo
- Department of Chemistry, University of California, Davis, Davis, CA, USA.
| |
Collapse
|
16
|
Karmakar S, Thakur S, Jain A. Can classical mechanics sense conical intersection? J Chem Phys 2024; 160:124110. [PMID: 38526110 DOI: 10.1063/5.0197381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2024] [Accepted: 03/05/2024] [Indexed: 03/26/2024] Open
Abstract
Conical intersection (CI) leads to fast electronic energy transfer. However, Hamm and Stock [Phys. Rev. Lett. 109, 173201 (2012)] showed the existence of a vibrational CI and its role in vibrational energy relaxation. In this paper, we further investigate the vibrational energy relaxation using an isolated model Hamiltonian system of four vibrational modes with two distinctively different timescales (two fast modes and two slow modes). We show that the excitation of the slow modes plays a crucial role in the energy relaxation mechanism. We also analyze the system from a mixed quantum-classical (surface hopping method) and a completely classical point of view. Notably, surface hopping and even classical simulations also capture fast energy relaxation, which is a signature of CI's existence.
Collapse
Affiliation(s)
- Sourav Karmakar
- Department of Chemistry, Indian Institute of Technology Bombay, Bombay, India
| | - Saumya Thakur
- Department of Chemistry, Indian Institute of Technology Bombay, Bombay, India
| | - Amber Jain
- Department of Chemistry, Indian Institute of Technology Bombay, Bombay, India
| |
Collapse
|
17
|
Tao Z, Qiu T, Bhati M, Bian X, Duston T, Rawlinson J, Littlejohn RG, Subotnik JE. Practical phase-space electronic Hamiltonians for ab initio dynamics. J Chem Phys 2024; 160:124101. [PMID: 38526114 DOI: 10.1063/5.0192084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 02/14/2024] [Indexed: 03/26/2024] Open
Abstract
Modern electronic structure theory is built around the Born-Oppenheimer approximation and the construction of an electronic Hamiltonian Ĥel(X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the well-known theory of electron translation (Γ') and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phase-space electronic Hamiltonian that depends on both nuclear position and momentum, ĤPS(X,P). While classical Born-Oppenheimer dynamics that run along the eigensurfaces of the operator Ĥel(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of ĤPS(X,P) can better capture both nuclear and electronic properties (including the elusive electronic momentum studied by Nafie). Moreover, only the latter (as opposed to the former) conserves the total linear and angular momentum in general.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Tian Qiu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Mansi Bhati
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Titouan Duston
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Jonathan Rawlinson
- Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
18
|
Dai D, Agrawal S, Prezhdo OV, Long R. Impact of large A-site cations on electron-vibrational interactions in 2D halide perovskites: Ab initio quantum dynamics. J Chem Phys 2024; 160:114704. [PMID: 38506296 DOI: 10.1063/5.0202251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Accepted: 03/03/2024] [Indexed: 03/21/2024] Open
Abstract
Using ab initio nonadiabatic molecular dynamics, we study the effect of large A-site cations on nonradiative electron-hole recombination in two-dimensional Ruddlesden-Popper perovskites HA2APb2I7, HA = n-hexylammonium, A = methylammonium (MA), or guanidinium (GA). The steric hindrance created by large GA cations distorts and stiffens the inorganic Pb-I lattice, reduces thermal structural fluctuations, and maintains the delocalization of electrons and holes at ambient and elevated temperatures. The delocalized charges interact more strongly in the GA system than in the MA system, and the charge recombination is accelerated. In contrast, replacement of only some MA cations with GA enhances disorder and increases charge lifetime, as seen in three-dimensional perovskites. This study highlights the key influence of structural fluctuations and disorder on the properties of charge carriers in metal halide perovskites, providing guidance for tuning materials' optoelectronic performance.
Collapse
Affiliation(s)
- Dandan Dai
- College of Chemistry, Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Sraddha Agrawal
- Department of Chemistry, University of Southern California, Los Angeles, California 90007, USA
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90007, USA
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People's Republic of China
| |
Collapse
|
19
|
Kim SS, Rhee YM. Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water. Phys Chem Chem Phys 2024; 26:9021-9036. [PMID: 38440829 DOI: 10.1039/d3cp05643k] [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
Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.
Collapse
Affiliation(s)
- Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
20
|
Ghosh A, Kumar S, Sarkar P. Point defect-mediated hot carrier relaxation dynamics of lead-free FASnI 3 perovskites. NANOSCALE 2024; 16:4737-4744. [PMID: 38299671 DOI: 10.1039/d3nr04039a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
In search of a promising optoelectronic performance, we herein investigated the hot carrier relaxation dynamics of a lead-free cubic phased bulk formamidinium tin triiodide (FASnI3) perovskite. To gain detailed theoretical insights, we should estimate the carrier relaxation dynamics of this pristine perovskite. To control the dynamics, point defects like central tin (Sn), iodine(I) anions, and formamidinium (FA) cations were introduced. With the iodine vacancy in the FASnI3 perovskite, the system seems to be unstable at room temperature, whereas the other three types of FASnI3 perovskites (pristine, Sn vacancy, and FA vacancy) are significantly stable at 300 K having semiconducting nature and excellent optical absorption in the UV-visible range. The computed electron-hole recombination time for the pristine system is 3.9 nanoseconds, which is in good agreement with the experimental investigation. The exciton relaxation processes in Sn and FA vacancy perovskites require 2.8 and 4.8 nanoseconds, respectively. These variations in the hot carrier relaxation dynamics processes are caused by the generation of significant changes in non-adiabatic coupling between energy levels, electron-phonon coupling, and quantum decoherence in different point defect analogous systems. The results presented here offer deeper insight into the temperature-dependent carrier relaxation dynamics of FASnI3 perovskites and thus open up opportunities for future exploration of their optoelectronic properties.
Collapse
Affiliation(s)
- Atish Ghosh
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India.
| | - Subhash Kumar
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India.
| | - Pranab Sarkar
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India.
| |
Collapse
|
21
|
Dupuy L, Rikus A, Maitra NT. Exact-Factorization-Based Surface Hopping without Velocity Adjustment. J Phys Chem Lett 2024:2643-2649. [PMID: 38422391 DOI: 10.1021/acs.jpclett.4c00115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2024]
Abstract
While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. The resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.
Collapse
Affiliation(s)
- Lucien Dupuy
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| | - Anton Rikus
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
- University of Münster, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, 48149 Münster, Germany
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
22
|
Xu J, Carney TE, Zhou R, Shepard C, Kanai Y. Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics. J Am Chem Soc 2024; 146:5011-5029. [PMID: 38362887 DOI: 10.1021/jacs.3c08226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. In this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.
Collapse
Affiliation(s)
- Jianhang Xu
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Thomas E Carney
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
- Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
23
|
Wang L, Liu X, Wan L, Gao Y, Wang X, Liu J, Tan S, Guo Q, Zhao W, Hu W, Li Q, Yang J. Excitonic Effects of the Excited-State Photocatalytic Reaction at the Molecule/Metal Oxide Interface. J Phys Chem Lett 2024; 15:2096-2104. [PMID: 38358755 DOI: 10.1021/acs.jpclett.3c03075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2024]
Abstract
Excitonic effects caused by the Coulomb interaction between electrons and holes play a crucial role in photocatalysis at the molecule/metal oxide interface. As an ideal model for investigating the excitonic effect, coadsorption and photodissociation of water and methanol molecules on titanium dioxide involve complex ground-state thermalcatalytic and excited-state photocatalytic reaction processes. Herein, we systemically investigate the excited-state electronic structures of the coadsorption of H2O and CH3OH molecules on a rutile TiO2(110) surface by linear-response time-dependent density functional theory calculations and probe the reaction path for generating HCOOH or CO2, from ground-state and excited-state perspectives. The reaction barriers in excited-state calculations are significantly different from those in ground-state calculations during three processes, with the largest decrease being 0.94 eV for the Ti5c-O-CH2-O-Ti5c formation process.
Collapse
Affiliation(s)
- Lei Wang
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiaofeng Liu
- School of Physics, Hefei University of Technology, Hefei, Anhui 230009, China
| | - Lingyun Wan
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yunzhi Gao
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiaoning Wang
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jie Liu
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Shijing Tan
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Qing Guo
- Department of Chemistry, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Wenhui Zhao
- Department of Physics, Ningbo University, Ningbo, Zhejiang 315211, China
| | - Wei Hu
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Qunxiang Li
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jinlong Yang
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
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
|
Tao Z, Bian X, Wu Y, Rawlinson J, Littlejohn RG, Subotnik JE. Total angular momentum conservation in Ehrenfest dynamics with a truncated basis of adiabatic states. J Chem Phys 2024; 160:054104. [PMID: 38310474 DOI: 10.1063/5.0177778] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2023] [Accepted: 12/18/2023] [Indexed: 02/05/2024] Open
Abstract
We show that standard Ehrenfest dynamics does not conserve linear and angular momentum when using a basis of truncated adiabatic states. However, we also show that previously proposed effective Ehrenfest equations of motion [M. Amano and K. Takatsuka, "Quantum fluctuation of electronic wave-packet dynamics coupled with classical nuclear motions," J. Chem. Phys. 122, 084113 (2005) and V. Krishna, "Path integral formulation for quantum nonadiabatic dynamics and the mixed quantum classical limit," J. Chem. Phys. 126, 134107 (2007)] involving the non-Abelian Berry force do maintain momentum conservation. As a numerical example, we investigate the Kramers doublet of the methoxy radical using generalized Hartree-Fock with spin-orbit coupling and confirm that angular momentum is conserved with the proper equations of motion. Our work makes clear some of the limitations of the Born-Oppenheimer approximation when using ab initio electronic structure theory to treat systems with unpaired electronic spin degrees of freedom, and we demonstrate that Ehrenfest dynamics can offer much improved, qualitatively correct results.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yanze Wu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Jonathan Rawlinson
- Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
26
|
Wu B, He X, Liu J. Nonadiabatic Field on Quantum Phase Space: A Century after Ehrenfest. J Phys Chem Lett 2024; 15:644-658. [PMID: 38205956 DOI: 10.1021/acs.jpclett.3c03385] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Nonadiabatic transition dynamics lies at the core of many electron/hole transfer, photoactivated, and vacuum field-coupled processes. About a century after Ehrenfest proposed "Phasenraum" and the Ehrenfest theorem, we report a conceptually novel trajectory-based nonadiabatic dynamics approach, nonadiabatic field (NAF), based on a generalized exact coordinate-momentum phase space formulation of quantum mechanics. It does not employ the conventional Born-Oppenheimer or Ehrenfest trajectory in the nonadiabatic coupling region. Instead, in NAF the equations of motion of the independent trajectory involve a nonadiabatic nuclear force term in addition to an adiabatic nuclear force term of a single electronic state. A few benchmark tests for gas phase and condensed phase systems indicate that NAF offers a practical tool to capture the correct correlation of electronic and nuclear dynamics for processes where the states remain coupled all the time as well as for the asymptotic region where the coupling of electronic states vanishes.
Collapse
Affiliation(s)
- Baihua Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Xin He
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
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
|
Akimov AV. Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2023; 14:11673-11683. [PMID: 38109379 DOI: 10.1021/acs.jpclett.3c03029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| |
Collapse
|
29
|
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
|
30
|
Dillon AD, Gieseking RLM. Convergence of Time-Derivative Nonadiabatic Couplings in Plane-Wave DFT Calculations. J Phys Chem A 2023; 127:9612-9620. [PMID: 37924298 DOI: 10.1021/acs.jpca.3c04858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2023]
Abstract
Accurate prediction of charge carrier relaxation rates is essential to design molecules and materials with the desired photochemical properties for applications like photocatalysis and solar energy conversion. Nonadiabatic molecular dynamics allows one to simulate the relaxation process of excited charge carriers. Plane-wave density functional theory (DFT) calculations make the time-derivative nonadiabatic couplings (TNACs) simple to compute because the basis is independent of the atomic positions. However, the effect of the kinetic energy cutoff for the plane-wave basis on the accuracy of the dynamics has not been studied. Here, we examine the effect of the kinetic energy cutoff on the TNACs and decay time scales for the prototypical model system of tetracene. These calculations show that the choice of kinetic energy cutoff can change the relaxation time by up to 30%. The relaxation times of states that have small TNACs to other states or are far from degenerate are more sensitive to the kinetic energy cutoff than those of states with large TNACs or near degeneracies. A kinetic energy cutoff of 60 Ry is sufficient for all states to reach semiquantitative agreement (absolute error <10%) with the decay times of our 110 Ry reference data, and a cutoff of 80 Ry is required for all states to reach quantitative agreement (absolute error <2%).
Collapse
Affiliation(s)
- Alva D Dillon
- Department of Chemistry, Brandeis University, 415 South Street, Waltham, Massachusetts 02453, United States
| | - Rebecca L M Gieseking
- Department of Chemistry, Brandeis University, 415 South Street, Waltham, Massachusetts 02453, United States
| |
Collapse
|
31
|
Polonius S, Zhuravel O, Bachmair B, Mai S. LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution. J Chem Theory Comput 2023; 19:7171-7186. [PMID: 37788824 PMCID: PMC10601485 DOI: 10.1021/acs.jctc.3c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Indexed: 10/05/2023]
Abstract
We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - Oleksandra Zhuravel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
32
|
Villaseco Arribas E, Vindel-Zandbergen P, Roy S, Maitra NT. Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics. Phys Chem Chem Phys 2023; 25:26380-26395. [PMID: 37750820 DOI: 10.1039/d3cp03464j] [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/27/2023]
Abstract
The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.
Collapse
Affiliation(s)
| | - Patricia Vindel-Zandbergen
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Saswata Roy
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
| |
Collapse
|
33
|
Wang B, Wu Y, Liu D, Vasenko AS, Casanova D, Prezhdo OV. Efficient Modeling of Quantum Dynamics of Charge Carriers in Materials Using Short Nonequilibrium Molecular Dynamics. J Phys Chem Lett 2023; 14:8289-8295. [PMID: 37681642 PMCID: PMC10518862 DOI: 10.1021/acs.jpclett.3c02187] [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/05/2023] [Accepted: 09/05/2023] [Indexed: 09/09/2023]
Abstract
Nonadiabatic molecular dynamics provides essential insights into excited-state processes, but it is computationally intense and simplifications are needed. The classical path approximation provides critical savings. Still, long heating and equilibration steps are required. We demonstrate that practical results can be obtained with short, partially equilibrated ab initio trajectories. Once the system's structure is adequate and essential fluctuations are sampled, the nonadiabatic Hamiltonian can be constructed. Local structures require only 1-2 ps trajectories, as demonstrated with point defects in metal halide perovskites. Short trajectories represent anharmonic motions common in defective structures, an essential improvement over the harmonic approximation around the optimized geometry. Glassy systems, such as grain boundaries, require different simulation protocols, e.g., involving machine learning force fields. 10-fold shorter trajectories generate 10-20% time scale errors, which are acceptable, given experimental uncertainties and other approximations. The practical NAMD protocol enables fast screening of excited-state dynamics for rapid exploration of new materials.
Collapse
Affiliation(s)
- Bipeng Wang
- Department
of Chemical Engineering, University of Southern
California, Los Angeles, California 90089, United States
| | - Yifan Wu
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | | | - Andrey S. Vasenko
- HSE
University, 101000 Moscow, Russia
- Donostia
International Physics Center (DIPC), 20018 San Sebastián-Donostia, Euskadi, Spain
| | - David Casanova
- Donostia
International Physics Center (DIPC), 20018 San Sebastián-Donostia, Euskadi, Spain
- IKERBASQUE, Basque Foundation for Science, 48009 Bilbao, Euskadi, Spain
| | - Oleg V. Prezhdo
- Department
of Chemical Engineering, University of Southern
California, Los Angeles, California 90089, United States
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
34
|
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
|
35
|
Ghosh A, Goswami B, Pal S, Sarkar P. How the Stacking Pattern Influences the Charge Transfer Dynamics of van der Waals Heterostructures: An Answer from a Time-Domain Ab Initio Study. J Phys Chem Lett 2023; 14:7672-7679. [PMID: 37603897 DOI: 10.1021/acs.jpclett.3c01626] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Here, we perform a time domain density functional study in conjunction with a non-adiabatic molecular dynamics (NAMD) simulation to investigate the charge carrier dynamics in a series of van der Waals heterostructures made of two-dimensional (2D) SnX2 (X = S or Se)-supported ZrS2, ZrSe2, and ZrSSe monolayers. Results from NAMD simulation reveal delayed electron-hole recombination (in the range of 0.53-2.13 ns) and ultrafast electron/hole transfer processes (electron transfer within 108.3-321.5 fs and hole transfer between 107.6 and 258.8 fs). The most interesting finding of our study is that switching from AB to AA stacking in the heterostructures extends the carrier lifespan by a significant amount. The delayed electron-hole recombination because of the switching stacking pattern can be rationalized by weak electron-phonon coupling, lower non-adiabatic coupling (NAC), and fast decoherence time. Thus, these insightful NAMD studies of excited charge carriers reveal that the stacking pattern variation is an effective tool to develop efficient photovoltaic devices based on 2D van der Waals heterostructures.
Collapse
Affiliation(s)
- Atish Ghosh
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| | - Biplab Goswami
- Department of Physics, Sreegopal Banerjee College, Bagati, Hoogly 712148, India
| | - Sougata Pal
- Department of Chemistry, University of Gour Banga, Malda 732103, India
| | - Pranab Sarkar
- Department of Chemistry, Visva-Bharati University, Santiniketan 731235, India
| |
Collapse
|
36
|
Gumber S, Prezhdo OV. Zeno and Anti-Zeno Effects in Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2023; 14:7274-7282. [PMID: 37556319 PMCID: PMC10440816 DOI: 10.1021/acs.jpclett.3c01831] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Accepted: 08/01/2023] [Indexed: 08/11/2023]
Abstract
Decoherence plays an important role in nonadiabatic (NA) molecular dynamics (MD) simulations because it provides a physical mechanism for trajectory hopping and can alter transition rates by orders of magnitude. Generally, decoherence effects slow quantum transitions, as exemplified by the quantum Zeno effect: in the limit of infinitely fast decoherence, the transitions stop. If the measurements are not sufficiently frequent, an opposite quantum anti-Zeno effect occurs, in which the transitions are accelerated with faster decoherence. Using two common NA-MD approaches, fewest switches surface hopping and decoherence-induced surface hopping, combined with analytic examination, we demonstrate that including decoherence into NA-MD slows down NA transitions; however, many realistic systems operate in the anti-Zeno regime. Therefore, it is important that NA-MD methods describe both Zeno and anti-Zeno effects. Numerical simulations of charge trapping and relaxation in graphitic carbon nitride suggest that time-dependent NA Hamiltonians encountered in realistic systems produce robust results with respect to errors in the decoherence time, a favorable feature for NA-MD simulations.
Collapse
Affiliation(s)
- Shriya Gumber
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
- Department
of Physics and Astronomy, University of
Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
37
|
Wang B, Winkler L, Wu Y, Müller KR, Sauceda HE, Prezhdo OV. Interpolating Nonadiabatic Molecular Dynamics Hamiltonian with Bidirectional Long Short-Term Memory Networks. J Phys Chem Lett 2023; 14:7092-7099. [PMID: 37530451 PMCID: PMC10424239 DOI: 10.1021/acs.jpclett.3c01723] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2023] [Accepted: 07/21/2023] [Indexed: 08/03/2023]
Abstract
Essential for understanding far-from-equilibrium processes, nonadiabatic (NA) molecular dynamics (MD) requires expensive calculations of the excitation energies and NA couplings. Machine learning (ML) can simplify computation; however, the NA Hamiltonian requires complex ML models due to its intricate relationship to atomic geometry. Working directly in the time domain, we employ bidirectional long short-term memory networks (Bi-LSTM) to interpolate the Hamiltonian. Applying this multiscale approach to three metal-halide perovskite systems, we achieve two orders of magnitude computational savings compared to direct ab initio calculation. Reasonable charge trapping and recombination times are obtained with NA Hamiltonian sampling every half a picosecond. The Bi-LSTM-NAMD method outperforms earlier models and captures both slow and fast time scales. In combination with ML force fields, the methodology extends NAMD simulation times from picoseconds to nanoseconds, comparable to charge carrier lifetimes in many materials. Nanosecond sampling is particularly important in systems containing defects, boundaries, interfaces, etc. that can undergo slow rearrangements.
Collapse
Affiliation(s)
- Bipeng Wang
- Department
of Chemical Engineering, University of Southern
California, Los Angeles, California 90089, United States
| | - Ludwig Winkler
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Yifan Wu
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- BIFOLD
- Berlin Institute for the Foundations of Learning and Data, 10587 Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul 136-713, Korea
- Max Planck
Institute for Informatics, Stuhlsatzenhausweg, 66123 Saarbrücken, Germany
- Google
Deepmind, 10587 Berlin, Germany
| | - Huziel E. Sauceda
- BASLEARN,
BASF-TU joint Lab, Technische Universität
Berlin, 10587 Berlin, Germany
- Departamento
de Materia Condensada, Instituto de Física, Universidad Nacional Autónoma de México, Apartado Postal 20-346, 01000 México, D.F., México
| | - Oleg V. Prezhdo
- Department
of Chemical Engineering, University of Southern
California, Los Angeles, California 90089, United States
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
38
|
Giri A, Walton SG, Tomko J, Bhatt N, Johnson MJ, Boris DR, Lu G, Caldwell JD, Prezhdo OV, Hopkins PE. Ultrafast and Nanoscale Energy Transduction Mechanisms and Coupled Thermal Transport across Interfaces. ACS NANO 2023; 17:14253-14282. [PMID: 37459320 PMCID: PMC10416573 DOI: 10.1021/acsnano.3c02417] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Accepted: 06/06/2023] [Indexed: 08/09/2023]
Abstract
The coupled interactions among the fundamental carriers of charge, heat, and electromagnetic fields at interfaces and boundaries give rise to energetic processes that enable a wide array of technologies. The energy transduction among these coupled carriers results in thermal dissipation at these surfaces, often quantified by the thermal boundary resistance, thus driving the functionalities of the modern nanotechnologies that are continuing to provide transformational benefits in computing, communication, health care, clean energy, power recycling, sensing, and manufacturing, to name a few. It is the purpose of this Review to summarize recent works that have been reported on ultrafast and nanoscale energy transduction and heat transfer mechanisms across interfaces when different thermal carriers couple near or across interfaces. We review coupled heat transfer mechanisms at interfaces of solids, liquids, gasses, and plasmas that drive the resulting interfacial heat transfer and temperature gradients due to energy and momentum coupling among various combinations of electrons, vibrons, photons, polaritons (plasmon polaritons and phonon polaritons), and molecules. These interfacial thermal transport processes with coupled energy carriers involve relatively recent research, and thus, several opportunities exist to further develop these nascent fields, which we comment on throughout the course of this Review.
Collapse
Affiliation(s)
- Ashutosh Giri
- Department
of Mechanical, Industrial and Systems Engineering, University of Rhode Island, Kingston, Rhode Island 02881, United States
| | - Scott G. Walton
- Plasma
Physics Division, Naval Research Laboratory, Washington, DC 22032, United States
| | - John Tomko
- Department
of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, Virginia 22904, United States
| | - Niraj Bhatt
- Department
of Mechanical, Industrial and Systems Engineering, University of Rhode Island, Kingston, Rhode Island 02881, United States
| | - Michael J. Johnson
- Plasma
Physics Division, Naval Research Laboratory, Washington, DC 22032, United States
| | - David R. Boris
- Plasma
Physics Division, Naval Research Laboratory, Washington, DC 22032, United States
| | - Guanyu Lu
- Department
of Mechanical Engineering, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Joshua D. Caldwell
- Department
of Mechanical Engineering, Vanderbilt University, Nashville, Tennessee 37235, United States
- Interdisciplinary
Materials Science, Vanderbilt University, Nashville, Tennessee 37235, United States
- Vanderbilt
Institute of Nanoscale Science and Engineering, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
- Department
of Physics and Astronomy, University of
Southern California, Los Angeles, California 90089, United States
| | - Patrick E. Hopkins
- Department
of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, Virginia 22904, United States
- Department
of Materials Science and Engineering, University
of Virginia, Charlottesville, Virginia 22904, United States
- Department
of Physics, University of Virginia, Charlottesville, Virginia 22904, United States
| |
Collapse
|
39
|
Chen CG, Giustini M, D'Abramo M, Amadei A. Unveiling the Excited State Dynamics of Indole in Solution. J Chem Theory Comput 2023. [PMID: 37329333 DOI: 10.1021/acs.jctc.3c00221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
In this paper, we reconstruct in detail the dynamics of the emitting electronic excited state of aqueous indole, investigating its relaxation mechanism and kinetics to be related to the time-dependent fluorescence signal. Taking advantage of the results shown in a very recent paper, we were able to model the relaxation process in solution in terms of the transitions between two gas-phase singlet electronic states (1La and 1Lb), subsequently irreversibly relaxing to the gas-phase singlet dark state (1πσ*). A comparison of the results with the available experimental data shows that the relaxation mechanism we obtain by our theoretical-computational model is reliable, reproducing rather accurately all the experimental observables.
Collapse
Affiliation(s)
| | - Mauro Giustini
- Department of Chemistry, Sapienza University of Rome, Rome 00185, Italy
| | - Marco D'Abramo
- Department of Chemistry, Sapienza University of Rome, Rome 00185, Italy
| | - Andrea Amadei
- Department of Technological and Chemical Sciences, Tor Vergata University of Rome, Rome 00133, Italy
| |
Collapse
|
40
|
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
|
41
|
Qiu T, Climent C, Subtonik JE. A Practical Approach to Wave Function Propagation, Hopping Probabilities, and Time Steps in Surface Hopping Calculations. J Chem Theory Comput 2023; 19:2744-2757. [PMID: 37130302 DOI: 10.1021/acs.jctc.3c00126] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
We compare several established approaches for propagating wave functions and calculating hopping probabilities within the fewest switches surface hopping (FSSH) algorithm for difficult cases with many electronic states and many trivial crossings. If only a single time step (Δtc) is employed, we find that no published approach can accurately capture the dynamics correctly unless Δtc → 0 (which is not computationally feasible). If multiple time steps are employed, for a fixed classical time step (Δtc), a robust scheme can be found for dynamically choosing quantum time steps (δtq1 and δtq2) and calculating hopping probabilities so that one can systematically reduce all errors and achieve maximally efficient accuracy; scattering calculations confirm that one can choose a fairly large classical time step. The robust scheme presented here uses both the "local diabatic" and adiabatic interpolation and thus borrows elements from both the Granucci/Persico and Meek/Levine algorithms. Our findings should be broadly applicable in the future.
Collapse
Affiliation(s)
- Tian Qiu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Clàudia Climent
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E Subtonik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
42
|
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
|
43
|
Lu D, Guo H. Quantum and Semiclassical Dynamics of Nonadiabatic Electronic Excitation of C( 3P) to C 1D) by Hyperthermal Collisions with N 2. J Phys Chem A 2023; 127:3190-3199. [PMID: 36989004 DOI: 10.1021/acs.jpca.3c00893] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
The dynamics and kinetics of nonadiabatic excitation of C(3P) to C(1D) induced by hyperthermal collisions with N2 molecules are investigated using a quantum mechanical and two semiclassical nonadiabatic methods. The full-dimensional interaction potential energy surfaces and spin-orbit coupling, which facilitates the spin-forbidden process, are represented by a recently constructed diabatic potential energy matrix. The multistate quantum dynamics for selected partial waves found small transition probabilities due to the weak spin-orbit coupling. The spin-flip transition is the most favored near the threshold due to effective curve crossing. Strong oscillations are also found in the probabilities, which are attributable to resonances supported by the deep well in the singlet-state potential. Vibrational state-specified rate coefficients are reported from J-shifted quantum dynamics calculations, and they follow the Arrhenius form. Vibrational excitation in the N2 collision partner is found to increase the excitation rate at low temperatures, but the trend is reversed at high temperatures. The two semiclassical methods qualitatively reproduce the quantum rate coefficients.
Collapse
Affiliation(s)
- Dandan Lu
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
44
|
Reiner M, Bachmair B, Tiefenbacher MX, Mai S, González L, Marquetand P, Dellago C. Nonadiabatic Forward Flux Sampling for Excited-State Rare Events. J Chem Theory Comput 2023; 19:1657-1671. [PMID: 36856706 PMCID: PMC10061683 DOI: 10.1021/acs.jctc.2c01088] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Indexed: 03/02/2023]
Abstract
We present a rare event sampling scheme applicable to coupled electronic excited states. In particular, we extend the forward flux sampling (FFS) method for rare event sampling to a nonadiabatic version (NAFFS) that uses the trajectory surface hopping (TSH) method for nonadiabatic dynamics. NAFFS is applied to two dynamically relevant excited-state models that feature an avoided crossing and a conical intersection with tunable parameters. We investigate how nonadiabatic couplings, temperature, and reaction barriers affect transition rate constants in regimes that cannot be otherwise obtained with plain, traditional TSH. The comparison with reference brute-force TSH simulations for limiting cases of rareness shows that NAFFS can be several orders of magnitude cheaper than conventional TSH and thus represents a conceptually novel tool to extend excited-state dynamics to time scales that are able to capture rare nonadiabatic events.
Collapse
Affiliation(s)
- Madlen
Maria Reiner
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Maximilian Xaver Tiefenbacher
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Leticia González
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Christoph Dellago
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Faculty
of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| |
Collapse
|
45
|
Toldo JM, do Casal MT, Ventura E, do Monte SA, Barbatti M. Surface hopping modeling of charge and energy transfer in active environments. Phys Chem Chem Phys 2023; 25:8293-8316. [PMID: 36916738 PMCID: PMC10034598 DOI: 10.1039/d3cp00247k] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
Abstract
An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.
Collapse
Affiliation(s)
| | | | - Elizete Ventura
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Silmar A do Monte
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille, France.
- Institut Universitaire de France, 75231, Paris, France
| |
Collapse
|
46
|
Dergachev ID, Dergachev VD, Rooein M, Mirzanejad A, Varganov SA. Predicting Kinetics and Dynamics of Spin-Dependent Processes. Acc Chem Res 2023; 56:856-866. [PMID: 36926853 DOI: 10.1021/acs.accounts.2c00843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
Abstract
ConspectusPredicting mechanisms and rates of nonadiabatic spin-dependent processes including photoinduced intersystem crossings, thermally activated spin-forbidden reactions, and spin crossovers in metal centers is a very active field of research. These processes play critical roles in transition-metal-based and metalloenzymatic catalysis, molecular magnets, light-harvesting materials, organic light-emitting diodes, photosensitizers for photodynamic therapy, and many other applications. Therefore, accurate modeling of spin-dependent processes in complex systems and on different time scales is important for many problems in chemistry, biochemistry, and materials sciences.Nonadiabatic statistical theory (NAST) and nonadiabatic molecular dynamics (NAMD) are two complementary approaches to modeling the kinetics and dynamics of spin-dependent processes. NAST predicts the probabilities and rate constants of nonradiative transitions between electronic states with different spin multiplicities using molecular properties at only few critical points on the potential energy surfaces (PESs), including the reactant minimum and the minimum energy crossing point (MECP) between two spin states. This makes it possible to obtain molecular properties for NAST calculations using accurate but often computationally expensive electronic structure methods, which is critical for predicting the rate constants of spin-dependent processes. Alternatively, NAST can be used to study spin-dependent processes in very large complex molecular systems using less computationally expensive electronic structure methods. The nuclear quantum effects, such as zero-point vibrational energy, tunneling, and interference between reaction paths can be easily incorporated. However, the statistical and local nature of NAST makes it more suitable for large systems and slow kinetics. In contrast, NAMD explores entire PESs of interacting electronic states, making it ideal for modeling fast barrierless spin-dependent processes. Because the knowledge of large portions of PESs is often needed, the simulations require a very large number of electronic structure calculations, which limits the NAMD applicability to relatively small molecular systems and ultrafast kinetics.In this Account, we discuss our contribution to the development of the NAST and NAMD approaches for predicting the rates and mechanism of spin-dependent processes. First, we briefly describe our NAST and NAMD implementations. The NAST implementation is an extension of the transition state theory to the processes involving two crossing potential energy surfaces of different spin multiplicities. The NAMD approach includes the trajectory surface hopping (TSH) and ab initio multiple spawning (AIMS) methods. Second, we discuss several applications of NAST and NAMD to model spin-dependent processes in different systems. The NAST applicability to large complex systems is demonstrated by the studies of the spin-forbidden isomerization of the active sites of metal-sulfur proteins. Our implementation of the MECP search algorithm within the fully ab initio fragment molecular orbital method allows applying NAST to systems with thousands of atoms, such as the solvated protein rubredoxin. Applications of NAMD to ultrafast spin-dependent processes are represented by the generalized AIMS simulations utilizing the fast GPU-based TeraChem electronic structure program to gain insight into the complex photoexcited state relaxation in 2-cyclopentenone.
Collapse
Affiliation(s)
- Ilya D Dergachev
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Vsevolod D Dergachev
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Mitra Rooein
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Amir Mirzanejad
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
47
|
Abstract
We present a nonadiabatic classical-trajectory approach that offers the best of both worlds between fewest-switches surface hopping (FSSH) and quasiclassical mapping dynamics. This mapping approach to surface hopping (MASH) propagates the nuclei on the active adiabatic potential-energy surface, such as in FSSH. However, unlike in FSSH, transitions between active surfaces are deterministic and occur when the electronic mapping variables evolve between specified regions of the electronic phase space. This guarantees internal consistency between the active surface and the electronic degrees of freedom throughout the dynamics. MASH is rigorously derivable from exact quantum mechanics as a limit of the quantum-classical Liouville equation (QCLE), leading to a unique prescription for momentum rescaling and frustrated hops. Hence, a quantum-jump procedure can, in principle, be used to systematically converge the accuracy of the results to that of the QCLE. This jump procedure also provides a rigorous framework for deriving approximate decoherence corrections similar to those proposed for FSSH. We apply MASH to simulate the nonadiabatic dynamics in various model systems and show that it consistently produces more accurate results than FSSH at a comparable computational cost.
Collapse
|
48
|
Li W, Xue T, Mora-Perez C, Prezhdo OV. Ab initio quantum dynamics of plasmonic charge carriers. TRENDS IN CHEMISTRY 2023. [DOI: 10.1016/j.trechm.2023.02.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023]
|
49
|
Cui P, Xue Y. Sulfication-induced non-radiative electron-hole recombination dynamics in graphene quantum dots for tuning photocatalytic performance. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 287:122117. [PMID: 36403541 DOI: 10.1016/j.saa.2022.122117] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Revised: 10/16/2022] [Accepted: 11/08/2022] [Indexed: 06/16/2023]
Abstract
GQDs, or graphene quantum dots, are promising materials for energy-related applications. Their optoelectronic properties can be modified by adding heteroatoms, making them good candidates for photocatalysts. However, the structure-property relationship of these materials still needs to be investigated to control their properties better. In particular, photocatalysis of GQDs is hindered by non-radiative electron-hole recombination. In this study, density functional theory (DFT) calculations were performed to investigate the electronic structures and optical properties of GQDs doped with three distinct sulfur functional groups, i.e., sulfur oxide (O3HS), sulfhydryl (SH), and thiophene (C4H4S), respectively. The results suggest that sulfur doping decreases the GQD bandgap. In particular, the asymmetric capping of the GQD edges with the C4H4S groups led to additional peaks at low excitation energies, whereas for GQDs functionalized with O3HS or SH groups, only a shift in the main absorption peak or a change in the absorption intensity was observed. SH functionalization drastically increased electronic coupling, while C4H4S functionalization induced more charge-relaxation channels in the GQDs. Thus, the results shed light on the mechanisms governing the photocatalytic efficiency of GQDs.
Collapse
Affiliation(s)
- Peng Cui
- School of New Materials and Shoes & Clothing Engineering, Liming Vocational University, Quanzhou 362000, Fujian Province, PR China; Nanotechnology Research Laboratory, Jiangnan University, Wuxi 214122, Jiangsu Province, PR China.
| | - Yuan Xue
- Nanotechnology Research Laboratory, Jiangnan University, Wuxi 214122, Jiangsu Province, PR China
| |
Collapse
|
50
|
Bondarenko AS, Tempelaar R. Overcoming positivity violations for density matrices in surface hopping. J Chem Phys 2023; 158:054117. [PMID: 36754802 DOI: 10.1063/5.0135456] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space.
Collapse
Affiliation(s)
- Anna S Bondarenko
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| | - Roel Tempelaar
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| |
Collapse
|