1
|
Barbatti M, Bondanza M, Crespo-Otero R, Demoulin B, Dral PO, Granucci G, Kossoski F, Lischka H, Mennucci B, Mukherjee S, Pederzoli M, Persico M, Pinheiro M, Pittner J, Plasser F, Sangiogo Gil E, Stojanovic L. Newton-X Platform: New Software Developments for Surface Hopping and Nuclear Ensembles. J Chem Theory Comput 2022; 18:6851-6865. [PMID: 36194696 PMCID: PMC9648185 DOI: 10.1021/acs.jctc.2c00804] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
![]()
Newton-X is an open-source computational platform to
perform nonadiabatic
molecular dynamics based on surface hopping and spectrum simulations
using the nuclear ensemble approach. Both are among the most common
methodologies in computational chemistry for photophysical and photochemical
investigations. This paper describes the main features of these methods
and how they are implemented in Newton-X. It emphasizes the newest
developments, including zero-point-energy leakage correction, dynamics
on complex-valued potential energy surfaces, dynamics induced by incoherent
light, dynamics based on machine-learning potentials, exciton dynamics
of multiple chromophores, and supervised and unsupervised machine
learning techniques. Newton-X is interfaced with several third-party
quantum-chemistry programs, spanning a broad spectrum of electronic
structure methods.
Collapse
Affiliation(s)
- Mario Barbatti
- Aix Marseille University, CNRS, ICR, 13013Marseille, France.,Institut Universitaire de France, 75231Paris, France
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Rachel Crespo-Otero
- Department of Chemistry, Queen Mary University of London, Mile End Road, E1 4NSLondon, U.K
| | | | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, 361005Xiamen, China
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000Toulouse, France
| | - Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas79409, United States
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | | | - Marek Pederzoli
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223Prague 8, Czech Republic
| | - Maurizio Persico
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Max Pinheiro
- Aix Marseille University, CNRS, ICR, 13013Marseille, France
| | - Jiří Pittner
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223Prague 8, Czech Republic
| | - Felix Plasser
- Department of Chemistry, Loughborough University, LE11 3TULoughborough, U.K
| | - Eduarda Sangiogo Gil
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Ljiljana Stojanovic
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BTLondon, U.K
| |
Collapse
|
2
|
Ibele LM, Curchod BFE, Agostini F. A Photochemical Reaction in Different Theoretical Representations. J Phys Chem A 2022; 126:1263-1281. [PMID: 35157450 PMCID: PMC8883471 DOI: 10.1021/acs.jpca.1c09604] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
The Born–Oppenheimer
picture has forged our representation
and interpretation of photochemical processes, from photoexcitation
down to the passage through a conical intersection, a funnel connecting
different electronic states. In this work, we analyze a full in silico
photochemical experiment, from the explicit electronic excitation
by a laser pulse to the formation of photoproducts following a nonradiative
decay through a conical intersection, by contrasting the picture offered
by Born–Oppenheimer and that proposed by the exact factorization.
The exact factorization offers an alternative understanding of photochemistry
that does not rely on concepts such as electronic states, nonadiabatic
couplings, and conical intersections. On the basis of nonadiabatic
quantum dynamics performed for a two-state 2D model system, this work
allows us to compare Born–Oppenheimer and exact factorization
for (i) an explicit photoexcitation with and without the Condon approximation,
(ii) the passage of a nuclear wavepacket through a conical intersection,
(iii) the formation of excited stationary states in the Franck–Condon
region, and (iv) the use of classical and quantum trajectories in
the exact factorization picture to capture nonadiabatic processes
triggered by a laser pulse.
Collapse
Affiliation(s)
- Lea M Ibele
- Department of Chemistry, Durham University, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Department of Chemistry, Durham University, Durham DH1 3LE, United Kingdom
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
3
|
Mandal A, Hunt KLC. Quantum transition probabilities due to overlapping electromagnetic pulses: Persistent differences between Dirac's form and nonadiabatic perturbation theory. J Chem Phys 2021; 154:024116. [PMID: 33445917 DOI: 10.1063/5.0020169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The probability of transition to an excited state of a quantum system in a time-dependent electromagnetic field determines the energy uptake from the field. The standard expression for the transition probability has been given by Dirac. Landau and Lifshitz suggested, instead, that the adiabatic effects of a perturbation should be excluded from the transition probability, leaving an expression in terms of the nonadiabatic response. In our previous work, we have found that these two approaches yield different results while a perturbing field is acting on the system. Here, we prove, for the first time, that differences between the two approaches may persist after the perturbing fields have been completely turned off. We have designed a pair of overlapping pulses in order to establish the possibility of lasting differences, in a case with dephasing. Our work goes beyond the analysis presented by Landau and Lifshitz, since they considered only linear response and required that a constant perturbation must remain as t → ∞. First, a "plateau" pulse populates an excited rotational state and produces coherences between the ground and excited states. Then, an infrared pulse acts while the electric field of the first pulse is constant, but after dephasing has occurred. The nonadiabatic perturbation theory permits dephasing, but dephasing of the perturbed part of the wave function cannot occur within Dirac's method. When the frequencies in both pulses are on resonance, the lasting differences in the calculated transition probabilities may exceed 35%. The predicted differences are larger for off-resonant perturbations.
Collapse
Affiliation(s)
- Anirban Mandal
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Katharine L C Hunt
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
4
|
Farfan CA, Turner DB. A systematic model study quantifying how conical intersection topography modulates photochemical reactions. Phys Chem Chem Phys 2020; 22:20265-20283. [PMID: 32966428 DOI: 10.1039/d0cp03464a] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.
Collapse
Affiliation(s)
- Camille A Farfan
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Daniel B Turner
- Department of Chemistry, New York University, New York, NY 10003, USA
| |
Collapse
|
5
|
Yaghoubi Jouybari M, Liu Y, Improta R, Santoro F. Ultrafast Dynamics of the Two Lowest Bright Excited States of Cytosine and 1-Methylcytosine: A Quantum Dynamical Study. J Chem Theory Comput 2020; 16:5792-5808. [PMID: 32687360 DOI: 10.1021/acs.jctc.0c00455] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The nonadiabatic quantum dynamics (QD) of cytosine and 1-methylcytosine in the gas phase is simulated for 250 fs after a photoexcitation to one of the first two bright states. The nuclear wavepacket is propagated on the coupled diabatic potential energy surfaces of the lowest seven excited states, including ππ*, nπ*, and Rydberg states along all the vibrational degrees of freedom. We focus in particular on the interplay between the bright and the dark nπ* states, not considering the decay to the ground electronic state. To run these simulations, we implemented an automatic general procedure to parametrize linear vibronic coupling (LVC) models with time-dependent density functional theory (DFT) computations and interfaced it with Gaussian package. The wavepacket was propagated with the multilayer version of the multiconfigurational time dependent Hartree method. Two different density functionals, PBE0 and CAM-B3LYP, which provide a different description of the relative stability of the lowest energy dark states, were used to parametrize the LVC Hamiltonian. Part of the photoexcited population on lowest HOMO-LUMO transition (πHπL*) decays within less than 100 fs to a nπ* state which mainly involves a promotion of an electron from the oxygen lone pair to the LUMO (nOπL*). The population of the second ππ* state decays almost completely, in <100 fs, not only to πHπL* and to nOπL* states but also to another nπL* state involving the nitrogen lone pair. The efficiency of the adopted protocol allowed us to check the accuracy of the predictions by repeating the QD simulations with different LVC Hamiltonians parametrized either at the ground-state minimum or at stationary structures of different relevant excited states.
Collapse
Affiliation(s)
- Martha Yaghoubi Jouybari
- 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
| | - Yanli Liu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Roberto Improta
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - 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
| |
Collapse
|
6
|
Photoinduced Double Proton Transfer in the Glyoxal-Methanol Complex Revisited: The Role of the Excited States. J Chem Theory Comput 2020; 16:3273-3286. [PMID: 32275423 DOI: 10.1021/acs.jctc.0c00007] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Under irradiation in the visible range, the glyoxal-methanol complex in a cryogenic argon matrix undergoes a double proton transfer (DPT) reaction through which the glyoxal molecule isomerizes into hydroxyketene. In this work, we employ electronic structure simulations in order to shed more light on the underlying mechanism. Rewardingly, we find that the lowest singlet excited state (S1) of the complex acts as a gateway to two previously unknown isomerization pathways, of which one takes place entirely in the singlet manifold and the other also involves the lowest triplet state (T1). Both of these pathways are fully compatible with the available experimental data, implying that either or both are operative under experimental conditions. In either pathway, the methanol molecule acts as a proton shuttle between the proton-donating and proton-accepting sites of glyoxal, resulting in a dramatic lowering of the potential energy barrier to isomerization with respect to the case of isolated glyoxal. The occurrence of DPT in the singlet manifold is demonstrated directly with the use of nonadiabatic molecular dynamics simulations at the spin-flip time-dependent density functional theory level.
Collapse
|
7
|
Titov E, Humeniuk A, Mitrić R. Comparison of moving and fixed basis sets for nonadiabatic quantum dynamics at conical intersections. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2019.110526] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
8
|
Tsuru S, Vidal ML, Pápai M, Krylov AI, Møller KB, Coriani S. Time-resolved near-edge X-ray absorption fine structure of pyrazine from electronic structure and nuclear wave packet dynamics simulations. J Chem Phys 2019; 151:124114. [PMID: 31575192 DOI: 10.1063/1.5115154] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
As a demonstration of the analysis of the electronic structure and the nuclear dynamics from time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS), we present the TR-NEXAFS spectra of pyrazine following the excitation to the 1B2u(ππ*) state. The spectra are calculated combining the frozen-core/core-valence separated equation-of-motion coupled cluster singles and doubles approach for the spectral signatures and the multiconfiguration time-dependent Hartree method for the wave packet propagation. The population decay from the 1B2u(ππ*) state to the 1B3u(nπ*) and 1Au(nπ*) states, followed by oscillatory flow of population between the 1B3u(nπ*) and 1Au(nπ*) states, is interpreted by means of visualization of the potential energy curves and the reduced nuclear densities. By examining the electronic structure of the three valence-excited states and the final core-excited states, we observe that the population dynamics is explicitly reflected in the TR-NEXAFS spectra, especially when the heteroatoms are selected as the X-ray absorption sites. This work illustrates the feasibility of extracting fine details of molecular photophysical processes from TR-NEXAFS spectra by using currently available theoretical methods.
Collapse
Affiliation(s)
- Shota Tsuru
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Marta L Vidal
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Mátyás Pápai
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Klaus B Møller
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Sonia Coriani
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
9
|
Peng J, Xie Y, Hu D, Lan Z. Performance of trajectory surface hopping method in the treatment of ultrafast intersystem crossing dynamics. J Chem Phys 2019; 150:164126. [PMID: 31042919 DOI: 10.1063/1.5079426] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We carried out extensive studies to examine the performance of the fewest-switches surface hopping method in the description of the ultrafast intersystem crossing dynamic of various singlet-triplet (S-T) models by comparison with the results of the exact full quantum dynamics. Different implementation details and some derivative approaches were examined. As expected, it is better to perform the trajectory surface hopping calculations in the spin-adiabatic representation or by the local diabatization approach, instead of in the spin-diabatic representation. The surface hopping method provides reasonable results for the short-time dynamics in the S-T model with weak spin-orbital coupling (diabatic coupling), although it does not perform well in the models with strong spin-orbital coupling (diabatic coupling). When the system accesses the S-T potential energy crossing with rather high kinetic energy, the trajectory surface hopping method tends to produce a good description of the nonadiabatic intersystem crossing dynamics. The impact of the decoherence correction on the performance of the trajectory surface hopping is system dependent. It improves the result accuracy in many cases, while its influence may also be minor for other cases.
Collapse
Affiliation(s)
- Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
10
|
Xie W, Sapunar M, Došlić N, Sala M, Domcke W. Assessing the performance of trajectory surface hopping methods: Ultrafast internal conversion in pyrazine. J Chem Phys 2019; 150:154119. [DOI: 10.1063/1.5084961] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Weiwei Xie
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Marin Sapunar
- Department of Physical Chemistry, Ruđer Bošković Institute, HR-10000 Zagreb, Croatia
| | - Nađa Došlić
- Department of Physical Chemistry, Ruđer Bošković Institute, HR-10000 Zagreb, Croatia
| | - Matthieu Sala
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS, Université de Bourgogne, BP 47870, F-21078 Dijon, France and Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstr. 40, D-24098 Kiel, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
| |
Collapse
|
11
|
Wibowo M, Persico M, Granucci G. Nonadiabatic dynamics simulations of singlet fission in 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene crystals. Phys Chem Chem Phys 2019; 21:692-701. [PMID: 30543235 DOI: 10.1039/c8cp05474f] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present simulations of the singlet fission dynamics in 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene (ThBF), a thienoquinoid compound recently investigated experimentally by Kawata et al. The simulation model consisted of two ThBF molecules embedded in their crystal environment. The aim was to understand the singlet fission mechanism, and to predict the excited state lifetimes and the singlet fission quantum yield, hitherto unknown. The simulations were performed by the trajectory surface hopping approach with on-the-fly calculations of the electronic wave functions and energies by the semiempirical FOMO-CI method. We found that the initially photogenerated excitonic bright state decays to the lower dark state with a biexponential behaviour, essentially due to transitions to other close lying states. The dark state in turn decays with a lifetime of about 1 ps to the double triplet 1TT state, which is long-lived, as ascertained by performing a simulation with inclusion of the spin-orbit coupling. The singlet fission quantum yield is predicted to be close to the theoretical maximum of 200%. In view of using this thienoquinoid compound in photovoltaic devices, a major drawback is the low energy of the T1 state at its equilibrium geometry.
Collapse
Affiliation(s)
- Meilani Wibowo
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, v. G. Moruzzi 13, I-56124 Pisa, Italy.
| | | | | |
Collapse
|
12
|
Crespo-Otero R, Barbatti M. Recent Advances and Perspectives on Nonadiabatic Mixed Quantum–Classical Dynamics. Chem Rev 2018; 118:7026-7068. [DOI: 10.1021/acs.chemrev.7b00577] [Citation(s) in RCA: 301] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | | |
Collapse
|
13
|
Granucci G, Melani G, Persico M, Van Leuven P. Energy Selection in Nonadiabatic Transitions. J Phys Chem A 2018; 122:678-689. [DOI: 10.1021/acs.jpca.7b10146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Giovanni Granucci
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi
13, I-56124 Pisa, Italy
| | - Giacomo Melani
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi
13, I-56124 Pisa, Italy
| | - Maurizio Persico
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi
13, I-56124 Pisa, Italy
| | - Piet Van Leuven
- Departement
Natuurkunde, Vrije Universiteit Brussel, Pleinlaan 2, B-1050 Brussel, Belgium
| |
Collapse
|
14
|
Xie W, Domcke W. Accuracy of trajectory surface-hopping methods: Test for a two-dimensional model of the photodissociation of phenol. J Chem Phys 2017; 147:184114. [DOI: 10.1063/1.5006788] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Weiwei Xie
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
| |
Collapse
|
15
|
Heaps CW, Mazziotti DA. Accurate non-adiabatic quantum dynamics from pseudospectral sampling of time-dependent Gaussian basis sets. J Chem Phys 2016. [DOI: 10.1063/1.4959872] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Charles W. Heaps
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - David A. Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
16
|
Humeniuk A, Mitrić R. Non-adiabatic dynamics around a conical intersection with surface-hopping coupled coherent states. J Chem Phys 2016; 144:234108. [DOI: 10.1063/1.4954189] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Alexander Humeniuk
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany
| |
Collapse
|
17
|
Kinzel D, Zilberg S, González L. Origin of the Regioselectivity in the Gas-Phase Aniline+CH 3+Electrophilic Aromatic Substitution. Chemphyschem 2015; 16:2366-74. [DOI: 10.1002/cphc.201500256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2015] [Indexed: 11/09/2022]
|
18
|
Kapral R. Quantum dynamics in open quantum-classical systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:073201. [PMID: 25634784 DOI: 10.1088/0953-8984/27/7/073201] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Often quantum systems are not isolated and interactions with their environments must be taken into account. In such open quantum systems these environmental interactions can lead to decoherence and dissipation, which have a marked influence on the properties of the quantum system. In many instances the environment is well-approximated by classical mechanics, so that one is led to consider the dynamics of open quantum-classical systems. Since a full quantum dynamical description of large many-body systems is not currently feasible, mixed quantum-classical methods can provide accurate and computationally tractable ways to follow the dynamics of both the system and its environment. This review focuses on quantum-classical Liouville dynamics, one of several quantum-classical descriptions, and discusses the problems that arise when one attempts to combine quantum and classical mechanics, coherence and decoherence in quantum-classical systems, nonadiabatic dynamics, surface-hopping and mean-field theories and their relation to quantum-classical Liouville dynamics, as well as methods for simulating the dynamics.
Collapse
Affiliation(s)
- Raymond Kapral
- Department of Chemistry, Chemical Physics Theory Group, University of Toronto, Toronto, ON, M5S 3H6 Canada
| |
Collapse
|
19
|
Sissa C, Delchiaro F, Di Maiolo F, Terenziani F, Painelli A. Vibrational coherences in charge-transfer dyes: A non-adiabatic picture. J Chem Phys 2014; 141:164317. [DOI: 10.1063/1.4898710] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Cristina Sissa
- Dipartimento di Chimica, Università degli studi di Parma, Parco Area delle Scienze 17/A, 43124 Parma, Italy
| | - Francesca Delchiaro
- Dipartimento di Chimica, Università degli studi di Parma, Parco Area delle Scienze 17/A, 43124 Parma, Italy
| | - Francesco Di Maiolo
- Dipartimento di Chimica, Università degli studi di Parma, Parco Area delle Scienze 17/A, 43124 Parma, Italy
| | - Francesca Terenziani
- Dipartimento di Chimica, Università degli studi di Parma, Parco Area delle Scienze 17/A, 43124 Parma, Italy
- National Interuniversity Consortium of Materials Science and Technology, INSTM-UdR Parma, Italy
| | - Anna Painelli
- Dipartimento di Chimica, Università degli studi di Parma, Parco Area delle Scienze 17/A, 43124 Parma, Italy
- National Interuniversity Consortium of Materials Science and Technology, INSTM-UdR Parma, Italy
| |
Collapse
|
20
|
Ryabinkin IG, Hsieh CY, Kapral R, Izmaylov AF. Analysis of geometric phase effects in the quantum-classical Liouville formalism. J Chem Phys 2014; 140:084104. [PMID: 24588145 DOI: 10.1063/1.4866366] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. As a consequence, GP effects in WA-QCL can be associated with a dynamical term in the corresponding equation of motion. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.
Collapse
Affiliation(s)
- Ilya G Ryabinkin
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| | - Chang-Yu Hsieh
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Raymond Kapral
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| |
Collapse
|
21
|
An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1526-1] [Citation(s) in RCA: 102] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
22
|
Semiclassical Monte-Carlo approach for modelling non-adiabatic dynamics in extended molecules. Nat Commun 2014; 4:2144. [PMID: 23864100 PMCID: PMC3759042 DOI: 10.1038/ncomms3144] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2013] [Accepted: 06/14/2013] [Indexed: 12/24/2022] Open
Abstract
Modelling of non-adiabatic dynamics in extended molecular systems and solids is a next frontier of atomistic electronic structure theory. The underlying numerical algorithms should operate only with a few quantities (that can be efficiently obtained from quantum chemistry), provide a controlled approximation (which can be systematically improved) and capture important phenomena such as branching (multiple products), detailed balance and evolution of electronic coherences. Here we propose a new algorithm based on Monte-Carlo sampling of classical trajectories, which satisfies the above requirements and provides a general framework for existing surface hopping methods for non-adiabatic dynamics simulations. In particular, our algorithm can be viewed as a post-processing technique for analysing numerical results obtained from the conventional surface hopping approaches. Presented numerical tests for several model problems demonstrate efficiency and accuracy of the new method. Many interesting chemical problems like photosynthesis and photovoltaics involve non-adiabatic dynamical phenomena, which are difficult to predict theoretically. Here, the authors develop a new numerical method capable of recovering quantum interferences that are neglected by conventional methods.
Collapse
|
23
|
Joubert-Doriol L, Ryabinkin IG, Izmaylov AF. Geometric phase effects in low-energy dynamics near conical intersections: A study of the multidimensional linear vibronic coupling model. J Chem Phys 2013; 139:234103. [DOI: 10.1063/1.4844095] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Yonehara T, Takatsuka K. Path-Branching Representation for Nonadiabatic Electron Dynamics in Conical Intersection. J Phys Chem A 2013; 117:8599-608. [DOI: 10.1021/jp402655q] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Takehiro Yonehara
- Department of Basic
Sciences, Graduate School of Arts and Sciences, The University of Tokyo, Komaba, 153-8902 Tokyo, Japan
| | - Kazuo Takatsuka
- Department of Basic
Sciences, Graduate School of Arts and Sciences, The University of Tokyo, Komaba, 153-8902 Tokyo, Japan
| |
Collapse
|
25
|
Hsieh CY, Kapral R. Analysis of the forward-backward trajectory solution for the mixed quantum-classical Liouville equation. J Chem Phys 2013; 138:134110. [DOI: 10.1063/1.4798221] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Borrelli R, Peluso A. Elementary electron transfer reactions: from basic concepts to recent computational advances. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1147] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
27
|
Kelly A, van Zon R, Schofield J, Kapral R. Mapping quantum-classical Liouville equation: Projectors and trajectories. J Chem Phys 2012; 136:084101. [DOI: 10.1063/1.3685420] [Citation(s) in RCA: 82] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Yonehara T, Hanasaki K, Takatsuka K. Fundamental Approaches to Nonadiabaticity: Toward a Chemical Theory beyond the Born–Oppenheimer Paradigm. Chem Rev 2011; 112:499-542. [DOI: 10.1021/cr200096s] [Citation(s) in RCA: 154] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Takehiro Yonehara
- Department of Basic Sciences, Graduate School of Arts and Sciences, The University of Tokyo Komaba, 153-8902 Tokyo, Japan
| | - Kota Hanasaki
- Department of Basic Sciences, Graduate School of Arts and Sciences, The University of Tokyo Komaba, 153-8902 Tokyo, Japan
| | - Kazuo Takatsuka
- Department of Basic Sciences, Graduate School of Arts and Sciences, The University of Tokyo Komaba, 153-8902 Tokyo, Japan
| |
Collapse
|
29
|
Kettunen JA, Niskanen J, Huttula M, Vapa M, Urpelainen S, Aksela H. Electron-ion coincidence study of photofragmentation of the CdCl(2) molecule. JOURNAL OF MASS SPECTROMETRY : JMS 2011; 46:901-907. [PMID: 21915954 DOI: 10.1002/jms.1967] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
In this work, the photofragmentation subsequent to valence and Cd4d photoionization of cadmium dichloride (CdCl(2)) were studied using He I and synchrotron excitation. The measurements were performed with a photoelectron-photoion coincidence (PEPICO) setup, and the connection between the singly ionized electronic states and cationic fragments was investigated. The valence-ionized states were found to lead to CdCl(2)(+), Cd(+) and CdCl(+). The Cd4d(- 1) states were found to lead only to Cl(+) ions. The observed charge transfer effect between Cd and Cl was concluded to take place due to internal conversion or fluorescence decay to dissociating valence states either directly or through consecutive fragmentation. The fragmentation energetics were investigated with molecular ab initio calculations, and the calculated energies were found to agree with the detected fragment appearances.
Collapse
Affiliation(s)
- J A Kettunen
- Department of Physics, University of Oulu, P.O. Box 3000, FIN-90014, University of Oulu, Finland.
| | | | | | | | | | | |
Collapse
|
30
|
Park K, Engelkemier J, Persico M, Manikandan P, Hase WL. Algorithms for Sampling a Quantum Microcanonical Ensemble of Harmonic Oscillators at Potential Minima and Conical Intersections. J Phys Chem A 2011; 115:6603-9. [DOI: 10.1021/jp110799m] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Kyoyeon Park
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| | - Joshua Engelkemier
- Department of Chemistry, Iowa State University, Ames, Iowa 50011-3111, United States
| | - Maurizio Persico
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| | - Paranjothy Manikandan
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| | - William L. Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| |
Collapse
|
31
|
Nonadiabatic dynamics with trajectory surface hopping method. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.64] [Citation(s) in RCA: 279] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
32
|
Picconi D, Barone V, Lami A, Santoro F, Improta R. The Interplay between ππ*/nπ* Excited States in Gas-Phase Thymine: A Quantum Dynamical Study. Chemphyschem 2011; 12:1957-68. [DOI: 10.1002/cphc.201001080] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2010] [Revised: 03/01/2011] [Indexed: 11/10/2022]
|
33
|
Vazdar M, Eckert-Maksić M, Barbatti M, Lischka H. Excited-state non-adiabatic dynamics simulations of pyrrole. Mol Phys 2010. [DOI: 10.1080/00268970802665639] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
34
|
Kelly A, Kapral R. Quantum-classical description of environmental effects on electronic dynamics at conical intersections. J Chem Phys 2010; 133:084502. [DOI: 10.1063/1.3475773] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Improta R, Barone V, Lami A, Santoro F. Quantum Dynamics of the Ultrafast ππ*/nπ* Population Transfer in Uracil and 5-Fluoro-Uracil in Water and Acetonitrile. J Phys Chem B 2009; 113:14491-503. [DOI: 10.1021/jp906524p] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Roberto Improta
- Dipartimento di Chimica and INSTM-Village, Universitá Federico II, Complesso Monte S. Angelo, via Cintia, I-80126 Napoli, Italy, Scuola Normale Superiore di Pisa, P.zza dei Cavalieri 7, I-56126 Pisa, Italy, and Istituto per i Processi Chimico-Fisici - CNR, Area della Ricerca del CNR Via Moruzzi,1 I-56124 Pisa, Italy
| | - Vincenzo Barone
- Dipartimento di Chimica and INSTM-Village, Universitá Federico II, Complesso Monte S. Angelo, via Cintia, I-80126 Napoli, Italy, Scuola Normale Superiore di Pisa, P.zza dei Cavalieri 7, I-56126 Pisa, Italy, and Istituto per i Processi Chimico-Fisici - CNR, Area della Ricerca del CNR Via Moruzzi,1 I-56124 Pisa, Italy
| | - Alessandro Lami
- Dipartimento di Chimica and INSTM-Village, Universitá Federico II, Complesso Monte S. Angelo, via Cintia, I-80126 Napoli, Italy, Scuola Normale Superiore di Pisa, P.zza dei Cavalieri 7, I-56126 Pisa, Italy, and Istituto per i Processi Chimico-Fisici - CNR, Area della Ricerca del CNR Via Moruzzi,1 I-56124 Pisa, Italy
| | - Fabrizio Santoro
- Dipartimento di Chimica and INSTM-Village, Universitá Federico II, Complesso Monte S. Angelo, via Cintia, I-80126 Napoli, Italy, Scuola Normale Superiore di Pisa, P.zza dei Cavalieri 7, I-56126 Pisa, Italy, and Istituto per i Processi Chimico-Fisici - CNR, Area della Ricerca del CNR Via Moruzzi,1 I-56124 Pisa, Italy
| |
Collapse
|
36
|
Yang S, Coe JD, Kaduk B, Martínez TJ. An “optimal” spawning algorithm for adaptive basis set expansion in nonadiabatic dynamics. J Chem Phys 2009; 130:134113. [DOI: 10.1063/1.3103930] [Citation(s) in RCA: 80] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
37
|
Santoro F, Improta R, Barone V. Three-dimensional diabatic models for the ππ* → nπ* excited-state decay of uracil derivatives in solution. Theor Chem Acc 2009. [DOI: 10.1007/s00214-009-0527-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
38
|
Villani G. A time-dependent quantum dynamics investigation of the guanine-cytosine system: A six-dimensional model. J Chem Phys 2008; 128:114306. [DOI: 10.1063/1.2890040] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
|
39
|
Barbatti M, Granucci G, Persico M, Ruckenbauer M, Vazdar M, Eckert-Maksić M, Lischka H. The on-the-fly surface-hopping program system Newton-X: Application to ab initio simulation of the nonadiabatic photodynamics of benchmark systems. J Photochem Photobiol A Chem 2007. [DOI: 10.1016/j.jphotochem.2006.12.008] [Citation(s) in RCA: 380] [Impact Index Per Article: 22.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
40
|
Villani G. Theoretical investigation of the coupling between the hydrogen transfer and the base pair opening in the adenine–thymine system. Chem Phys 2007. [DOI: 10.1016/j.chemphys.2007.06.005] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
41
|
Granucci G, Persico M. Critical appraisal of the fewest switches algorithm for surface hopping. J Chem Phys 2007; 126:134114. [PMID: 17430023 DOI: 10.1063/1.2715585] [Citation(s) in RCA: 470] [Impact Index Per Article: 27.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper the authors address the problem of internal consistency in trajectory surface hopping methods, i.e., the requirement that the fraction of trajectories running on each electronic state equals the probabilities computed by the electronic time-dependent Schrodinger equation, after averaging over all trajectories. They derive a formula for the hopping probability in Tully's "fewest switches" spirit that would yield a rigorously consistent treatment. They show the relationship of Tully's widely used surface hopping algorithm with the "exact" prescription that cannot be applied when running each trajectory independently. They also bring out the connection of the consistency problem with the coherent propagation of the electronic wave function and the artifacts caused by coherent Rabi-type oscillations of the state probabilities in weak coupling regimes. A real molecule (azobenzene) and two ad hoc models serve as examples to illustrate the above theoretical arguments. Following a proposal by Truhlar's group [Zhu et al., J. Chem. Phys. 121, 7658 (2004) Zhu et al., J. Chem. Theory Comput. 1, 527 (2005)], they apply a decoherence correction to the state probabilities, in conjunction with Tully's algorithm, and they obtain satisfactory results in terms of internal consistency and of agreement with the outcomes of quantum wave packet calculations.
Collapse
Affiliation(s)
- Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento 35, 56126 Pisa, Italy
| | | |
Collapse
|
42
|
Barbatti M, Lischka H. Can the Nonadiabatic Photodynamics of Aminopyrimidine Be a Model for the Ultrafast Deactivation of Adenine? J Phys Chem A 2007; 111:2852-8. [PMID: 17388405 DOI: 10.1021/jp070089w] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The reaction path for the ultrafast deactivation of 6-aminopyrimidine (6AP) has been investigated by means of ab initio surface-hopping dynamics. The electronic vertical excitation spectrum, excited-state S1 minima, and minima on the crossing seam of 6AP resemble well those found for adenine. The deactivation from the S1 to the S0 state takes place at the ultrafast time scale of 400 fs and involves the out-of-plane ring deformation of the C4 atom, a position that is sterically restricted in adenine by the imidazole ring. Mechanical restrictions have been used to simulate in a simple way the role of the imidazole group in adenine. As a result, deactivation via out-of-plane ring deformation of the C2 and N3 atoms are observed in good agreement with predictions for adenine. These dynamics results show that the previously suggested ring puckering deactivation paths really exist at a time scale, which is compatible with experimentally observed life times. The electronic structure of the crossing seam has been shown to have the same nature as those of simple biradicaloid systems, a feature which seems to be valid for any cyclic system.
Collapse
Affiliation(s)
- Mario Barbatti
- Institute for Theoretical Chemistry, University of Vienna, Waehringerstrasse 17, A-1090, Vienna, Austria.
| | | |
Collapse
|
43
|
Domcke W, Stock G. Theory of Ultrafast Nonadiabatic Excited-State Processes and their Spectroscopic Detection in Real Time. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470141595.ch1] [Citation(s) in RCA: 250] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
44
|
Abstract
An approach to the quantum-classical mechanics of phase space dependent operators, which has been proposed recently, is remodeled as a formalism for wave fields. Such wave fields obey a system of coupled nonlinear equations that can be written by means of a suitable non-Hamiltonian bracket. As an example, the theory is applied to the relaxation dynamics of the spin-boson model. In the adiabatic limit, a good agreement with calculations performed by the operator approach is obtained. Moreover, the theory proposed in this paper can take nonadiabatic effects into account without resorting to surface-hopping approximations. Hence, the results obtained follow qualitatively those of previous surface-hopping calculations and increase by a factor of (at least) 2, the time length over which nonadiabatic dynamics can be propagated with small statistical errors. Moreover, it is worth to note that the dynamics of quantum-classical wave fields proposed here is a straightforward non-Hamiltonian generalization of the formalism for nonlinear quantum mechanics that Weinberg introduced recently.
Collapse
Affiliation(s)
- Alessandro Sergi
- Dipartimento di Fisica, Universitá degli Studi di Messina, Contrada Papardo 98166 Messina, Italy.
| |
Collapse
|
45
|
Fernandez-Ramos A, Miller JA, Klippenstein SJ, Truhlar DG. Modeling the kinetics of bimolecular reactions. Chem Rev 2007; 106:4518-84. [PMID: 17091928 DOI: 10.1021/cr050205w] [Citation(s) in RCA: 393] [Impact Index Per Article: 23.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Antonio Fernandez-Ramos
- Departamento de Quimica Fisica, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | | | | | | |
Collapse
|
46
|
Complex excited dynamics around a plateau on a retinal-like potential surface: chaos, multi-exponential decays and quantum/classical differences. Theor Chem Acc 2007. [DOI: 10.1007/s00214-006-0220-3] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
47
|
|
48
|
|
49
|
Burghardt I, Hynes JT. Excited-State Charge Transfer at a Conical Intersection: Effects of an Environment. J Phys Chem A 2006; 110:11411-23. [PMID: 17020251 DOI: 10.1021/jp057569c] [Citation(s) in RCA: 78] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The influence of a polar and polarizable environment on charge transfer processes at a conical intersection (CI) can be described by a diabatic free energy model yielding coupled surfaces as a function of both molecular coordinates and a solvent coordinate. We extend and apply this model for the S1-S0 CI in protonated Schiff bases, representing a model for retinal isomerization (Faraday Discuss. 2004, 127, 395, 2004). A dielectric continuum description of the solvent is combined with a minimal, two-electron-two-orbital electronic structure model according to Bonacić-Koutecký, Koutecký, and Michl (Angew. Chem. 1987, 26, 170), which characterizes the charge translocation effects at the CI. The model predicts that the nonequilibrium solvent state resulting from the S0-->S1 Franck-Condon transition can entail the disappearance of the CI, such that solvent motion is necessary to reach the CI seam. The concerted evolution of the intramolecular coordinates and the solvent coordinate is illustrated by an excited-state minimum energy path.
Collapse
Affiliation(s)
- Irene Burghardt
- Département de Chimie, Ecole Normale Supérieure, 24 rue Lhomond, F-75231 Paris Cedex 05, France
| | | |
Collapse
|
50
|
Villani G. Theoretical investigation of hydrogen transfer mechanism in the guanine–cytosine base pair. Chem Phys 2006. [DOI: 10.1016/j.chemphys.2005.11.006] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|