1
|
Zhu Y, Peng J, Xu C, Lan Z. Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation. J Phys Chem Lett 2024; 15:9601-9619. [PMID: 39270134 DOI: 10.1021/acs.jpclett.4c01751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2024]
Abstract
The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.
Collapse
Affiliation(s)
- Yifei Zhu
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
2
|
Duston T, Tao Z, Bian X, Bhati M, Rawlinson J, Littlejohn RG, Pei Z, Shao Y, Subotnik JE. A Phase-Space Electronic Hamiltonian For Vibrational Circular Dichroism. J Chem Theory Comput 2024. [PMID: 39226223 DOI: 10.1021/acs.jctc.4c00662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
We show empirically that a phase-space non-Born-Oppenheimer electronic Hamiltonian approach to quantum chemistry (where the electronic Hamiltonian is parametrized by both nuclear position and momentum, ĤPS(R,P)) is both a practical and accurate means to recover vibrational circular dichroism spectra. We further hypothesize that such a phase-space approach may lead to very new dynamical physics beyond spectroscopic circular dichroism, with potential implications for understanding chiral induced spin selectivity (CISS), noting that classical phase-space approaches conserve the total nuclear plus electronic momentum, whereas classical Born-Oppenheimer approaches do not (they conserve only the nuclear momentum).
Collapse
Affiliation(s)
- Titouan Duston
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Mansi Bhati
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Jonathan Rawlinson
- Department of Mathematics, Nottingham Trent University, Nottingham NG1 4FQ, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, United States
| | - Zheng Pei
- Department of Chemistry, The University of Oklahoma, Norman, Oklahoma 73104, United States
| | - Yihan Shao
- Department of Chemistry, The University of Oklahoma, Norman, Oklahoma 73104, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
3
|
Zhang L, Pios SV, Martyka M, Ge F, Hou YF, Chen Y, Chen L, Jankowska J, Barbatti M, Dral PO. MLatom Software Ecosystem for Surface Hopping Dynamics in Python with Quantum Mechanical and Machine Learning Methods. J Chem Theory Comput 2024; 20:5043-5057. [PMID: 38836623 DOI: 10.1021/acs.jctc.4c00468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2024]
Abstract
We present an open-source MLatom@XACS software ecosystem for on-the-fly surface hopping nonadiabatic dynamics based on the Landau-Zener-Belyaev-Lebedev algorithm. The dynamics can be performed via Python API with a wide range of quantum mechanical (QM) and machine learning (ML) methods, including ab initio QM (CASSCF and ADC(2)), semiempirical QM methods (e.g., AM1, PM3, OMx, and ODMx), and many types of ML potentials (e.g., KREG, ANI, and MACE). Combinations of QM and ML methods can also be used. While the user can build their own combinations, we provide AIQM1, which is based on Δ-learning and can be used out-of-the-box. We showcase how AIQM1 reproduces the isomerization quantum yield of trans-azobenzene at a low cost. We provide example scripts that, in dozens of lines, enable the user to obtain the final population plots by simply providing the initial geometry of a molecule. Thus, those scripts perform geometry optimization, normal mode calculations, initial condition sampling, parallel trajectories propagation, population analysis, and final result plotting. Given the capabilities of MLatom to be used for training different ML models, this ecosystem can be seamlessly integrated into the protocols building ML models for nonadiabatic dynamics. In the future, a deeper and more efficient integration of MLatom with Newton-X will enable a vast range of functionalities for surface hopping dynamics, such as fewest-switches surface hopping, to facilitate similar workflows via the Python API.
Collapse
Affiliation(s)
- Lina Zhang
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Sebastian V Pios
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Mikołaj Martyka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Fuchun Ge
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yi-Fan Hou
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yuxinxin Chen
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Lipeng Chen
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, Marseille 13397, France
- Institut Universitaire de France, Paris 75231, France
| | - Pavlo O Dral
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| |
Collapse
|
4
|
Wada S, Tsutsumi T, Saita K, Taketsugu T. Ab initio molecular dynamics study of intersystem crossing dynamics for MH 2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces. J Comput Chem 2024; 45:552-562. [PMID: 38009451 DOI: 10.1002/jcc.27271] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 10/18/2023] [Accepted: 11/14/2023] [Indexed: 11/28/2023]
Abstract
Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.
Collapse
Affiliation(s)
- Satoi Wada
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo, Japan
| | - Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
| | - Kenichiro Saita
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Japan
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Lei Y, Zheng Z, Vasquez L, Zhao J, Ma J, Ma H. Enhanced Electron Transfer and Spin Flip through Spin-Orbital Couplings in Organic/Inorganic Heterojunctions: A Nonadiabatic Surface Hopping Simulation. J Phys Chem Lett 2022; 13:4840-4848. [PMID: 35616399 DOI: 10.1021/acs.jpclett.2c01177] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The circumstances of transferred electrons across organic/inorganic interfaces have attracted intensive interest because of the distinctive electronic structure properties of those two components. Leveraging ab initio nonadiabatic molecular dynamics methods in conjunction with spin dynamics induced by spin-orbital couplings (SOCs), this study reports two competitive channels during photoinduced dynamical processes in the prototypical ZnPc/monolayer MoS2 heterojunction. Interestingly, the electron-transfer and relaxation processes occur simultaneously because of the enhancement of electron-phonon couplings and expansion of dynamical pathways by SOCs, suggesting that the electron-transfer rate and relaxation processes can be tuned by SOCs, hence yielding the performance promotion of photovoltaic and photocatalytic devices. Additionally, approximately half of the transferred electrons flip their spin within 1.6 ps because of strong SOCs in MoS2, achieving great agreement with experimental measurements. This investigation provides instructive perspectives for designing novel devices and applications based on organic/inorganic heterojunctions, demonstrating the importance of spin dynamics simulations in exploring sophisticated photoinduced processes in materials.
Collapse
Affiliation(s)
- Yuli Lei
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Zhenfa Zheng
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Luis Vasquez
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Jin Zhao
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jing Ma
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
8
|
Liu J, Lan Z, Yang J. An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics. Phys Chem Chem Phys 2021; 23:22313-22323. [PMID: 34591049 DOI: 10.1039/d1cp03477d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.
Collapse
Affiliation(s)
- Jie Liu
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| |
Collapse
|
9
|
Abstract
When nonadiabatic dynamics are described on the basis of trajectories, severe trajectory branching occurs when the nuclear wave packets on some potential energy surfaces are reflected while those on the remaining surfaces are not. As a result, the traditional Ehrenfest mean field (EMF) approximation breaks down. In this study, two versions of the branching corrected mean field (BCMF) method are proposed. Namely, when trajectory branching is identified, BCMF stochastically selects either the reflected or the nonreflected group to build the new mean field trajectory or splits the mean field trajectory into two new trajectories with the corresponding weights. As benchmarked in six standard model systems and an extensive model base with two hundred diverse scattering models, BCMF significantly improves the accuracy while retaining the high efficiency of the traditional EMF. In fact, BCMF closely reproduces the exact quantum dynamics in all investigated systems, thus highlighting the essential role of branching correction in nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
- Jiabo Xu
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
10
|
Smith B, Akimov AV. Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:073001. [PMID: 31661681 DOI: 10.1088/1361-648x/ab5246] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States of America
| | | |
Collapse
|