1
|
Karashima S, Humeniuk A, Suzuki T. Vibrational Motions in Ultrafast Electronic Relaxation of Pyrazine. J Am Chem Soc 2024. [PMID: 38621175 DOI: 10.1021/jacs.4c02886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2024]
Abstract
Ultrafast internal conversion via a conical intersection is ubiquitous in highly efficient photochemical reactions. Internal conversion from the 1ππ* to the 1nπ* state of pyrazine is the paradigm for this phenomenon; however, the relaxation occurs in such a short time (<20 fs) that the nuclear motion is difficult to observe in real time. The present study precisely measures the vibrational coherence transferred from the 1ππ* state to the 1nπ* state using time-resolved photoelectron spectroscopy with an unprecedented time resolution of 13.3 fs and reveals the key nuclear motions that drive the internal conversion.
Collapse
Affiliation(s)
- Shutaro Karashima
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| | - Alexander Humeniuk
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| | - Toshinori Suzuki
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| |
Collapse
|
2
|
Taylor JT, Tozer DJ, Curchod BFE. On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2). J Chem Phys 2023; 159:214115. [PMID: 38059547 DOI: 10.1063/5.0176140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Accepted: 11/14/2023] [Indexed: 12/08/2023] Open
Abstract
Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.
Collapse
Affiliation(s)
- Jack T Taylor
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David J Tozer
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
3
|
Xu Q, Liu Y, Wang M, Cerezo J, Improta R, Santoro F. The Resonance Raman Spectrum of Cytosine in Water: Analysis of the Effect of Specific Solute-Solvent Interactions and Non-Adiabatic Couplings. Molecules 2023; 28:molecules28052286. [PMID: 36903532 PMCID: PMC10005559 DOI: 10.3390/molecules28052286] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Revised: 02/21/2023] [Accepted: 02/22/2023] [Indexed: 03/06/2023] Open
Abstract
In this contribution, we report a computational study of the vibrational Resonance Raman (vRR) spectra of cytosine in water, on the grounds of potential energy surfaces (PES) computed by time-dependent density functional theory (TD-DFT) and CAM-B3LYP and PBE0 functionals. Cytosine is interesting because it is characterized by several close-lying and coupled electronic states, challenging the approach commonly used to compute the vRR for systems where the excitation frequency is in quasi-resonance with a single state. We adopt two recently developed time-dependent approaches, based either on quantum dynamical numerical propagations of vibronic wavepackets on coupled PES or on analytical correlation functions for cases in which inter-state couplings were neglected. In this way, we compute the vRR spectra, considering the quasi-resonance with the eight lowest-energy excited states, disentangling the role of their inter-state couplings from the mere interference of their different contributions to the transition polarizability. We show that these effects are only moderate in the excitation energy range explored by experiments, where the spectral patterns can be rationalized from the simple analysis of displacements of the equilibrium positions along the different states. Conversely, at higher energies, interference and inter-state couplings play a major role, and the adoption of a fully non-adiabatic approach is strongly recommended. We also investigate the effect of specific solute-solvent interactions on the vRR spectra, by considering a cluster of cytosine, hydrogen-bonded by six water molecules, and embedded in a polarizable continuum. We show that their inclusion remarkably improves the agreement with the experiments, mainly altering the composition of the normal modes, in terms of internal valence coordinates. We also document cases, mostly for low-frequency modes, in which a cluster model is not sufficient, and more elaborate mixed quantum classical approaches, in explicit solvent models, need to be applied.
Collapse
Affiliation(s)
- Qiushuang Xu
- School of Physics Engineering, Qufu Normal University, Qufu 273165, China
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China
- 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, Yantai 264025, China
| | - Meishan Wang
- School of Physics Engineering, Qufu Normal University, Qufu 273165, China
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China
| | - Javier Cerezo
- 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
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95, I-80145 Napoli, Italy
- Correspondence: (R.I.); (F.S.)
| | - Fabrizio Santoro
- 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
- Correspondence: (R.I.); (F.S.)
| |
Collapse
|
4
|
Komarova K. Density Matrix via Few Dominant Observables for the Ultrafast Non-Radiative Decay in Pyrazine. J Chem Theory Comput 2023; 19:746-757. [PMID: 36657738 PMCID: PMC11137821 DOI: 10.1021/acs.jctc.2c01211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Indexed: 01/21/2023]
Abstract
Unraveling the density matrix of a non-stationary quantum state as an explicit function of a few observables provides a complementary view of quantum dynamics. We have recently developed a practical way to identify the minimal set of the dominant observables that govern the quantal dynamics even in the case of strong non-adiabatic effects and large anharmonicity [Komarova et al., J. Chem. Phys. 155, 204110 (2021)]. Fast convergence in the number of the dominant contributions is achieved when instead of the density matrix we describe the time-evolution of the surprisal, the logarithm of the density operator. In the present work, we illustrate the efficiency of the proposed approach using an example of the early time dynamics in pyrazine in a Hilbert space accounting for up to four vibrational normal modes, {Q10a, Q6a, Q1, and Q9a}, and two coupled electronic states, the optically dark B 1 3 u ( n π * ) and the bright B 1 2 u ( π π * ) states. Dynamics in four-dimensional (4D) configurational space involve 19,600 vibronic eigenstates. Our results reveal that the rate of the ultrafast population decay as well as the shape of the nuclear wave packets in 2D, accounting only for {Q10a,Q6a} normal modes, are accurately captured with only six dominant time-independent observables in the surprisal. Extension of the dynamics to 3D and 4D vibrational subspace requires only five additional constraints. The time-evolution of a quantum state in 4D vibrational space on two electronic states is thus compacted to only 11 time-dependent coefficients of these observables.
Collapse
Affiliation(s)
- Ksenia Komarova
- The Fritz Haber Center for Molecular
Dynamics and Institute of Chemistry, The
Hebrew University of Jerusalem, Jerusalem91904, Israel
| |
Collapse
|
5
|
Aranda D, García-González F, Avila Ferrer FJ, López-Tocón I, Soto J, Otero JC. Computational Model for Electrochemical Surface-Enhanced Raman Scattering: Key Role of the Surface Charges and Synergy between Electromagnetic and Charge-Transfer Enhancement Mechanisms. J Chem Theory Comput 2022; 18:6802-6815. [PMID: 36222738 DOI: 10.1021/acs.jctc.2c00633] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a computational model for electrochemical surface-enhanced Raman scattering (EC-SERS). The surface excess of charge induced by the electrode potential (Vel) was introduced by applying an external electric field to a set of clusters [Agn]q with (n, q) of (19, ±1) or (20, 0) on which a molecule adsorbs. Using DFT/TD-DFT calculations, these metal-molecule complexes were classified by the adsorbate partial charge, and the main Vel-dependent properties were simultaneously studied with the aid of vibronic resonance Raman computations, namely, changes on the vibrational wavenumbers, relative intensities, and enhancement factors (EFs) for all SERS mechanisms: chemical or nonresonant, and resonance Raman with bright states of the adsorbate, charge-transfer (CT) states, and plasmon-like excitations on the metal cluster. We selected two molecules to test our model, pyridine, for which Vel has a remarkable effect, and 9,10-bis((E)-2-(pyridin-4-yl)vinyl)anthracene, which is almost insensitive to the applied bias. The results nicely reproduced most of the experimental observations, while the limitations of our approach were critically evaluated. We detected that accounting explicitly for the surface charges is key for EC-SERS models and that the highest calculated EFs, up to 107 to 108, are obtained by interstate coupling of bright local excitations of the metal cluster and CT states. These results highlight the importance of nonadiabatic effects in SERS and the capabilities of EC-SERS as a technique with potential to study excited-state coupling by tuning the CT and plasmon-like states by manipulating Vel.
Collapse
Affiliation(s)
- Daniel Aranda
- Instituto de Ciencia Molecular (ICMol), Universidad de Valencia, Catedrático J. Beltrán 2, 46980Paterna, Valencia, Spain.,Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Francisco García-González
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Francisco José Avila Ferrer
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Isabel López-Tocón
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Juan Soto
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Juan Carlos Otero
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| |
Collapse
|
6
|
Xu Q, Aranda D, Yaghoubi Jouybari M, Liu Y, Wang M, Cerezo J, Improta R, Santoro F. Nonadiabatic Vibrational Resonance Raman Spectra from Quantum Dynamics Propagations with LVC Models. Application to Thymine. J Phys Chem A 2022; 126:7468-7479. [PMID: 36099554 PMCID: PMC9596142 DOI: 10.1021/acs.jpca.2c05271] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We present a viable protocol to compute vibrational resonance
Raman
(vRR) spectra for systems with several close-lying and potentially
coupled electronic states. It is based on the parametrization of linear
vibronic coupling (LVC) models from time-dependent density functional
theory (TD-DFT) calculations and quantum dynamics propagations of
vibronic wavepackets with the multilayer version of the multiconfiguration
time-dependent Hartree (ML-MCTDH) method. Our approach is applied
to thymine considering seven coupled electronic states, comprising
the three lowest bright states, and all vibrational coordinates. Computed
vRR at different excitation wavelengths are in good agreement with
the available experimental data. Up to 250 nm the signal is dominated
by the lowest HOMO → LUMO transition, whereas at 233 nm, in
the valley between the two lowest energy absorption bands, the contributions
of all the three bright states, and their interferences and couplings,
are important. Inclusion of solvent (water) effects improves the agreement
with experiment, reproducing the coalescence of vibrational bands
due to CC and C=O stretchings. With our approach we disentangle
and assess the effect of interferences between the contribution of
different quasi-resonant states to the transition polarizability and
the effect of interstate couplings. Our findings strongly suggest
that in cases of close-lying and potentially coupled states a simple
inclusion of interference effects is not sufficient, and a fully nonadiabatic
computation should instead be performed. We also document that for
systems with strong couplings and quasi-degenerate states, the use
of HT perturbative approach, not designed for these cases, may lead
to large artifacts.
Collapse
Affiliation(s)
- Qiushuang Xu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China.,School of Physics Engineering, Qufu Normal University, 2673100 Qufu, Shandong, PR China.,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
| | - Daniel Aranda
- Instituto de Ciencia Molecular (ICMol)., Universidad de Valencia, c/Catedrático José Beltrán, 2, 46980 Paterna, Spain
| | - Martha Yaghoubi Jouybari
- 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
| | - Meishan Wang
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Javier Cerezo
- 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.,Departamento de Química, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95, I-80145 Napoli, Italy.,DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Fabrizio Santoro
- 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
|
7
|
Majumdar A, Jansen TLC. Quantum-Classical Simulation of Molecular Motors Driven Only by Light. J Phys Chem Lett 2021; 12:5512-5518. [PMID: 34096725 PMCID: PMC8279737 DOI: 10.1021/acs.jpclett.1c00951] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2021] [Accepted: 05/29/2021] [Indexed: 06/12/2023]
Abstract
Molecular motors that exhibit controlled unidirectional rotation provide great prospects for many types of applications, including nanorobotics. Existing rotational motors have two key components: photoisomerization around a π-bond followed by a thermally activated helical inversion, the latter being the rate-determining step. We propose an alternative molecular system in which the rotation is caused by the electric coupling of chromophores. This is used to engineer the excited state energy surface and achieve unidirectional rotation using light as the only input and avoid the slow thermally activated step, potentially leading to much faster operational speeds. To test the working principle, we employ quantum-classical calculations to study the dynamics of such a system. We estimate that motors built on this principle should be able to work on a subnanosecond time scale for such a full rotation. We explore the parameter space of our model to guide the design of a molecule that can act as such a motor.
Collapse
|
8
|
Choi S, Vaníček J. A time-reversible integrator for the time-dependent Schrödinger equation on an adaptive grid. J Chem Phys 2019; 151:234102. [PMID: 31864241 DOI: 10.1063/1.5127856] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
One of the most accurate methods for solving the time-dependent Schrödinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid points, we let the grid move together with the wavepacket but find that the naïve algorithm based on an alternate evolution of the wavefunction and grid destroys the time reversibility of the exact evolution. Yet, we show that the time reversibility is recovered if the wavefunction and grid are evolved simultaneously during each kinetic or potential step; this is achieved by using the Ehrenfest theorem together with the splitting method. The proposed algorithm is conditionally stable, symmetric, and time-reversible and conserves the norm of the wavefunction. The preservation of these geometric properties is shown analytically and demonstrated numerically on a three-dimensional harmonic model and collinear model of He-H2 scattering. We also show that the proposed algorithm can be symmetrically composed to obtain time-reversible integrators of an arbitrary even order. We observed 10 000-fold speedup by using the tenth-order instead of the second-order method to obtain a solution with a time discretization error below 10-9. Moreover, using the adaptive grid instead of the fixed grid resulted in a 64-fold reduction in the required number of grid points in the harmonic system and made it possible to simulate the He-H2 scattering for six times longer while maintaining reasonable accuracy. Applicability of the algorithm to high-dimensional quantum dynamics is demonstrated using the strongly anharmonic eight-dimensional Hénon-Heiles model.
Collapse
Affiliation(s)
- Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
9
|
Roulet J, Choi S, Vaníček J. Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation. J Chem Phys 2019; 150:204113. [PMID: 31153180 DOI: 10.1063/1.5094046] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Exact nonadiabatic quantum evolution preserves many geometric properties of the molecular Hilbert space. In the first paper of this series ["Paper I," S. Choi and J. Vaníček, J. Chem. Phys. 150, 204112 (2019)], we presented numerical integrators of arbitrary-order of accuracy that preserve these geometric properties exactly even in the adiabatic representation, in which the molecular Hamiltonian is not separable into kinetic and potential terms. Here, we focus on the separable Hamiltonian in diabatic representation, where the split-operator algorithm provides a popular alternative because it is explicit and easy to implement, while preserving most geometric invariants. Whereas the standard version has only second-order accuracy, we implemented, in an automated fashion, its recursive symmetric compositions, using the same schemes as in Paper I, and obtained integrators of arbitrary even order that still preserve the geometric properties exactly. Because the automatically generated splitting coefficients are redundant, we reduce the computational cost by pruning these coefficients and lower memory requirements by identifying unique coefficients. The order of convergence and preservation of geometric properties are justified analytically and confirmed numerically on a one-dimensional two-surface model of NaI and a three-dimensional three-surface model of pyrazine. As for efficiency, we find that to reach a convergence error of 10-10, a 600-fold speedup in the case of NaI and a 900-fold speedup in the case of pyrazine are obtained with the higher-order compositions instead of the second-order split-operator algorithm. The pyrazine results suggest that the efficiency gain survives in higher dimensions.
Collapse
Affiliation(s)
- Julien Roulet
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | - Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| |
Collapse
|
10
|
Choi S, Vaníček J. Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation. J Chem Phys 2019; 150:204112. [PMID: 31153205 DOI: 10.1063/1.5092611] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Geometric integrators of the Schrödinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement but, unfortunately, is restricted to systems whose Hamiltonian is separable into kinetic and potential terms. Here, we describe several implicit geometric integrators applicable to both separable and nonseparable Hamiltonians and, in particular, to the nonadiabatic molecular Hamiltonian in the adiabatic representation. These integrators combine the dynamic Fourier method with the recursive symmetric composition of the trapezoidal rule or implicit midpoint method, which results in an arbitrary order of accuracy in the time step. Moreover, these integrators are exactly unitary, symplectic, symmetric, time-reversible, and stable and, in contrast to the split-operator algorithm, conserve energy exactly, regardless of the accuracy of the solution. The order of convergence and conservation of geometric properties are proven analytically and demonstrated numerically on a two-surface NaI model in the adiabatic representation. Although each step of the higher order integrators is more costly, these algorithms become the most efficient ones if higher accuracy is desired; a thousand-fold speedup compared to the second-order trapezoidal rule (the Crank-Nicolson method) was observed for a wavefunction convergence error of 10-10. In a companion paper [J. Roulet, S. Choi, and J. Vaníček, J. Chem. Phys. 150, 204113 (2019)], we discuss analogous, arbitrary-order compositions of the split-operator algorithm and apply both types of geometric integrators to a higher-dimensional system in the diabatic representation.
Collapse
Affiliation(s)
- Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
11
|
Chen L, Gelin MF, Domcke W. Multimode quantum dynamics with multiple Davydov D2 trial states: Application to a 24-dimensional conical intersection model. J Chem Phys 2019; 150:024101. [DOI: 10.1063/1.5066022] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lipeng Chen
- Department of Chemistry, Technische Universität München, D-85747, Garching, Germany
| | - Maxim F. Gelin
- Department of Chemistry, Technische Universität München, D-85747, Garching, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technische Universität München, D-85747, Garching, Germany
| |
Collapse
|
12
|
Mignolet B, Kanno M, Shimakura N, Koseki S, Remacle F, Kono H, Fujimura Y. Ultrafast nonradiative transition pathways in photo-excited pyrazine: Ab initio analysis of time-resolved vacuum ultraviolet photoelectron spectrum. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.07.015] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
13
|
Antipov SV, Bhattacharyya S, El Hage K, Xu ZH, Meuwly M, Rothlisberger U, Vaníček J. Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2017; 4:061509. [PMID: 29376107 PMCID: PMC5758379 DOI: 10.1063/1.4996559] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Accepted: 10/11/2017] [Indexed: 06/07/2023]
Abstract
Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.
Collapse
Affiliation(s)
- Sergey V Antipov
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Swarnendu Bhattacharyya
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Krystel El Hage
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Zhen-Hao Xu
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
14
|
Greene SM, Batista VS. Tensor-Train Split-Operator Fourier Transform (TT-SOFT) Method: Multidimensional Nonadiabatic Quantum Dynamics. J Chem Theory Comput 2017; 13:4034-4042. [DOI: 10.1021/acs.jctc.7b00608] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Samuel M. Greene
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
- Energy Sciences Institute, Yale University, P.O. Box 27394, West Haven, Connecticut 06516-7394, United States
| | - Victor S. Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
- Energy Sciences Institute, Yale University, P.O. Box 27394, West Haven, Connecticut 06516-7394, United States
| |
Collapse
|
15
|
Farag MH, Jansen TLC, Knoester J. Probing the Interstate Coupling near a Conical Intersection by Optical Spectroscopy. J Phys Chem Lett 2016; 7:3328-3334. [PMID: 27509384 DOI: 10.1021/acs.jpclett.6b01463] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Conical intersections are points where adiabatic potential energy surfaces cross. The interstate coupling between the potential energy surfaces plays a crucial role in many processes associated with conical intersections. Still no method exists to measure this coupling driving the chemical reactions between the potential energy surfaces involved. In this Letter, using a generic model for photoisomerization, we propose a novel experimental approach to estimate the coupling that mixes the electronic states near a conical intersection. The approach is based on analyzing the vibrational wavepacket of the reactant in the adiabatic ground and excited electronic states. The nuclear wavepacket dynamics are extracted from linear absorption and two-dimensional electronic spectroscopy. Comparing the frequencies of the coupling mode in the adiabatic ground and excited states from models with and without coupling between the potential energy surfaces suggests an experimental tool to determine the interstate coupling.
Collapse
Affiliation(s)
- Marwa H Farag
- Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Jasper Knoester
- Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
16
|
Reddy CS, Prasad MD. A Gaussian Wave Packet Propagation Approach to Vibrationally Resolved Optical Spectra at Non-Zero Temperatures. J Phys Chem A 2016; 120:2583-90. [PMID: 27035861 DOI: 10.1021/acs.jpca.6b00308] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An effective time dependent approach based on a method that is similar to the Gaussian wave packet propagation (GWP) technique of Heller is developed for the computation of vibrationally resolved electronic spectra at finite temperatures in the harmonic, Franck-Condon/Hertzberg-Teller approximations. Since the vibrational thermal density matrix of the ground electronic surface and the time evolution operator on that surface commute, it is possible to write the spectrum generating correlation function as a trace of the time evolved doorway state. In the stated approximations, the doorway state is a superposition of the harmonic oscillator zero and one quantum eigenfunctions and thus can be propagated by the GWP. The algorithm has an O(N(3)) dependence on the number of vibrational modes. An application to pyrene absorption spectrum at two temperatures is presented as a proof of the concept.
Collapse
Affiliation(s)
- Ch Sridhar Reddy
- School of Chemistry, University of Hyderabad , Hyderabad 500 046, India
| | - M Durga Prasad
- School of Chemistry, University of Hyderabad , Hyderabad 500 046, India
| |
Collapse
|
17
|
Chen L, Gelin MF, Chernyak VY, Domcke W, Zhao Y. Dissipative dynamics at conical intersections: simulations with the hierarchy equations of motion method. Faraday Discuss 2016; 194:61-80. [DOI: 10.1039/c6fd00088f] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The effect of a dissipative environment on the ultrafast nonadiabatic dynamics at conical intersections is analyzed for a two-state two-mode model chosen to represent the S2(ππ*)–S1(nπ*) conical intersection in pyrazine (the system) which is bilinearly coupled to infinitely many harmonic oscillators in thermal equilibrium (the bath). The system–bath coupling is modeled by the Drude spectral function. The equation of motion for the reduced density matrix of the system is solved numerically exactly with the hierarchy equation of motion method using graphics-processor-unit (GPU) technology. The simulations are valid for arbitrary strength of the system–bath coupling and arbitrary bath memory relaxation time. The present computational studies overcome the limitations of weak system–bath coupling and short memory relaxation time inherent in previous simulations based on multi-level Redfield theory [A. Kühl and W. Domcke, J. Chem. Phys. 2002, 116, 263]. Time evolutions of electronic state populations and time-dependent reduced probability densities of the coupling and tuning modes of the conical intersection have been obtained. It is found that even weak coupling to the bath effectively suppresses the irregular fluctuations of the electronic populations of the isolated two-mode conical intersection. While the population of the upper adiabatic electronic state (S2) is very efficiently quenched by the system–bath coupling, the population of the diabatic ππ* electronic state exhibits long-lived oscillations driven by coherent motion of the tuning mode. Counterintuitively, the coupling to the bath can lead to an enhanced lifetime of the coherence of the tuning mode as a result of effective damping of the highly excited coupling mode, which reduces the strong mode–mode coupling inherent to the conical intersection. The present results extend previous studies of the dissipative dynamics at conical intersections to the nonperturbative regime of system–bath coupling. They pave the way for future first-principles simulations of femtosecond time-resolved four-wave-mixing spectra of chromophores in condensed phases which are nonperturbative in the system dynamics, the system–bath coupling as well as the field-matter coupling.
Collapse
Affiliation(s)
- Lipeng Chen
- Division of Materials Science
- Nanyang Technological University
- Singapore 639798
| | - Maxim F. Gelin
- Department of Chemistry
- Technische Universität München
- Garching D-85747
- Germany
| | | | - Wolfgang Domcke
- Department of Chemistry
- Technische Universität München
- Garching D-85747
- Germany
| | - Yang Zhao
- Division of Materials Science
- Nanyang Technological University
- Singapore 639798
| |
Collapse
|
18
|
Krčmář J, Gelin MF, Domcke W. Simulation of femtosecond two-dimensional electronic spectra of conical intersections. J Chem Phys 2015; 143:074308. [PMID: 26298135 DOI: 10.1063/1.4928685] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have simulated femtosecond two-dimensional (2D) electronic spectra for an excited-state conical intersection using the wave-function version of the equation-of-motion phase-matching approach. We show that 2D spectra at fixed values of the waiting time provide information on the structure of the vibronic eigenstates of the conical intersection, while the evolution of the spectra with the waiting time reveals predominantly ground-state wave-packet dynamics. The results show that 2D spectra of conical intersection systems differ significantly from those obtained for chromophores with well separated excited-state potential-energy surfaces. The spectral signatures which can be attributed to conical intersections are discussed.
Collapse
Affiliation(s)
- Jindřich Krčmář
- Department of Chemistry, Technische Universität München, D-85747 Garching, Germany
| | - Maxim F Gelin
- Department of Chemistry, Technische Universität München, D-85747 Garching, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technische Universität München, D-85747 Garching, Germany
| |
Collapse
|
19
|
Reddy CS, Prasad MD. Finite temperature vibronic spectra of harmonic surfaces: a time-dependent coupled cluster approach. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1070928] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
20
|
Zimmermann T, Vaníček J. Efficient on-the-fly ab initio semiclassical method for computing time-resolved nonadiabatic electronic spectra with surface hopping or Ehrenfest dynamics. J Chem Phys 2015; 141:134102. [PMID: 25296779 DOI: 10.1063/1.4896735] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive a somewhat crude, yet very efficient semiclassical approximation for computing nonadiabatic spectra. The resulting method, which is a generalization of the multiple-surface dephasing representation, includes quantum effects through interference of mixed quantum-classical trajectories and through quantum treatment of the collective electronic degree of freedom. The method requires very little computational effort beyond the fewest-switches surface hopping or Ehrenfest locally mean-field dynamics and is very easy to implement. The proposed approximation is tested by computing the absorption and time-resolved stimulated emission spectra of pyrazine using the four-dimensional three-surface model which allows for comparison with the numerically exact quantum spectra. As expected, the multiple-surface dephasing representation is not suitable for high-resolution linear spectra, yet it seems to capture all the important features of pump-probe spectra. Finally, the method is combined with on-the-fly ab initio evaluation of the electronic structure (i.e., energies, forces, electric-dipole, and nonadiabatic couplings) in order to compute fully dimensional nonadiabatic spectra of pyrazine without approximations inherent to analytical, including vibronic-coupling models. The Appendix provides derivations of perturbative expressions for linear and pump-probe spectra of arbitrary mixed states and for arbitrary laser pulse shapes.
Collapse
Affiliation(s)
- Tomáš Zimmermann
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
21
|
Sala M, Lasorne B, Gatti F, Guérin S. The role of the low-lying dark nπ* states in the photophysics of pyrazine: a quantum dynamics study. Phys Chem Chem Phys 2015; 16:15957-67. [PMID: 24964033 DOI: 10.1039/c4cp02165g] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The excited state dynamics of pyrazine has attracted considerable attention in the last three decades. It has long been recognized that after UV excitation, the dynamics of the molecule is impacted by strong non-adiabatic effects due to the existence of a conical intersection between the B2u(ππ*) and B3u(nπ*) electronic states. However, a recent study based on trajectory surface hopping dynamics simulations suggested the participation of the Au(nπ*) and B2g(nπ*) low-lying dark electronic states in the ultrafast radiationless decay of the molecule after excitation to the B2u(ππ*) state. The purpose of this work was to pursue the investigation of the role of the Au(nπ*) and B2g(nπ*) states in the photophysics of pyrazine. A linear vibronic coupling model hamiltonian including the four lowest excited electronic states and the sixteen most relevant vibrational degrees of freedom was constructed using high level XMCQDPT2 electronic structure calculations. Wavepacket propagations using the MCTDH method were then performed and used to simulate the absorption spectrum and the electronic state population dynamics of the system. Our results show that the Au(nπ*) state plays an important role in the photophysics of pyrazine.
Collapse
Affiliation(s)
- Matthieu Sala
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS, Université de Bourgogne, BP 47870, F-21078 Dijon, France.
| | | | | | | |
Collapse
|
22
|
Kanno M, Ito Y, Shimakura N, Koseki S, Kono H, Fujimura Y. Ab initio quantum dynamical analysis of ultrafast nonradiative transitions via conical intersections in pyrazine. Phys Chem Chem Phys 2015; 17:2012-24. [DOI: 10.1039/c4cp04807e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Ab initio quantum dynamical calculations reveal the contributions of the optically dark S1 (1B3u), S3 (1Au), and S4 (1B2g) states to ultrafast nonradiative transition from the optically bright S2 (1B2u) state in pyrazine and suggest a much faster transition rate than previously believed.
Collapse
Affiliation(s)
- Manabu Kanno
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Yuta Ito
- Department of Chemistry
- Niigata University
- Niigata 950-2181
- Japan
| | | | - Shiro Koseki
- Department of Chemistry
- Graduate School of Science, and Research Institute for Molecular Electronic Devices (RIMED)
- Osaka Prefecture University
- Sakai 599-8531
- Japan
| | - Hirohiko Kono
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Yuichi Fujimura
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| |
Collapse
|
23
|
Sala M, Guérin S, Gatti F. Quantum dynamics of the photostability of pyrazine. Phys Chem Chem Phys 2015; 17:29518-30. [DOI: 10.1039/c5cp04605j] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
We propose a new mechanism for the radiationless decay of photoexcited pyrazine to its ground electronic state involving a conical intersection between the dark Au(nπ) state and the ground state.
Collapse
Affiliation(s)
- Matthieu Sala
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS
- Université de Bourgogne Franche-Comté
- F-21078 Dijon
- France
| | - Stéphane Guérin
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS
- Université de Bourgogne Franche-Comté
- F-21078 Dijon
- France
| | - Fabien Gatti
- CTMM
- Institut Charles Gerhardt UMR 5253 CNRS
- F-34095 Montpellier
- France
| |
Collapse
|
24
|
Wehrle M, Šulc M, Vaníček J. On-the-fly ab initio semiclassical dynamics: Identifying degrees of freedom essential for emission spectra of oligothiophenes. J Chem Phys 2014; 140:244114. [DOI: 10.1063/1.4884718] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Marius Wehrle
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Miroslav Šulc
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
25
|
Sala M, Saab M, Lasorne B, Gatti F, Guérin S. Laser control of the radiationless decay in pyrazine using the dynamic Stark effect. J Chem Phys 2014; 140:194309. [DOI: 10.1063/1.4875736] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Šulc M, Hernández H, Martínez TJ, Vaníček J. Relation of exact Gaussian basis methods to the dephasing representation: theory and application to time-resolved electronic spectra. J Chem Phys 2014; 139:034112. [PMID: 23883015 DOI: 10.1063/1.4813124] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We recently showed that the dephasing representation (DR) provides an efficient tool for computing ultrafast electronic spectra and that further acceleration is possible with cellularization [M. Šulc and J. Vaníček, Mol. Phys. 110, 945 (2012)]. Here, we focus on increasing the accuracy of this approximation by first implementing an exact Gaussian basis method, which benefits from the accuracy of quantum dynamics and efficiency of classical dynamics. Starting from this exact method, the DR is derived together with ten other methods for computing time-resolved spectra with intermediate accuracy and efficiency. These methods include the Gaussian DR, an exact generalization of the DR, in which trajectories are replaced by communicating frozen Gaussian basis functions evolving classically with an average Hamiltonian. The newly obtained methods are tested numerically on time correlation functions and time-resolved stimulated emission spectra in the harmonic potential, pyrazine S0∕S1 model, and quartic oscillator. Numerical results confirm that both the Gaussian basis method and the Gaussian DR increase the accuracy of the DR. Surprisingly, in chaotic systems the Gaussian DR can outperform the presumably more accurate Gaussian basis method, in which the two bases are evolved separately.
Collapse
Affiliation(s)
- Miroslav Šulc
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | | | | | |
Collapse
|
27
|
Zambrano E, Šulc M, Vaníček J. Improving the accuracy and efficiency of time-resolved electronic spectra calculations: Cellular dephasing representation with a prefactor. J Chem Phys 2013; 139:054109. [DOI: 10.1063/1.4817005] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Avila Ferrer FJ, Barone V, Cappelli C, Santoro F. Duschinsky, Herzberg–Teller, and Multiple Electronic Resonance Interferential Effects in Resonance Raman Spectra and Excitation Profiles. The Case of Pyrene. J Chem Theory Comput 2013; 9:3597-611. [DOI: 10.1021/ct400197y] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Francisco J. Avila Ferrer
- Consiglio Nazionale delle Ricerche
− CNR, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), UOS di Pisa, Area della Ricerca, Via G. Moruzzi 1,
I-56124 Pisa, Italy
- Physical Chemistry, Faculty
of Science, University of Málaga, Málaga 29071, Spain
| | - Vincenzo Barone
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa,
Italy
| | - Chiara Cappelli
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa,
Italy
- Dipartimento di Chimica
e Chimica
Industriale dell’Università di Pisa, Via Risorgimento 35, I-56126 Pisa, Italy
| | - Fabrizio Santoro
- Consiglio Nazionale delle Ricerche
− CNR, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), UOS di Pisa, Area della Ricerca, Via G. Moruzzi 1,
I-56124 Pisa, Italy
| |
Collapse
|
29
|
Shiozaki T, Woywod C, Werner HJ. Pyrazine excited states revisited using the extended multi-state complete active space second-order perturbation method. Phys Chem Chem Phys 2013; 15:262-9. [DOI: 10.1039/c2cp43381h] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
30
|
Nikoobakht B, Köppel H, Gindensperger E, Cederbaum LS. Efficient computation of adiabatic electronic populations in multi-mode vibronic systems: theory, implementation, and application. J Chem Phys 2012; 137:114110. [PMID: 22998252 DOI: 10.1063/1.4753966] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The effective-mode formalism developed earlier is applied to efficiently compute adiabatic electronic populations of multi-mode vibronic systems. Two different versions of the formalism are utilized. In the first one, the effective modes are used as new vibrational basis, and the time-dependent wave function as well as populations are calculated numerically exactly. In the second variant, the hierarchy-of-modes formalism is applied as an approximation scheme which leads to accurate results when including typically 7-10 members of the hierarchy. While in the first version, the propagation of the wave packet becomes numerically tedious, the computation of the adiabatic populations is rendered possible at all. Through the hierarchy-of-modes formalism, both the propagation as well as computation of adiabatic populations are speeded up by 1-3 orders of magnitude in typical cases. The formalism is applied to study the dynamics of pyrazine-type models for varying vibronic coupling strength and a (very) large number of bath modes.
Collapse
Affiliation(s)
- Behnam Nikoobakht
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, INF 229, D-69120 Heidelberg, Germany
| | | | | | | |
Collapse
|
31
|
Thanopulos I, Li X, Brumer P, Shapiro M. Time-dependent partitioning theory of the control of radiationless transitions in 24-mode pyrazine. J Chem Phys 2012; 137:064111. [DOI: 10.1063/1.4739275] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
|
32
|
Silverstein DW, Jensen L. Vibronic coupling simulations for linear and nonlinear optical processes: Theory. J Chem Phys 2012; 136:064111. [DOI: 10.1063/1.3684236] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Lin CK, Niu Y, Zhu C, Shuai Z, Lin SH. The Role of the nπ* 1Au State in the Photoabsorption and Relaxation of Pyrazine. Chem Asian J 2011; 6:2977-85. [PMID: 21922672 DOI: 10.1002/asia.201100472] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2011] [Indexed: 11/10/2022]
Affiliation(s)
- Chih-Kai Lin
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu 30010, Taiwan.
| | | | | | | | | |
Collapse
|
34
|
Thanopulos I, Brumer P, Shapiro M. Intramolecular energy transfer dynamics in 24-mode pyrazine by partitioning technique: A time-dependent approach. J Chem Phys 2010; 133:154111. [DOI: 10.1063/1.3495953] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
|
35
|
Suzuki YI, Fuji T, Horio T, Suzuki T. Time-resolved photoelectron imaging of ultrafast S2→S1 internal conversion through conical intersection in pyrazine. J Chem Phys 2010; 132:174302. [PMID: 20459164 DOI: 10.1063/1.3395206] [Citation(s) in RCA: 80] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yoshi-Ichi Suzuki
- CREST, Japan Science and Technology Agency, Sanbancho, Chiyoda-ku, Tokyo 102-0075, Japan
| | | | | | | |
Collapse
|
36
|
Woywod C, Papp A, Halász GJ, Vibók Á. Theoretical investigation of the electronic spectrum of pyrazine. Theor Chem Acc 2009. [DOI: 10.1007/s00214-009-0678-x] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
37
|
Ab initio studies of excited electronic state S2 of pyrazine and Franck–Condon simulation of its absorption spectrum. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2009.05.043] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
38
|
Islampour R, Miralinaghi M. Theoretical Study of Internal Conversion Decay Rates Associated with the Three Lowest Singlet Electronic States in Pyrazine. J Phys Chem A 2009; 113:2340-9. [DOI: 10.1021/jp8072206] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Reza Islampour
- Department of Chemistry, Tarbiat Moallem University, 49 Mofateh Avenue, Tehran, Iran
| | - Mahsa Miralinaghi
- Department of Chemistry, Tarbiat Moallem University, 49 Mofateh Avenue, Tehran, Iran
| |
Collapse
|
39
|
Burghardt I, Giri K, Worth GA. Multimode quantum dynamics using Gaussian wavepackets: The Gaussian-based multiconfiguration time-dependent Hartree (G-MCTDH) method applied to the absorption spectrum of pyrazine. J Chem Phys 2008; 129:174104. [DOI: 10.1063/1.2996349] [Citation(s) in RCA: 145] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
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]
|
41
|
Puzari P, Sarkar B, Adhikari S. A quantum-classical approach to the molecular dynamics of pyrazine with a realistic model Hamiltonian. J Chem Phys 2006; 125:194316. [PMID: 17129112 DOI: 10.1063/1.2393228] [Citation(s) in RCA: 37] [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 investigate the molecular dynamics of pyrazine after excitation to the S2 electronic state by using the time-dependent discrete variable representation (TDDVR) method. The investigation has been carried out with a realistic 24-mode model Hamiltonian consisting of all the vibrational degrees of freedom of pyrazine molecule. First, we perform the simulation on a basic four-mode model, and then by including additional eight important modes and finally, by introducing 20 bath modes on the basic model. This sequential inclusion of bath modes demonstrates the effect of weak modes on the subsystem, where the calculations of energy and population transfer from basic model to the bath quantify the same effect. The spectral profile obtained by using TDDVR approach shows reasonably good agreement with the results calculated by quantum mechanical approach. It appears that the TDDVR approach for those large systems where quantum mechanical description is needed in a restricted region is a good compromise between accuracy and speed.
Collapse
Affiliation(s)
- Panchanan Puzari
- Department of Chemistry, Indian Institute of Technology, Guwahati, North Guwahati, Guwahati 781 039, India
| | | | | |
Collapse
|
42
|
Chen X, Batista VS. Matching-pursuit/split-operator-Fourier-transform simulations of excited-state nonadiabatic quantum dynamics in pyrazine. J Chem Phys 2006; 125:124313. [PMID: 17014180 DOI: 10.1063/1.2356477] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A simple approach for numerically exact simulations of nonadiabatic quantum dynamics in multidimensional systems is introduced and applied to the description of the photoabsorption spectroscopy of pyrazine. The propagation scheme generalizes the recently developed matching-pursuit/split-operator-Fourier-transform (MP/SOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] to simulations of nonadiabatic quantum dynamics. The time-evolution operator is applied, as defined by the Trotter expansion to second order accuracy, in dynamically adaptive coherent-state expansions. These representations are obtained by combining the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting computational approach are demonstrated in calculations of time-dependent survival amplitudes and photoabsorption cross sections, using a model Hamiltonian that allows for direct comparisons with benchmark calculations. Simulations in full-dimensional potential energy surfaces involve the propagation of a 24-dimensional wave packet to describe the S(1)S(2) interconversion of pyrazine after S(0)-->S(2) photoexcitation. The reported results show that the generalized MP/SOFT method is a practical and accurate approach to model nonadiabatic reaction dynamics in polyatomic systems.
Collapse
Affiliation(s)
- Xin Chen
- Department of Chemistry, Yale University, New Haven, CT 06520-8107, USA
| | | |
Collapse
|
43
|
Neugebauer J, Hess BA. Resonance Raman spectra of uracil based on Kramers–Kronig relations using time-dependent density functional calculations and multireference perturbation theory. J Chem Phys 2004; 120:11564-77. [PMID: 15268191 DOI: 10.1063/1.1697371] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The use of time-dependent density functional calculations for the optimization of excited-state structures and the subsequent calculation of resonance Raman intensities within the transform-theory framework is compared to calculations of Hartree-Fock/configuration interaction singles-type (CIS). The transform theory of resonance Raman scattering is based on Kramers-Kronig relations between polarizability tensor components and the optical absorption. Stationary points for the two lowest excited singlet states of uracil are optimized and characterized by means of numerical differentiation of analytical excited-state gradients. It is shown that the effect of electron correlation leads to substantial modifications of the relative intensities. Calculations of vibrational frequencies for ground and excited states are carried out, which show that the neglect of Duschinsky mixing and the assumption of equal wave numbers for ground and excited state are not in all cases good approximations. We also compare the transform-theory resonance Raman intensities with those obtained within a simple approximation from excited-state gradients at the ground-state equilibrium position, and find that they are in qualitative agreement in the case of CIS, but show some important differences in calculations based on density functional theory. Since the results from CIS calculations are in better agreement with experiment, we also present approximate resonance Raman spectra obtained using excited-state gradients from multireference perturbation theory calculations, which confirm the CIS gradients.
Collapse
Affiliation(s)
- Johannes Neugebauer
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nurnberg, Egerlandstrasse 3, D-91058 Erlangen, Germany.
| | | |
Collapse
|
44
|
Thoss M, Wang H. SEMICLASSICAL DESCRIPTION OF MOLECULAR DYNAMICS BASED ON INITIAL-VALUE REPRESENTATION METHODS. Annu Rev Phys Chem 2004; 55:299-332. [PMID: 15117255 DOI: 10.1146/annurev.physchem.55.091602.094429] [Citation(s) in RCA: 137] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Recent progress in the development of semiclassical methods to describe quantum effects in molecular dynamics is reviewed. Focusing on rigorous semiclassical methods that are based on the initial-value representation of the semiclassical propagator, we discuss several promising schemes that have been developed in the past few years to extend the applicability of semiclassical approaches to complex molecular systems. In particular, integral-filtering techniques and forward-backward methods are surveyed. Furthermore, recently proposed approaches that allow the semiclassical description of nonadiabatic molecular dynamics are discussed. The potential and efficiency of these methods is illustrated by selected applications.
Collapse
Affiliation(s)
- Michael Thoss
- Theoretische Chemie, Technische Universitat Munchen, D-85747 Garching, Germany.
| | | |
Collapse
|
45
|
Suzuki YI, Stener M, Seideman T. Multidimensional calculation of time-resolved photoelectron angular distributions: The internal conversion dynamics of pyrazine. J Chem Phys 2003. [DOI: 10.1063/1.1536981] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Bonella S, Coker DF. Semiclassical implementation of the mapping Hamiltonian approach for nonadiabatic dynamics using focused initial distribution sampling. J Chem Phys 2003. [DOI: 10.1063/1.1542883] [Citation(s) in RCA: 69] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
|
47
|
Suzuki YI, Stener M, Seideman T. Theory of time-resolved photoelectron imaging: nonperturbative calculation for an internally converting polyatomic molecule. PHYSICAL REVIEW LETTERS 2002; 89:233002. [PMID: 12485003 DOI: 10.1103/physrevlett.89.233002] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2002] [Indexed: 05/24/2023]
Abstract
We present the first calculation of time-resolved photoelectron angular distributions for a polyatomic system. Our method takes rotations into exact account, treats the laser field nonperturbatively, and computes the electronic dynamics from first principles. Our results point to the information content of time-resolved photoelectron imaging observables and illustrate the role played by the field intensity.
Collapse
Affiliation(s)
- Yoshi-Ichi Suzuki
- Steacie Institute for Molecular Science, National Research Council of Canada, Ottawa, Ontario, Canada K1A OR6
| | | | | |
Collapse
|
48
|
|
49
|
Tran-Thi TH, Prayer C, Millié P, Uznanski P, Hynes JT. Substituent and Solvent Effects on the Nature of the Transitions of Pyrenol and Pyranine. Identification of an Intermediate in the Excited-State Proton-Transfer Reaction†. J Phys Chem A 2002. [DOI: 10.1021/jp0125606] [Citation(s) in RCA: 83] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
50
|
Talaga DS, Zink JI. Symmetry and Local Mode Coupling in Absorption and Raman Spectroscopy of Intervalence Electronic Transitions. J Phys Chem A 2001. [DOI: 10.1021/jp003285a] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David S. Talaga
- Department of Chemistry and Biochemistry,University of California, Los Angeles, Los Angeles, California 90095-1569
| | - Jeffrey I. Zink
- Department of Chemistry and Biochemistry,University of California, Los Angeles, Los Angeles, California 90095-1569
| |
Collapse
|