1
|
Anglada JM, Martins-Costa MTC, Francisco JS, Ruiz-López MF. Triplet State Radical Chemistry: Significance of the Reaction of 3SO 2 with HCOOH and HNO 3. J Am Chem Soc 2024; 146:14297-14306. [PMID: 38722613 PMCID: PMC11117184 DOI: 10.1021/jacs.4c03938] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2024] [Revised: 04/25/2024] [Accepted: 04/29/2024] [Indexed: 05/23/2024]
Abstract
The triplet excited states of sulfur dioxide can be accessed in the UV region and have a lifetime large enough that they can react with atmospheric trace gases. In this work, we report high level ab initio calculations for the reaction of the a3B1 and b3A2 excited states of SO2 with weak and strong acidic species such as HCOOH and HNO3, aimed to extend the chemistry reported in previous studies with nonacidic H atoms (water and alkanes). The reactions investigated in this work are very versatile and follow different kinds of mechanisms, namely, proton-coupled electron transfer (pcet) and conventional hydrogen atom transfer (hat) mechanisms. The study provides new insights into a general and very important class of excited-state-promoted reactions, opening up interesting chemical perspectives for technological applications of photoinduced H-transfer reactions. It also reveals that atmospheric triplet chemistry is more significant than previously thought.
Collapse
Affiliation(s)
- Josep M. Anglada
- Departament
de Química Biològica (IQAC − CSIC), c/Jordi Girona 18, Barcelona E-08034, Spain
| | - Marilia T. C. Martins-Costa
- Laboratoire
de Physique et Chimie Théoriques, UMR CNRS 7019, University of Lorraine, CNRS, BP 70239, Vandoeuvre-lès-Nancy 54506, France
| | - Joseph S. Francisco
- Department
of Earth and Environmental Science and Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6316, United States
| | - Manuel F. Ruiz-López
- Laboratoire
de Physique et Chimie Théoriques, UMR CNRS 7019, University of Lorraine, CNRS, BP 70239, Vandoeuvre-lès-Nancy 54506, France
| |
Collapse
|
2
|
Wu Y, Sun J, Li Z, Zhang Z, Luo Z, Chang Y, Wu G, Zhang W, Yu S, Yuan K, Yang X. Photodissociation dynamics of SO2 via the G̃1B1 state: The O(1D2) and O(1S0) product channels. J Chem Phys 2024; 160:164311. [PMID: 38661196 DOI: 10.1063/5.0208090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2024] [Accepted: 04/09/2024] [Indexed: 04/26/2024] Open
Abstract
Produced by both nature and human activities, sulfur dioxide (SO2) is an important species in the earth's atmosphere. SO2 has also been found in the atmospheres of other planets and satellites in the solar system. The photoabsorption cross sections and photodissociation of SO2 have been studied for several decades. In this paper, we reported the experimental results for photodissociation dynamics of SO2 via the G̃1B1 state. By analyzing the images from the time-sliced velocity map ion imaging method, the vibrational state population distributions and anisotropy parameters were obtained for the O(1D2) + SO(X3Σ-, a1Δ, b1Σ+) and O(1S0) + SO(X3Σ-) channels, and the branching ratios for the channels O(1D2) + SO(X3Σ-), O(1D2) + SO(a1Δ), and O(1D2) + SO(b1Σ+) were determined to be ∼0.3, ∼0.6, and ∼0.1, respectively. The SO products were dominant in electronically and rovibrationally excited states, which may have yet unrecognized roles in the upper planetary atmosphere.
Collapse
Affiliation(s)
- Yucheng Wu
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jitao Sun
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zhenxing Li
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
| | - Zhaoxue Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
| | - Zijie Luo
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
| | - Yao Chang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
| | - Guorong Wu
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
| | - Weiqing Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
| | - Shengrui Yu
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou, Zhejiang 311231, China
| | - Kaijun Yuan
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Hefei National Laboratory, Hefei 230088, China
| | - Xueming Yang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, China
- Hefei National Laboratory, Hefei 230088, China
- Department of Chemistry and Center for Advanced Light Source Research, College of Science, Southern University of Science and Technology, Shenzhen 518055, China
| |
Collapse
|
3
|
Li Z, Fu YL, Luo Z, Yang S, Wu Y, Wu H, Wu G, Zhang W, Fu B, Yuan K, Zhang D, Yang X. Roaming in highly excited states: The central atom elimination of triatomic molecule decomposition. Science 2024; 383:746-750. [PMID: 38359138 DOI: 10.1126/science.adn3357] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Accepted: 01/18/2024] [Indexed: 02/17/2024]
Abstract
Chemical reactions are generally assumed to proceed from reactants to products along the minimum energy path (MEP). However, straying from the MEP-roaming-has been recognized as an unconventional reaction mechanism and found to occur in both the ground and first excited states. Its existence in highly excited states is however not yet established. We report a dissociation channel to produce electronically excited fragments, S(1D)+O2(a1Δg), from SO2 photodissociation in highly excited states. The results revealed two dissociation pathways: One proceeds through the MEP to produce vibrationally colder O2(a1Δg) and the other yields vibrationally hotter O2(a1Δg) by means of a roaming pathway involving an intramolecular O abstraction during reorientation motion. Such roaming dynamics may well be the rule rather than the exception for molecular photodissociation through highly excited states.
Collapse
Affiliation(s)
- Zhenxing Li
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
| | - Yan-Lin Fu
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
| | - Zijie Luo
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
- Marine Engineering College, Dalian Maritime University, Liaoning 116026, P. R. China
| | - Shuaikang Yang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
| | - Yucheng Wu
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Hao Wu
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Guorong Wu
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
| | - Weiqing Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| | - Kaijun Yuan
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| | - Donghui Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| | - Xueming Yang
- State Key Laboratory of Molecular Reaction Dynamics and Dalian Coherent Light Source, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Road, Dalian 116023, P. R. China
- Department of Chemistry and Center for Advanced Light Source Research, College of Science, Southern University of Science and Technology, Shenzhen 518055, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| |
Collapse
|
4
|
Akimov AV. Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2023; 14:11673-11683. [PMID: 38109379 DOI: 10.1021/acs.jpclett.3c03029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| |
Collapse
|
5
|
Domenianni LI, Bauer M, Schmidt-Räntsch T, Lindner J, Schneider S, Vöhringer P. Photoinduced Metallonitrene Formation by N 2 Elimination from Azide Diradical Ligands. Angew Chem Int Ed Engl 2023; 62:e202309618. [PMID: 37549374 DOI: 10.1002/anie.202309618] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Revised: 08/04/2023] [Accepted: 08/07/2023] [Indexed: 08/09/2023]
Abstract
Transition-metal nitrides/nitrenes are highly promising reagents for catalytic nitrogen-atom-transfer reactivity. They are typically prepared in situ upon optically induced N2 elimination from azido precursors. A full exploitation of their catalytic potential, however, requires in-depth knowledge of the primary photo-induced processes and the structural/electronic factors mediating the N2 loss with birth of the terminal metal-nitrogen core. Using femtosecond infrared spectroscopy, we elucidate here the primary molecular-level mechanisms responsible for the formation of a unique platinum(II) nitrene with a triplet ground state from a closed-shell platinum(II) azide precursor. The spectroscopic data in combination with quantum-chemical calculations provide compelling evidence that product formation requires the initial occupation of a singlet excited state with an anionic azide diradical ligand that is bound to a low-spin d8 -configured PtII ion. Subsequent intersystem crossing generates the Pt-bound triplet azide diradical, which smoothly evolves into the triplet nitrene via N2 loss in a near barrierless adiabatic dissociation. Our data highlight the importance of the productive, N2 -releasing state possessing azide ππ* character as a design principle for accessing efficient N-atom-transfer catalysts.
Collapse
Affiliation(s)
- Luis I Domenianni
- Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115, Bonn, Germany
| | - Markus Bauer
- Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115, Bonn, Germany
| | - Till Schmidt-Räntsch
- Georg-August-Universität Göttingen, Institut für Anorganische Chemie, Tammannstraße 4, 37077, Göttingen, Germany
| | - Jörg Lindner
- Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115, Bonn, Germany
| | - Sven Schneider
- Georg-August-Universität Göttingen, Institut für Anorganische Chemie, Tammannstraße 4, 37077, Göttingen, Germany
| | - Peter Vöhringer
- Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115, Bonn, Germany
| |
Collapse
|
6
|
Pourestarabadi S, Dehestani M. Non-adiabatic coupling in the potential energy surfaces of SO 2 molecule. Phys Chem Chem Phys 2023; 25:24526-24538. [PMID: 37661660 DOI: 10.1039/d3cp02127k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2023]
Abstract
To investigate the potential energy surfaces and the coupling between the adiabatic states of SO2 molecules, it is necessary to consider the non-adiabatic coupling terms (NACTs), where the Born-Oppenheimer approximation breaks down. In this work, we analyze the conical intersections between 1 1A1 and 1 1B2 states (the A' states in Cs symmetry) and 1 1A2 and 1 1B1 states (the A'' states in Cs symmetry) using NACTs and adiabatic-to-diabatic transformation (ADT) angles. Our results confirm reasonable interaction between 1 1A1 and 1 1B2 states and strong interaction between 1 1A2 and 1 1B1 states.
Collapse
Affiliation(s)
- Sedigheh Pourestarabadi
- Department of Chemistry, Shahid Bahonar University of Kerman, Kerman, Iran.
- Young Researchers Society, Shahid Bahonar University of Kerman, Kerman, Iran
| | - Maryam Dehestani
- Department of Chemistry, Shahid Bahonar University of Kerman, Kerman, Iran.
| |
Collapse
|
7
|
Schnack-Petersen AK, Pápai M, Coriani S, Møller KB. A theoretical study of the time-resolved x-ray absorption spectrum of the photoionized BT-1T cation. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2023; 10:034102. [PMID: 37250952 PMCID: PMC10224778 DOI: 10.1063/4.0000183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 05/01/2023] [Indexed: 05/31/2023]
Abstract
The time-resolved x-ray absorption spectrum of the BT-1T cation (BT-1T+) is theoretically simulated in order to investigate the charge transfer reaction of the system. We employ both trajectory surface hopping and quantum dynamics to simulate the structural evolution over time and the changes in the state populations. To compute the static x-ray absorption spectra (XAS) of the ground and excited states, we apply both the time-dependent density functional theory and the coupled cluster singles and doubles method. The results obtained are in good agreement between the methods. It is, furthermore, found that the small structural changes that occur during the reaction have little effect on the static XAS. Hence, the tr-XAS can be computed based on the state populations determined from a nuclear dynamics simulation and one set of static XAS calculations, utilizing the ground state optimized geometry. This approach can save considerable computational resources, as the static spectra need not to be calculated for all geometries. As BT-1T is a relatively rigid molecule, the outlined approach should only be considered when investigating non-radiative decay processes in the vicinity of the Franck-Condon point.
Collapse
Affiliation(s)
| | | | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Klaus Braagaard Møller
- Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
8
|
Shu Y, Zhang L, Wu D, Chen X, Sun S, Truhlar DG. New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States. J Chem Theory Comput 2023; 19:2419-2429. [PMID: 37079755 DOI: 10.1021/acs.jctc.2c01173] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/22/2023]
Abstract
It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Dihua Wu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Xiye Chen
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
9
|
Penfold TJ, Eng J. Mind the GAP: quantifying the breakdown of the linear vibronic coupling Hamiltonian. Phys Chem Chem Phys 2023; 25:7195-7204. [PMID: 36820783 DOI: 10.1039/d2cp05576g] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
Excited state dynamics play a critical role across a broad range of scientific fields. Importantly, the highly non-equilibrium nature of the states generated by photoexcitation means that excited state simulations should usually include an accurate description of the coupled electronic-nuclear motion, which often requires solving the time-dependent Schrödinger equation (TDSE). One of the biggest challenges for these simulations is the requirement to calculate the PES over which the nuclei evolve. An effective approach for addressing this challenge is to use the approximate linear vibronic coupling (LVC) Hamiltonian, which enables a model potential to be parameterised using relatively few quantum chemistry calculations. However, this approach is only valid provided there are no large amplitude motions in the excited state dynamics. In this paper we introduce and deploy a metric, the global anharmonicity parameter (GAP), which can be used to assess the accuracy of an LVC potential. Following its derivation, we illustrate its utility by applying it to three molecules exhibiting different rigidity in their excited states.
Collapse
Affiliation(s)
- Thomas J Penfold
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK.
| | - Julien Eng
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK.
| |
Collapse
|
10
|
Zobel JP, Heindl M, Plasser F, Mai S, González L. Surface Hopping Dynamics on Vibronic Coupling Models. Acc Chem Res 2021; 54:3760-3771. [PMID: 34570472 PMCID: PMC8529708 DOI: 10.1021/acs.accounts.1c00485] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The simulation of photoinduced non-adiabatic dynamics is of great
relevance in many scientific disciplines, ranging from physics and
materials science to chemistry and biology. Upon light irradiation,
different relaxation processes take place in which electronic and
nuclear motion are intimately coupled. These are best described by
the time-dependent molecular Schrödinger equation, but its
solution poses fundamental practical challenges to contemporary theoretical
chemistry. Two widely used and complementary approaches to this problem
are multiconfigurational time-dependent Hartree (MCTDH) and trajectory
surface hopping (SH). MCTDH is an accurate fully quantum-mechanical
technique but often is feasible only in reduced dimensionality, in
combination with approximate vibronic coupling (VC) Hamiltonians,
or both (i.e., reduced-dimensional VC potentials). In contrast, SH
is a quantum–classical technique that neglects most nuclear
quantum effects but allows nuclear dynamics in full dimensionality
by calculating potential energy surfaces on the fly. If nuclear quantum
effects do not play a central role and a linear VC (LVC) Hamiltonian
is appropriate—e.g., for stiff molecules that generally keep
their conformation in the excited state—then it seems advantageous
to combine the efficient LVC and SH techniques. In this Account, we
describe how surface hopping based on an LVC Hamiltonian (SH/LVC)—as
recently implemented in the SHARC surface hopping package—can
provide an economical and automated approach to simulate non-adiabatic
dynamics. First, we illustrate the potential of SH/LVC in a number
of showcases, including intersystem crossing in SO2, intra-Rydberg
dynamics in acetone, and several photophysical studies on large transition-metal
complexes, which would be much more demanding or impossible to perform
with other methods. While all of the applications provide very useful
insights into light-induced phenomena, they also hint at difficulties
faced by the SH/LVC methodology that need to be addressed in the future.
Second, we contend that the SH/LVC approach can be useful to benchmark
SH itself. By the use of the same (LVC) potentials as MCTDH calculations
have employed for decades and by relying on the efficiency of SH/LVC,
it is possible to directly compare multiple SH test calculations with
a MCTDH reference and ponder the accuracy of various correction algorithms
behind the SH methodology, such as decoherence corrections or momentum
rescaling schemes. Third, we demonstrate how the efficiency of SH/LVC
can also be exploited to identify essential nuclear and electronic
degrees of freedom to be employed in more accurate MCTDH calculations.
Lastly, we show that SH/LVC is able to advance the development of
SH protocols that can describe nuclear dynamics including explicit
laser fields—a very challenging endeavor for trajectory-based
schemes. To end, this Account compiles the typical costs of contemporary
SH simulations, evidencing the great advantages of using parametrized
potentials. The LVC model is a sleeping beauty that, kissed by SH,
is fueling the field of excited-state molecular dynamics. We hope
that this Account will stimulate future research in this direction,
leveraging the advantages of the SH/VC schemes to larger extents and
extending their applicability to uncharted territories.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| | - Moritz Heindl
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| | - Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währingerstr. 19, 1090 Vienna, Austria
| |
Collapse
|
11
|
M J C, Nag P, Vennapusa SR. Surface hopping dynamics reveal ultrafast triplet generation promoted by S 1-T 2-T 1 spin-vibronic coupling in 2-mercaptobenzothiazole. Phys Chem Chem Phys 2021; 23:20183-20192. [PMID: 34473155 DOI: 10.1039/d1cp02587b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the T1 formation upon populating the optically "bright" S2 in 2-mercaptobenzothiazole to interpret the underlying relaxation pathways associated with the experimental decay constants reported by D. Koyama and A. J. Orr-Ewing, Phys. Chem. Chem. Phys., 2016, 18, 26224-26235. Energetics, electronic populations and geometries of various stationary points of low-lying electronic states are computed using the semi-classical ab initio surface hopping dynamics simulations. Estimated decay constants of S2-S1 internal conversion (IC) and S1-T2 intersystem crossing (ISC) are in excellent agreement with the experiment. The observed ultrafast ISC is analyzed based on the S1-T2-T1 spin-vibronic coupling mechanism. In contrast to the previous assignment of 6 ps to the T2-T1 IC, our findings enable us to attribute this decay constant to the combined events of T2-T1 IC followed by relaxation of vibrationally hot T1.
Collapse
Affiliation(s)
- Chithra M J
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Maruthamala PO, Vithura, Thiruvananthapuram 695551, Kerala, India.
| | - Probal Nag
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Maruthamala PO, Vithura, Thiruvananthapuram 695551, Kerala, India.
| | - Sivaranjana Reddy Vennapusa
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Maruthamala PO, Vithura, Thiruvananthapuram 695551, Kerala, India.
| |
Collapse
|
12
|
Guan Y, Xie C, Guo H, Yarkony DR. Enabling a Unified Description of Both Internal Conversion and Intersystem Crossing in Formaldehyde: A Global Coupled Quasi-Diabatic Hamiltonian for Its S 0, S 1, and T 1 States. J Chem Theory Comput 2021; 17:4157-4168. [PMID: 34132545 DOI: 10.1021/acs.jctc.1c00370] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
In our recent work, a diabatic Hamiltonian that couples the S0 and S1 states of formaldehyde was constructed using a robust fitting-and-diabatizing procedure with artificial neural networks, which is capable of representing adiabatic energies, energy gradients, and derivative couplings over a wide range of geometries including seams of conical intersection. In this work, based on the diabatization of S0 and S1, the spin-orbit couplings between singlet states (S0, S1) and triplet state T1 are also determined in the same diabatic representation. The diabatized spin-orbit couplings are then fit with a symmetrized neural-network functional form. The ab initio spin-orbit couplings are well reproduced in large configuration space. Together with the neural-network-based potential energy surface for T1, the full quasi-diabatic Hamiltonian for the S0, S1, and T1 states is completed, enabling a unified description of both internal conversion and intersystem crossing in formaldehyde. The vibrational levels on the three adiabatic states are found to be in good agreement with known experimental band origins.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, People's Republic of China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
13
|
Karamatskos ET, Yarlagadda S, Patchkovskii S, Vrakking MJJ, Welsch R, Küpper J, Rouzée A. Time-resolving the UV-initiated photodissociation dynamics of OCS. Faraday Discuss 2021; 228:413-431. [PMID: 33570531 DOI: 10.1039/d0fd00119h] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We present a time-resolved study of the photodissociation dynamics of OCS after UV-photoexcitation at λ = 237 nm. OCS molecules (X1Σ+) were primarily excited to the 11A'' and the 21A' Renner-Teller components of the 1Σ- and 1Δ states. Dissociation into CO and S fragments was observed through time-delayed strong-field ionisation and imaging of the kinetic energy of the resulting CO+ and S+ fragments by intense 790 nm laser pulses. Surprisingly, fast oscillations with a period of ∼100 fs were observed in the S+ channel of the UV dissociation. Based on wavepacket-dynamics simulations coupled with a simple electrostatic-interaction model, these oscillations do not correspond to the known highly-excited rotational motion of the leaving CO(X1Σ+, J ≫ 0) fragments, which has a timescale of ∼140 fs. Instead, we suggest to assign the observed oscillations to the excitation of vibrational wavepackets in the 23A'' or 21A'' states of the molecule that predissociate to form S(3PJ) photoproducts.
Collapse
Affiliation(s)
- Evangelos T Karamatskos
- Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, 22607 Hamburg, Germany. and Department of Physics, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | | | | | | | - Ralph Welsch
- Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, 22607 Hamburg, Germany. and Center for Ultrafast Imaging, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Jochen Küpper
- Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, 22607 Hamburg, Germany. and Department of Physics, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany and Center for Ultrafast Imaging, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Arnaud Rouzée
- Max Born Institute, Max-Born-Straße 2a, 12489 Berlin, Germany.
| |
Collapse
|
14
|
Abstract
In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Institut de Chimie Radicalaire, CNRS 7273, Aix-Marseille University, 13013 Marseille, France;
| | - Dmitry A Fedorov
- Oak Ridge Associated Universities, Oak Ridge, Tennessee 37830, USA;
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, Nevada 89557-0216, USA;
| |
Collapse
|
15
|
Heindl M, González L. Validating fewest-switches surface hopping in the presence of laser fields. J Chem Phys 2021; 154:144102. [PMID: 33858152 DOI: 10.1063/5.0044807] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The capability of fewest-switches surface hopping (FSSH) to describe non-adiabatic dynamics under explicit excitation with external fields is evaluated. Different FSSH parameters are benchmarked against multi-configurational time dependent Hartree (MCTDH) reference calculations using SO2 and 2-thiocytosine as model, yet realistic, molecular systems. Qualitatively, FSSH is able to reproduce the trends in the MCTDH dynamics with (also without) an explicit external field; however, no set of FSSH parameters is ideal. The adequate treatment of the overcoherence in FSSH is revealed as the driving factor to improve the description of the excitation process with respect to the MCTDH reference. Here, two corrections were tested: the augmented-FSSH (AFSSH) correction and the energy-based decoherence correction. A dependence on the employed basis is detected in AFSSH, performing better when spin-orbit and external laser field couplings are treated as off-diagonal elements instead of projecting them onto the diagonal of the Hamilton operator. In the presence of an electric field, the excited state dynamics was found to depend strongly on the vector used to rescale the kinetic energy along after a transition between surfaces. For SO2, recurrence of the excited wave packet throughout the duration of the applied laser pulse is observed for laser pulses (>100 fs), resulting in additional interferences missed by FSSH and only visible in variational multi-configurational Gaussian when utilizing a large number of Gaussian basis functions. This feature vanishes when going toward larger molecules, such as 2-thiocytosine, where this effect is barely visible in a laser pulse 200 fs long.
Collapse
Affiliation(s)
- Moritz Heindl
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| |
Collapse
|
16
|
Moitra T, Karak P, Chakraborty S, Ruud K, Chakrabarti S. Behind the scenes of spin-forbidden decay pathways in transition metal complexes. Phys Chem Chem Phys 2021; 23:59-81. [PMID: 33319894 DOI: 10.1039/d0cp05108j] [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/19/2022]
Abstract
The interpretation of the ultrafast photophysics of transition metal complexes following photo-absorption is quite involved as the heavy metal center leads to a complicated and entangled singlet-triplet manifold. This opens up multiple pathways for deactivation, often with competitive rates. As a result, intersystem crossing (ISC) and phosphorescence are commonly observed in transition metal complexes. A detailed understanding of such an excited-state structure and dynamics calls for state-of-the-art experimental and theoretical methodologies. In this review, we delve into the inability of non-relativistic quantum theory to describe spin-forbidden transitions, which can be overcome by taking into account spin-orbit coupling, whose importance grows with increasing atomic number. We present the quantum chemical theory of phosphorescence and ISC together with illustrative examples. Finally, a few applications are highlighted, bridging the gap between theoretical studies and experimental applications, such as photofunctional materials.
Collapse
Affiliation(s)
- Torsha Moitra
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | | | | | | | | |
Collapse
|
17
|
Menger MFS, Ehrmaier J, Faraji S. PySurf: A Framework for Database Accelerated Direct Dynamics. J Chem Theory Comput 2020; 16:7681-7689. [PMID: 33231447 PMCID: PMC7726901 DOI: 10.1021/acs.jctc.0c00825] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Indexed: 11/28/2022]
Abstract
The greatest restriction to the theoretical study of the dynamics of photoinduced processes is computationally expensive electronic structure calculations. Machine learning algorithms have the potential to reduce the number of these computations significantly. Here, PySurf is introduced as an innovative code framework, which is specifically designed for rapid prototyping and development tasks for data science applications in computational chemistry. It comes with powerful Plugin and Workflow engines, which allows intuitive customization for individual tasks. Data is automatically stored through the database framework, which enables additional interpolation of properties in previously evaluated regions of the conformational space. To illustrate the potential of the framework, a code for nonadiabatic surface hopping simulations based on the Landau-Zener algorithm is presented here. Deriving gradients from the interpolated potential energy surfaces allows for full-dimensional nonadiabatic surface hopping simulations using only adiabatic energies (energy only). Simulations of a pyrazine model and ab initio-based calculations of the SO2 molecule show that energy-only calculations with PySurf are able to correctly predict the nonadiabatic dynamics of these prototype systems. The results reveal the degree of sophistication, which can be achieved by the database accelerated energy-only surface hopping simulations being competitive to commonly used semiclassical approaches.
Collapse
Affiliation(s)
- Maximilian F. S.
J. Menger
- Zernike Institute
for Advanced
Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Johannes Ehrmaier
- Zernike Institute
for Advanced
Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Shirin Faraji
- Zernike Institute
for Advanced
Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| |
Collapse
|
18
|
Richings GW, Habershon S. Direct Grid-Based Nonadiabatic Dynamics on Machine-Learned Potential Energy Surfaces: Application to Spin-Forbidden Processes. J Phys Chem A 2020; 124:9299-9313. [DOI: 10.1021/acs.jpca.0c06125] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
19
|
Talotta F, Morisset S, Rougeau N, Lauvergnat D, Agostini F. Internal Conversion and Intersystem Crossing with the Exact Factorization. J Chem Theory Comput 2020; 16:4833-4848. [PMID: 32633509 DOI: 10.1021/acs.jctc.0c00493] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a detailed derivation of the generalized coupled-trajectory mixed quantum-classical (G-CT-MQC) algorithm based on the exact-factorization equations. The ultimate goal is to propose an algorithm that can be employed for molecular dynamics simulations of nonradiative phenomena, as the spin-allowed internal conversions and the spin-forbidden intersystem crossings. Internal conversions are nonadiabatic processes driven by the kinetic coupling between electronic states, whereas intersystem crossings are mediated by the spin-orbit coupling. In this paper, we discuss computational issues related to the suitable representation for electronic dynamics and the different natures of kinetic and spin-orbit coupling. Numerical studies on model systems allow us to test the performance of the G-CT-MQC algorithm in different situations.
Collapse
Affiliation(s)
- Francesco Talotta
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France.,Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Sabine Morisset
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Nathalie Rougeau
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - David Lauvergnat
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Federica Agostini
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
20
|
Frandsen BN, Farahani S, Vogt E, Lane JR, Kjaergaard HG. Spectroscopy of OSSO and Other Sulfur Compounds Thought to be Present in the Venus Atmosphere. J Phys Chem A 2020; 124:7047-7059. [DOI: 10.1021/acs.jpca.0c04388] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Benjamin N. Frandsen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen Ø, Denmark
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - Sara Farahani
- School of Science, University of Waikato, Private Bag 3105, Hamilton 3240, New Zealand
| | - Emil Vogt
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen Ø, Denmark
| | - Joseph R. Lane
- School of Science, University of Waikato, Private Bag 3105, Hamilton 3240, New Zealand
| | - Henrik G. Kjaergaard
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen Ø, Denmark
| |
Collapse
|
21
|
Lischka H, Shepard R, Müller T, Szalay PG, Pitzer RM, Aquino AJA, Araújo do Nascimento MM, Barbatti M, Belcher LT, Blaudeau JP, Borges I, Brozell SR, Carter EA, Das A, Gidofalvi G, González L, Hase WL, Kedziora G, Kertesz M, Kossoski F, Machado FBC, Matsika S, do Monte SA, Nachtigallová D, Nieman R, Oppel M, Parish CA, Plasser F, Spada RFK, Stahlberg EA, Ventura E, Yarkony DR, Zhang Z. The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry. J Chem Phys 2020; 152:134110. [PMID: 32268762 DOI: 10.1063/1.5144267] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.
Collapse
Affiliation(s)
- Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Ron Shepard
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Thomas Müller
- Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, Jülich 52428, Germany
| | - Péter G Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Budapest, Hungary
| | - Russell M Pitzer
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adelia J A Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | | | | | - Lachlan T Belcher
- Laser and Optics Research Center, Department of Physics, US Air Force Academy, Colorado 80840, USA
| | | | - Itamar Borges
- Departamento de Química, Instituto Militar de Engenharia, Rio de Janeiro, RJ 22290-270, Brazil
| | - Scott R Brozell
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Emily A Carter
- Office of the Chancellor and Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Box 951405, Los Angeles, California 90095-1405, USA
| | - Anita Das
- Indian Institute of Engineering Science and Technology, Shibpur, Howrah, India
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - William L Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Gary Kedziora
- Air Force Research Laboratory, Wright-Patterson Air Force Base, Ohio 45433, USA
| | - Miklos Kertesz
- Department of Chemistry, Georgetown University, 37th and O Streets, NW, Washington, DC 20057-1227, USA
| | | | - Francisco B C Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, 1901 N. 13th St., Philadelphia, Pennsylvania 19122, USA
| | | | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 160610 Prague 6, Czech Republic
| | - Reed Nieman
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Markus Oppel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Carol A Parish
- Department of Chemistry, Gottwald Center for the Sciences, University of Richmond, Richmond, Virginia 23173, USA
| | - Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Rene F K Spada
- Departamento de Física, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | - Eric A Stahlberg
- Biomedical Informatics and Data Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, USA
| | - Elizete Ventura
- Universidade Federal da Paraíba, 58059-900 João Pessoa, PB, Brazil
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, 3400 N. Charles Street, Baltimore, Maryland 21218, USA
| | - Zhiyong Zhang
- Stanford Research Computing Center, Stanford University, 255 Panama Street, Stanford, California 94305, USA
| |
Collapse
|
22
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
23
|
Neville SP, Seidu I, Schuurman MS. Propagative block diagonalization diabatization of DFT/MRCI electronic states. J Chem Phys 2020; 152:114110. [DOI: 10.1063/1.5143126] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Simon P. Neville
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Issaka Seidu
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Michael S. Schuurman
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
- National Research Council of Canada, 100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
24
|
Talotta F, Morisset S, Rougeau N, Lauvergnat D, Agostini F. Spin-Orbit Interactions in Ultrafast Molecular Processes. PHYSICAL REVIEW LETTERS 2020; 124:033001. [PMID: 32031839 DOI: 10.1103/physrevlett.124.033001] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Indexed: 06/10/2023]
Abstract
We investigate spin-orbit interactions in ultrafast molecular processes employing the exact factorization of the electron-nuclear wave function. We revisit the original derivation by including spin-orbit coupling, and show how the dynamics driven by the time-dependent potential energy surface alleviates inconsistencies arising from different electronic representations. We propose a novel trajectory-based scheme to simulate spin-forbidden non-radiative processes, and we show its performance in the treatment of excited-state dynamics where spin-orbit effects couple different spin multiplets.
Collapse
Affiliation(s)
- Francesco Talotta
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Sabine Morisset
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Nathalie Rougeau
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - David Lauvergnat
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Federica Agostini
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
25
|
Plasser F, Mai S, Fumanal M, Gindensperger E, Daniel C, González L. Strong Influence of Decoherence Corrections and Momentum Rescaling in Surface Hopping Dynamics of Transition Metal Complexes. J Chem Theory Comput 2019; 15:5031-5045. [DOI: 10.1021/acs.jctc.9b00525] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Felix Plasser
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K
| | - Sebastian Mai
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| | - Maria Fumanal
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Etienne Gindensperger
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Chantal Daniel
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Leticia González
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| |
Collapse
|
26
|
Mai S, Atkins AJ, Plasser F, González L. The Influence of the Electronic Structure Method on Intersystem Crossing Dynamics. The Case of Thioformaldehyde. J Chem Theory Comput 2019; 15:3470-3480. [DOI: 10.1021/acs.jctc.9b00282] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| | - Andrew J. Atkins
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| |
Collapse
|
27
|
Kumar P, Kłos J, Poirier B, Alexander MH, Guo H. Accurate characterization of the lowest triplet potential energy surface of SO 2 with a coupled cluster method. J Chem Phys 2019; 150:144303. [DOI: 10.1063/1.5088959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Praveen Kumar
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Jacek Kłos
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, USA
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Millard H. Alexander
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, USA
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
28
|
Anglada JM, Martins-Costa MTC, Francisco JS, Ruiz-López MF. Triplet state promoted reaction of SO2 with H2O by competition between proton coupled electron transfer (pcet) and hydrogen atom transfer (hat) processes. Phys Chem Chem Phys 2019; 21:9779-9784. [DOI: 10.1039/c9cp01105f] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The excited triplet electronic state of SO2 (a3B1) reacts with water through a proton coupled electron transfer (pcet) mechanism rather than via a conventional hydrogen atom transfer (hat) process.
Collapse
Affiliation(s)
- Josep M. Anglada
- Departament de Química Biològica (IQAC – CSIC)
- E-08034 Barcelona
- Spain
| | - Marilia T. C. Martins-Costa
- Laboratoire de Physique et Chimie Théoriques
- UMR CNRS 7019
- University of Lorraine
- CNRS
- 54506 Vandoeuvre-lès-Nancy
| | - Joseph S. Francisco
- Department of Earth and Environmental Science and Department of Chemistry
- University of Pennsylvania
- Philadelphia
- USA
| | - Manuel F. Ruiz-López
- Laboratoire de Physique et Chimie Théoriques
- UMR CNRS 7019
- University of Lorraine
- CNRS
- 54506 Vandoeuvre-lès-Nancy
| |
Collapse
|
29
|
Plasser F, Gómez S, Menger MFSJ, Mai S, González L. Highly efficient surface hopping dynamics using a linear vibronic coupling model. Phys Chem Chem Phys 2018; 21:57-69. [PMID: 30306987 DOI: 10.1039/c8cp05662e] [Citation(s) in RCA: 72] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We report an implementation of the linear vibronic coupling (LVC) model within the surface hopping dynamics approach and present utilities for parameterizing this model in a blackbox fashion. This results in an extremely efficient method to obtain qualitative and even semi-quantitative information about the photodynamical behavior of a molecule, and provides a new route toward benchmarking the results of surface hopping computations. The merits and applicability of the method are demonstrated in a number of applications. First, the method is applied to the SO2 molecule showing that it is possible to compute its absorption spectrum beyond the Condon approximation, and that all the main features and timescales of previous on-the-fly dynamics simulations of intersystem crossing are reproduced while reducing the computational effort by three orders of magnitude. The dynamics results are benchmarked against exact wavepacket propagations on the same LVC potentials and against a variation of the electronic structure level. Four additional test cases are presented to exemplify the broader applicability of the model. The photodynamics of the isomeric adenine and 2-aminopurine molecules are studied and it is shown that the LVC model correctly predicts ultrafast decay in the former and an extended excited-state lifetime in the latter. Futhermore, the method correctly predicts ultrafast intersystem crossing in the modified nucleobase 2-thiocytosine and its absence in 5-azacytosine while it fails to describe the ultrafast internal conversion to the ground state in the latter.
Collapse
Affiliation(s)
- Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough, LE11 3TU, UK.
| | | | | | | | | |
Collapse
|
30
|
Mai S, Marquetand P, González L. Nonadiabatic dynamics: The SHARC approach. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2018; 8:e1370. [PMID: 30450129 PMCID: PMC6220962 DOI: 10.1002/wcms.1370] [Citation(s) in RCA: 224] [Impact Index Per Article: 37.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/14/2018] [Revised: 04/05/2018] [Accepted: 04/06/2018] [Indexed: 12/12/2022]
Abstract
We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.
Collapse
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry University of Vienna Vienna Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry University of Vienna Vienna Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry University of Vienna Vienna Austria
| |
Collapse
|
31
|
Mai S, Plasser F, Pabst M, Neese F, Köhn A, González L. Surface hopping dynamics including intersystem crossing using the algebraic diagrammatic construction method. J Chem Phys 2018; 147:184109. [PMID: 29141436 DOI: 10.1063/1.4999687] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
We report an implementation for employing the algebraic diagrammatic construction to second order [ADC(2)] ab initio electronic structure level of theory in nonadiabatic dynamics simulations in the framework of the SHARC (surface hopping including arbitrary couplings) dynamics method. The implementation is intended to enable computationally efficient, reliable, and easy-to-use nonadiabatic dynamics simulations of intersystem crossing in organic molecules. The methodology is evaluated for the 2-thiouracil molecule. It is shown that ADC(2) yields reliable excited-state energies, wave functions, and spin-orbit coupling terms for this molecule. Dynamics simulations are compared to previously reported results using high-level multi-state complete active space perturbation theory, showing favorable agreement.
Collapse
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Mathias Pabst
- Institute of Physical Chemistry, University of Mainz, Duesbergweg 10, D-55099 Mainz, Germany
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Andreas Köhn
- Institute of Physical Chemistry, University of Mainz, Duesbergweg 10, D-55099 Mainz, Germany
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
32
|
Hirshberg B, Gerber RB, Krylov AI. Autocorrelation of electronic wave-functions: a new approach for describing the evolution of electronic structure in the course of dynamics. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1464675] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Barak Hirshberg
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry , Jerusalem, Israel
| | - R. Benny Gerber
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry , Jerusalem, Israel
- Department of Chemistry, University of California , Irvine, CA, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California , Los Angeles, CA, USA
| |
Collapse
|
33
|
Penfold TJ, Gindensperger E, Daniel C, Marian CM. Spin-Vibronic Mechanism for Intersystem Crossing. Chem Rev 2018; 118:6975-7025. [DOI: 10.1021/acs.chemrev.7b00617] [Citation(s) in RCA: 401] [Impact Index Per Article: 66.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Thomas J. Penfold
- Chemistry - School of Natural and Environmental Sciences, Newcastle University, Newcastle Upon-Tyne NE1 7RU, United Kingdom
| | - Etienne Gindensperger
- Laboratoire de Chimie Quantique, Institut de Chimie UMR-7177, CNRS - Université de Strasbourg, 1 Rue Blaise Pascal 67008 Strasbourg, France
| | - Chantal Daniel
- Laboratoire de Chimie Quantique, Institut de Chimie UMR-7177, CNRS - Université de Strasbourg, 1 Rue Blaise Pascal 67008 Strasbourg, France
| | - Christel M. Marian
- Institut für Theoretische Chemie, Heinrich-Heine-Universität Düsseldorf, Universitätsstrasse 1, 40225 Düsseldorf, Germany
| |
Collapse
|
34
|
Bie L, Liu F, Li Y, Dong T, Gao J, Du L, Yuan Q. Spin crossover dynamics studies on the thermally activated molecular oxygen binding mechanism on a model copper complex. Phys Chem Chem Phys 2018; 20:15852-15862. [DOI: 10.1039/c8cp02482k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The theoretical description of the primary dioxygen (O2) binding and activation step in many copper or iron enzymes, suffers from the instrinsically electronic non-adiabaticity of the spin flip events of the triplet dioxygen molecule (3O2), mediated by spin–orbit couplings.
Collapse
Affiliation(s)
- Lihua Bie
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Fang Liu
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Yanwei Li
- Environment Research Institute
- Shandong University
- Jinan
- P. R. China
| | - Tiange Dong
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Likai Du
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Qiaoxia Yuan
- College of Engineering
- Huazhong Agricultural University
- Wuhan 430070
- P. R. China
| |
Collapse
|
35
|
Duan JX, Zhou Y, Xie ZZ, Sun TL, Cao J. Incorporating spin–orbit effects into surface hopping dynamics using the diagonal representation: a linear-response time-dependent density functional theory implementation with applications to 2-thiouracil. Phys Chem Chem Phys 2018; 20:15445-15454. [DOI: 10.1039/c8cp01852a] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Evaluation of SOC values employs Casida's wave functions and the Breit–Pauli spin–orbit Hamiltonian with effective charge approximation.
Collapse
Affiliation(s)
- Jun-Xin Duan
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing
- Wuhan University of Technology
- Wuhan
- P. R. China
| | - Yun Zhou
- Guizhou Provincial Key Laboratory of Computational Nano-material Science
- Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology
- Guizhou Normal College
- Guiyang
- China
| | - Zhi-Zhong Xie
- Department of Chemistry
- School of Chemistry
- Chemical Engineering and Life Sciences
- Wuhan University of Technology
- Wuhan 430070
| | - Tao-Lei Sun
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing
- Wuhan University of Technology
- Wuhan
- P. R. China
| | - Jun Cao
- Guizhou Provincial Key Laboratory of Computational Nano-material Science
- Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology
- Guizhou Normal College
- Guiyang
- China
| |
Collapse
|
36
|
Atkins AJ, González L. Trajectory Surface-Hopping Dynamics Including Intersystem Crossing in [Ru(bpy) 3] 2. J Phys Chem Lett 2017; 8:3840-3845. [PMID: 28766339 DOI: 10.1021/acs.jpclett.7b01479] [Citation(s) in RCA: 90] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Surface-hopping dynamics coupled to linear response TDDFT and explicit nonadiabatic and spin-orbit couplings have been used to model the ultrafast intersystem crossing (ISC) dynamics in [Ru(bpy)3]2+. Simulations using an ensemble of trajectories starting from the singlet metal-to-ligand charge transfer (1MLCT) band show that the manifold of 3MLCT triplet states is first populated from high-lying singlet states within 26 ± 3 fs. ISC competes with an intricate internal conversion relaxation process within the singlet manifold to the lowest singlet state. Normal-mode analysis and principal component analysis, combined with further dynamical simulations where the nuclei are frozen, unequivocally demonstrate that it is not only the high density of states and the large spin-orbit couplings of the system that promote ISC. Instead, geometrical relaxation involving the nitrogen atoms is required to allow for state mixing and efficient triplet population transfer.
Collapse
Affiliation(s)
- Andrew J Atkins
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna , Währinger Straße 17, A-1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna , Währinger Straße 17, A-1090 Vienna, Austria
| |
Collapse
|
37
|
Marquetand P, Nogueira JJ, Mai S, Plasser F, González L. Challenges in Simulating Light-Induced Processes in DNA. Molecules 2016. [PMCID: PMC6155660 DOI: 10.3390/molecules22010049] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.
Collapse
|
38
|
Cao J, Xie ZZ, Yu X. Excited-state dynamics of oxazole: A combined electronic structure calculations and dynamic simulations study. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.05.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
39
|
Zeidabadinejad L, Dehestani M. Effects of Displacement–Distortion of Potential Energy Surfaces on Nonadiabatic Electron Transfers via Conical Intersections: Application to SO 2 and trans-1,3,5-Hexatriene. J Phys Chem A 2016; 120:4431-46. [DOI: 10.1021/acs.jpca.6b01849] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Leila Zeidabadinejad
- Department of Chemistry, Shahid Bahonar University of Kerman, Pajoohesh Sq., Kerman, 76169-14111 Iran,
P.O. Box: 76169-133
| | - Maryam Dehestani
- Department of Chemistry, Shahid Bahonar University of Kerman, Pajoohesh Sq., Kerman, 76169-14111 Iran,
P.O. Box: 76169-133
| |
Collapse
|
40
|
Kłos J, Alexander MH, Kumar P, Poirier B, Jiang B, Guo H. New ab initio adiabatic potential energy surfaces and bound state calculations for the singlet ground X̃1A1 and excited C̃1B2(21A′) states of SO2. J Chem Phys 2016; 144:174301. [DOI: 10.1063/1.4947526] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jacek Kłos
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, USA
| | - Millard H. Alexander
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - Praveen Kumar
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Bin Jiang
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
41
|
Fedorov DA, Pruitt SR, Keipert K, Gordon MS, Varganov SA. Ab Initio Multiple Spawning Method for Intersystem Crossing Dynamics: Spin-Forbidden Transitions between 3B1 and 1A1 States of GeH2. J Phys Chem A 2016; 120:2911-9. [DOI: 10.1021/acs.jpca.6b01406] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Dmitry A. Fedorov
- Department
of Chemistry, University of Nevada, Reno, 1664 North Virginia Street, Reno, Nevada 89557-0216, United States
| | - Spencer R. Pruitt
- Argonne
Leadership Computing Facility, Argonne National Laboratory, 9700 South
Cass Avenue, Lemont, Illinois 60439, United States
| | - Kristopher Keipert
- Ames
Laboratory, Department of Chemistry, Iowa State University, Ames, Iowa 50010, United States
| | - Mark S. Gordon
- Ames
Laboratory, Department of Chemistry, Iowa State University, Ames, Iowa 50010, United States
| | - Sergey A. Varganov
- Department
of Chemistry, University of Nevada, Reno, 1664 North Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
42
|
Plasser F, Ruckenbauer M, Mai S, Oppel M, Marquetand P, González L. Efficient and Flexible Computation of Many-Electron Wave Function Overlaps. J Chem Theory Comput 2016; 12:1207-19. [PMID: 26854874 PMCID: PMC4785508 DOI: 10.1021/acs.jctc.5b01148] [Citation(s) in RCA: 117] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2015] [Indexed: 11/29/2022]
Abstract
A new algorithm for the computation of the overlap between many-electron wave functions is described. This algorithm allows for the extensive use of recurring intermediates and thus provides high computational efficiency. Because of the general formalism employed, overlaps can be computed for varying wave function types, molecular orbitals, basis sets, and molecular geometries. This paves the way for efficiently computing nonadiabatic interaction terms for dynamics simulations. In addition, other application areas can be envisaged, such as the comparison of wave functions constructed at different levels of theory. Aside from explaining the algorithm and evaluating the performance, a detailed analysis of the numerical stability of wave function overlaps is carried out, and strategies for overcoming potential severe pitfalls due to displaced atoms and truncated wave functions are presented.
Collapse
Affiliation(s)
- Felix Plasser
- Institute for Theoretical
Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Matthias Ruckenbauer
- Institute for Theoretical
Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute for Theoretical
Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Markus Oppel
- Institute for Theoretical
Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute for Theoretical
Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Leticia González
- Institute for Theoretical
Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| |
Collapse
|
43
|
Cao J, Xie ZZ. Internal conversion and intersystem crossing in α,β-enones: a combination of electronic structure calculations and dynamics simulations. Phys Chem Chem Phys 2016; 18:6931-45. [DOI: 10.1039/c5cp06833a] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The geometrical constraint of the ring gives rise to a smaller spin–orbital coupling in the singlet–triplet crossing region, resulting in a lower intersystem crossing rate.
Collapse
Affiliation(s)
- Jun Cao
- Guizhou Provincial Key Laboratory of Computational Nano-material Science
- Guizhou Education University
- Guiyang
- China
- Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology
| | - Zhi-Zhong Xie
- Department of Chemistry
- School of Chemistry
- Chemical Engineering and Life Sciences
- Wuhan University of Technology
- Wuhan 430070
| |
Collapse
|
44
|
Franco de Carvalho F, Tavernelli I. Nonadiabatic dynamics with intersystem crossings: A time-dependent density functional theory implementation. J Chem Phys 2015; 143:224105. [DOI: 10.1063/1.4936864] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- F. Franco de Carvalho
- Centre Européen de Calcul Atomique et Moléculaire, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - I. Tavernelli
- IBM Research GmbH, Zurich Research Laboratory, 8803 Ruschlikon, Switzerland
| |
Collapse
|
45
|
Xia SH, Liu XY, Fang Q, Cui G. Photodissociation dynamics of CH3C(O)SH in argon matrix: A QM/MM nonadiabatic dynamics simulation. J Chem Phys 2015; 143:194303. [DOI: 10.1063/1.4935598] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
46
|
Richings GW, Worth GA. A Practical Diabatisation Scheme for Use with the Direct-Dynamics Variational Multi-Configuration Gaussian Method. J Phys Chem A 2015; 119:12457-70. [DOI: 10.1021/acs.jpca.5b07921] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gareth W. Richings
- School
of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, U.K
| | - Graham A. Worth
- School
of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, U.K
| |
Collapse
|
47
|
|
48
|
Ouyang W, Dou W, Subotnik JE. Surface hopping with a manifold of electronic states. I. Incorporating surface-leaking to capture lifetimes. J Chem Phys 2015; 142:084109. [PMID: 25725714 DOI: 10.1063/1.4908032] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the incorporation of the surface-leaking (SL) algorithm into Tully's fewest-switches surface hopping (FSSH) algorithm to simulate some electronic relaxation induced by an electronic bath in conjunction with some electronic transitions between discrete states. The resulting SL-FSSH algorithm is benchmarked against exact quantum scattering calculations for three one-dimensional model problems. The results show excellent agreement between SL-FSSH and exact quantum dynamics in the wide band limit, suggesting the potential for a SL-FSSH algorithm. Discrepancies and failures are investigated in detail to understand the factors that will limit the reliability of SL-FSSH, especially the wide band approximation. Considering the easiness of implementation and the low computational cost, we expect this method to be useful in studying processes involving both a continuum of electronic states (where electronic dynamics are probabilistic) and processes involving only a few electronic states (where non-adiabatic processes cannot ignore short-time coherence).
Collapse
Affiliation(s)
- Wenjun Ouyang
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
49
|
White AJ, Gorshkov VN, Tretiak S, Mozyrsky D. Non-adiabatic molecular dynamics by accelerated semiclassical Monte Carlo. J Chem Phys 2015; 143:014115. [DOI: 10.1063/1.4923473] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexander J. White
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Center for Nonlinear Studies (CNLS), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Center for Nonlinear Studies (CNLS), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Dmitry Mozyrsky
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
50
|
Cui G, Thiel W. Generalized trajectory surface-hopping method for internal conversion and intersystem crossing. J Chem Phys 2015; 141:124101. [PMID: 25273406 DOI: 10.1063/1.4894849] [Citation(s) in RCA: 117] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Trajectory-based fewest-switches surface-hopping (FSSH) dynamics simulations have become a popular and reliable theoretical tool to simulate nonadiabatic photophysical and photochemical processes. Most available FSSH methods model internal conversion. We present a generalized trajectory surface-hopping (GTSH) method for simulating both internal conversion and intersystem crossing processes on an equal footing. We consider hops between adiabatic eigenstates of the non-relativistic electronic Hamiltonian (pure spin states), which is appropriate for sufficiently small spin-orbit coupling. This choice allows us to make maximum use of existing electronic structure programs and to minimize the changes to available implementations of the traditional FSSH method. The GTSH method is formulated within the quantum mechanics (QM)/molecular mechanics framework, but can of course also be applied at the pure QM level. The algorithm implemented in the GTSH code is specified step by step. As an initial GTSH application, we report simulations of the nonadiabatic processes in the lowest four electronic states (S0, S1, T1, and T2) of acrolein both in vacuo and in acetonitrile solution, in which the acrolein molecule is treated at the ab initio complete-active-space self-consistent-field level. These dynamics simulations provide detailed mechanistic insight by identifying and characterizing two nonadiabatic routes to the lowest triplet state, namely, direct S1 → T1 hopping as major pathway and sequential S1 → T2 → T1 hopping as minor pathway, with the T2 state acting as a relay state. They illustrate the potential of the GTSH approach to explore photoinduced processes in complex systems, in which intersystem crossing plays an important role.
Collapse
Affiliation(s)
- Ganglong Cui
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|