151
|
Song C, Martínez TJ, Neaton JB. A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods. J Chem Phys 2021; 155:024108. [PMID: 34266268 DOI: 10.1063/5.0055914] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Physics, University of California Berkeley, Berkeley, California 94720, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Jeffrey B Neaton
- Department of Physics, University of California Berkeley, Berkeley, California 94720, USA
| |
Collapse
|
152
|
Chen ZH, Wang Y, Xu RX, Yan Y. Correlated vibration-solvent effects on the non-Condon exciton spectroscopy. J Chem Phys 2021; 154:244105. [PMID: 34241336 DOI: 10.1063/5.0053169] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Excitation energy transfer is crucially involved in a variety of systems. During the process, the non-Condon vibronic coupling and the surrounding solvent interaction may synergetically play important roles. In this work, we study the correlated vibration-solvent influences on the non-Condon exciton spectroscopy. Statistical analysis is elaborated for the overall vibration-plus-solvent environmental effects. Analytic solutions are derived for the linear absorption of monomer systems. General simulations are accurately carried out via the dissipaton-equation-of-motion approach. The resulted spectra in either the linear absorption or strong field regime clearly demonstrate the coherence enhancement due to the synergetic vibration-solvent correlation.
Collapse
Affiliation(s)
- Zi-Hao Chen
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yao Wang
- Hefei National Laboratory for Physical Sciences at the Microscale and iChEM and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Rui-Xue Xu
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
153
|
Song H, Nam Y, Keefer D, Garavelli M, Mukamel S, Tretiak S. Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene. J Phys Chem Lett 2021; 12:5716-5722. [PMID: 34128675 DOI: 10.1021/acs.jpclett.1c01397] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S2/S3(D2d) → S1(D8h) → S0(D4h) → S0(D2d) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. The obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.
Collapse
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
| | - Yeonsig Nam
- Department of Chemistry, University of California, Irvine 92697, California, United States
| | - Daniel Keefer
- Department of Chemistry, University of California, Irvine 92697, California, United States
| | - Marco Garavelli
- Department of Industrial Chemistry, "T. Montanari", Università degli Studi di Bologna, Viale del Risorgimento, 4, Bologna 40136, Italy
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine 92697, California, United States
- Department of Physics and Astronomy, University of California, Irvine 92697, California, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
| |
Collapse
|
154
|
Li W, She Y, Vasenko AS, Prezhdo OV. Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites. NANOSCALE 2021; 13:10239-10265. [PMID: 34031683 DOI: 10.1039/d1nr01990b] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.
Collapse
Affiliation(s)
- Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China.
| | | | | | | |
Collapse
|
155
|
Shao C, Xu J, Wang L. Branching and phase corrected surface hopping: A benchmark of nonadiabatic dynamics in multilevel systems. J Chem Phys 2021; 154:234109. [PMID: 34241240 DOI: 10.1063/5.0056224] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Since the seminal work of Tully [J. Chem. Phys. 93, 1061 (1990)], two-level scattering models have been extensively adopted as the standard benchmark systems to assess the performance of different trajectory surface hopping methods for nonadiabatic dynamics simulations. Here, we extend the branching and phase corrections to multilevel systems and combine them with both the traditional fewest switches surface hopping (FSSH) and its variant global flux surface hopping (GFSH) algorithms. To get a comprehensive evaluation of the proposed methods, we construct a series of more challenging and diverse three-level and four-level scattering models and use exact quantum solutions as references. Encouragingly, both FSSH and GFSH with the branching and phase corrections produce excellent and nearly identical results in all investigated systems, indicating that the new surface hopping methods are robust to describe multilevel problems and the reliability is insensitive to the definition of self-consistent hopping probabilities in the adiabatic representation. Furthermore, the branching correction is found to be especially important when dealing with strongly repulsive potential energy surfaces, which are common in realistic systems, thus promising for general applications.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
156
|
Krotz A, Provazza J, Tempelaar R. A reciprocal-space formulation of mixed quantum-classical dynamics. J Chem Phys 2021; 154:224101. [PMID: 34241207 DOI: 10.1063/5.0053177] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We derive a formulation of mixed quantum-classical dynamics for modeling electronic carriers interacting with phonons in reciprocal space. For dispersionless phonons, we start by expressing the real-space classical coordinates in terms of complex variables. Taking these variables as a Fourier series then yields the reciprocal-space coordinates. Evaluating the electron-phonon interaction term through Ehrenfest's theorem, we arrive at a reciprocal-space formalism that is equivalent to mean-field mixed quantum-classical dynamics in real space. This equivalence is numerically verified for the Holstein and Peierls models, for which we find the reciprocal-space Hellmann-Feynman forces to involve momentum-derivative contributions in addition to the position-derivative terms commonly seen in real space. To illustrate the advantage of the reciprocal-space formulation, we present a proof of concept for the inexpensive modeling of low-momentum carriers interacting with phonons using a truncated reciprocal-space basis, which is not possible within a real-space formulation.
Collapse
Affiliation(s)
- Alex Krotz
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| | - Justin Provazza
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| | - Roel Tempelaar
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| |
Collapse
|
157
|
Gonzalvez Perez I, Barford W. Ultrafast Fluorescence Depolarization in Conjugated Polymers. J Phys Chem Lett 2021; 12:5344-5348. [PMID: 34076446 DOI: 10.1021/acs.jpclett.1c01354] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We report on large-scale simulations of intrachain exciton dynamics in poly(para-phenylenevinylene). Our theoretical model describes Frenkel exciton coupling to both fast, quantized C-C bond vibrations and slow, classical torsional modes. We also incorporate system-bath interactions. The dynamics is simulated using the time evolution block decimation method, which avoids the failures of the Ehrenfest approximation to describe decoherence processes and nonadiabatic interstate conversion. System-bath interactions are modeled using quantum trajectories and Lindblad quantum jump operators. We find that following photoexcitation, the quantum mechanical entanglement of the exciton and C-C bond phonons causes exciton-site decoherence. Next, system-bath interactions cause the stochastic collapse of high-energy delocalized excitons into chromophores. Finally, torsional relaxation causes additional exciton-density localization. We relate these dynamical processes to the predicted fluorescence depolarization, extract the time scales corresponding to them, and thus interpret the observed sub-ps fluorescence depolarization.
Collapse
Affiliation(s)
- Isabel Gonzalvez Perez
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - William Barford
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
158
|
Alfonso-Hernandez L, Oldani N, Athanasopoulos S, Lupton JM, Tretiak S, Fernandez-Alberti S. Photoinduced Energy Transfer in Linear Guest-Host Chromophores: A Computational Study. J Phys Chem A 2021; 125:5303-5313. [PMID: 34106721 DOI: 10.1021/acs.jpca.1c02644] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Polymer-based guest-host systems represent a promising class of materials for efficient light-emitting diodes. The energy transfer from the polymer host to the guest is the key process in light generation. Therefore, microscopic descriptions of the different mechanisms involved in the energy transfer can contribute to enlighten the basis of the highly efficient light harvesting observed in this kind of materials. Herein, the nature of intramolecular energy transfer in a dye-end-capped conjugated polymer is explored by using atomistic nonadiabatic excited-state molecular dynamics. Linear perylene end-capped (PEC) polyindenofluorenes (PIF), consisting of n (n = 2, 4, and 6) repeat units, i.e., PEC-PIFn oligomers, are considered as model systems. After photoexcitation at the oligomer absorption maximum, an initial exciton becomes self-trapped on one of the monomer units (donors). Thereafter, an efficient ultrafast through-space energy transfer from this unit to the perylene acceptor takes place. We observe that this energy transfer occurs equally well from any monomer unit on the chain. Effective specific vibronic couplings between each monomer and the acceptor are identified. These oligomer → end-cap energy transfer steps do not match with the rates predicted by Förster-type energy transfer. The through-space and through-bond mechanisms are two distinct channels of energy transfer. The former dominates the overall process, whereas the through-bond energy transfer between indenofluorene monomer units along the oligomer backbone only makes a minor contribution.
Collapse
Affiliation(s)
- L Alfonso-Hernandez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - N Oldani
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - S Athanasopoulos
- Departamento de Física, Universidad Carlos III de Madrid, Avenida Universidad 30, 28911 Leganés, Madrid, Spain
| | - J M Lupton
- Institut für Angewandte und Experimentelle Physik, Universität Regensburg, Universitätsstrasse 31, 93053 Regensburg, Germany
| | - S Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - S Fernandez-Alberti
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| |
Collapse
|
159
|
Song H, Freixas VM, Fernandez-Alberti S, White AJ, Zhang Y, Mukamel S, Govind N, Tretiak S. An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem. J Chem Theory Comput 2021; 17:3629-3643. [PMID: 34014085 DOI: 10.1021/acs.jctc.1c00131] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.
Collapse
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD, Bernal, Argentina
| | | | - Alexander J 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
| | - Shaul Mukamel
- Departments of Chemistry, Physics, and Astronomy, University of California, Irvine, California 92697, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.,Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
160
|
Abstract
Theoretical simulations of electronic excitations and associated processes in molecules are indispensable for fundamental research and technological innovations. However, such simulations are notoriously challenging to perform with quantum mechanical methods. Advances in machine learning open many new avenues for assisting molecular excited-state simulations. In this Review, we track such progress, assess the current state of the art and highlight the critical issues to solve in the future. We overview a broad range of machine learning applications in excited-state research, which include the prediction of molecular properties, improvements of quantum mechanical methods for the calculations of excited-state properties and the search for new materials. Machine learning approaches can help us understand hidden factors that influence photo-processes, leading to a better control of such processes and new rules for the design of materials for optoelectronic applications.
Collapse
|
161
|
Nematiaram T, Padula D, Troisi A. Bright Frenkel Excitons in Molecular Crystals: A Survey. CHEMISTRY OF MATERIALS : A PUBLICATION OF THE AMERICAN CHEMICAL SOCIETY 2021; 33:3368-3378. [PMID: 34526736 PMCID: PMC8432684 DOI: 10.1021/acs.chemmater.1c00645] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/23/2021] [Revised: 04/14/2021] [Indexed: 05/12/2023]
Abstract
We computed the optical properties of a large set of molecular crystals (∼2200 structures) composed of molecules whose lowest excited states are strongly coupled and generate wide excitonic bands. Such bands are classified in terms of their dimensionality (1-, 2-, and 3-dimensional), the position of the optically allowed state in relation with the excitonic density of states, and the presence of Davydov splitting. The survey confirms that one-dimensional aggregates are rare in molecular crystals highlighting the need to go beyond the simple low-dimensional models. Furthermore, this large set of data is used to search for technologically interesting and less common properties. For instance, we considered the largest excitonic bandwidth that is achievable within known molecular crystals and identified materials with strong super-radiant states. Finally, we explored the possibility that strong excitonic coupling can be used to generate emissive states in the near-infrared region in materials formed by molecules with bright visible absorption and we could identify the maximum allowable red shift in this material class. These insights with the associated searchable database provide practical guidelines for designing materials with interesting optical properties.
Collapse
Affiliation(s)
- Tahereh Nematiaram
- Department
of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, U.K.
| | - Daniele Padula
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, via A. Moro 2, Siena 53100, Italy
| | - Alessandro Troisi
- Department
of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, U.K.
| |
Collapse
|
162
|
Barbatti M. Velocity Adjustment in Surface Hopping: Ethylene as a Case Study of the Maximum Error Caused by Direction Choice. J Chem Theory Comput 2021; 17:3010-3018. [PMID: 33844922 DOI: 10.1021/acs.jctc.1c00012] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The most common surface hopping dynamics algorithms require velocity adjustment after hopping to ensure total-energy conservation. Based on the semiclassical analysis, this adjustment must be made parallel to the nonadiabatic coupling vector's direction. Nevertheless, this direction is not always known, and the common practice has been to adjust the velocity in either the linear momentum or velocity directions. This paper benchmarks surface hopping dynamics of photoexcited ethylene with velocity adjustment in several directions, including those of the nonadiabatic coupling vector, the momentum, and the energy gradient difference. It is shown that differences in time constants and structural evolution fall within the statistical uncertainty of the method considering up to 500 trajectories in each dynamics set, rendering the three approaches statistically equivalent. For larger ensembles beyond 1000 trajectories, significant differences between the results arise, limiting the validity of adjustment in alternative directions. Other possible adjustment directions (velocity, single-state gradients, angular momentum) are evaluated as well. Given the small size of ethylene, the results reported in this paper should be considered an upper limit for the error caused by the choice of the velocity-adjustment direction on surface hopping dynamics.
Collapse
|
163
|
Cao X, Tian P. "Dividing and Conquering" and "Caching" in Molecular Modeling. Int J Mol Sci 2021; 22:5053. [PMID: 34068835 PMCID: PMC8126232 DOI: 10.3390/ijms22095053] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Revised: 04/26/2021] [Accepted: 04/27/2021] [Indexed: 11/17/2022] Open
Abstract
Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, most important methodological advancements in more than half century of molecular modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but results are not transferable. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science. This framework is based on a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.
Collapse
Affiliation(s)
- Xiaoyong Cao
- School of Life Sciences, Jilin University, Changchun 130012, China;
| | - Pu Tian
- School of Life Sciences, Jilin University, Changchun 130012, China;
- School of Artificial Intelligence, Jilin University, Changchun 130012, China
| |
Collapse
|
164
|
Madjet MEA, Ali E, Carignano M, Vendrell O, Chakraborty HS. Ultrafast Transfer and Transient Entrapment of Photoexcited Mg Electron in Mg@C_{60}. PHYSICAL REVIEW LETTERS 2021; 126:183002. [PMID: 34018762 DOI: 10.1103/physrevlett.126.183002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Revised: 02/12/2021] [Accepted: 03/26/2021] [Indexed: 06/12/2023]
Abstract
Electron relaxation is studied in endofullerene Mg@C_{60} after an initial localized photoexcitation in Mg by nonadiabatic molecular dynamics simulations. Two approaches to the electronic structure of the excited electronic states are used: (i) an independent particle approximation based on a density-functional theory description of molecular orbitals and (ii) a configuration-interaction description of the many-body effects. Both methods exhibit similar relaxation times, leading to an ultrafast decay and charge transfer from Mg to C_{60} within tens of femtoseconds. Method (i) further elicits a transient trap of the transferred electron that can delay the electron-hole recombination. Results shall motivate experiments to probe these ultrafast processes by two-photon transient absorption or photoelectron spectroscopy in gas phase, in solution, or as thin films.
Collapse
Affiliation(s)
- Mohamed El-Amine Madjet
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Straße 38, 01187 Dresden, Germany
- Department of Natural Sciences, D. L. Hubbard Center for Innovation, Northwest Missouri State University, Maryville, Missouri 64468, USA
| | - Esam Ali
- Department of Natural Sciences, D. L. Hubbard Center for Innovation, Northwest Missouri State University, Maryville, Missouri 64468, USA
| | - Marcelo Carignano
- Department of Biomedical Engineering, Northwestern University, Evanston, Illinois 60208, USA
| | - Oriol Vendrell
- Theoretical Chemistry, Institute of Physical Chemistry and Centre for Advanced Materials, Heidelberg University, Im Neuenheimer Feld 229 and 225, 69120 Heidelberg, Germany
| | - Himadri S Chakraborty
- Department of Natural Sciences, D. L. Hubbard Center for Innovation, Northwest Missouri State University, Maryville, Missouri 64468, USA
| |
Collapse
|
165
|
Fortino M, Cozza C, Bonomi M, Pietropaolo A. Multi-replica biased sampling for photoswitchable π-conjugated polymers. J Chem Phys 2021; 154:174108. [PMID: 34241080 DOI: 10.1063/5.0045944] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In recent years, π-conjugated polymers are attracting considerable interest in view of their light-dependent torsional reorganization around the π-conjugated backbone, which determines peculiar light-emitting properties. Motivated by the interest in designing conjugated polymers with tunable photoswitchable pathways, we devised a computational framework to enhance the sampling of the torsional conformational space and, at the same time, estimate ground- to excited-state free-energy differences. This scheme is based on a combination of Hamiltonian Replica Exchange Method (REM), parallel bias metadynamics, and free-energy perturbation theory. In our scheme, each REM samples an intermediate unphysical state between the ground and the first two excited states, which are characterized by time-dependent density functional theory simulations at the B3LYP/6-31G* level of theory. We applied the method to a 5-mer of 9,9-dioctylfluorene and found that upon irradiation, this system can undergo a dihedral inversion from -155° to 155°, crossing a barrier that decreases from 0.1 eV in the ground state (S0) to 0.05 eV and 0.04 eV in the first (S1) and second (S2) excited states. Furthermore, S1 and even more S2 were predicted to stabilize coplanar dihedrals, with a local free-energy minimum located at ±44°. The presence of a free-energy barrier of 0.08 eV for the S1 state and 0.12 eV for the S2 state can trap this conformation in a basin far from the global free-energy minimum located at 155°. The simulation results were compared with the experimental emission spectrum, showing a quantitative agreement with the predictions provided by our framework.
Collapse
Affiliation(s)
- Mariagrazia Fortino
- Dipartimento di Scienze Della Salute, Università di Catanzaro, Viale Europa, 88100 Catanzaro, Italy
| | - Concetta Cozza
- Dipartimento di Scienze Della Salute, Università di Catanzaro, Viale Europa, 88100 Catanzaro, Italy
| | - Massimiliano Bonomi
- Structural Bioinformatics Unit, Department of Structural Biology and Chemistry; CNRS UMR 3528; C3BI, CNRS USR 3756; Institut Pasteur, Paris, France
| | - Adriana Pietropaolo
- Dipartimento di Scienze Della Salute, Università di Catanzaro, Viale Europa, 88100 Catanzaro, Italy
| |
Collapse
|
166
|
Lykhin AO, Truhlar DG, Gagliardi L. Role of Triplet States in the Photodynamics of Aniline. J Am Chem Soc 2021; 143:5878-5889. [PMID: 33843225 DOI: 10.1021/jacs.1c00989] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The dynamics of excited heteroaromatic molecules is a key to understanding the photoprotective properties of many biologically relevant chromophores that dissipate their excitation energy nonreactively and thereby prevent the detrimental effects of ultraviolet radiation. Despite their structural variability, most substituted aromatic compounds share a common feature of a repulsive 1πσ* potential energy surface. This surface can lead to photoproducts, and it can also facilitate the population transfer back to the ground electronic state by means of a 1πσ*/S0 conical intersection. Here, we explore a hidden relaxation route involving the triplet electronic state of aniline, which has recently been discovered by means of time-selected photofragment translational spectroscopy [J. Chem. Phys. 2019, 151, 141101]. By using the recently available analytical gradients for multiconfiguration pair-density functional theory, it is now possible to locate the minimum-energy crossing points between states of different spin and therefore compute the intersystem crossing rates with a multireference method, rather than with the less reliable single-reference methods. Using such calculations, we demonstrate that the population loss of aniline in the T1(3ππ*) state is dominated by C6H5NH2 → C6H5NH· + H· dissociation, and we explain the long nonradiative lifetimes of the T1(3ππ*) state at the excitation wavelengths of 294-264 nm.
Collapse
Affiliation(s)
- Aleksandr O Lykhin
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
167
|
Nottoli M, Cupellini L, Lipparini F, Granucci G, Mennucci B. Multiscale Models for Light-Driven Processes. Annu Rev Phys Chem 2021; 72:489-513. [PMID: 33561359 DOI: 10.1146/annurev-physchem-090419-104031] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| |
Collapse
|
168
|
Di Maiolo F, Worth GA, Burghardt I. Multi-layer Gaussian-based multi-configuration time-dependent Hartree (ML-GMCTDH) simulations of ultrafast charge separation in a donor-acceptor complex. J Chem Phys 2021; 154:144106. [PMID: 33858146 DOI: 10.1063/5.0046933] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report on first applications of the Multi-Layer Gaussian-based Multi-Configuration Time-Dependent Hartree (ML-GMCTDH) method [Römer et al., J. Chem. Phys. 138, 064106 (2013)] beyond its basic two-layer variant. The ML-GMCTDH scheme provides an embedding of a variationally evolving Gaussian wavepacket basis into a hierarchical tensor representation of the wavefunction. A first-principles parameterized model Hamiltonian for ultrafast non-adiabatic dynamics in an oligothiophene-fullerene charge transfer complex is employed, relying on a two-state linear vibronic coupling model that combines a distribution of tuning type modes with an intermolecular coordinate that also modulates the electronic coupling. Efficient ML-GMCTDH simulations are carried out for up to 300 vibrational modes using an implementation within the QUANTICS program. Excellent agreement with reference ML-MCTDH calculations is obtained.
Collapse
Affiliation(s)
- Francesco Di Maiolo
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Graham A Worth
- Department of Chemistry, University College London, 20 Gordon St., London WC1H 0AJ, United Kingdom
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| |
Collapse
|
169
|
Santoro F, Green JA, Martinez-Fernandez L, Cerezo J, Improta R. Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives. Phys Chem Chem Phys 2021; 23:8181-8199. [PMID: 33875988 DOI: 10.1039/d0cp05907b] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.
Collapse
Affiliation(s)
- Fabrizio Santoro
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - James A Green
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy.
| | - Lara Martinez-Fernandez
- Departamento de Química, Facultad de Ciencias and Institute for Advanced Research in Chemistry (IADCHEM), Universidad Autónoma de Madrid, Campus de Excelencia UAM-CSIC, 28049 Madrid, Spain
| | - Javier Cerezo
- Departamento de Química, Facultad de Ciencias and Institute for Advanced Research in Chemistry (IADCHEM), Universidad Autónoma de Madrid, Campus de Excelencia UAM-CSIC, 28049 Madrid, Spain
| | - Roberto Improta
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy.
| |
Collapse
|
170
|
Niehaus TA. Ground-to-excited derivative couplings for the density functional-based tight-binding method: semi-local and long-range corrected formulations. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02735-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
171
|
Freixas VM, White AJ, Nelson T, Song H, Makhov DV, Shalashilin D, Fernandez-Alberti S, Tretiak S. Nonadiabatic Excited-State Molecular Dynamics Methodologies: Comparison and Convergence. J Phys Chem Lett 2021; 12:2970-2982. [PMID: 33730495 DOI: 10.1021/acs.jpclett.1c00266] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum-classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.
Collapse
Affiliation(s)
- Victor M Freixas
- Universidad Nacional de Quilmes, Roque Saénz Peña 352, B1876BXD Bernal, Argentina
| | - Alexander J White
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie Nelson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dmitry V Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, U.K
- School of Mathematics, University of Bristol, Bristol BS8 1TW, U.K
| | | | | | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
172
|
Gelin MF, Huang X, Xie W, Chen L, Došlić NA, Domcke W. Ab Initio Surface-Hopping Simulation of Femtosecond Transient-Absorption Pump-Probe Signals of Nonadiabatic Excited-State Dynamics Using the Doorway-Window Representation. J Chem Theory Comput 2021; 17:2394-2408. [PMID: 33755464 DOI: 10.1021/acs.jctc.1c00109] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Xiang Huang
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| | - Weiwei Xie
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems, D-01187 Dresden, Germany
| | - Nad A Došlić
- Department of Physical Chemistry, Ruder Boscovic Institute, HR-10000 Zagreb, Croatia
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| |
Collapse
|
173
|
Uratani H, Yoshikawa T, Nakai H. Trajectory Surface Hopping Approach to Condensed-Phase Nonradiative Relaxation Dynamics Using Divide-and-Conquer Spin-Flip Time-Dependent Density-Functional Tight Binding. J Chem Theory Comput 2021; 17:1290-1300. [PMID: 33577323 DOI: 10.1021/acs.jctc.0c01155] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nonradiative relaxation of excited molecules is central to many crucial issues in photochemistry. Condensed phases are typical contexts in which such problems are considered, and the nonradiative relaxation dynamics are expected to be significantly affected by interactions with the environment, for example, a solvent. We developed a nonadiabatic molecular dynamics simulation technique that can treat the nonradiative relaxation and explicitly include the environment in the calculations without a heavy computational burden. Specifically, we combined trajectory surface hopping with Tully's fewest-switches algorithm, a tight-binding approximated version of spin-flip time-dependent density-functional theory, and divide-and-conquer (DC) spatial fragmentation scheme. Numerical results showed that this method can treat systems with thousands of atoms within reasonable computational resources, and the error arising from DC fragmentation is negligibly small. Using this method, we obtained molecular insights into the solvent dependence of the photoexcited-state dynamics of trans-azobenzene, which demonstrate the importance of the environment for condensed-phase nonradiative relaxation.
Collapse
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
174
|
Fortino M, Collini E, Bloino J, Pedone A. Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations. J Chem Phys 2021; 154:094110. [PMID: 33685164 DOI: 10.1063/5.0039949] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The non-radiative relaxation process within the Q-bands of chlorophylls represents a crucial preliminary step during the photosynthetic mechanism. Despite several experimental and theoretical efforts performed in order to clarify the complex dynamics characterizing this stage, a complete understanding of this mechanism is still far to be reached. In this study, non-adiabatic excited-state molecular dynamic simulations have been performed to model the non-radiative process within the Q-bands for a model system of chlorophylls. This system has been considered in the gas phase and then, to have a more representative picture of the environment, with implicit and mixed implicit-explicit solvation models. In the first part of this analysis, absorption spectra have been simulated for each model in order to guide the setup for the non-adiabatic excited-state molecular dynamic simulations. Then, non-adiabatic excited-state molecular dynamic simulations have been performed on a large set of independent trajectories and the population of the Qx and Qy states has been computed as the average of all the trajectories, estimating the rate constant for the process. Finally, with the aim of investigating the possible role played by the solvent in the Qx-Qy crossing mechanism, an essential dynamic analysis has been performed on the generated data, allowing one to find the most important motions during the simulated dynamics.
Collapse
Affiliation(s)
| | | | | | - Alfonso Pedone
- Università di Modena e Reggio Emilia, Modena 45125, Italy
| |
Collapse
|
175
|
Popp W, Brey D, Binder R, Burghardt I. Quantum Dynamics of Exciton Transport and Dissociation in Multichromophoric Systems. Annu Rev Phys Chem 2021; 72:591-616. [PMID: 33636997 DOI: 10.1146/annurev-physchem-090419-040306] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Due to the subtle interplay of site-to-site electronic couplings, exciton delocalization, nonadiabatic effects, and vibronic couplings, quantum dynamical studies are needed to elucidate the details of ultrafast photoinduced energy and charge transfer events in organic multichromophoric systems. In this vein, we review an approach that combines first-principles parameterized lattice Hamiltonians with accurate quantum dynamical simulations using advanced multiconfigurational methods. Focusing on the elementary transfer steps in organic functional materials, we address coherent exciton migration and creation of charge transfer excitons in homopolymers, notably representative of the poly(3-hexylthiophene) material, as well as exciton dissociation at polymer:fullerene heterojunctions. We emphasize the role of coherent transfer, trapping effects due to high-frequency phonon modes, and thermal activation due to low-frequency soft modes that drive a diffusive dynamics.
Collapse
Affiliation(s)
- Wjatscheslaw Popp
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Dominik Brey
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| |
Collapse
|
176
|
Mališ M, Luber S. ΔSCF with Subsystem Density Embedding for Efficient Nonadiabatic Molecular Dynamics in Condensed-Phase Systems. J Chem Theory Comput 2021; 17:1653-1661. [DOI: 10.1021/acs.jctc.0c01200] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
177
|
Darghouth AAMHM, Casida ME, Zhu X, Natarajan B, Su H, Humeniuk A, Titov E, Miao X, Mitrić R. Effect of varying the TD-lc-DFTB range-separation parameter on charge and energy transfer in a model pentacene/buckminsterfullerene heterojunction. J Chem Phys 2021; 154:054102. [PMID: 33557554 DOI: 10.1063/5.0024559] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Atomistic modeling of energy and charge transfer at the heterojunction of organic solar cells is an active field with many remaining outstanding questions owing, in part, to the difficulties in performing reliable photodynamics calculations on very large systems. One approach to being able to overcome these difficulties is to design and apply an appropriate simplified method. Density-functional tight binding (DFTB) has become a popular form of approximate density-functional theory based on a minimal valence basis set and neglect of all but two center integrals. We report the results of our tests of a recent long-range correction (lc) [A. Humeniuk and R. Mitrić, J. Chem. Phys. 143, 134120 (2015)] for time-dependent (TD) lc-DFTB by carrying out TD-lc-DFTB fewest switches surface hopping calculations of energy and charge transfer times using the relatively new DFTBABY [A. Humeniuk and R. Mitrić, Comput. Phys. Commun. 221, 174 (2017)] program. An advantage of this method is the ability to run enough trajectories to get meaningful ensemble averages. Our interest in the present work is less in determining exact energy and charge transfer rates than in understanding how the results of these calculations vary with the value of the range-separation parameter (Rlc = 1/μ) for a model organic solar cell heterojunction consisting of a gas-phase van der Waals complex P/F made up of a single pentacene (P) molecule together with a single buckminsterfullerene (F) molecule. The default value of Rlc = 3.03 a0 is found to be much too small as neither energy nor charge transfer is observed until Rlc ≈ 10 a0. Tests at a single geometry show that the best agreement with high-quality ab initio spectra is obtained in the limit of no lc (i.e., very large Rlc). A plot of energy and charge transfer rates as a function of Rlc is provided, which suggests that a value of Rlc ≈ 15 a0 yields the typical literature (condensed-phase) charge transfer time of about 100 fs. However, energy and charge transfer times become as high as ∼300 fs for Rlc ≈ 25 a0. A closer examination of the charge transfer process P*/F → P+/F- shows that the initial electron transfer is accompanied by a partial delocalization of the P hole onto F, which then relocalizes back onto P, consistent with a polaron-like picture in which the nuclei relax to stabilize the resultant redistribution of charges.
Collapse
Affiliation(s)
| | - Mark E Casida
- Laboratoire de Spectrométrie, Interactions et Chimie Théorique (SITh), Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Grenoble-Alpes, 301 rue de la Chimie, CS 40700, 38058 Grenoble Cedex 9, France
| | - Xi Zhu
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Bhaarathi Natarajan
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Haibin Su
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Alexander Humeniuk
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Evgenii Titov
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Xincheng Miao
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| |
Collapse
|
178
|
Smith B, Shakiba M, Akimov AV. Nonadiabatic Dynamics in Si and CdSe Nanoclusters: Many-Body vs Single-Particle Treatment of Excited States. J Chem Theory Comput 2021; 17:678-693. [DOI: 10.1021/acs.jctc.0c01009] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| | - Mohammad Shakiba
- Department of Materials Science and Engineering, Shahid Bahonar University of Kerman, Kerman 76169-14111, Iran
| | - Alexey V. Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| |
Collapse
|
179
|
Rodríguez-Hernández B, Nelson T, Oldani N, Martínez-Mesa A, Uranga-Piña L, Segawa Y, Tretiak S, Itami K, Fernandez-Alberti S. Exciton Spatial Dynamics and Self-Trapping in Carbon Nanocages. J Phys Chem Lett 2021; 12:224-231. [PMID: 33326240 DOI: 10.1021/acs.jpclett.0c03364] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Three-dimensional cage-shaped molecules formed from chainlike structures hold potential as unique optoelectronic materials and host compounds. Their optical, structural, and dynamical features are tunable by changes in shape and size. We perform a comparison of these properties for three sizes of strained conjugated [n.n.n]carbon nanocages composed of three paraphenylene chains (bridges) of length n = 4, 5, or 6. The exciton intramolecular redistribution occurring during nonradiative relaxation has been explored using nonadiabatic excited-state molecular dynamics. Our results provide atomistic insight into the conformational features associated with the observed red- and blue-shift trends in the absorption and fluorescence spectra, respectively, with increasing nanocage size. Their internal conversion processes involve intramolecular energy transfer that leads to exciton self-trapping on a few phenylene units at the center of a single bridge. The dependence of these dynamical features on the size of the nanocage can be used to tune their host-guest chemical properties and their use for organic electronics and catenane-like applications.
Collapse
Affiliation(s)
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nicolas Oldani
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Aliezer Martínez-Mesa
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
- DynAMoS (Dynamical processes in Atomic and Molecular Systems), Facultad de Física, Universidad de La Habana, San Lázaro y L, La Habana 10400, Cuba
| | - Llinersy Uranga-Piña
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
- DynAMoS (Dynamical processes in Atomic and Molecular Systems), Facultad de Física, Universidad de La Habana, San Lázaro y L, La Habana 10400, Cuba
| | - Yasutomo Segawa
- Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- JST, ERATO, Itami Molecular Nanocarbon Project, Nagoya University, Nagoya 464-8602, Japan
- Institute for Molecular Science, Myodaiji, Okazaki 444-8787, Japan
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Myodaiji, Okazaki 444-8787, Japan
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Kenichiro Itami
- Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- JST, ERATO, Itami Molecular Nanocarbon Project, Nagoya University, Nagoya 464-8602, Japan
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Nagoya 464-8602, Japan
| | | |
Collapse
|
180
|
Abstract
Coordination compounds, characterized by fascinating and tunable electronic properties, are capable of binding easily to proteins, polymers, wires and DNA. Upon irradiation, these molecular systems develop functions finding applications in solar cells, photocatalysis, luminescent and conformational probes, electron transfer triggers and diagnostic or therapeutic tools. The control of these functions is activated by the light wavelength, the metal/ligand cooperation and the environment within the first picoseconds (ps). After a brief summary of the theoretical background, this perspective reviews case studies, from 1st row to 3rd row transition metal complexes, that illustrate how spin-orbit, vibronic coupling and quantum effects drive the photophysics of this class of molecules at the early stage of the photoinduced elementary processes within the fs-ps time scale range.
Collapse
Affiliation(s)
- Chantal Daniel
- Laboratoire de Chimie Quantique, Université de Strasbourg, CNRS UMR7177, Institut Le Bel, 4 Rue Blaise Pascal, 67000 Strasbourg, France.
| |
Collapse
|
181
|
De Sio A, Sommer E, Nguyen XT, Groß L, Popović D, Nebgen BT, Fernandez-Alberti S, Pittalis S, Rozzi CA, Molinari E, Mena-Osteritz E, Bäuerle P, Frauenheim T, Tretiak S, Lienau C. Intermolecular conical intersections in molecular aggregates. NATURE NANOTECHNOLOGY 2021; 16:63-68. [PMID: 33199882 DOI: 10.1038/s41565-020-00791-2] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Accepted: 10/07/2020] [Indexed: 06/11/2023]
Abstract
Conical intersections (CoIns) of multidimensional potential energy surfaces are ubiquitous in nature and control pathways and yields of many photo-initiated intramolecular processes. Such topologies can be potentially involved in the energy transport in aggregated molecules or polymers but are yet to be uncovered. Here, using ultrafast two-dimensional electronic spectroscopy (2DES), we reveal the existence of intermolecular CoIns in molecular aggregates relevant for photovoltaics. Ultrafast, sub-10-fs 2DES tracks the coherent motion of a vibrational wave packet on an optically bright state and its abrupt transition into a dark state via a CoIn after only 40 fs. Non-adiabatic dynamics simulations identify an intermolecular CoIn as the source of these unusual dynamics. Our results indicate that intermolecular CoIns may effectively steer energy pathways in functional nanostructures for optoelectronics.
Collapse
Affiliation(s)
- Antonietta De Sio
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany.
| | - Ephraim Sommer
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany
| | - Xuan Trung Nguyen
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany
| | - Lynn Groß
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - Duško Popović
- Institut für Organische Chemie II und Neue Materialien, Universität Ulm, Ulm, Germany
| | | | - Sebastian Fernandez-Alberti
- National University of Quilmes/CONICET, Department of Science and Technology, Bernal (B1876BXD), Buenos Aires Province, Argentina
| | | | | | - Elisa Molinari
- Istituto Nanoscienze-CNR, Modena, Italy
- Università di Modena e Reggio Emilia, Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Modena, Italy
| | - Elena Mena-Osteritz
- Institut für Organische Chemie II und Neue Materialien, Universität Ulm, Ulm, Germany
| | - Peter Bäuerle
- Institut für Organische Chemie II und Neue Materialien, Universität Ulm, Ulm, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
- Computational Science Research Center, Beijing and Computational Science and Applied Research Institute Shenzhen, Shenzhen, China
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Christoph Lienau
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany
| |
Collapse
|
182
|
Ollitrault PJ, Mazzola G, Tavernelli I. Nonadiabatic Molecular Quantum Dynamics with Quantum Computers. PHYSICAL REVIEW LETTERS 2020; 125:260511. [PMID: 33449795 DOI: 10.1103/physrevlett.125.260511] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 12/10/2020] [Indexed: 05/11/2023]
Abstract
The theoretical investigation of nonadiabatic processes is hampered by the complexity of the coupled electron-nuclear dynamics beyond the Born-Oppenheimer approximation. Classically, the simulation of such reactions is limited by the unfavorable scaling of the computational resources as a function of the system size. While quantum computing exhibits proven quantum advantage for the simulation of real-time dynamics, the study of quantum algorithms for the description of nonadiabatic phenomena is still unexplored. In this Letter, we propose a quantum algorithm for the simulation of fast nonadiabatic chemical processes together with an initialization scheme for quantum hardware calculations. In particular, we introduce a first-quantization method for the time evolution of a wave packet on two coupled harmonic potential energy surfaces (Marcus model). In our approach, the computational resources scale polynomially in the system dimensions, opening up new avenues for the study of photophysical processes that are classically intractable.
Collapse
Affiliation(s)
- Pauline J Ollitrault
- IBM Quantum, IBM Research-Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Guglielmo Mazzola
- IBM Quantum, IBM Research-Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| | - Ivano Tavernelli
- IBM Quantum, IBM Research-Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| |
Collapse
|
183
|
Mukazhanova A, Malone W, Negrin-Yuvero H, Fernandez-Alberti S, Tretiak S, Sharifzadeh S. Photoexcitation dynamics in perylene diimide dimers. J Chem Phys 2020; 153:244117. [PMID: 33380092 DOI: 10.1063/5.0031485] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We utilize first-principles theory to investigate photo-induced excited-state dynamics of functionalized perylene diimide. This class of materials is highly suitable for solar energy conversion because of the strong optical absorbance, efficient energy transfer, and chemical tunability. We couple time-dependent density functional theory to a recently developed time-resolved non-adiabatic dynamics approach based on a semi-empirical description. By studying the monomer and dimer, we focus on the role stacking plays on the time-scales associated with excited-state non-radiative relaxation from a high excitonic state to the lowest energy exciton. We predict that the time-scale for energy conversion in the dimer is significantly faster than that in the monomer when equivalent excited states are accounted for. Additionally, for the dimer, the decay from the second to the nearly degenerate lowest energy excited-state involves two time-scales: a rapid decay on the order of ∼10 fs followed by a slower decay of ∼100 fs. Analysis of the spatial localization of the electronic transition density during the internal conversion process points out the existence of localized states on individual monomers, indicating that the strength of thermal fluctuations exceeds electronic couplings between the states such that the exciton hops between localized states throughout the simulation.
Collapse
Affiliation(s)
- Aliya Mukazhanova
- Division of Materials Science and Engineering, Boston University, Boston, Massachusetts 02215, USA
| | - Walter Malone
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Hassiel Negrin-Yuvero
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | | | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Sahar Sharifzadeh
- Division of Materials Science and Engineering, Boston University, Boston, Massachusetts 02215, USA
| |
Collapse
|
184
|
Freixas VM, Nelson T, Ondarse-Alvarez D, Nijjar P, Mikhailovsky A, Zhou C, Fernandez-Alberti S, Bazan GC, Tretiak S. Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful? J Chem Phys 2020; 153:244114. [PMID: 33380074 DOI: 10.1063/5.0028126] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born-Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. Here, we address this issue by modeling the electronic energy transfer in a donor-acceptor-donor (D-A-D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are "sufficiently good" to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.
Collapse
Affiliation(s)
- Victor M Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | | - Parmeet Nijjar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA
| | - Alexander Mikhailovsky
- Department of Chemistry and Biochemistry, Center for Polymers and Organic Solids, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Cheng Zhou
- Department of Chemistry, National University of Singapore, Singapore 117543, Singapore
| | | | - Guillermo C Bazan
- Department of Chemistry and Biochemistry, Center for Polymers and Organic Solids, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
185
|
Negrin-Yuvero H, Freixas VM, Rodriguez-Hernandez B, Rojas-Lorenzo G, Tretiak S, Bastida A, Fernandez-Alberti S. Photoinduced Dynamics with Constrained Vibrational Motion: FrozeNM Algorithm. J Chem Theory Comput 2020; 16:7289-7298. [PMID: 33201709 DOI: 10.1021/acs.jctc.0c00930] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Ab initio molecular dynamics (AIMD) simulation, analyzed in terms of vibrational normal modes, is a widely used technique that facilitates understanding of complex structural motions and coupling between electronic and nuclear degrees of freedom. Usually, only a subset of vibrations is directly involved in the process of interest. The impact of these vibrations can be evaluated by performing AIMD simulations by selectively freezing certain motions. Herein, we present frozen normal mode (FrozeNM), a new algorithm to apply normal-mode constraints in AIMD simulations, as implemented in the nonadiabatic excited state molecular dynamics code. We further illustrate its capacity by analyzing the impact of normal-mode constraints on the photoinduced energy transfer between polyphenylene ethynylene dendrimer building blocks. Our results show that the electronic relaxation can be significantly slowed down by freezing a well-selected small subset of active normal modes characterized by their contributions in the direction of energy transfer. The application of these constraints reduces the nonadiabatic coupling between electronic excited states during the entire dynamical simulations. Furthermore, we validate reduced dimensionality models by freezing all the vibrations, except a few active modes. Altogether, we consider FrozeNM as a useful tool that can be broadly used to underpin the role of vibrational motion in a studied process and to formulate reduced models that describe essential physical phenomena.
Collapse
Affiliation(s)
- H Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - V M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - B Rodriguez-Hernandez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - G Rojas-Lorenzo
- Departamento de Física Atómica y Molecular, Instituto Superior de Tecnologías y Ciencias Aplicadas, Universidad de La Habana , La Habana, Cuba
| | - S Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - A Bastida
- Departamento de Química Física, Universidad de Murcia, Murcia 30100, Spain
| | - S Fernandez-Alberti
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| |
Collapse
|
186
|
Kim HW, Rhee YM. Two-oscillator mapping modification of the Poisson bracket mapping equation formulation of the quantum-classical Liouville equation. J Chem Phys 2020; 153:214103. [PMID: 33291890 DOI: 10.1063/5.0027799] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Mapping basis solutions provide efficient ways for simulating mixed quantum-classical (MQC) dynamics in complex systems by matching multiple quantum states of interest to some fictitious physical states. Recently, various MQC methods were devised such that two harmonic oscillators are employed to represent each electronic state, showing improvements over one-oscillator-based methods. Here, we introduce and analyze newly modified mapping approximations of the quantum-classical Liouville equation (QCLE) using two oscillators for each electronic state. We design two separate mapping relations that we can adopt toward simulating dynamics and computing expectation values. Through the process, two MQC methods can be constructed, one of which actually reproduces the population dynamics of the forward and backward trajectory solution of QCLE. By applying the methods to spin-boson systems with a range of parameters, we find out that the choice of mapping relations greatly affects the simulation results. We also show that further improvement is possible through using modified identity operator formulations. Our findings may be helpful in constructing improved MQC methods in the future.
Collapse
Affiliation(s)
- Hyun Woo Kim
- Center for Molecular Modeling and Simulation, Korea Research Institute of Chemical Technology (KRICT), Daejeon 34114, South Korea
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| |
Collapse
|
187
|
Kuan KY, Singleton DA. Vibrationally Hot and Cold Triplets. Sensitizer-Dependent Dynamics and Localized Vibrational Promotion of a Di-π-methane Rearrangement. J Am Chem Soc 2020; 142:19885-19888. [PMID: 33179917 PMCID: PMC8522498 DOI: 10.1021/jacs.0c10468] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Large intramolecular 13C kinetic isotope effects (KIEs) for the di-π-methane rearrangement of benzobarrelene fit with statistical expectations from heavy-atom tunneling when a low-energy sensitizer is employed, but much lower KIEs are observed with higher-energy sensitizers. These results in combination with trajectory studies suggest that the excess vibrational energy available from triplet energy transfer leads to hot and nonstatistical dynamics in the rearrangement.
Collapse
Affiliation(s)
- Kai-Yuan Kuan
- Department of Chemistry, Texas A&M University, P.O. Box 30012, College Station, Texas 77842, United States
| | - Daniel A Singleton
- Department of Chemistry, Texas A&M University, P.O. Box 30012, College Station, Texas 77842, United States
| |
Collapse
|
188
|
Uratani H, Morioka T, Yoshikawa T, Nakai H. Fast Nonadiabatic Molecular Dynamics via Spin-Flip Time-Dependent Density-Functional Tight-Binding Approach: Application to Nonradiative Relaxation of Tetraphenylethylene with Locked Aromatic Rings. J Chem Theory Comput 2020; 16:7299-7313. [PMID: 33197192 DOI: 10.1021/acs.jctc.0c00936] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Nonadiabatic dynamics around conical intersections between ground and excited states are crucial to understand excited-state phenomena in complex chemical systems. With this background in mind, we present an approach combining fewest-switches trajectory surface hopping and spin-flip (SF) time-dependent (TD) density-functional tight binding (DFTB), which is a simplified version of SF-TD density functional theory (DFT) with semiempirical parametrizations, for computationally efficient nonadiabatic molecular dynamics simulations. The estimated computational time of the SF-TD-DFTB approach is several orders of magnitude lower than that of SF-TD-DFT. In addition, the proposed method reproduces the time scales and quantum yields in photoisomerization reactions of azobenzene at a level comparable with conventional ab initio approaches, demonstrating reasonable accuracy. Finally, we report a practical application of the developed technique to explore the nonradiative relaxation processes of tetraphenylethylene and its derivative with torsionally locked aromatic rings and discuss the effect of locking the rings on the excited-state lifetime.
Collapse
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Toshiki Morioka
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
189
|
Di Maiolo F, Brey D, Binder R, Burghardt I. Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature. J Chem Phys 2020; 153:184107. [PMID: 33187420 DOI: 10.1063/5.0027588] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.
Collapse
Affiliation(s)
- Francesco Di Maiolo
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Dominik Brey
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
190
|
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
|
191
|
Wang S, Kind T, Tantillo DJ, Fiehn O. Predicting in silico electron ionization mass spectra using quantum chemistry. J Cheminform 2020; 12:63. [PMID: 33372633 PMCID: PMC7576811 DOI: 10.1186/s13321-020-00470-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Accepted: 10/13/2020] [Indexed: 12/20/2022] Open
Abstract
Compound identification by mass spectrometry needs reference mass spectra. While there are over 102 million compounds in PubChem, less than 300,000 curated electron ionization (EI) mass spectra are available from NIST or MoNA mass spectral databases. Here, we test quantum chemistry methods (QCEIMS) to generate in silico EI mass spectra (MS) by combining molecular dynamics (MD) with statistical methods. To test the accuracy of predictions, in silico mass spectra of 451 small molecules were generated and compared to experimental spectra from the NIST 17 mass spectral library. The compounds covered 43 chemical classes, ranging up to 358 Da. Organic oxygen compounds had a lower matching accuracy, while computation time exponentially increased with molecular size. The parameter space was probed to increase prediction accuracy including initial temperatures, the number of MD trajectories and impact excess energy (IEE). Conformational flexibility was not correlated to the accuracy of predictions. Overall, QCEIMS can predict 70 eV electron ionization spectra of chemicals from first principles. Improved methods to calculate potential energy surfaces (PES) are still needed before QCEIMS mass spectra of novel molecules can be generated at large scale.
Collapse
Affiliation(s)
- Shunyang Wang
- West Coast Metabolomics Center, UC Davis Genome Center, University of California, 451 Health Sciences Drive, Davis, CA, 95616, USA.,Department of Chemistry, University of California, 1 Shields Ave, Davis, CA, 95616, USA
| | - Tobias Kind
- West Coast Metabolomics Center, UC Davis Genome Center, University of California, 451 Health Sciences Drive, Davis, CA, 95616, USA
| | - Dean J Tantillo
- Department of Chemistry, University of California, 1 Shields Ave, Davis, CA, 95616, USA
| | - Oliver Fiehn
- West Coast Metabolomics Center, UC Davis Genome Center, University of California, 451 Health Sciences Drive, Davis, CA, 95616, USA.
| |
Collapse
|
192
|
Abstract
When nonadiabatic dynamics are described on the basis of trajectories, severe trajectory branching occurs when the nuclear wave packets on some potential energy surfaces are reflected while those on the remaining surfaces are not. As a result, the traditional Ehrenfest mean field (EMF) approximation breaks down. In this study, two versions of the branching corrected mean field (BCMF) method are proposed. Namely, when trajectory branching is identified, BCMF stochastically selects either the reflected or the nonreflected group to build the new mean field trajectory or splits the mean field trajectory into two new trajectories with the corresponding weights. As benchmarked in six standard model systems and an extensive model base with two hundred diverse scattering models, BCMF significantly improves the accuracy while retaining the high efficiency of the traditional EMF. In fact, BCMF closely reproduces the exact quantum dynamics in all investigated systems, thus highlighting the essential role of branching correction in nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
- Jiabo Xu
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
193
|
Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models. APPLIED SCIENCES-BASEL 2020. [DOI: 10.3390/app10196821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In biological and life science applications, photosynthesis is an important process that involves the absorption and transformation of sunlight into chemical energy. During the photosynthesis process, the light photons are captured by the green chlorophyll pigments in their photosynthetic antennae and further funneled to the reaction center. One of the most important light harvesting complexes that are highly important in the study of photosynthesis is the membrane-attached Fenna–Matthews–Olson (FMO) complex found in the green sulfur bacteria. In this review, we discuss the mathematical formulations and computational modeling of some of the light harvesting complexes including FMO. The most recent research developments in the photosynthetic light harvesting complexes are thoroughly discussed. The theoretical background related to the spectral density, quantum coherence and density functional theory has been elaborated. Furthermore, details about the transfer and excitation of energy in different sites of the FMO complex along with other vital photosynthetic light harvesting complexes have also been provided. Finally, we conclude this review by providing the current and potential applications in environmental science, energy, health and medicine, where such mathematical and computational studies of the photosynthesis and the light harvesting complexes can be readily integrated.
Collapse
|
194
|
Strong SE, Hestand NJ. Modeling nonlocal electron-phonon coupling in organic crystals using interpolative maps: The spectroscopy of crystalline pentacene and 7,8,15,16-tetraazaterrylene. J Chem Phys 2020; 153:124113. [PMID: 33003728 DOI: 10.1063/5.0021731] [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
Electron-phonon coupling plays a central role in the transport properties and photophysics of organic crystals. Successful models describing charge- and energy-transport in these systems routinely include these effects. Most models for describing photophysics, on the other hand, only incorporate local electron-phonon coupling to intramolecular vibrational modes, while nonlocal electron-phonon coupling is neglected. One might expect nonlocal coupling to have an important effect on the photophysics of organic crystals because it gives rise to large fluctuation in the charge-transfer couplings, and charge-transfer couplings play an important role in the spectroscopy of many organic crystals. Here, we study the effects of nonlocal coupling on the absorption spectrum of crystalline pentacene and 7,8,15,16-tetraazaterrylene. To this end, we develop a new mixed quantum-classical approach for including nonlocal coupling into spectroscopic and transport models for organic crystals. Importantly, our approach does not assume that the nonlocal coupling is linear, in contrast to most modern charge-transport models. We find that the nonlocal coupling broadens the absorption spectrum non-uniformly across the absorption line shape. In pentacene, for example, our model predicts that the lower Davydov component broadens considerably more than the upper Davydov component, explaining the origin of this experimental observation for the first time. By studying a simple dimer model, we are able to attribute this selective broadening to correlations between the fluctuations of the charge-transfer couplings. Overall, our method incorporates nonlocal electron-phonon coupling into spectroscopic and transport models with computational efficiency, generalizability to a wide range of organic crystals, and without any assumption of linearity.
Collapse
Affiliation(s)
- Steven E Strong
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| | - Nicholas J Hestand
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
195
|
Chen WC, Cheng YC. Elucidating the Magnitude of Internal Reorganization Energy of Molecular Excited States from the Perspective of Transition Density. J Phys Chem A 2020; 124:7644-7657. [PMID: 32864966 DOI: 10.1021/acs.jpca.0c06482] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Quantifying vibronic couplings in molecular excited states is crucial for the elucidation of a broad range of photophysical phenomena. In this study, we compare different theoretical approaches for the calculation of reorganization energy, a measure of vibronic coupling strength, and provide a rigorous derivation to show that molecular transition density characterizing electron-hole excitation could be used to quantify the magnitude of reorganization energy. The theory enables a descriptor based on molecular-orbital coefficients and atomic transition densities to quantify the magnitude of reorganization energies in molecular excited states. Applying the approach to low-lying excited states of polyacenes, we demonstrate that transition density distribution explains the difference in the magnitude of the reorganization energy of different excited states. Furthermore, to clarify the applicability of the transition density descriptor in molecular design for small-reorganization energy molecules, we investigate a broad range of molecular chromophores to show the effectiveness of the proposed theory. With this perspective on the relationship between reorganization energy and transition density, we successfully provide a quantitative rule to identify π-conjugated systems with small reorganization energy in the excited state, which should be useful for the development of novel optoelectronic materials.
Collapse
Affiliation(s)
- Wei-Chih Chen
- Department of Chemistry and Center for Quantum Science and Engineering, National Taiwan University, Taipei City 106, Taiwan, R.O.C
| | - Yuan-Chung Cheng
- Department of Chemistry and Center for Quantum Science and Engineering, National Taiwan University, Taipei City 106, Taiwan, R.O.C
| |
Collapse
|
196
|
Esch MP, Levine BG. Decoherence-corrected Ehrenfest molecular dynamics on many electronic states. J Chem Phys 2020; 153:114104. [DOI: 10.1063/5.0022529] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
197
|
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]
|
198
|
Conti I, Cerullo G, Nenov A, Garavelli M. Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going. J Am Chem Soc 2020; 142:16117-16139. [PMID: 32841559 PMCID: PMC7901644 DOI: 10.1021/jacs.0c04952] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Computational spectroscopy is becoming a mandatory tool for the interpretation of the
complex, and often congested, spectral maps delivered by modern non-linear multi-pulse
techniques. The fields of Electronic Structure Methods,
Non-Adiabatic Molecular Dynamics, and Theoretical
Spectroscopy represent the three pillars of the virtual ultrafast
optical spectrometer, able to deliver transient spectra in
silico from first principles. A successful simulation strategy requires a
synergistic approach that balances between the three fields, each one having its very
own challenges and bottlenecks. The aim of this Perspective is to demonstrate that,
despite these challenges, an impressive agreement between theory and experiment is
achievable now regarding the modeling of ultrafast photoinduced processes in complex
molecular architectures. Beyond that, some key recent developments in the three fields
are presented that we believe will have major impacts on spectroscopic simulations in
the very near future. Potential directions of development, pending challenges, and
rising opportunities are illustrated.
Collapse
Affiliation(s)
- Irene Conti
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Giulio Cerullo
- Dipartimento di Fisica, Politecnico di Milano, IFN-CNR, Piazza Leonardo da Vinci 32, I-20133 Milano, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| |
Collapse
|
199
|
Kossoski F, Barbatti M. Nonadiabatic dynamics in multidimensional complex potential energy surfaces. Chem Sci 2020; 11:9827-9835. [PMID: 34094243 PMCID: PMC8162122 DOI: 10.1039/d0sc04197a] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Despite the continuous development of theoretical methodologies for describing nonadiabatic dynamics of molecular systems, there is a lack of approaches for processes where the norm of the wave function is not conserved, i.e., when an imaginary potential accounts for some irreversible decaying mechanism. Current approaches rely on building potential energy surfaces of reduced dimensionality, which is not optimal for more involving and realistic multidimensional problems. Here, we present a novel methodology for describing the dynamics of complex-valued molecular Hamiltonians, which is a generalisation of the trajectory surface hopping method. As a first application, the complex surface fewest switches surface hopping (CS-FSSH) method was employed to survey the relaxation mechanisms of the shape resonant anions of iodoethene. We have provided the first detailed and dynamical picture of the π*/σ* mechanism of dissociative electron attachment in halogenated unsaturated compounds, which is believed to underlie electron-induced reactions of several molecules of interest. Electron capture into the π* orbital promotes C
Created by potrace 1.16, written by Peter Selinger 2001-2019
]]>
C stretching and out-of-plane vibrations, followed by charge transfer from the double bond into the σ* orbital at the C–I bond, and, finally, release of the iodine ion, all within only 15 fs. On-the-fly dynamics simulations of a vast class of processes can be envisioned with the CS-FSSH methodology, including autoionisation from transient anions, core-ionised and superexcited states, Auger and interatomic coulombic decay, and time-dependent luminescence. Despite the continuous development of methods for describing nonadiabatic dynamics, there is a lack of multidimensional approaches for processes where the wave function norm is not conserved. A new surface hopping variant closes this knowledge gap.![]()
Collapse
|
200
|
Ichiryu K, Naito T. Steady-state Observation of Interacting Photochemical Species in Photoexcited Solid States. CHEM LETT 2020. [DOI: 10.1246/cl.200335] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Kohsuke Ichiryu
- Graduate School of Science and Engineering, Ehime University, Matsuyama, Ehime 790-8577, Japan
| | - Toshio Naito
- Graduate School of Science and Engineering, Ehime University, Matsuyama, Ehime 790-8577, Japan
- Geodynamics Research Center (GRC), Ehime University, Matsuyama, Ehime 790-8577, Japan
- Research Unit for Development of Organic Superconductors, Ehime University, Matsuyama, Ehime 790-8577, Japan
| |
Collapse
|