151
|
Ding B, Yu Y, Geng S, Liu B, Hao Y, Liang G. Computational Methods for the Interaction between Cyclodextrins and Natural Compounds: Technology, Benefits, Limitations, and Trends. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2022; 70:2466-2482. [PMID: 35170315 DOI: 10.1021/acs.jafc.1c07018] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Cyclodextrins (CDs) have a hollow structure with a hydrophobic interior and hydrophilic exterior. Forming inclusion complexes with CDs will maximize the bioavailability of natural compounds and enable active components to be processed into functional foods, medicines, additives, and so forth. However, experimental methods cannot explain CD-guest binding at the atomic level. Different models have been recently developed to simulate the interaction between CDs and guests to study the binding conformation and analyze noncovalent forces. This review paper summarizes modeling methods of CD-natural compound complexes. The methods include quantitative structure-activity relationships, molecular docking, molecular dynamics simulations, and quantum-chemical calculations. The applications of these methods to enhance the solubility and bioactivities of guest molecules, assist material transportation, and promote compound extraction are also discussed. The purpose of this review is to explore interaction mechanisms of CDs and guests and to help expand new applications of CDs.
Collapse
Affiliation(s)
- Botian Ding
- Key Laboratory of Biorheological Science and Technology, Ministry of Education, Bioengineering College, Chongqing University, Chongqing 400044, China
| | - Yuandong Yu
- Key Laboratory of Biorheological Science and Technology, Ministry of Education, Bioengineering College, Chongqing University, Chongqing 400044, China
| | - Sheng Geng
- Key Laboratory of Biorheological Science and Technology, Ministry of Education, Bioengineering College, Chongqing University, Chongqing 400044, China
| | - Benguo Liu
- School of Food Science, Henan Institute of Science and Technology, Xinxiang 453003, China
| | - Youjin Hao
- College of Life Sciences, Chongqing Normal University, Chongqing, 401331, China
| | - Guizhao Liang
- Key Laboratory of Biorheological Science and Technology, Ministry of Education, Bioengineering College, Chongqing University, Chongqing 400044, China
| |
Collapse
|
152
|
NAST: Nonadiabatic Statistical Theory Package for Predicting Kinetics of Spin-Dependent Processes. Top Curr Chem (Cham) 2022; 380:15. [PMID: 35201520 DOI: 10.1007/s41061-022-00366-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 01/15/2022] [Indexed: 10/19/2022]
Abstract
We present a nonadiabatic statistical theory (NAST) package for predicting kinetics of spin-dependent processes, such as intersystem crossings, spin-forbidden unimolecular reactions, and spin crossovers. The NAST package can calculate the probabilities and rates of transitions between the electronic states of different spin multiplicities. Both the microcanonical (energy-dependent) and canonical (temperature-dependent) rate constants can be obtained. Quantum effects, including tunneling, zero-point vibrational energy, and reaction path interference, can be accounted for. In the limit of an adiabatic unimolecular reaction proceeding on a single electronic state, NAST reduces to the traditional transition state theory. Because NAST requires molecular properties at only a few points on potential energy surfaces, it can be applied to large molecular systems, used with accurate high-level electronic structure methods, and employed to study slow nonadiabatic processes. The essential NAST input data include the nuclear Hessian at the reactant minimum, as well as the nuclear Hessians, energy gradients, and spin-orbit coupling at the minimum energy crossing point (MECP) between two states. The additional computational tools included in the NAST package can be used to extract the required input data from the output files of electronic structure packages, calculate the effective Hessian at the MECP, and fit the reaction coordinate for more advanced NAST calculations. We describe the theory, its implementation, and three examples of application to different molecular systems.
Collapse
|
153
|
Ibele LM, Curchod BFE, Agostini F. A Photochemical Reaction in Different Theoretical Representations. J Phys Chem A 2022; 126:1263-1281. [PMID: 35157450 PMCID: PMC8883471 DOI: 10.1021/acs.jpca.1c09604] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
The Born–Oppenheimer
picture has forged our representation
and interpretation of photochemical processes, from photoexcitation
down to the passage through a conical intersection, a funnel connecting
different electronic states. In this work, we analyze a full in silico
photochemical experiment, from the explicit electronic excitation
by a laser pulse to the formation of photoproducts following a nonradiative
decay through a conical intersection, by contrasting the picture offered
by Born–Oppenheimer and that proposed by the exact factorization.
The exact factorization offers an alternative understanding of photochemistry
that does not rely on concepts such as electronic states, nonadiabatic
couplings, and conical intersections. On the basis of nonadiabatic
quantum dynamics performed for a two-state 2D model system, this work
allows us to compare Born–Oppenheimer and exact factorization
for (i) an explicit photoexcitation with and without the Condon approximation,
(ii) the passage of a nuclear wavepacket through a conical intersection,
(iii) the formation of excited stationary states in the Franck–Condon
region, and (iv) the use of classical and quantum trajectories in
the exact factorization picture to capture nonadiabatic processes
triggered by a laser pulse.
Collapse
Affiliation(s)
- Lea M Ibele
- Department of Chemistry, Durham University, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Department of Chemistry, Durham University, Durham DH1 3LE, United Kingdom
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
154
|
Vindel-Zandbergen P, Matsika S, Maitra NT. Exact-Factorization-Based Surface Hopping for Multistate Dynamics. J Phys Chem Lett 2022; 13:1785-1790. [PMID: 35170972 DOI: 10.1021/acs.jpclett.1c04132] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
A surface-hopping algorithm recently derived from the exact factorization approach, SHXF [Ha et al. J. Phys. Chem. Lett. 2018, 9, 1097], introduces an additional term in the electronic equation of surface hopping that couples electronic states through the quantum momentum. This term not only provides a first-principles description of decoherence, but here we show it is crucial to accurately capture nonadiabatic dynamics when more than two states are occupied at any given time. Using a vibronic coupling model of the uracil cation, we show that the lack of this term in traditional surface-hopping methods, including those with decoherence corrections, leads to failure to predict the dynamics through a three-state intersection, while SHXF performs similarly to the multiconfiguration time-dependent Hartree quantum dynamics benchmark.
Collapse
Affiliation(s)
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
155
|
Shu Y, Varga Z, Kanchanakungwankul S, Zhang L, Truhlar DG. Diabatic States of Molecules. J Phys Chem A 2022; 126:992-1018. [PMID: 35138102 DOI: 10.1021/acs.jpca.1c10583] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
156
|
Sasikumar D, Vinod K, Sunny J, Hariharan M. Exciton interactions in helical crystals of a hydrogen-bonded eumelanin monomer. Chem Sci 2022; 13:2331-2338. [PMID: 35310511 PMCID: PMC8864807 DOI: 10.1039/d1sc06755a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Accepted: 01/19/2022] [Indexed: 12/04/2022] Open
Abstract
Eumelanin, a naturally occurring group of heterogeneous polymers/aggregates providing photoprotection to living organisms, consist of 5,6-dihydroxyindole (DHI) and 5,6-dihydroxyindole-2-carboxylic acid (DHICA) building blocks. Despite their prevalence in the animal world, the structure and therefore the mechanism behind the photoprotective broadband absorption and non-radiative decay of eumelanin remain largely unknown. As a small step towards solving the incessant mystery, DHI is crystallized in a non-protic solvent environment to obtain DHI crystals having a helical packing motif. The present approach reflects the solitary directional effect of hydrogen bonds between the DHI chromophores for generating the crystalline assembly and filters out any involvement of the surrounding solvent environment. The DHI single crystals having an atypical chiral packing motif (P212121 Sohncke space group) incorporate enantiomeric zig-zag helical stacks arranged in a herringbone fashion with respect to each other. Each of the zig-zag helical stacks originates from a bifurcated hydrogen bonding interaction between the hydroxyl substituents in adjacent DHI chromophores which act as the backbone structure for the helical assembly. Fragment-based excited state analysis performed on the DHI crystalline assembly demonstrates exciton delocalization along the DHI units that connect each enantiomeric helical stack while, within each stack, the excitons remain localized. Fascinatingly, over the time evolution for generation of single-crystals of the DHI-monomer, mesoscopic double-helical crystals are formed, possibly attributed to the presence of covalently connected DHI trimers in chloroform solution. The oligomeric DHI (in line with the chemical disorder model) along with the characteristic crystalline packing observed for DHI provides insights into the broadband absorption feature exhibited by the chromophore. Single crystals of DHI monomer, a eumelanin precursor, adopt an atypical chiral packing arrangement incorporating enantiomeric zig-zag helical stacks while its covalently connected DHI trimer forms double-helical crystals in the mesoscopic scale.![]()
Collapse
Affiliation(s)
- Devika Sasikumar
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Maruthamala P.O., Vithura Thiruvananthapuram Kerala 695551 India
| | - Kavya Vinod
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Maruthamala P.O., Vithura Thiruvananthapuram Kerala 695551 India
| | - Jeswin Sunny
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Maruthamala P.O., Vithura Thiruvananthapuram Kerala 695551 India
| | - Mahesh Hariharan
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Maruthamala P.O., Vithura Thiruvananthapuram Kerala 695551 India
| |
Collapse
|
157
|
Talbot JJ, Head-Gordon M, Miller WH, Cotton SJ. Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation. Phys Chem Chem Phys 2022; 24:4820-4831. [PMID: 35156112 DOI: 10.1039/d1cp04090a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Using a simulated pump-probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer-Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D1Σ+ adiabatic state that is embedded in the manifold of near-dissociation C1Σ+ vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.
Collapse
Affiliation(s)
- Justin J Talbot
- Department of Chemistry, University of California, Berkeley, California 94720, USA.
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, 94720, USA
| | - William H Miller
- Department of Chemistry, University of California, Berkeley, California 94720, USA.
| | - Stephen J Cotton
- Department of Chemistry, University of California, Berkeley, California 94720, USA.
| |
Collapse
|
158
|
Poštulka J, Slavíček P, Pysanenko A, Poterya V, Fárník M. Bimolecular reactions on sticky and slippery clusters: Electron-induced reactions of hydrogen peroxide. J Chem Phys 2022; 156:054306. [DOI: 10.1063/5.0079283] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
- Jan Poštulka
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Andriy Pysanenko
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| | - Viktoriya Poterya
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| | - Michal Fárník
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| |
Collapse
|
159
|
Shu Y, Zhang L, Chen X, Sun S, Huang Y, Truhlar DG. Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping. J Chem Theory Comput 2022; 18:1320-1328. [PMID: 35104136 DOI: 10.1021/acs.jctc.1c01080] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Direct dynamics by mixed quantum-classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive time steps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature-driven coherent switching with decay of mixing (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computation time and accuracy as compared to previous mixed quantum-classical nonadiabatic dynamics algorithms. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.,School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Xiye Chen
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.,School of Chemistry and Chemical 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
| | - Yudong Huang
- School of Chemistry and Chemical 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, USA
| |
Collapse
|
160
|
Mukherjee S, Ravi S, Dutta J, Sardar S, Adhikari S. Beyond Born-Oppenheimer based diabatic surfaces of 1,3,5-C 6H 3F 3+ to generate the photoelectron spectra using time-dependent discrete variable representation approach. Phys Chem Chem Phys 2022; 24:2185-2202. [PMID: 35006221 DOI: 10.1039/d1cp04733g] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this article, Beyond Born-Oppenheimer (BBO) treatment is implemented to construct diabatic potential energy surfaces (PESs) of 1,3,5-C6H3F3+ over a series [eighteen (18)] of two-dimensional (2D) nuclear planes constituted with eleven normal modes (Q2, Q9x, Q9y, Q13x, Q13y, Q18x, Q18y, Q10x, Q10y, Q12x and Q12y) to include all possible nonadiabatic interactions among six coupled electronic states (X̃2E'', , B̃2E' and ). We had formulated explicit expressions of adiabatic to diabatic transformation (ADT) equations [S. Mukherjee, J. Dutta, B. Mukherjee, S. Sardar and S. Adhikari, J. Chem. Phys., 2019, 150, 064308] for the same system forming six state sub-Hilbert space and at present, these ADT equations are solved by incorporating MRCI level ab initio adiabatic PESs and CP-MCSCF calculated nonadiabatic coupling terms (NACTs) to derive diabatic PESs and couplings. Such single-valued, smooth, symmetric and continuous diabatic surface matrices are utilized to carry out multi-state multi-mode nuclear dynamics with the help of time-dependent discrete variable representation (TDDVR) methodology to compute the photoelectron (PE) spectra of 1,3,5-C6H3F3. Our theoretically calculated spectra for X̃2E'', and states using BBO treatment and TDDVR dynamics show peak by peak correspondence with the experimental results as well as better than the findings of the multi-configuration time-dependent Hartree (MCTDH) method.
Collapse
Affiliation(s)
- Soumya Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India.
| | - Satyam Ravi
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India. .,School of Advance Science and Languages, VIT Bhopal University, Bhopal - 466114, India
| | - Joy Dutta
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India.
| | - Subhankar Sardar
- Department of Chemistry, Bhatter College, Dantan, Paschim Medinipur - 721426, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India.
| |
Collapse
|
161
|
Ananth N. Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics. Annu Rev Phys Chem 2022; 73:299-322. [PMID: 35081325 DOI: 10.1146/annurev-physchem-082620-021809] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, USA;
| |
Collapse
|
162
|
Baudin P, Mouvet F, Rothlisberger U. A multiple time step algorithm for trajectory surface hopping simulations. J Chem Phys 2022; 156:034107. [DOI: 10.1063/5.0065728] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Pablo Baudin
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - François Mouvet
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
163
|
Scheidegger A, Vaníček J, Golubev NV. Search for long-lasting electronic coherence using on-the-fly ab initio semiclassical dynamics. J Chem Phys 2022; 156:034104. [DOI: 10.1063/5.0076609] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Alan Scheidegger
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Nikolay V. Golubev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
164
|
Xu C, Lin K, Hu D, Gu FL, Gelin MF, Lan Z. Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods. J Phys Chem Lett 2022; 13:661-668. [PMID: 35023755 DOI: 10.1021/acs.jpclett.1c03373] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.
Collapse
Affiliation(s)
- Chao Xu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou, 510006, P. R. China
| | - Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou, 510006, P. R. China
| | - Deping Hu
- 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
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou, 510006, P. R. China
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou, 310018, P. R. 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
| |
Collapse
|
165
|
List NH, Jones CM, Martínez TJ. Internal conversion of the anionic GFP chromophore: in and out of the I-twisted S 1/S 0 conical intersection seam. Chem Sci 2022; 13:373-385. [PMID: 35126970 PMCID: PMC8729814 DOI: 10.1039/d1sc05849e] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2021] [Accepted: 11/11/2021] [Indexed: 11/21/2022] Open
Abstract
The functional diversity of the green fluorescent protein (GFP) family is intimately connected to the interplay between competing photo-induced transformations of the chromophore motif, anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI-). Its ability to undergo Z/E-isomerization is of particular importance for super-resolution microscopy and emerging opportunities in optogenetics. Yet, key dynamical features of the underlying internal conversion process in the native HBDI- chromophore remain largely elusive. We investigate the intrinsic excited-state behavior of isolated HBDI- to resolve competing decay pathways and map out the factors governing efficiency and the stereochemical outcome of photoisomerization. Based on non-adiabatic dynamics simulations, we demonstrate that non-selective progress along the two bridge-torsional (i.e., phenolate, P, or imidazolinone, I) pathways accounts for the three decay constants reported experimentally, leading to competing ultrafast relaxation primarily along the I-twisted pathway and S1 trapping along the P-torsion. The majority of the population (∼70%) is transferred to S0 in the vicinity of two approximately enantiomeric minima on the I-twisted intersection seam (MECI-Is). Despite their sloped, reactant-biased topographies (suggesting low photoproduct yields), we find that decay through these intersections leads to products with a surprisingly high quantum yield of ∼30%. This demonstrates that E-isomer generation results at least in part from direct isomerization on the excited state. A photoisomerization committor analysis reveals a difference in intrinsic photoreactivity of the two MECI-Is and that the observed photoisomerization is the combined result of two effects: early, non-statistical dynamics around the less reactive intersection followed by later, near-statistical behavior around the more reactive MECI-I. Our work offers new insight into internal conversion of HBDI- that both establishes the intrinsic properties of the chromophore and enlightens principles for the design of chromophore derivatives and protein variants with improved photoswitching properties.
Collapse
Affiliation(s)
- Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Chey M Jones
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| |
Collapse
|
166
|
Zhu M, Zheng C. Post-spin crossing dynamics determine the regioselectivity in open-shell singlet biradical recombination. Org Chem Front 2022. [DOI: 10.1039/d1qo01757h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Comprehensive computational studies reveal unique dynamic effects in a multi-spin-state reaction that determine the regioselectivity of a biradical recombination process.
Collapse
Affiliation(s)
- Min Zhu
- State Key Laboratory of Organometallic Chemistry, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, 345 Lingling Lu, Shanghai 200032, China
- School of Physical Science and Technology, ShanghaiTech University, 100 Haike Road, Shanghai 201210, China
| | - Chao Zheng
- State Key Laboratory of Organometallic Chemistry, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, 345 Lingling Lu, Shanghai 200032, China
| |
Collapse
|
167
|
Zhao Y, Sun K, Chen L, Gelin M. The hierarchy of Davydov's Ansätze and its applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1589] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Yang Zhao
- Division of Materials Science Nanyang Technological University Singapore Singapore
| | - Kewei Sun
- Division of Materials Science Nanyang Technological University Singapore Singapore
- School of Science, Hanghzhou Dianzi University Hangzhou China
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems Dresden Germany
| | - Maxim Gelin
- School of Science, Hanghzhou Dianzi University Hangzhou China
| |
Collapse
|
168
|
Akimov AV. Extending the Time Scales of Nonadiabatic Molecular Dynamics via Machine Learning in the Time Domain. J Phys Chem Lett 2021; 12:12119-12128. [PMID: 34913701 DOI: 10.1021/acs.jpclett.1c03823] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
A novel methodology for direct modeling of long-time scale nonadiabatic dynamics in extended nanoscale and solid-state systems is developed. The presented approach enables forecasting the vibronic Hamiltonians as a direct function of time via machine-learning models trained directly in the time domain. The use of periodic and aperiodic functions that transform time into effective input modes of the artificial neural network is demonstrated to be essential for such an approach to work for both abstract and atomistic models. The best strategies and possible limitations pertaining to the new methodology are explored and discussed. An exemplary direct simulation of unprecedentedly long 20 picosecond trajectories is conducted for a divacancy-containing monolayer black phosphorus system, and the importance of conducting such extended simulations is demonstrated. New insights into the excited states photophysics in this system are presented, including the role of decoherence and model definition.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
169
|
Runeson JE, Richardson JO. Quantum Entanglement from Classical Trajectories. PHYSICAL REVIEW LETTERS 2021; 127:250403. [PMID: 35029436 DOI: 10.1103/physrevlett.127.250403] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Revised: 08/20/2021] [Accepted: 11/05/2021] [Indexed: 06/14/2023]
Abstract
A long-standing challenge in mixed quantum-classical trajectory simulations is the treatment of entanglement between the classical and quantal degrees of freedom. We present a novel approach that describes the emergence of entangled states entirely in terms of independent and deterministic Ehrenfest-like classical trajectories. For a two-level quantum system in a classical environment, this is derived by mapping the quantum system onto a path-integral representation of a spin 1/2. We demonstrate that the method correctly accounts for coherence and decoherence and thus reproduces the splitting of a wave packet in a nonadiabatic scattering problem. This discovery opens up a new class of simulations as an alternative to stochastic surface-hopping, coupled-trajectory, or semiclassical approaches.
Collapse
Affiliation(s)
- Johan E Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
170
|
Wang S, Li C, Evangelista FA. Analytic Energy Gradients for the Driven Similarity Renormalization Group Multireference Second-Order Perturbation Theory. J Chem Theory Comput 2021; 17:7666-7681. [PMID: 34839660 DOI: 10.1021/acs.jctc.1c00980] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving for the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of the three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.
Collapse
Affiliation(s)
- Shuhe Wang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
171
|
Andrade X, Pemmaraju CD, Kartsev A, Xiao J, Lindenberg A, Rajpurohit S, Tan LZ, Ogitsu T, Correa AA. Inq, a Modern GPU-Accelerated Computational Framework for (Time-Dependent) Density Functional Theory. J Chem Theory Comput 2021; 17:7447-7467. [PMID: 34726888 DOI: 10.1021/acs.jctc.1c00562] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present inq, a new implementation of density functional theory (DFT) and time-dependent DFT (TDDFT) written from scratch to work on graphic processing units (GPUs). Besides GPU support, inq makes use of modern code design features and takes advantage of newly available hardware. By designing the code around algorithms, rather than against specific implementations and numerical libraries, we aim to provide a concise and modular code. The result is a fairly complete DFT/TDDFT implementation in roughly 12 000 lines of open-source C++ code representing a modular platform for community-driven application development on emerging high-performance computing architectures.
Collapse
Affiliation(s)
- Xavier Andrade
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94551, United States
| | - Chaitanya Das Pemmaraju
- Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Alexey Kartsev
- Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Jun Xiao
- Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Aaron Lindenberg
- Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Sangeeta Rajpurohit
- The Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Liang Z Tan
- The Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Tadashi Ogitsu
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94551, United States
| | - Alfredo A Correa
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94551, United States
| |
Collapse
|
172
|
Ollitrault PJ, Miessen A, Tavernelli I. Molecular Quantum Dynamics: A Quantum Computing Perspective. Acc Chem Res 2021; 54:4229-4238. [PMID: 34787398 DOI: 10.1021/acs.accounts.1c00514] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
ConspectusSimulating molecular dynamics (MD) within a comprehensive quantum framework has been a long-standing challenge in computational chemistry. An exponential scaling of computational cost renders solving the time dependent Schrödinger equation (TDSE) of a molecular Hamiltonian, including both electronic and nuclear degrees of freedom (DOFs), as well as their couplings, infeasible for more than a few DOFs. In the Born-Oppenheimer (BO), or adiabatic, picture, electronic and nuclear parts of the wave function are decoupled and treated separately. Within this framework, the nuclear wave function evolves along potential energy surfaces (PESs) computed as solutions to the electronic Schrödinger equation parametrized in the nuclear DOFs. This approximation, together with increasingly elaborate numerical approaches to solve the nuclear time dependent Schrödinger equation (TDSE), enabled the treatment of up to a few dozens of degrees of freedom (DOFs). However, for particular applications, such as photochemistry, the BO approximation breaks down. In this regime of non-adiabatic dynamics, solving the full molecular problem including electron-nuclear couplings becomes essential, further increasing the complexity of the numerical solution. Although valuable methods such as multiconfigurational time-dependent Hartree (MCTDH) have been proposed for the solution of the coupled electron-nuclear dynamics, they remain hampered by an exponential scaling in the number of nuclear DOFs and by the difficulty of finding universal variational forms.In this Account, we present a perspective on novel quantum computational algorithms, aiming to alleviate the exponential scaling inherent to the simulation of many-body quantum dynamics. In particular, we focus on the derivation and application of quantum algorithms for adiabatic and non-adiabatic quantum dynamics, which include efficient approaches for the calculation of the BO potential energy surfaces (PESs). Thereafter, we study the time-evolution of a model system consisting of two coupled PESs in first and second quantization. In a first application, we discuss a recently introduced quantum algorithm for the evolution of a wavepacket in first quantization and exploit the potential quantum advantage of mapping its spatial grid representation to logarithmically many qubits. For the second demonstration, we move to the second quantization framework and review the scaling properties of two alternative time-evolution algorithms, namely, a variational quantum algorithm (VQA) (based on the McLachlan variational principle) and conventional Trotter-type evolution (based on a Lie-Trotter-Suzuki formula). Both methods clearly demonstrate the potential of quantum algorithms and their favorable scaling compared to the available classical approaches. However, a clear demonstration of quantum advantage in the context of molecular quantum dynamics may require the implementation of these algorithms in fault-tolerant quantum computers, while their application in near-term, noisy quantum devices is still unclear and deserves further investigation.
Collapse
|
173
|
Abstract
Photochemical reactions are increasingly being used for chemical and materials synthesis, for example, in photoredox catalysis, and generally involve photoexcitation of molecular chromophores dissolved in a liquid solvent. The choice of solvent influences the outcomes of the photochemistry because solute-solvent interactions modify the energies of and crossings between electronic states of the chromophores, and they affect the evolving structures of the photoexcited molecules. Ultrafast laser spectroscopy methods with femtosecond to picosecond time resolution can resolve the dynamics of these photoexcited molecules as they undergo structural and electronic changes, relax back to the ground state, dissipate their excess internal energy to the surrounding solvent, or undergo photochemical reactions. In this Account, we illustrate how experimental studies using ultrafast lasers can reveal the influences that different solvents or cosolutes exert on the photoinduced nonadiabatic dynamics of internal conversion and intersystem crossing in nonradiative relaxation pathways. Although the environment surrounding a solute molecule is rapidly changing, with fluctuations in the coordination to neighboring solvent molecules occurring on femtosecond or picosecond time scales, we show that it is possible to photoexcite selectively only those molecular chromophores transiently experiencing specific solute-solvent interactions such as intermolecular hydrogen bonding.The effects of different solvation environments on the photodynamics are illustrated using four selected examples of photochemical processes in which the solvent has a marked effect on the outcomes. We first consider two aromatic carbonyl compounds, benzophenone and acetophenone, which are known to undergo fast intersystem crossing to populate the first excited triplet state on time scales of a few picoseconds. We show that the nonadiabatic excited-state dynamics are modified by transient hydrogen bonding of the carbonyl group to a protic solvent or by coordination to a metal cation cosolute. We then examine how different solvents modify the competition between two alternative relaxation pathways in a photoexcited UVA-sunscreen molecule, diethylamino hydroxybenzoyl hexyl benzoate (DHHB). This relaxation back to the ground electronic state is an essential part of the effective operation of the sunscreen compound, but the dynamics are sensitive to the surrounding environment. Finally, we consider how solvents of different polarity affect the energies and lifetimes of excited states with locally excited or charge-transfer character in heterocyclic organic compounds used as excited-state electron donors for photoredox catalysis. With these and other examples, we seek to develop a molecular level understanding of how the choice of solution environment might be used to control the outcomes of photochemical reactions.
Collapse
Affiliation(s)
- Ravi Kumar Venkatraman
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Andrew J. Orr-Ewing
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
174
|
Esch MP, Levine BG. An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics. J Chem Phys 2021; 155:214101. [PMID: 34879667 DOI: 10.1063/5.0070686] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.
Collapse
Affiliation(s)
- Michael P Esch
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
175
|
Zhang J, Peng J, Hu D, Lan Z. Investigation of nonadiabatic dynamics in the photolysis of methyl nitrate (CH 3ONO 2) by on-the-fly surface hopping simulation. Phys Chem Chem Phys 2021; 23:25597-25611. [PMID: 34546246 DOI: 10.1039/d1cp03226g] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The photolysis mechanism of methyl nitrate (CH3ONO2) was studied using the on-the-fly surface hopping dynamics at the XMS-CASPT2 level. Several critical geometries, including electronic state minima and conical intersections, were obtained, which play essential roles in the nonadiabatic dynamics of CH3ONO2. The ultrafast nonadiabatic decay dynamics to the ground state were simulated, which gives a proper explanation on the broad and structureless absorption spectra of CH3ONO2. The photodissociation channels, including CH3O + NO2, CH3O + NO + O, and others, as well as their branching ratios, were identified. When the dynamics starts from the lowest two electronic states (S1 and S2), the CH3O + NO2 channel is the dominant photolysis pathway, although we observed the minor contributions of other channels. In contrast, when the trajectories start from the third excited state S3, both CH3O + NO2 and CH3O + NO + O channels become important. Here the CH3O-NO2 bond dissociation takes place first, and then for some trajectories, the N-O bond of the NO2 part breaks successively. The quasi-degeneracy of electronic states may exist in the dissociation limits of both CH3O + NO2 and CH3O + NO + O channels. The current work provides valuable information in the understanding of experimental findings of the wavelength-dependent photolysis mechanism of CH3ONO2.
Collapse
Affiliation(s)
- Juanjuan Zhang
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
176
|
Guan Y, Xie C, Yarkony DR, Guo H. High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics. Phys Chem Chem Phys 2021; 23:24962-24983. [PMID: 34473156 DOI: 10.1039/d1cp03008f] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Nonadiabatic dynamics, which goes beyond the Born-Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China.
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131, USA.
| |
Collapse
|
177
|
Tao Z, Yu Q, Roy S, Hammes-Schiffer S. Direct Dynamics with Nuclear-Electronic Orbital Density Functional Theory. Acc Chem Res 2021; 54:4131-4141. [PMID: 34726895 DOI: 10.1021/acs.accounts.1c00516] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Direct dynamics simulations of chemical reactions typically require the selection of a method for generating the potential energy surfaces and a method for the dynamical propagation of the nuclei on these surfaces. The nuclear-electronic orbital (NEO) framework avoids this Born-Oppenheimer separation by treating specified nuclei on the same level as the electrons with wave function methods or density functional theory (DFT). The NEO approach is particularly applicable to proton, hydride, and proton-coupled electron transfer reactions, where the transferring proton(s) and all electrons are treated quantum mechanically. In this manner, the zero-point energy, density delocalization, and anharmonicity of the transferring protons are inherently and efficiently included in the energies, optimized geometries, and dynamics.This Account describes how various NEO methods can be used for direct dynamics simulations on electron-proton vibronic surfaces. The strengths and limitations of these approaches are discussed, and illustrative examples are presented. The NEO-DFT method can be used to simulate chemical reactions on the ground state vibronic surface, as illustrated by the application to hydride transfer in C4H9+. The NEO multistate DFT (NEO-MSDFT) method is useful for simulating ground state reactions in which the proton density becomes bilobal during the dynamics, a characteristic of hydrogen tunneling, as illustrated by proton transfer in malonaldehyde. The NEO time-dependent DFT (NEO-TDDFT) method produces excited electronic, vibrational, and vibronic surfaces. The application of linear-response NEO-TDDFT to H2 and H3+, as well as the partially and fully deuterated counterparts, shows that this approach produces accurate fundamental vibrational excitation energies when all nuclei and all electrons are treated quantum mechanically. Moreover, when only specified nuclei are treated quantum mechanically, this approach can be used to optimize geometries on excited state vibronic surfaces, as illustrated by photoinduced single and double proton transfer systems, and to conduct adiabatic dynamics on these surfaces. The real-time NEO-TDDFT method provides an alternative approach for simulating nonequilibrium nuclear-electronic dynamics of such systems. These various NEO methods can be combined with nonadiabatic dynamics methods such as Ehrenfest and surface hopping dynamics to include the nonadiabatic effects between the quantum and classical subsystems. The real-time NEO-TDDFT Ehrenfest dynamics simulation of excited state intramolecular proton transfer in o-hydroxybenzaldehyde illustrates the power of this type of combined approach. The field of multicomponent quantum chemistry is in the early stages, and the methods discussed herein provide the foundation for a wide range of promising future directions to be explored. An appealing future direction is the expansion of the real-time NEO-TDDFT method to describe the dynamics of all nuclei and electrons on the same level. Direct dynamics simulations using NEO wave function methods such as equation-of-motion coupled cluster or multiconfigurational approaches are also attractive but computationally expensive options. The further development of NEO direct dynamics methods will enable the simulation of the nuclear-electronic dynamics for a vast array of chemical and biological processes that extend beyond the Born-Oppenheimer approximation.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Saswata Roy
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
178
|
Ibele LM, Curchod BFE. Dynamics near a conical intersection-A diabolical compromise for the approximations of ab initio multiple spawning. J Chem Phys 2021; 155:174119. [PMID: 34742188 DOI: 10.1063/5.0071376] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Full multiple spawning (FMS) offers an exciting framework for the development of strategies to simulate the excited-state dynamics of molecular systems. FMS proposes to depict the dynamics of nuclear wavepackets by using a growing set of traveling multidimensional Gaussian functions called trajectory basis functions (TBFs). Perhaps the most recognized method emanating from FMS is the so-called ab initio multiple spawning (AIMS). In AIMS, the couplings between TBFs-in principle exact in FMS-are approximated to allow for the on-the-fly evaluation of required electronic-structure quantities. In addition, AIMS proposes to neglect the so-called second-order nonadiabatic couplings and the diagonal Born-Oppenheimer corrections. While AIMS has been applied successfully to simulate the nonadiabatic dynamics of numerous complex molecules, the direct influence of these missing or approximated terms on the nonadiabatic dynamics when approaching and crossing a conical intersection remains unknown to date. It is also unclear how AIMS could incorporate geometric-phase effects in the vicinity of a conical intersection. In this work, we assess the performance of AIMS in describing the nonadiabatic dynamics through a conical intersection for three two-dimensional, two-state systems that mimic the excited-state dynamics of bis(methylene)adamantyl, butatriene cation, and pyrazine. The population traces and nuclear density dynamics are compared with numerically exact quantum dynamics and trajectory surface hopping results. We find that AIMS offers a qualitatively correct description of the dynamics through a conical intersection for the three model systems. However, any attempt at improving the AIMS results by accounting for the originally neglected second-order nonadiabatic contributions appears to be stymied by the hermiticity requirement of the AIMS Hamiltonian and the independent first-generation approximation.
Collapse
Affiliation(s)
- Lea M Ibele
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
179
|
Mukherjee S, Varganov SA. Intersystem crossing and internal conversion dynamics with GAIMS-TeraChem: Excited state relaxation in 2-cyclopentenone. J Chem Phys 2021; 155:174107. [PMID: 34742200 DOI: 10.1063/5.0068040] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Excited states relaxation in complex molecules often involves two types of nonradiative transitions, internal conversion (IC) and intersystem crossing (ISC). In the situations when the timescales of IC and ISC are comparable, an interplay between these two types of transitions can lead to complex nonadiabatic dynamics on multiple electronic states of different characters and spin multiplicities. We demonstrate that the generalized ab initio multiple spawning (GAIMS) method interfaced with the fast graphics processing unit-based TeraChem electronic structure code can be used to model such nonadiabatic dynamics involving both the IC and ISC transitions in molecules of moderate size. We carried out 1500 fs GAIMS simulations leading to the creation of up to 2500 trajectory basis functions to study the excited states relaxation in 2-cyclopentenone. After a vertical excitation from the ground state to the bright S2 state, the molecule quickly relaxes to the S1 state via conical intersection. The following relaxation proceeds along two competing pathways: one involves IC to the ground state, and the other is dominated by ISC to the low-lying triplet states. The time constants describing the population transfer between the six lowest singlet and triplet states predicted by the GAIMS dynamics are in good agreement with the characteristic times of IC and ISC obtained from the analysis of the time-resolved photoelectron spectrum.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA
| |
Collapse
|
180
|
Wang C, Ulusoy IS, Aebersold LE, Wilson AK. Multi-configuration electron-nuclear dynamics: An open-shell approach. J Chem Phys 2021; 155:154103. [PMID: 34686063 DOI: 10.1063/5.0063478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.
Collapse
Affiliation(s)
- Cong Wang
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Lucas E Aebersold
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| |
Collapse
|
181
|
Lin K, Peng J, Gu FL, Lan Z. Simulation of Open Quantum Dynamics with Bootstrap-Based Long Short-Term Memory Recurrent Neural Network. J Phys Chem Lett 2021; 12:10225-10234. [PMID: 34647736 DOI: 10.1021/acs.jpclett.1c02672] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The recurrent neural network with the long short-term memory cell (LSTM-NN) is employed to simulate the long-time dynamics of open quantum systems. The bootstrap method is applied in the LSTM-NN construction and prediction, which provides a Monte Carlo estimation of a forecasting confidence interval. Within this approach, a large number of LSTM-NNs are constructed by resampling time-series sequences that were obtained from the early stage quantum evolution given by numerically exact multilayer multiconfigurational time-dependent Hartree method. The built LSTM-NN ensemble is used for the reliable propagation of the long-time quantum dynamics, and the simulated result is highly consistent with the exact evolution. The forecasting uncertainty that partially reflects the reliability of the LSTM-NN prediction is also given. This demonstrates the bootstrap-based LSTM-NN approach is a practical and powerful tool to propagate the long-time quantum dynamics of open systems with high accuracy and low computational cost.
Collapse
Affiliation(s)
- Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- 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
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. 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
| |
Collapse
|
182
|
Lee JWL, Tikhonov DS, Chopra P, Maclot S, Steber AL, Gruet S, Allum F, Boll R, Cheng X, Düsterer S, Erk B, Garg D, He L, Heathcote D, Johny M, Kazemi MM, Köckert H, Lahl J, Lemmens AK, Loru D, Mason R, Müller E, Mullins T, Olshin P, Passow C, Peschel J, Ramm D, Rompotis D, Schirmel N, Trippel S, Wiese J, Ziaee F, Bari S, Burt M, Küpper J, Rijs AM, Rolles D, Techert S, Eng-Johnsson P, Brouard M, Vallance C, Manschwetus B, Schnell M. Time-resolved relaxation and fragmentation of polycyclic aromatic hydrocarbons investigated in the ultrafast XUV-IR regime. Nat Commun 2021; 12:6107. [PMID: 34671016 PMCID: PMC8528970 DOI: 10.1038/s41467-021-26193-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2020] [Accepted: 09/17/2021] [Indexed: 11/18/2022] Open
Abstract
Polycyclic aromatic hydrocarbons (PAHs) play an important role in interstellar chemistry and are subject to high energy photons that can induce excitation, ionization, and fragmentation. Previous studies have demonstrated electronic relaxation of parent PAH monocations over 10-100 femtoseconds as a result of beyond-Born-Oppenheimer coupling between the electronic and nuclear dynamics. Here, we investigate three PAH molecules: fluorene, phenanthrene, and pyrene, using ultrafast XUV and IR laser pulses. Simultaneous measurements of the ion yields, ion momenta, and electron momenta as a function of laser pulse delay allow a detailed insight into the various molecular processes. We report relaxation times for the electronically excited PAH*, PAH+* and PAH2+* states, and show the time-dependent conversion between fragmentation pathways. Additionally, using recoil-frame covariance analysis between ion images, we demonstrate that the dissociation of the PAH2+ ions favors reaction pathways involving two-body breakup and/or loss of neutral fragments totaling an even number of carbon atoms.
Collapse
Affiliation(s)
- J. W. L. Lee
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - D. S. Tikhonov
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9764.c0000 0001 2153 9986Institute of Physical Chemistry, Christian-Albrechts-Universität zu Kiel, Kiel, Germany
| | - P. Chopra
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9764.c0000 0001 2153 9986Institute of Physical Chemistry, Christian-Albrechts-Universität zu Kiel, Kiel, Germany
| | - S. Maclot
- grid.4514.40000 0001 0930 2361Department of Physics, Lund University, Lund, Sweden ,grid.8761.80000 0000 9919 9582Physics Department, University of Gothenburg, Gothenburg, Sweden
| | - A. L. Steber
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9764.c0000 0001 2153 9986Institute of Physical Chemistry, Christian-Albrechts-Universität zu Kiel, Kiel, Germany ,grid.9026.d0000 0001 2287 2617Center for Ultrafast Imaging, Universität Hamburg, Hamburg, Germany
| | - S. Gruet
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - F. Allum
- grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - R. Boll
- grid.434729.f0000 0004 0590 2900European XFEL, Schenefeld, Germany
| | - X. Cheng
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - S. Düsterer
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - B. Erk
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - D. Garg
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9026.d0000 0001 2287 2617Department of Physics, Universität Hamburg, Hamburg, Germany
| | - L. He
- grid.7683.a0000 0004 0492 0453Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - D. Heathcote
- grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - M. Johny
- grid.7683.a0000 0004 0492 0453Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - M. M. Kazemi
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - H. Köckert
- grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - J. Lahl
- grid.4514.40000 0001 0930 2361Department of Physics, Lund University, Lund, Sweden
| | - A. K. Lemmens
- grid.5590.90000000122931605Radboud University, FELIX Laboratory, Nijmegen, The Netherlands ,grid.7177.60000000084992262Van’t Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - D. Loru
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9764.c0000 0001 2153 9986Institute of Physical Chemistry, Christian-Albrechts-Universität zu Kiel, Kiel, Germany
| | - R. Mason
- grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - E. Müller
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - T. Mullins
- grid.7683.a0000 0004 0492 0453Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - P. Olshin
- grid.15447.330000 0001 2289 6897Saint Petersburg State University, Saint Petersburg, Russia
| | - C. Passow
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - J. Peschel
- grid.4514.40000 0001 0930 2361Department of Physics, Lund University, Lund, Sweden
| | - D. Ramm
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - D. Rompotis
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.434729.f0000 0004 0590 2900European XFEL, Schenefeld, Germany
| | - N. Schirmel
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - S. Trippel
- grid.9026.d0000 0001 2287 2617Center for Ultrafast Imaging, Universität Hamburg, Hamburg, Germany ,grid.7683.a0000 0004 0492 0453Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - J. Wiese
- grid.7683.a0000 0004 0492 0453Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9026.d0000 0001 2287 2617Department of Chemistry, Universität Hamburg, Hamburg, Germany
| | - F. Ziaee
- grid.36567.310000 0001 0737 1259J.R. Macdonald Laboratory, Department of Physics, Kansas State University, Manhattan, KS USA
| | - S. Bari
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - M. Burt
- grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - J. Küpper
- grid.9026.d0000 0001 2287 2617Center for Ultrafast Imaging, Universität Hamburg, Hamburg, Germany ,grid.9026.d0000 0001 2287 2617Department of Physics, Universität Hamburg, Hamburg, Germany ,grid.7683.a0000 0004 0492 0453Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9026.d0000 0001 2287 2617Department of Chemistry, Universität Hamburg, Hamburg, Germany
| | - A. M. Rijs
- grid.5590.90000000122931605Radboud University, FELIX Laboratory, Nijmegen, The Netherlands ,grid.12380.380000 0004 1754 9227Division of BioAnalytical Chemistry, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - D. Rolles
- grid.36567.310000 0001 0737 1259J.R. Macdonald Laboratory, Department of Physics, Kansas State University, Manhattan, KS USA
| | - S. Techert
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.7450.60000 0001 2364 4210Institute for X-Ray Physics, Georg-August-Universität, Göttingen, Germany
| | - P. Eng-Johnsson
- grid.4514.40000 0001 0930 2361Department of Physics, Lund University, Lund, Sweden
| | - M. Brouard
- grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - C. Vallance
- grid.4991.50000 0004 1936 8948The Chemistry Research Laboratory, University of Oxford, Oxford, United Kingdom
| | - B. Manschwetus
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - M. Schnell
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9764.c0000 0001 2153 9986Institute of Physical Chemistry, Christian-Albrechts-Universität zu Kiel, Kiel, Germany
| |
Collapse
|
183
|
Zobel JP, Heindl M, Plasser F, Mai S, González L. Surface Hopping Dynamics on Vibronic Coupling Models. Acc Chem Res 2021; 54:3760-3771. [PMID: 34570472 PMCID: PMC8529708 DOI: 10.1021/acs.accounts.1c00485] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The simulation of photoinduced non-adiabatic dynamics is of great
relevance in many scientific disciplines, ranging from physics and
materials science to chemistry and biology. Upon light irradiation,
different relaxation processes take place in which electronic and
nuclear motion are intimately coupled. These are best described by
the time-dependent molecular Schrödinger equation, but its
solution poses fundamental practical challenges to contemporary theoretical
chemistry. Two widely used and complementary approaches to this problem
are multiconfigurational time-dependent Hartree (MCTDH) and trajectory
surface hopping (SH). MCTDH is an accurate fully quantum-mechanical
technique but often is feasible only in reduced dimensionality, in
combination with approximate vibronic coupling (VC) Hamiltonians,
or both (i.e., reduced-dimensional VC potentials). In contrast, SH
is a quantum–classical technique that neglects most nuclear
quantum effects but allows nuclear dynamics in full dimensionality
by calculating potential energy surfaces on the fly. If nuclear quantum
effects do not play a central role and a linear VC (LVC) Hamiltonian
is appropriate—e.g., for stiff molecules that generally keep
their conformation in the excited state—then it seems advantageous
to combine the efficient LVC and SH techniques. In this Account, we
describe how surface hopping based on an LVC Hamiltonian (SH/LVC)—as
recently implemented in the SHARC surface hopping package—can
provide an economical and automated approach to simulate non-adiabatic
dynamics. First, we illustrate the potential of SH/LVC in a number
of showcases, including intersystem crossing in SO2, intra-Rydberg
dynamics in acetone, and several photophysical studies on large transition-metal
complexes, which would be much more demanding or impossible to perform
with other methods. While all of the applications provide very useful
insights into light-induced phenomena, they also hint at difficulties
faced by the SH/LVC methodology that need to be addressed in the future.
Second, we contend that the SH/LVC approach can be useful to benchmark
SH itself. By the use of the same (LVC) potentials as MCTDH calculations
have employed for decades and by relying on the efficiency of SH/LVC,
it is possible to directly compare multiple SH test calculations with
a MCTDH reference and ponder the accuracy of various correction algorithms
behind the SH methodology, such as decoherence corrections or momentum
rescaling schemes. Third, we demonstrate how the efficiency of SH/LVC
can also be exploited to identify essential nuclear and electronic
degrees of freedom to be employed in more accurate MCTDH calculations.
Lastly, we show that SH/LVC is able to advance the development of
SH protocols that can describe nuclear dynamics including explicit
laser fields—a very challenging endeavor for trajectory-based
schemes. To end, this Account compiles the typical costs of contemporary
SH simulations, evidencing the great advantages of using parametrized
potentials. The LVC model is a sleeping beauty that, kissed by SH,
is fueling the field of excited-state molecular dynamics. We hope
that this Account will stimulate future research in this direction,
leveraging the advantages of the SH/VC schemes to larger extents and
extending their applicability to uncharted territories.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| | - Moritz Heindl
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| | - Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| |
Collapse
|
184
|
Dunnett AJ, Gowland D, Isborn CM, Chin AW, Zuehlsdorff TJ. Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue. J Chem Phys 2021; 155:144112. [PMID: 34654312 DOI: 10.1063/5.0062950] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute-solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute-solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.
Collapse
Affiliation(s)
- Angus J Dunnett
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Duncan Gowland
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
185
|
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
|
186
|
Abstract
In this paper, we discuss coupled-trajectory schemes for molecular-dynamics simulations of excited-state processes. New coupled-trajectory strategies to capture decoherence effects, revival of coherence and nonadiabatic interferences in long-time dynamics are proposed, and compared to independent-trajectory schemes. The working framework is provided by the exact factorization of the electron-nuclear wave function, and it exploits ideas emanating from various surface-hopping schemes. The new coupled-trajectory algorithms are tested on a one-dimensional two-state system using different model parameters which allow one to induce different dynamics. The benchmark is provided by the numerically exact solution of the time-dependent Schrödinger equation.
Collapse
Affiliation(s)
- Carlotta Pieroni
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France.,Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Federica Agostini
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
187
|
Sršeň Š, Slavíček P. Optimal Representation of the Nuclear Ensemble: Application to Electronic Spectroscopy. J Chem Theory Comput 2021; 17:6395-6404. [PMID: 34542278 DOI: 10.1021/acs.jctc.1c00749] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Nuclear densities are frequently represented by an ensemble of nuclear configurations or points in the phase space in various contexts of molecular simulations. The size of the ensemble directly affects the accuracy and computational cost of subsequent calculations of observable quantities. In the present work, we address the question of how many configurations do we need and how to select them most efficiently. We focus on the nuclear ensemble method in the context of electronic spectroscopy, where thousands of sampled configurations are usually needed for sufficiently converged spectra. The proposed representative sampling technique allows for a dramatic reduction of the sample size. By using an exploratory method, we model the density from a large sample in the space of transition properties. The representative subset of nuclear configurations is optimized by minimizing its Kullback-Leibler divergence to the full density with simulated annealing. High-level calculations are then performed only for the selected subset of configurations. We tested the algorithm on electronic absorption spectra of three molecules: (E)-azobenzene, the simplest Criegee intermediate, and hydrated nitrate anion. Typically, dozens of nuclear configurations provided sufficiently accurate spectra. A strongly forbidden transition of the nitrate anion presented the most challenging case due to rare geometries with disproportionately high transition intensities. This problematic case was easily diagnosed within the present approach. We also discuss various exploratory methods and a possible extension to dynamical simulations.
Collapse
Affiliation(s)
- Štěpán Sršeň
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 16628, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 16628, Czech Republic
| |
Collapse
|
188
|
Sun J, Sasmal S, Vendrell O. A bosonic perspective on the classical mapping of fermionic quantum dynamics. J Chem Phys 2021; 155:134110. [PMID: 34624967 DOI: 10.1063/5.0066740] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider the application of the original Meyer-Miller (MM) Hamiltonian to mapping fermionic quantum dynamics to classical equations of motion. Non-interacting fermionic and bosonic systems share the same one-body density dynamics when evolving from the same initial many-body state. The MM classical mapping is exact for non-interacting bosons, and therefore, it yields the exact time-dependent one-body density for non-interacting fermions as well. Starting from this observation, the MM mapping is compared to different mappings specific for fermionic systems, namely, the spin mapping with and without including a Jordan-Wigner transformation and the Li-Miller mapping (LMM). For non-interacting systems, the inclusion of fermionic anti-symmetry through the Jordan-Wigner transform does not lead to any improvement in the performance of the mappings, and instead, it worsens the classical description. For an interacting impurity model and for models of excitonic energy transfer, the MM and LMM mappings perform similarly, and in some cases, the former outperforms the latter when compared to a full quantum description. The classical mappings are able to capture interference effects, both constructive and destructive, that originate from equivalent energy transfer pathways in the models.
Collapse
Affiliation(s)
- Jing Sun
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | - Sudip Sasmal
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| |
Collapse
|
189
|
Hu D, Peng J, Chen L, Gelin MF, Lan Z. Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations. J Phys Chem Lett 2021; 12:9710-9719. [PMID: 34590858 DOI: 10.1021/acs.jpclett.1c02640] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.
Collapse
Affiliation(s)
- Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- School of Environment, South China Normal University, Guangzhou 510006, China
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
190
|
Lin K, Hu D, Peng J, Xu C, Gu FL, Lan Z. Prediction of the excited-state reaction channels in photo-induced processes of nitrofurantoin using first-principle calculations and dynamics simulations. CHEMOSPHERE 2021; 281:130831. [PMID: 34289597 DOI: 10.1016/j.chemosphere.2021.130831] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2021] [Revised: 03/25/2021] [Accepted: 05/06/2021] [Indexed: 06/13/2023]
Abstract
The understanding of the photochemistry of antibiotic compounds is important because it gives the direct information on the possible environmental pollution caused by them. Due to their large size, the theoretical studies of their excited-state reactions are rather challenging. In current work, we combined the on-the-fly trajectory surface-hopping dynamics, conical-intersection optimizations and excited-state pathway calculations to study the photochemistry of the trans-isomer of nitrofurantoin, a widely-used drug to treat the urinary tract infections. The dynamics-then-pathway approach was taken. First the trajectory surface hopping dynamics at the state-averaged complete-active-space self-consistent-field (SA-CASSCF) level with small active space and small basis sets were run. Second, the minimum-energy conical-intersection optimizations were performed. Finally the excited pathways from the Frank-Condon region to different reaction channels were built at the multi-state multi-reference second-order perturbation (MS-CASPT2) level with large active space and large basis set. Several possible channels responsible for the photo-induced reaction mechanism of the trans-nitrofurantoin were obtained, including the cleavage of the NO bond of the NO2 moiety, the photoisomerization at the central CN bond, and other internal conversion channels. Our findings give some preliminary explanations on available experimental observations. It is also demonstrates that the current theoretical approach is a powerful tool to explore the excited-state reactions in the photochemistry of media-sized or large-sized drug compounds.
Collapse
Affiliation(s)
- Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou, 510006, PR China
| | - Deping Hu
- 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, PR China
| | - Jiawei Peng
- 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, PR China
| | - Chao Xu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou, 510006, PR China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou, 510006, PR 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, PR China.
| |
Collapse
|
191
|
Xiong R, Hu R, Zhang Y, Yang X, Lin P, Wen C, Sa B, Sun Z. Computational discovery of PtS 2/GaSe van der Waals heterostructure for solar energy applications. Phys Chem Chem Phys 2021; 23:20163-20173. [PMID: 34551041 DOI: 10.1039/d1cp02436a] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
2D van der Waals (vdW) heterostructures as potential materials for solar energy-related applications have been brought to the forefront for researchers. Here, by employing first-principles calculations, we proposed that the PtS2/GaSe vdW heterostructure is a distinguished candidate for photocatalytic water splitting and solar cells. It is shown that the PtS2/GaSe heterostructure exhibits high thermal stability with an indirect band gap of 1.81 eV. We further highlighted the strain induced type-V to type-II band alignment transitions and band gap variations in PtS2/GaSe heterostructures. More importantly, the outstanding absorption coefficients in the visible light region and high carrier mobility further guarantee the photo energy conversion efficiency of PtS2/GaSe heterostructures. Interestingly, the natural type-V band alignments of PtS2/GaSe heterostructures are appropriate for the redox potential of water. On the other hand, the power conversion efficiency of ZnO/(PtS2/GaSe heterostructure)/CIGS (copper indium gallium diselenide) solar cells can achieve ∼17.4%, which can be further optimized up to ∼18.5% by increasing the CIGS thickness. Our present study paves the way for facilitating the potential application of vdW heterostructures as a promising photocatalyst for water splitting as well as the buffer layer for solar cells.
Collapse
Affiliation(s)
- Rui Xiong
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, P. R. China.
| | - Rong Hu
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, P. R. China.
| | - Yinggan Zhang
- College of Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen 361005, P. R. China
| | - Xuhui Yang
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, P. R. China.
| | - Peng Lin
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, P. R. China.
| | - Cuilian Wen
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, P. R. China.
| | - Baisheng Sa
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, P. R. China.
| | - Zhimei Sun
- School of Materials Science and Engineering, and Center for Integrated Computational Materials Science, International Research Institute for Multidisciplinary Science, Beihang University, Beijing 100191, P. R. China.
| |
Collapse
|
192
|
Chen Z, Yang J. Nucleus-electron correlation revising molecular bonding fingerprints from the exact wavefunction factorization. J Chem Phys 2021; 155:104111. [PMID: 34525813 DOI: 10.1063/5.0056773] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
We present a novel theory and implementation for computing coupled electronic and quantal nuclear subsystems on a single potential energy surface, moving beyond the standard Born-Oppenheimer (BO) separation of nuclei and electrons. We formulate an exact self-consistent nucleus-electron embedding potential from the single product molecular wavefunction and demonstrate that the fundamental behavior of the correlated nucleus-electron can be computed for mean-field electrons that are responsive to a quantal anharmonic vibration of selected nuclei in a discrete variable representation. Geometric gauge choices are discussed and necessary for formulating energy invariant biorthogonal electronic equations. Our method is further applied to characterize vibrationally averaged molecular bonding properties of molecular energetics, bond lengths, and protonic and electron densities. Moreover, post-Hartree-Fock electron correlation can be conveniently computed on the basis of nucleus-electron coupled molecular orbitals, as demonstrated for correlated models of second-order Møllet-Plesset perturbation and full configuration interaction theories. Our approach not only accurately quantifies non-classical nucleus-electron couplings for revising molecular bonding properties but also provides an alternative time-independent approach for deploying non-BO molecular quantum chemistry.
Collapse
Affiliation(s)
- Ziyong Chen
- Department of Chemistry, The University of Hong Kong, Hong Kong, People's Republic of China
| | - Jun Yang
- Department of Chemistry, The University of Hong Kong, Hong Kong, People's Republic of China
| |
Collapse
|
193
|
Merritt ICD, Jacquemin D, Vacher M. Attochemistry: Is Controlling Electrons the Future of Photochemistry? J Phys Chem Lett 2021; 12:8404-8415. [PMID: 34436903 DOI: 10.1021/acs.jpclett.1c02016] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Controlling matter with light has always been a great challenge, leading to the ever-expanding field of photochemistry. In addition, since the first generation of light pulses of attosecond (1 as = 10-18 s) duration, a great deal of effort has been devoted to observing and controlling electrons on their intrinsic time scale. Because of their short duration, attosecond pulses have a large spectral bandwidth populating several electronically excited states in a coherent manner, i.e., an electronic wavepacket. Because of interference, such a wavepacket has a new electronic distribution implying a potentially different and totally new reactivity as compared to traditional photochemistry, leading to the novel concept of "attochemistry". This nascent field requires the support of theory right from the start. In this Perspective, we discuss the opportunities offered by attochemistry, the related challenges, and the current and future state-of-the-art developments in theoretical chemistry needed to model it accurately.
Collapse
Affiliation(s)
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Morgane Vacher
- Université de Nantes, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| |
Collapse
|
194
|
Jones CM, List NH, Martínez TJ. Resolving the ultrafast dynamics of the anionic green fluorescent protein chromophore in water. Chem Sci 2021; 12:11347-11363. [PMID: 34667545 PMCID: PMC8447926 DOI: 10.1039/d1sc02508b] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 07/13/2021] [Indexed: 11/21/2022] Open
Abstract
The chromophore of the green fluorescent protein (GFP) is critical for probing environmental influences on fluorescent protein behavior. Using the aqueous system as a bridge between the unconfined vacuum system and a constricting protein scaffold, we investigate the steric and electronic effects of the environment on the photodynamical behavior of the chromophore. Specifically, we apply ab initio multiple spawning to simulate five picoseconds of nonadiabatic dynamics after photoexcitation, resolving the excited-state pathways responsible for internal conversion in the aqueous chromophore. We identify an ultrafast pathway that proceeds through a short-lived (sub-picosecond) imidazolinone-twisted (I-twisted) species and a slower (several picoseconds) channel that proceeds through a long-lived phenolate-twisted (P-twisted) intermediate. The molecule navigates the non-equilibrium energy landscape via an aborted hula-twist-like motion toward the one-bond-flip dominated conical intersection seams, as opposed to following the pure one-bond-flip paths proposed by the excited-state equilibrium picture. We interpret our simulations in the context of time-resolved fluorescence experiments, which use short- and long-time components to describe the fluorescence decay of the aqueous GFP chromophore. Our results suggest that the longer time component is caused by an energetically uphill approach to the P-twisted intersection seam rather than an excited-state barrier to reach the twisted intramolecular charge-transfer species. Irrespective of the location of the nonadiabatic population events, the twisted intersection seams are inefficient at facilitating isomerization in aqueous solution. The disordered and homogeneous nature of the aqueous solvent environment facilitates non-selective stabilization with respect to I- and P-twisted species, offering an important foundation for understanding the consequences of selective stabilization in heterogeneous and rigid protein environments.
Collapse
Affiliation(s)
- Chey M Jones
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| |
Collapse
|
195
|
Weight BM, Mandal A, Huo P. Ab initio symmetric quasi-classical approach to investigate molecular Tully models. J Chem Phys 2021; 155:084106. [PMID: 34470343 DOI: 10.1063/5.0061934] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.
Collapse
Affiliation(s)
- Braden M Weight
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Arkajit Mandal
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
196
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 162] [Impact Index Per Article: 54.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
197
|
Keith JA, Vassilev-Galindo V, Cheng B, Chmiela S, Gastegger M, Müller KR, Tkatchenko A. Combining Machine Learning and Computational Chemistry for Predictive Insights Into Chemical Systems. Chem Rev 2021; 121:9816-9872. [PMID: 34232033 PMCID: PMC8391798 DOI: 10.1021/acs.chemrev.1c00107] [Citation(s) in RCA: 190] [Impact Index Per Article: 63.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Indexed: 12/23/2022]
Abstract
Machine learning models are poised to make a transformative impact on chemical sciences by dramatically accelerating computational algorithms and amplifying insights available from computational chemistry methods. However, achieving this requires a confluence and coaction of expertise in computer science and physical sciences. This Review is written for new and experienced researchers working at the intersection of both fields. We first provide concise tutorials of computational chemistry and machine learning methods, showing how insights involving both can be achieved. We follow with a critical review of noteworthy applications that demonstrate how computational chemistry and machine learning can be used together to provide insightful (and useful) predictions in molecular and materials modeling, retrosyntheses, catalysis, and drug design.
Collapse
Affiliation(s)
- John A. Keith
- Department
of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Valentin Vassilev-Galindo
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Bingqing Cheng
- Accelerate
Programme for Scientific Discovery, Department
of Computer Science and Technology, 15 J. J. Thomson Avenue, Cambridge CB3 0FD, United Kingdom
| | - Stefan Chmiela
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Michael Gastegger
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587, Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea
- Max-Planck-Institut für Informatik, 66123 Saarbrücken, Germany
- Google Research, Brain Team, 10117 Berlin, Germany
| | - Alexandre Tkatchenko
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
198
|
Jankowska J, Sobolewski AL. Modern Theoretical Approaches to Modeling the Excited-State Intramolecular Proton Transfer: An Overview. Molecules 2021; 26:molecules26175140. [PMID: 34500574 PMCID: PMC8434569 DOI: 10.3390/molecules26175140] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Revised: 08/20/2021] [Accepted: 08/23/2021] [Indexed: 02/02/2023] Open
Abstract
The excited-state intramolecular proton transfer (ESIPT) phenomenon is nowadays widely acknowledged to play a crucial role in many photobiological and photochemical processes. It is an extremely fast transformation, often taking place at sub-100 fs timescales. While its experimental characterization can be highly challenging, a rich manifold of theoretical approaches at different levels is nowadays available to support and guide experimental investigations. In this perspective, we summarize the state-of-the-art quantum-chemical methods, as well as molecular- and quantum-dynamics tools successfully applied in ESIPT process studies, focusing on a critical comparison of their specific properties.
Collapse
Affiliation(s)
- Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
- Correspondence:
| | | |
Collapse
|
199
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
200
|
Hao D, Zhang Z, Ji Y. Responsive polymeric drug delivery systems for combination anticancer therapy: experimental design and computational insights. INT J POLYM MATER PO 2021. [DOI: 10.1080/00914037.2021.1960340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Dule Hao
- Jiangsu Province Hi-Tech Key Laboratory for Bio-medical Research, School of Chemistry and Chemical Engineering, Southeast University, Nanjing, China
| | - Zheng Zhang
- Jiangsu Province Hi-Tech Key Laboratory for Bio-medical Research, School of Chemistry and Chemical Engineering, Southeast University, Nanjing, China
| | - Yuanhui Ji
- Jiangsu Province Hi-Tech Key Laboratory for Bio-medical Research, School of Chemistry and Chemical Engineering, Southeast University, Nanjing, China
| |
Collapse
|