151
|
Heindl M, Hongyan J, Hua SA, Oelschlegel M, Meyer F, Schwarzer D, González L. Excited-State Dynamics of [Ru( S-Sbpy)(bpy) 2] 2+ to Form Long-Lived Localized Triplet States. Inorg Chem 2021; 60:1672-1682. [PMID: 33434007 PMCID: PMC7880568 DOI: 10.1021/acs.inorgchem.0c03163] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
![]()
The
novel photosensitizer [Ru(S–Sbpy)(bpy)2]2+ harbors two distinct sets of excited states
in the UV/Vis region of the absorption spectrum located on either
bpy or S–Sbpy ligands. Here, we address the question
of whether following excitation into these two types of states could
lead to the formation of different long-lived excited states from
where energy transfer to a reactive species could occur. Femtosecond
transient absorption spectroscopy identifies the formation of the
final state within 80 fs for both excitation wavelengths. The recorded
spectra hint at very similar dynamics following excitation toward
either the parent or sulfur-decorated bpy ligands, indicating ultrafast
interconversion into a unique excited-state species regardless of
the initial state. Non-adiabatic surface hopping dynamics simulations
show that ultrafast spin–orbit-mediated mixing of the states
within less than 50 fs strongly increases the localization of the
excited electron at the S–Sbpy ligand. Extensive
structural relaxation within this sulfurated ligand is possible, via
S–S bond cleavage that results in triplet state energies that
are lower than those in the analogue [Ru(bpy)3]2+. This structural relaxation upon localization of the charge on S–Sbpy is found to be the reason for the formation of
a single long-lived species independent of the excitation wavelength. We investigate selective excitation in
two types of ligands
present in [Ru(S−Sbpy)(bpy)2]2+. Femtosecond transient absorption spectroscopy identifies the formation
of the same final state within 80 fs for both excitation wavelengths.
Surface hopping simulations reveal ultrafast mixing of singlet and
triplet states within less than 50 fs. Energy lowering due to S−S
bond cleavage is identified as the driving factor for convergence
to a single final state.
Collapse
Affiliation(s)
- Moritz Heindl
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, A-1090 Vienna, Austria
| | - Jiang Hongyan
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, D-37077 Göttingen, Germany
| | - Shao-An Hua
- Institute of Inorganic Chemistry, University of Göttingen, D-37077 Göttingen, Germany
| | - Manuel Oelschlegel
- Institute of Inorganic Chemistry, University of Göttingen, D-37077 Göttingen, Germany
| | - Franc Meyer
- Institute of Inorganic Chemistry, University of Göttingen, D-37077 Göttingen, Germany
| | - Dirk Schwarzer
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, D-37077 Göttingen, Germany
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, A-1090 Vienna, Austria.,Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, A-1090 Vienna, Austria
| |
Collapse
|
152
|
Pieroni C, Marsili E, Lauvergnat D, Agostini F. Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies. J Chem Phys 2021; 154:034104. [PMID: 33499611 DOI: 10.1063/5.0036726] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.
Collapse
Affiliation(s)
- Carlotta Pieroni
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Emanuele Marsili
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
153
|
Menger MFS, Ehrmaier J, Faraji S. PySurf: A Framework for Database Accelerated Direct Dynamics. J Chem Theory Comput 2020; 16:7681-7689. [PMID: 33231447 PMCID: PMC7726901 DOI: 10.1021/acs.jctc.0c00825] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Indexed: 11/28/2022]
Abstract
The greatest restriction to the theoretical study of the dynamics of photoinduced processes is computationally expensive electronic structure calculations. Machine learning algorithms have the potential to reduce the number of these computations significantly. Here, PySurf is introduced as an innovative code framework, which is specifically designed for rapid prototyping and development tasks for data science applications in computational chemistry. It comes with powerful Plugin and Workflow engines, which allows intuitive customization for individual tasks. Data is automatically stored through the database framework, which enables additional interpolation of properties in previously evaluated regions of the conformational space. To illustrate the potential of the framework, a code for nonadiabatic surface hopping simulations based on the Landau-Zener algorithm is presented here. Deriving gradients from the interpolated potential energy surfaces allows for full-dimensional nonadiabatic surface hopping simulations using only adiabatic energies (energy only). Simulations of a pyrazine model and ab initio-based calculations of the SO2 molecule show that energy-only calculations with PySurf are able to correctly predict the nonadiabatic dynamics of these prototype systems. The results reveal the degree of sophistication, which can be achieved by the database accelerated energy-only surface hopping simulations being competitive to commonly used semiclassical approaches.
Collapse
Affiliation(s)
- Maximilian F. S.
J. Menger
- Zernike Institute
for Advanced
Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Johannes Ehrmaier
- Zernike Institute
for Advanced
Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Shirin Faraji
- Zernike Institute
for Advanced
Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| |
Collapse
|
154
|
Liu Y, Rozgonyi T, Marquetand P, Weinacht T. Excited-state dynamics of CH 2I 2 and CH 2IBr studied with UV-pump VUV-probe momentum-resolved photoion spectroscopy. J Chem Phys 2020; 153:184304. [PMID: 33187419 DOI: 10.1063/5.0026177] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We perform time-resolved ionization spectroscopy measurements of the excited state dynamics of CH2I2 and CH2IBr following photoexcitation in the deep UV. The fragment ions produced by ionization with a vacuum-ultraviolet probe pulse are measured with velocity map imaging, and the momentum resolved yields are compared with trajectory surface hopping calculations of the measurement observable. Together with recent time-resolved photoelectron spectroscopy measurements of the same dynamics, these results provide a detailed picture of the coupled electronic and nuclear dynamics involved. Our measurements highlight the non-adiabatic coupling between electronic states, which leads to notable differences in the dissociation dynamics for the two molecules.
Collapse
Affiliation(s)
- Yusong Liu
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, USA
| | - Tamás Rozgonyi
- Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| | - Philipp Marquetand
- University of Vienna, Faculty of Chemistry, Institute of Theoretical Chemistry, Währinger Straße 17, 1090 Wien, Austria
| | - Thomas Weinacht
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
155
|
Westermayr J, Marquetand P. Deep learning for UV absorption spectra with SchNarc: First steps toward transferability in chemical compound space. J Chem Phys 2020; 153:154112. [DOI: 10.1063/5.0021915] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- J. Westermayr
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - P. Marquetand
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Faculty of Chemistry, Data Science @ Uni Vienna, University of Vienna, Währinger Str. 29, 1090 Vienna, Austria
| |
Collapse
|
156
|
Zhao L, Zhou G, Jia B, Teng G, Zhan K, Zheng H, Luo J, Liu B. New insight into the ultrafast excited state deactivation mechanism of guanosine in the gas phase. J Photochem Photobiol A Chem 2020. [DOI: 10.1016/j.jphotochem.2020.112753] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
157
|
Li YZ, Yang JW, Makroni L, Wang WL, Liu FY. Photodynamics of methyl-vinyl Criegee intermediate: Different conical intersections govern the fates of syn/anti configurations. CHINESE J CHEM PHYS 2020. [DOI: 10.1063/1674-0068/cjcp2006088] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Ya-zhen Li
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an 710062, China
| | - Jia-wei Yang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an 710062, China
| | - Lily Makroni
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an 710062, China
| | - Wen-liang Wang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an 710062, China
| | - Feng-yi Liu
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an 710062, China
| |
Collapse
|
158
|
Mai S, González L. Molecular Photochemistry: Recent Developments in Theory. Angew Chem Int Ed Engl 2020; 59:16832-16846. [PMID: 32052547 PMCID: PMC7540682 DOI: 10.1002/anie.201916381] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 02/12/2020] [Indexed: 12/16/2022]
Abstract
Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.
Collapse
Affiliation(s)
- Sebastian Mai
- Photonics InstituteVienna University of TechnologyGusshausstrasse 27–291040ViennaAustria
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
159
|
Westermayr J, Marquetand P. Machine learning and excited-state molecular dynamics. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab9c3e] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
160
|
Ibele LM, Sánchez-Murcia PA, Mai S, Nogueira JJ, González L. Excimer Intermediates en Route to Long-Lived Charge-Transfer States in Single-Stranded Adenine DNA as Revealed by Nonadiabatic Dynamics. J Phys Chem Lett 2020; 11:7483-7488. [PMID: 32794719 PMCID: PMC7503858 DOI: 10.1021/acs.jpclett.0c02193] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 08/14/2020] [Indexed: 05/25/2023]
Abstract
The ultrafast time evolution of a single-stranded adenine DNA is studied using a hybrid multiscale quantum mechanics/molecular mechanics (QM/MM) scheme coupled to nonadiabatic surface hopping dynamics. As a model, we use (dA)20 where a stacked adenine tetramer is treated quantum chemically. The dynamical simulations combined with on-the-fly quantitative wave function analysis evidence the nature of the long-lived electronically excited states formed upon absorption of UV light. After a rapid decrease of the initially excited excitons, relaxation to monomer-like states and excimers occurs within 100 fs. The former monomeric states then relax into additional excimer states en route to forming stabilized charge-transfer states on a longer timescale of hundreds of femtoseconds. The different electronic-state characters is reflected on the spatial separation between the adenines: excimers and charge-transfer states show a much smaller spatial separation than the monomer-like states and the initially formed excitons.
Collapse
Affiliation(s)
- Lea M. Ibele
- Department
of Chemistry, Durham University, Durham DH1 3LE, U.K.
| | - Pedro A. Sánchez-Murcia
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Juan J. Nogueira
- Chemistry
Department, Universidad Autónoma
de Madridand IADCHEM, Institute for Advanced Research in Chemistry, Calle Francisco Tomás y Valiente
7, 28049 Madrid, Spain
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
161
|
Talotta F, Boggio‐Pasqua M, González L. Early Relaxation Dynamics in the Photoswitchable Complex trans-[RuCl(NO)(py) 4 ] 2. Chemistry 2020; 26:11522-11528. [PMID: 32281169 PMCID: PMC7539916 DOI: 10.1002/chem.202000507] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Revised: 03/16/2020] [Indexed: 11/22/2022]
Abstract
The design of photoswitchable transition metal complexes with tailored properties is one of the most important challenges in chemistry. Studies explaining the underlying mechanisms are, however, scarce. Herein, the early relaxation dynamics towards NO photoisomerization in trans-[RuCl(NO)(py)4 ]2+ is elucidated by means of non-adiabatic dynamics, which provided time-resolved information and branching ratios. Three deactivation mechanisms (I, II, III) in the ratio 3:2:4 were identified. Pathways I and III involve ultrafast intersystem crossing and internal conversion, whereas pathway II involves only internal conversion.
Collapse
Affiliation(s)
- Francesco Talotta
- lnstitute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
- Laboratoire de Chimie et Physique Quantiques, UMR 5626, lRSAMCCNRS et Université Toulouse 3118 route de Narbonne31062ToulouseFrance
- Present Address: Laboratoire de Chimie PhysiqueUMR 8000 CNRS/University Paris-SudUniversity Paris-Saclay, 91403 Orsay, andlnstitut de Sciences Moleculaires d'Orsay, UMR8214 CNRS/University Paris-Sud, University Paris-Saclay91403OrsayFrance
| | - Martial Boggio‐Pasqua
- Laboratoire de Chimie et Physique Quantiques, UMR 5626, lRSAMCCNRS et Université Toulouse 3118 route de Narbonne31062ToulouseFrance
| | - Leticia González
- lnstitute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
- Vienna Research Platform on Accelerating Photoreaction DiscoveryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
162
|
Chen HT, Zhou Z, Subotnik JE. On the proper derivation of the Floquet-based quantum classical Liouville equation and surface hopping describing a molecule or material subject to an external field. J Chem Phys 2020; 153:044116. [PMID: 32752688 DOI: 10.1063/5.0013873] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate different approaches to derive the proper Floquet-based quantum-classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis and then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form that is similar to the usual QCLE with two modifications: (1) The electronic degrees of freedom are expanded to infinite dimension and (2) the nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a fictitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., [J. Chem. Phys. 140, 084104 (2014)] who found that transforming to the adiabatic representation must always be the last operation applied, although now we have extended this result to a time-dependent Hamiltonian. This paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.
Collapse
Affiliation(s)
- Hsing-Ta Chen
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Zeyu Zhou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
163
|
Wang G, Li Y, Song T, Shang C, Yang J, Lily M, Fang Y, Liu F. Fluorescence Toggling Mechanism of Photochromic Phenylhydrazones: N–N Single Bond Rotation-Assisting E/Z Photoisomerization Differs from Imine. J Phys Chem A 2020; 124:6411-6419. [DOI: 10.1021/acs.jpca.0c03525] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gang Wang
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| | - Yazhen Li
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| | - Tingting Song
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| | - Congdi Shang
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| | - Jiawei Yang
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| | - Makroni Lily
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| | - Yu Fang
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| | - Fengyi Liu
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, Shaanxi 710119, People’s Republic of China
| |
Collapse
|
164
|
Prlj A, Ibele LM, Marsili E, Curchod BFE. On the Theoretical Determination of Photolysis Properties for Atmospheric Volatile Organic Compounds. J Phys Chem Lett 2020; 11:5418-5425. [PMID: 32543205 PMCID: PMC7372557 DOI: 10.1021/acs.jpclett.0c01439] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Volatile organic compounds (VOCs) are ubiquitous atmospheric molecules that generate a complex network of chemical reactions in the troposphere, often triggered by the absorption of sunlight. Understanding the VOC composition of the atmosphere relies on our ability to characterize all of their possible reaction pathways. When considering reactions of (transient) VOCs with sunlight, the availability of photolysis rate constants, utilized in general atmospheric models, is often out of experimental reach due to the unstable nature of these molecules. Here, we show how recent advances in computational photochemistry allow us to calculate in silico the different ingredients of a photolysis rate constant, namely, the photoabsorption cross-section and wavelength-dependent quantum yields. The rich photochemistry of tert-butyl hydroperoxide, for which experimental data are available, is employed to test our protocol and highlight the strengths and weaknesses of different levels of electronic structure and nonadiabatic molecular dynamics to study the photochemistry of (transient) VOCs.
Collapse
|
165
|
Shu Y, Zhang L, Sun S, Truhlar DG. Time-Derivative Couplings for Self-Consistent Electronically Nonadiabatic Dynamics. J Chem Theory Comput 2020; 16:4098-4106. [PMID: 32456433 DOI: 10.1021/acs.jctc.0c00409] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Electronically nonadiabatic dynamics methods based on a self-consistent potential, such as semiclassical Ehrenfest and coherent switching with decay of mixing, have a number of advantages but are computationally slower than approximations based on an unaveraged potential because they require evaluation of all components of the nonadiabatic coupling vector. Here we introduce a new approximation to the self-consistent potential that does not have this computational drawback. The new approximation uses time-derivative couplings evaluated by overlap integrals of electronic wave functions to approximate the nonadiabatic coupling terms in the equations of motion. We present a numerical test of the method for ethylene that shows there is little loss of accuracy in the ensemble-averaged results. This new approximation to the self-consistent potential makes direct dynamics calculations with self-consistent potentials more efficient for complex systems and makes them practically affordable for some cases where the cost was previously too high.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
166
|
Malone W, Nebgen B, White A, Zhang Y, Song H, Bjorgaard JA, Sifain AE, Rodriguez-Hernandez B, Freixas VM, Fernandez-Alberti S, Roitberg AE, Nelson TR, Tretiak S. NEXMD Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:5771-5783. [DOI: 10.1021/acs.jctc.0c00248] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Walter Malone
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Benjamin Nebgen
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Alexander White
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Josiah A. Bjorgaard
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Andrew E. Sifain
- U.S. Army Research Laboratory, Aberdeen Proving Ground, Aberdeen, Maryland 21005, United States
| | | | - Victor M. Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | | | - Adrian E. Roitberg
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Tammie R. Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
167
|
Gong Q, Li Y, Liu F. Trajectory surface-hopping dynamics of light-driven imine switches. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
168
|
Mai S, González L. Molekulare Photochemie: Moderne Entwicklungen in der theoretischen Chemie. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201916381] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Sebastian Mai
- Institut für Photonik Technische Universität Wien Gußhausstraße 27–29 1040 Wien Österreich
| | - Leticia González
- Institut für theoretische Chemie Fakultät für Chemie Universität Wien Währinger Straße 17 1090 Wien Österreich
| |
Collapse
|
169
|
Cao J, Chen DC. The excited-state relaxation mechanism of potential UVA-activated phototherapeutic molecules: trajectory surface hopping simulations of both 4-thiothymine and 2,4-dithiothymine. Phys Chem Chem Phys 2020; 22:10924-10933. [PMID: 32373808 DOI: 10.1039/d0cp01450h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Recent experimental investigations of the photochemical properties of a series of sulfur-substituted pyrimidine derivatives provide insights into the phototherapeutic potential of these nucleobase variants. Herein we elucidate the triplet formation mechanism of two prospective UVA-activated phototherapeutic molecules, 4-thiothymine and 2,4-dithiothymine, upon photo-excitation by applying the trajectory surface hopping dynamics at the LR-TDDFT level. Our simulations reasonably reproduce the experimental time constants and demonstrate the preferred triplet formation pathway which starts from the S1(nSπ*) state for both molecules. It is found that deactivation of the first bright state to the S1(nSπ*) state proceeds through a mechanism involving elongation of the C5-C6 and C4-S8 bond-lengths and C2-pyramidalization in 4-thiothymine and involving elongation of the C5-C6 and C2-S7 bond-lengths in 2,4-dithiothymine. The intersystem crossing of 2,4-dithiothymine occurs either at geometries characterized by elongated C5-C6 and C2-S7 bond-lengths or at geometries showing elongated C5-C6 and C4-S8 bond-lengths as seen in 4-thiothymine. The solvents are found to affect the S2 state decay of 4-thiothymine, leading to a competing pathway between S2→ S1 and S2→ T3. This study provides a molecular-level understanding of the underlying excited-state relaxation of the two UVA-activated thiopyrimidines, which may be linked to their potential applications in pharmacological science and also prove helpful for designing more effective phototherapeutic agents.
Collapse
Affiliation(s)
- Jun Cao
- School of Materials Science and Energy Engineering, Foshan University, Foshan, Guangdong 528000, P. R. China.
| | | |
Collapse
|
170
|
Westermayr J, Gastegger M, Marquetand P. Combining SchNet and SHARC: The SchNarc Machine Learning Approach for Excited-State Dynamics. J Phys Chem Lett 2020; 11:3828-3834. [PMID: 32311258 PMCID: PMC7246974 DOI: 10.1021/acs.jpclett.0c00527] [Citation(s) in RCA: 86] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 04/20/2020] [Indexed: 05/26/2023]
Abstract
In recent years, deep learning has become a part of our everyday life and is revolutionizing quantum chemistry as well. In this work, we show how deep learning can be used to advance the research field of photochemistry by learning all important properties-multiple energies, forces, and different couplings-for photodynamics simulations. We simplify such simulations substantially by (i) a phase-free training skipping costly preprocessing of raw quantum chemistry data; (ii) rotationally covariant nonadiabatic couplings, which can either be trained or (iii) alternatively be approximated from only ML potentials, their gradients, and Hessians; and (iv) incorporating spin-orbit couplings. As the deep-learning method, we employ SchNet with its automatically determined representation of molecular structures and extend it for multiple electronic states. In combination with the molecular dynamics program SHARC, our approach termed SchNarc is tested on two polyatomic molecules and paves the way toward efficient photodynamics simulations of complex systems.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Michael Gastegger
- Machine
Learning Group, Technical University of
Berlin, 10587 Berlin, Germany
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Str. 29, 1090 Vienna, Austria
| |
Collapse
|
171
|
Peters LM, Kussmann J, Ochsenfeld C. Combining Graphics Processing Units, Simplified Time-Dependent Density Functional Theory, and Finite-Difference Couplings to Accelerate Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2020; 11:3955-3961. [PMID: 32374606 PMCID: PMC7304892 DOI: 10.1021/acs.jpclett.0c00320] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Accepted: 04/23/2020] [Indexed: 06/11/2023]
Abstract
Starting from our recently published implementation of nonadiabatic molecular dynamics (NAMD) on graphics processing units (GPUs), we explore further approaches to accelerate ab initio NAMD calculations at the time-dependent density functional theory (TDDFT) level of theory. We employ (1) the simplified TDDFT schemes of Grimme et al. and (2) the Hammes-Schiffer-Tully approach to obtain nonadiabatic couplings from finite-difference calculations. The resulting scheme delivers an accurate physical picture while virtually eliminating the two computationally most demanding steps of the algorithm. Combined with our GPU-based integral routines for SCF, TDDFT, and TDDFT derivative calculations, NAMD simulations of systems of a few hundreds of atoms at a reasonable time scale become accessible on a single compute node. To demonstrate this and to present a first, illustrative example, we perform TDDFT/MM-NAMD simulations of the rhodopsin protein.
Collapse
Affiliation(s)
- Laurens
D. M. Peters
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Jörg Kussmann
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
172
|
Westermayr J, Faber FA, Christensen AS, von Lilienfeld OA, Marquetand P. Neural networks and kernel ridge regression for excited states dynamics of CH2NH$_2^+$: From single-state to multi-state representations and multi-property machine learning models. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab88d0] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
173
|
Balasubramani SG, Chen GP, Coriani S, Diedenhofen M, Frank MS, Franzke YJ, Furche F, Grotjahn R, Harding ME, Hättig C, Hellweg A, Helmich-Paris B, Holzer C, Huniar U, Kaupp M, Marefat Khah A, Karbalaei Khani S, Müller T, Mack F, Nguyen BD, Parker SM, Perlt E, Rappoport D, Reiter K, Roy S, Rückert M, Schmitz G, Sierka M, Tapavicza E, Tew DP, van Wüllen C, Voora VK, Weigend F, Wodyński A, Yu JM. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J Chem Phys 2020; 152:184107. [PMID: 32414256 PMCID: PMC7228783 DOI: 10.1063/5.0004635] [Citation(s) in RCA: 523] [Impact Index Per Article: 130.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Accepted: 04/07/2020] [Indexed: 01/30/2023] Open
Abstract
TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.
Collapse
Affiliation(s)
- Sree Ganesh Balasubramani
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Guo P Chen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Build. 207, DK-2800 Kongens Lyngby, Denmark
| | - Michael Diedenhofen
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Yannick J Franzke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | | | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Arnim Hellweg
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Uwe Huniar
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Alireza Marefat Khah
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | | | - Thomas Müller
- Forschungszentrum Jülich, Jülich Supercomputer Centre, Wilhelm-Jonen Straße, 52425 Jülich, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Brian D Nguyen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, USA
| | - Eva Perlt
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Kevin Reiter
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Saswata Roy
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Matthias Rückert
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, Langelandsgade 140, DK-8000 Aarhus, Denmark
| | - Marek Sierka
- TURBOMOLE GmbH, Litzenhardtstraße 19, 76135 Karlsruhe, Germany
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, USA
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstaße 1, 70569 Stuttgart, Germany
| | - Christoph van Wüllen
- Fachbereich Chemie and Forschungszentrum OPTIMAS, Technische Universität Kaiserslautern, Erwin-Schrödinger-Staße 52, 67663 Kaiserslautern, Germany
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Florian Weigend
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Artur Wodyński
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
174
|
Straub S, Stubbe J, Lindner J, Sarkar B, Vöhringer P. Vibrational Relaxation Dynamics of an Azido–Cobalt(II) Complex from Femtosecond UV-Pump/MIR-Probe Spectroscopy and Model Simulations with Ab Initio Anharmonic Couplings. Inorg Chem 2020; 59:14629-14642. [DOI: 10.1021/acs.inorgchem.0c00553] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Steffen Straub
- Lehrstuhl für Molekulare Physikalische Chemie, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115 Bonn, Germany
| | - Jessica Stubbe
- Institut für Chemie und Biochemie, Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/34, 14195 Berlin, Germany
| | - Jörg Lindner
- Lehrstuhl für Molekulare Physikalische Chemie, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115 Bonn, Germany
| | - Biprajit Sarkar
- Institut für Chemie und Biochemie, Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/34, 14195 Berlin, Germany
- Lehrstuhl für Anorganische Koordinationschemie, Institut für Anorganische Chemie Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Peter Vöhringer
- Lehrstuhl für Molekulare Physikalische Chemie, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115 Bonn, Germany
| |
Collapse
|
175
|
Lischka H, Shepard R, Müller T, Szalay PG, Pitzer RM, Aquino AJA, Araújo do Nascimento MM, Barbatti M, Belcher LT, Blaudeau JP, Borges I, Brozell SR, Carter EA, Das A, Gidofalvi G, González L, Hase WL, Kedziora G, Kertesz M, Kossoski F, Machado FBC, Matsika S, do Monte SA, Nachtigallová D, Nieman R, Oppel M, Parish CA, Plasser F, Spada RFK, Stahlberg EA, Ventura E, Yarkony DR, Zhang Z. The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry. J Chem Phys 2020; 152:134110. [PMID: 32268762 DOI: 10.1063/1.5144267] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.
Collapse
Affiliation(s)
- Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Ron Shepard
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Thomas Müller
- Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, Jülich 52428, Germany
| | - Péter G Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Budapest, Hungary
| | - Russell M Pitzer
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adelia J A Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | | | | | - Lachlan T Belcher
- Laser and Optics Research Center, Department of Physics, US Air Force Academy, Colorado 80840, USA
| | | | - Itamar Borges
- Departamento de Química, Instituto Militar de Engenharia, Rio de Janeiro, RJ 22290-270, Brazil
| | - Scott R Brozell
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Emily A Carter
- Office of the Chancellor and Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Box 951405, Los Angeles, California 90095-1405, USA
| | - Anita Das
- Indian Institute of Engineering Science and Technology, Shibpur, Howrah, India
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - William L Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Gary Kedziora
- Air Force Research Laboratory, Wright-Patterson Air Force Base, Ohio 45433, USA
| | - Miklos Kertesz
- Department of Chemistry, Georgetown University, 37th and O Streets, NW, Washington, DC 20057-1227, USA
| | | | - Francisco B C Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, 1901 N. 13th St., Philadelphia, Pennsylvania 19122, USA
| | | | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 160610 Prague 6, Czech Republic
| | - Reed Nieman
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Markus Oppel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Carol A Parish
- Department of Chemistry, Gottwald Center for the Sciences, University of Richmond, Richmond, Virginia 23173, USA
| | - Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Rene F K Spada
- Departamento de Física, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | - Eric A Stahlberg
- Biomedical Informatics and Data Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, USA
| | - Elizete Ventura
- Universidade Federal da Paraíba, 58059-900 João Pessoa, PB, Brazil
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, 3400 N. Charles Street, Baltimore, Maryland 21218, USA
| | - Zhiyong Zhang
- Stanford Research Computing Center, Stanford University, 255 Panama Street, Stanford, California 94305, USA
| |
Collapse
|
176
|
Shu Y, Zhang L, Mai S, Sun S, González L, Truhlar DG. Implementation of Coherent Switching with Decay of Mixing into the SHARC Program. J Chem Theory Comput 2020; 16:3464-3475. [DOI: 10.1021/acs.jctc.0c00112] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry 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, People’s Republic of China
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090, Vienna, Austria
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People’s Republic of China
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090, Vienna, Austria
| | - Donald G. Truhlar
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
177
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
178
|
Mai S, Menger MFSJ, Marazzi M, Stolba DL, Monari A, González L. Competing ultrafast photoinduced electron transfer and intersystem crossing of [Re(CO) 3 (Dmp)(His124)(Trp122)] + in Pseudomonas aeruginosa azurin: a nonadiabatic dynamics study. Theor Chem Acc 2020; 139:65. [PMID: 32214889 PMCID: PMC7078154 DOI: 10.1007/s00214-020-2555-6] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2019] [Accepted: 01/20/2020] [Indexed: 12/28/2022]
Abstract
We present a computational study of sub-picosecond nonadiabatic dynamics in a rhenium complex coupled electronically to a tryptophan (Trp) side chain of Pseudomonas aeruginosa azurin, a prototypical protein used in the study of electron transfer in proteins. To gain a comprehensive understanding of the photoinduced processes in this system, we have carried out vertical excitation calculations at the TDDFT level of theory as well as nonadiabatic dynamics simulations using the surface hopping including arbitrary couplings (SHARC) method coupled to potential energy surfaces represented with a linear vibronic coupling model. The results show that the initial photoexcitation populates both singlet metal-to-ligand charge transfer (MLCT) and singlet charge-separated (CS) states, where in the latter an electron was transferred from the Trp amino acid to the complex. Subsequently, a complex mechanism of simultaneous intersystem crossing and electron transfer leads to the sub-picosecond population of triplet MLCT and triplet CS states. These results confirm the assignment of the sub-ps time constants of previous experimental studies and constitute the first computational evidence for the ultrafast formation of the charge-separated states in Re-sensitized azurin.
Collapse
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Present Address: Photonics Institute, Vienna University of Technology, Gußhausstr. 27–29, 1040 Vienna, Austria
| | - Maximilian F. S. J. Menger
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Present Address: Zernike Institute for Advanced Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Marco Marazzi
- Department of Analytical Chemistry, Physical Chemistry and Chemical Engineering, Universidad de Alcalá, Ctra. Madrid-Barcelona Km. 33,600, 28871 Alcalá de Henares, Madrid Spain
- Chemical Research Institute “Andrés M. del Río” (IQAR), Universidad de Alcalá, 28871 Alcalá de Henares, Madrid Spain
| | - Dario L. Stolba
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Antonio Monari
- Université de Lorraine and CNRS, LPTC UMR, 7019 Nancy, France
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
179
|
Guan Y, Yarkony DR. Accurate Neural Network Representation of the Ab Initio Determined Spin-Orbit Interaction in the Diabatic Representation Including the Effects of Conical Intersections. J Phys Chem Lett 2020; 11:1848-1858. [PMID: 32062966 DOI: 10.1021/acs.jpclett.0c00074] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
A method for fitting ab initio determined spin-orbit coupling interactions, in the Breit-Pauli approximation, based on quasidiabatic representations using neural network fits is reported. The algorithm generalizes our recently reported neural network approach for representing the dipole interaction. The S0, S1, and T1 states of formaldehyde are used as an example. First, the two singlet states S0 and S1 are diabatized with a modified Boys Localization diabatization method. Second, the spin-orbit coupling between singlet and triplet states is transformed to the diabatic representation. This removes the discontinuities in the adiabatic representation. The diabatized spin-orbit couplings are then fit with smooth neural network functions. The analytic representation of spin-orbit coupling interactions in a diabatic basis by neural networks will make accurate full-dimensional quantum dynamical treatment of both internal conversion and intersystem crossing possible, which will help us to gain better understanding of both processes.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
180
|
Norell J, Odelius M, Vacher M. Ultrafast dynamics of photo-excited 2-thiopyridone: Theoretical insights into triplet state population and proton transfer pathways. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2020; 7:024101. [PMID: 32206689 PMCID: PMC7078009 DOI: 10.1063/1.5143228] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Accepted: 02/29/2020] [Indexed: 06/10/2023]
Abstract
Ultrafast non-adiabatic dynamics of the small heteroaromatic compound 2-thiopyridone has been studied with surface hopping simulations based on multi-configurational quantum chemistry. Initial excitation of the bright S 2 ( π , π * ) state is found to promptly relax to S 1 ( n , π * ) through in-plane motion. The subsequent dynamics are oppositely driven by out-of-plane motion, which results in both complex population transfers among all the available states and intersystem crossing predominantly through the "El-Sayed forbidden" S 1 ( n , π * ) to T 2 ( n , π * ) channel, through significant mixing of electronic excitation characters. Despite this complexity, the femto- to picosecond triplet population, expected from several spectroscopic measurements, is well described as a simple exponential decay of the singlet state manifold. No proton transfer is found in the reported trajectories, but two mechanisms for its possible mediation in previously reported experiments are proposed based on the observed structural dynamics: (i) ultrafast intra-molecular transfer driven by the initially coherent in-plane motion and (ii) inter-molecular solvent-mediated transfer driven by the out-of-plane modes that dominate the later motion.
Collapse
Affiliation(s)
- Jesper Norell
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Michael Odelius
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | | |
Collapse
|
181
|
Heim P, Mai S, Thaler B, Cesnik S, Avagliano D, Bella-Velidou D, Ernst WE, González L, Koch M. Revealing Ultrafast Population Transfer between Nearly Degenerate Electronic States. J Phys Chem Lett 2020; 11:1443-1449. [PMID: 31918552 PMCID: PMC7052817 DOI: 10.1021/acs.jpclett.9b03462] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2019] [Accepted: 01/10/2020] [Indexed: 06/10/2023]
Abstract
The response of a molecule to photoexcitation is governed by the coupling of its electronic states. However, if the energetic spacing between the electronically excited states at the Franck-Condon window becomes sufficiently small, it is infeasible to selectively excite and monitor individual states with conventional time-resolved spectroscopy, preventing insight into the energy transfer and relaxation dynamics of the molecule. Here, we demonstrate how the combination of time-resolved spectroscopy and extensive surface hopping dynamics simulations with a global fit approach on individually excited ensembles overcomes this limitation and resolves the dynamics in the n3p Rydberg states in acetone. Photoelectron transients of the three closely spaced states n3px, n3py, and n3pz are used to validate the theoretical results, which in turn allow retrieving a comprehensive kinetic model describing the mutual interactions of these states for the first time.
Collapse
Affiliation(s)
- Pascal Heim
- Institute
of Experimental Physics, Graz University
of Technology, Petersgasse 16, A-8010 Graz, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, A-1090 Vienna, Austria
| | - Bernhard Thaler
- Institute
of Experimental Physics, Graz University
of Technology, Petersgasse 16, A-8010 Graz, Austria
| | - Stefan Cesnik
- Institute
of Experimental Physics, Graz University
of Technology, Petersgasse 16, A-8010 Graz, Austria
| | - Davide Avagliano
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, A-1090 Vienna, Austria
| | - Dimitra Bella-Velidou
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, A-1090 Vienna, Austria
| | - Wolfgang E. Ernst
- Institute
of Experimental Physics, Graz University
of Technology, Petersgasse 16, A-8010 Graz, Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, A-1090 Vienna, Austria
| | - Markus Koch
- Institute
of Experimental Physics, Graz University
of Technology, Petersgasse 16, A-8010 Graz, Austria
| |
Collapse
|
182
|
Smith B, Akimov AV. Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:073001. [PMID: 31661681 DOI: 10.1088/1361-648x/ab5246] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States of America
| | | |
Collapse
|
183
|
Zhou Z, Chen HT, Nitzan A, Subotnik JE. Nonadiabatic Dynamics in a Laser Field: Using Floquet Fewest Switches Surface Hopping To Calculate Electronic Populations for Slow Nuclear Velocities. J Chem Theory Comput 2020; 16:821-834. [PMID: 31951404 DOI: 10.1021/acs.jctc.9b00950] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We investigate two well-known approaches for extending the fewest switches surface hopping (FSSH) algorithm to periodic time-dependent couplings. The first formalism acts as if the instantaneous adiabatic electronic states were standard adiabatic states, which just happen to evolve in time. The second formalism replaces the role of the usual adiabatic states by the time-independent adiabatic Floquet states. For a set of modified Tully model problems, the Floquet FSSH (F-FSSH) formalism gives a better estimate for both transmission and reflection probabilities than the instantaneous adiabatic FSSH (IA-FSSH) formalism, especially for slow nuclear velocities. More importantly, only F-FSSH predicts the correct final scattering momentum. Finally, in order to use Floquet theory accurately, we find that it is crucial to account for the interference between wavepackets on different Floquet states. Our results should be of interest to all those interested in laser-induced molecular dynamics.
Collapse
Affiliation(s)
- Zeyu Zhou
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| | - Hsing-Ta Chen
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| | - Abraham Nitzan
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| | - Joseph Eli Subotnik
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| |
Collapse
|
184
|
Nelson TR, White AJ, Bjorgaard JA, Sifain AE, Zhang Y, Nebgen B, Fernandez-Alberti S, Mozyrsky D, Roitberg AE, Tretiak S. Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials. Chem Rev 2020; 120:2215-2287. [PMID: 32040312 DOI: 10.1021/acs.chemrev.9b00447] [Citation(s) in RCA: 198] [Impact Index Per Article: 49.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Josiah A Bjorgaard
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Andrew E Sifain
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States.,U.S. Army Research Laboratory , Aberdeen Proving Ground , Maryland 21005 , United States
| | - Yu Zhang
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Benjamin Nebgen
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Adrian E Roitberg
- Department of Chemistry , University of Florida , Gainesville , Florida 32611 , United States
| | - Sergei Tretiak
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|
185
|
Bellonzi N, Alguire E, Fatehi S, Shao Y, Subotnik JE. TD-DFT spin-adiabats with analytic nonadiabatic derivative couplings. J Chem Phys 2020; 152:044112. [PMID: 32007078 PMCID: PMC7043850 DOI: 10.1063/1.5126440] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Accepted: 12/22/2019] [Indexed: 11/14/2022] Open
Abstract
We present an algorithm for efficient calculation of analytic nonadiabatic derivative couplings between spin-adiabatic, time-dependent density functional theory states within the Tamm-Dancoff approximation. Our derivation is based on the direct differentiation of the Kohn-Sham pseudowavefunction using the framework of Ou et al. Our implementation is limited to the case of a system with an even number of electrons in a closed shell ground state, and we validate our algorithm against finite difference at an S1/T2 crossing of benzaldehyde. Through the introduction of a magnetic field spin-coupling operator, we break time-reversal symmetry to generate complex valued nonadiabatic derivative couplings. Although the nonadiabatic derivative couplings are complex valued, we find that a phase rotation can generate an almost entirely real-valued derivative coupling vector for the case of benzaldehyde.
Collapse
Affiliation(s)
- Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | | | - Shervin Fatehi
- Department of Chemistry, The University of Texas Rio Grande Valley, Edinburg, Texas 78539, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
186
|
Francés‐Monerris A, Carmona‐García J, Acuña AU, Dávalos JZ, Cuevas CA, Kinnison DE, Francisco JS, Saiz‐Lopez A, Roca‐Sanjuán D. Photodissociation Mechanisms of Major Mercury(II) Species in the Atmospheric Chemical Cycle of Mercury. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201915656] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Antonio Francés‐Monerris
- Université de LorraineCNRS, LPCT 54000 Nancy France
- Departamento de Química FísicaUniversitat de València 46100 Burjassot Spain
| | | | - A. Ulises Acuña
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | - Juan Z. Dávalos
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | - Carlos A. Cuevas
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | | | - Joseph S. Francisco
- Department of Earth and Environmental Sciences and Department of ChemistryUniversity of Pennsylvania Philadelphia PA 19104 USA
| | - Alfonso Saiz‐Lopez
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | | |
Collapse
|
187
|
Francés‐Monerris A, Carmona‐García J, Acuña AU, Dávalos JZ, Cuevas CA, Kinnison DE, Francisco JS, Saiz‐Lopez A, Roca‐Sanjuán D. Photodissociation Mechanisms of Major Mercury(II) Species in the Atmospheric Chemical Cycle of Mercury. Angew Chem Int Ed Engl 2020; 59:7605-7610. [DOI: 10.1002/anie.201915656] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2019] [Indexed: 12/23/2022]
Affiliation(s)
- Antonio Francés‐Monerris
- Université de LorraineCNRS, LPCT 54000 Nancy France
- Departamento de Química FísicaUniversitat de València 46100 Burjassot Spain
| | | | - A. Ulises Acuña
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | - Juan Z. Dávalos
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | - Carlos A. Cuevas
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | | | - Joseph S. Francisco
- Department of Earth and Environmental Sciences and Department of ChemistryUniversity of Pennsylvania Philadelphia PA 19104 USA
| | - Alfonso Saiz‐Lopez
- Department of Atmospheric Chemistry and ClimateInstitute of Physical Chemistry RocasolanoCSIC 28006 Madrid Spain
| | | |
Collapse
|
188
|
Heller ER, Richardson JO. Instanton formulation of Fermi’s golden rule in the Marcus inverted regime. J Chem Phys 2020; 152:034106. [DOI: 10.1063/1.5137823] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Eric R. Heller
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
189
|
Abstract
The conical intersection dynamics of thiophenol is studied theoretically using the stimulated X-ray Raman imaging (SXRI) technique. SXRI employs a hard X-ray narrowband/broadband hybrid probe field and provides a real-time and real-space image of the passage through conical intersections. The signal, calculated using the minimal-coupling radiation/matter Hamiltonian, carries the phase information, and the real-space image of the transition charge density can be reconstructed by its Fourier transform. The two conical intersections (S2/S1 (11ππ*/1πσ*) and S1/S0 (1πσ*/S0)) can be distinguished and identified by the diffraction patterns in the level crossing regimes.
Collapse
Affiliation(s)
- Daeheum Cho
- Department of Chemistry , University of California , Irvine , California 92697-2025 , United States
| | - Shaul Mukamel
- Department of Chemistry , University of California , Irvine , California 92697-2025 , United States
| |
Collapse
|
190
|
|
191
|
Ibele LM, Curchod BFE. A molecular perspective on Tully models for nonadiabatic dynamics. Phys Chem Chem Phys 2020; 22:15183-15196. [DOI: 10.1039/d0cp01353f] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
We present a series of standardized molecular tests for nonadiabatic dynamics, reminiscent of the one-dimensional Tully models proposed in 1990.
Collapse
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry
- Durham University
- Durham DH1 3LE
- UK
| | | |
Collapse
|
192
|
Mai S, González L. Identification of important normal modes in nonadiabatic dynamics simulations by coherence, correlation, and frequency analyses. J Chem Phys 2019; 151:244115. [DOI: 10.1063/1.5129335] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
193
|
Peters LM, Kussmann J, Ochsenfeld C. Nonadiabatic Molecular Dynamics on Graphics Processing Units: Performance and Application to Rotary Molecular Motors. J Chem Theory Comput 2019; 15:6647-6659. [PMID: 31763834 PMCID: PMC6909237 DOI: 10.1021/acs.jctc.9b00859] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2019] [Indexed: 11/29/2022]
Abstract
Nonadiabatic molecular dynamics (NAMD) simulations of molecular systems require the efficient evaluation of excited-state properties, such as energies, gradients, and nonadiabatic coupling vectors. Here, we investigate the use of graphics processing units (GPUs) in addition to central processing units (CPUs) to efficiently calculate these properties at the time-dependent density functional theory (TDDFT) level of theory. Our implementation in the FermiONs++ program package uses the J-engine and a preselective screening procedure for the calculation of Coulomb and exchange kernels, respectively. We observe good speed-ups for small and large molecular systems (comparable to those observed in ground-state calculations) and reduced (down to sublinear) scaling behavior with respect to the system size (depending on the spatial locality of the investigated excitation). As a first illustrative application, we present efficient NAMD simulations of a series of newly designed light-driven rotary molecular motors and compare their S1 lifetimes. Although all four rotors show different S1 excitation energies, their ability to rotate upon excitation is conserved, making the series an interesting starting point for rotary molecular motors with tunable excitation energies.
Collapse
Affiliation(s)
- Laurens
D. M. Peters
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Jörg Kussmann
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
194
|
Park JW. Analytical Gradient Theory for Quasidegenerate N-Electron Valence State Perturbation Theory (QD-NEVPT2). J Chem Theory Comput 2019; 16:326-339. [DOI: 10.1021/acs.jctc.9b00919] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
195
|
Mai S, González L. Unconventional two-step spin relaxation dynamics of [Re(CO) 3(im)(phen)] + in aqueous solution. Chem Sci 2019; 10:10405-10411. [PMID: 32110331 PMCID: PMC6988600 DOI: 10.1039/c9sc03671g] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Accepted: 09/19/2019] [Indexed: 12/12/2022] Open
Abstract
Changes of molecular spin are ubiquitous in chemistry and biology. Among spin flip processes, one of the fastest is intersystem crossing (ISC) in transition metal complexes. Here, we investigate the spin relaxation dynamics and emission spectrum of [Re(CO)3(im)(phen)]+ (im = imidazole, phen = phenanthroline) using extensive full-dimensional excited-state dynamics simulations in explicit aqueous solution. Contrary to what has been observed in other transition metal complexes, the transition from the singlet to triplet states occurs via a two-step process, with clearly separable electronic and nuclear-driven components with two different time scales. The initially excited electronic wave function is a "molecular spin-orbit wave packet" that evolves almost instantaneously, with an 8 fs time constant, into an approximate 25 : 75 singlet-to-triplet equilibrium. Surprisingly, this ISC process is an order of magnitude faster than it was previously documented for this and other rhenium(i) carbonyl diimine complexes from emission spectra. Simulations including explicit laser field interactions evidence that few-cycle UV laser pulses are required to follow the creation and evolution of such molecular spin-orbit wave packets. The analysis of the dynamics also reveals a retarded ISC component, with a time constant of 420 fs, which can be explained invoking intramolecular vibrational energy redistribution. The emission spectrum is shown to be characterized by ISC convoluted with internal conversion and vibrational relaxation. These results provide fundamental understanding of ultrafast intersystem crossing in transition metal complexes.
Collapse
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry , Faculty of Chemistry , University of Vienna , Währinger Straße 17 , 1090 Vienna , Austria .
| | - Leticia González
- Institute of Theoretical Chemistry , Faculty of Chemistry , University of Vienna , Währinger Straße 17 , 1090 Vienna , Austria .
| |
Collapse
|
196
|
Zhou W, Mandal A, Huo P. Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations. J Phys Chem Lett 2019; 10:7062-7070. [PMID: 31665889 DOI: 10.1021/acs.jpclett.9b02747] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.
Collapse
Affiliation(s)
- Wanghuai Zhou
- Advanced Functional Material and Photoelectric Technology Research Institution, School of Science , Hubei University of Automotive Technology , Shiyan , Hubei 442002 , People's Republic of China
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Arkajit Mandal
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Pengfei Huo
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| |
Collapse
|
197
|
Guan Y, Guo H, Yarkony DR. Extending the Representation of Multistate Coupled Potential Energy Surfaces To Include Properties Operators Using Neural Networks: Application to the 1,21A States of Ammonia. J Chem Theory Comput 2019; 16:302-313. [DOI: 10.1021/acs.jctc.9b00898] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
198
|
Schmidt B, Klein R, Cancissu Araujo L. WavePacket: A Matlab package for numerical quantum dynamics. III. Quantum-classical simulations and surface hopping trajectories. J Comput Chem 2019; 40:2677-2688. [PMID: 31411345 DOI: 10.1002/jcc.26045] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2019] [Revised: 07/24/2019] [Accepted: 07/25/2019] [Indexed: 11/06/2022]
Abstract
WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I (Schmidt and Lorenz, Comput. Phys. Commun. 2017, 213, 223] and Part II (Schmidt and Hartmann, Comput. Phys. Commun. 2018, 228, 229] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagation techniques to WavePacket. There classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces, trajectories may switch between them. To model these transitions, two classes of stochastic algorithms have been implemented: (1) Tully's fewest switches surface hopping and (2) Landau-Zener-based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring nonadiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.1.0, which is essentially an object-oriented rewrite of previous versions, allowing to perform fully classical, quantum-classical and quantum-mechanical simulations on an equal footing, that is, for the same physical system described by the same WavePacket input. The software package is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics are available. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Burkhard Schmidt
- Institut für Mathematik, Freie Universität Berlin, Arnimallee 6, D-14195, Berlin, Germany
| | - Rupert Klein
- Institut für Mathematik, Freie Universität Berlin, Arnimallee 6, D-14195, Berlin, Germany
| | | |
Collapse
|
199
|
Taffet EJ, Lee BG, Toa ZSD, Pace N, Rumbles G, Southall J, Cogdell RJ, Scholes GD. Carotenoid Nuclear Reorganization and Interplay of Bright and Dark Excited States. J Phys Chem B 2019; 123:8628-8643. [PMID: 31553605 DOI: 10.1021/acs.jpcb.9b04027] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report quantum chemical calculations using multireference perturbation theory (MRPT) with the density matrix renormalization group (DMRG) plus photothermal deflection spectroscopy measurements to investigate the manifold of carotenoid excited states and establish their energies relative to the bright state (S2) as a function of nuclear reorganization. We conclude that the primary photophysics and function of carotenoids are determined by interplay of only the bright (S2) and lowest-energy dark (S1) states. The lowest-lying dark state, far from being energetically distinguishable from the lowest-lying bright state along the entire excited-state nuclear reorganization pathway, is instead computed to be either the second or first excited state depending on what equilibrium geometry is considered. This result suggests that, rather than there being a dark intermediate excited state bridging a non-negligible energy gap from the lowest-lying dark state to the lowest-lying bright state, there is in fact no appreciable energy gap to bridge following photoexcitation. Instead, excited-state nuclear reorganization constitutes the bridge from S2 to S1, in the sense that these two states attain energetic degeneracy along this pathway.
Collapse
Affiliation(s)
- Elliot J Taffet
- Department of Chemistry , Princeton University , Princeton , New Jersey 08544 , United States
| | - Benjamin G Lee
- Chemical and Materials Science Center , National Renewable Energy Laboratory , Golden , Colorado 80401 , United States
| | - Zi S D Toa
- Department of Chemistry , Princeton University , Princeton , New Jersey 08544 , United States
| | - Natalie Pace
- Chemical and Materials Science Center , National Renewable Energy Laboratory , Golden , Colorado 80401 , United States
| | - Garry Rumbles
- Chemical and Materials Science Center , National Renewable Energy Laboratory , Golden , Colorado 80401 , United States
| | - June Southall
- Institute of Molecular, Cell and Systems Biology, College of Medical Veterinary and Life Sciences , University of Glasgow , University Avenue, Glasgow G12 8QQ , U.K
| | - Richard J Cogdell
- Institute of Molecular, Cell and Systems Biology, College of Medical Veterinary and Life Sciences , University of Glasgow , University Avenue, Glasgow G12 8QQ , U.K
| | - Gregory D Scholes
- Department of Chemistry , Princeton University , Princeton , New Jersey 08544 , United States
| |
Collapse
|
200
|
Wang H, Odelius M, Prendergast D. A combined multi-reference pump-probe simulation method with application to XUV signatures of ultrafast methyl iodide photodissociation. J Chem Phys 2019; 151:124106. [PMID: 31575206 DOI: 10.1063/1.5116816] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
UV pump-XUV/X-ray probe measurements have been successfully applied in the study of photo-induced chemical reactions. Although rich element-specific electronic structure information is accessible within XUV/X-ray (inner-shell) absorption spectra, it can be difficult to interpret the chemistry directly from the spectrum without supporting theoretical simulations. A multireference method to completely simulate UV pump-XUV/X-ray probe measurement has been developed and applied to study the methyl iodide photodissociation process. Multireference, fewest-switches surface hopping (FSSH) trajectories were used to explore the coupled electronic and ionic dynamics upon photoexcitation of methyl iodide. Interpretation of previous measurements is provided by associated multireference, restricted active space, inner-shell spectral simulations. This combination of multireference FSSH trajectories and XUV spectra provides an interpretation of transient features appearing in previous measurements within the first 100 fs after photoexcitation and validates the significant branching ratio in the final excited-state population. This methodology should prove useful for interpretation of the increasing number of inner-shell probe studies of molecular excited states or for directing new experiments toward interesting regions of the potential energy landscape.
Collapse
Affiliation(s)
- Han Wang
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Michael Odelius
- Department of Physics, Stockholm University, AlbaNova University Center, SE-106 91 Stockholm, Sweden
| | - David Prendergast
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|