1
|
Huang X, Xie W, Došlić N, Gelin MF, Domcke W. Ab Initio Quasiclassical Simulation of Femtosecond Time-Resolved Two-Dimensional Electronic Spectra of Pyrazine. J Phys Chem Lett 2021; 12:11736-11744. [PMID: 34851116 DOI: 10.1021/acs.jpclett.1c03589] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.
Collapse
Affiliation(s)
- Xiang Huang
- Department of Chemistry, Technical University of Munich, Garching, D-85747, Germany
| | - Weiwei Xie
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, 76131, Germany
| | - Nađa Došlić
- Department of Physical Chemistry, Ruder Boscovic Institute, Zagreb, HR-10000, Croatia
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou, 310018, China
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Garching, D-85747, Germany
| |
Collapse
|
2
|
The excited state effective dipole moment of 2,3-benzofuran from thermochromic shifts in absorption and emission spectra. J Photochem Photobiol A Chem 2021. [DOI: 10.1016/j.jphotochem.2021.113476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
3
|
Segatta F, Nenov A, Nascimento DR, Govind N, Mukamel S, Garavelli M. iSPECTRON: A simulation interface for linear and nonlinear spectra with ab-initio quantum chemistry software. J Comput Chem 2021; 42:644-659. [PMID: 33556195 DOI: 10.1002/jcc.26485] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2020] [Revised: 01/05/2021] [Accepted: 01/11/2021] [Indexed: 12/18/2022]
Abstract
We introduce iSPECTRON, a program that parses data from common quantum chemistry software (NWChem, OpenMolcas, Gaussian, Cobramm, etc.), produces the input files for the simulation of linear and nonlinear spectroscopy of molecules with the Spectron code, and analyzes the spectra with a broad range of tools. Vibronic spectra are expressed in term of the electronic eigenstates, obtained from quantum chemistry computations, and vibrational/bath effects are incorporated in the framework of the displaced harmonic oscillator model, where all required quantities are computed at the Franck-Condon point. The program capabilities are illustrated by simulating linear absorption, transient absorption and two dimensional electronic spectra of the pyrene molecule. Calculations at two levels of electronic structure theory, time-dependent density functional theory (with NWChem) and RASSCF/RASPT2 (with OpenMolcas) are presented and compared where possible. The iSPECTRON program is available online at https://github.com/ispectrongit/iSPECTRON/ and distributed open source under the terms of the Educational Community License version 2.0 (ECL 2.0).
Collapse
Affiliation(s)
- Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| | - Daniel R Nascimento
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington, USA
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington, USA
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California, USA
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| |
Collapse
|
4
|
Aleotti F, Aranda D, Yaghoubi Jouybari M, Garavelli M, Nenov A, Santoro F. Parameterization of a linear vibronic coupling model with multiconfigurational electronic structure methods to study the quantum dynamics of photoexcited pyrene. J Chem Phys 2021; 154:104106. [PMID: 33722019 DOI: 10.1063/5.0044693] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
With this work, we present a protocol for the parameterization of a Linear Vibronic Coupling (LVC) Hamiltonian for quantum dynamics using highly accurate multiconfigurational electronic structure methods such as RASPT2/RASSCF, combined with a maximum-overlap diabatization technique. Our approach is fully portable and can be applied to many medium-size rigid molecules whose excited state dynamics requires a quantum description. We present our model and discuss the details of the electronic structure calculations needed for the parameterization, analyzing critical situations that could arise in the case of strongly interacting excited states. The protocol was applied to the simulation of the excited state dynamics of the pyrene molecule, starting from either the first or the second bright state (S2 or S5). The LVC model was benchmarked against state-of-the-art quantum mechanical calculations with optimizations and energy scans and turned out to be very accurate. The dynamics simulations, performed including all active normal coordinates with the multilayer multiconfigurational time-dependent Hartree method, show good agreement with the available experimental data, endorsing prediction of the excited state mechanism, especially for S5, whose ultrafast deactivation mechanism was not yet clearly understood.
Collapse
Affiliation(s)
- Flavia Aleotti
- Dipartimento di Chimica Industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Daniel Aranda
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Martha Yaghoubi Jouybari
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
5
|
Abou-Hatab S, Carnevale V, Matsika S. Modeling solvation effects on absorption and fluorescence spectra of indole in aqueous solution. J Chem Phys 2021; 154:064104. [PMID: 33588532 PMCID: PMC7878019 DOI: 10.1063/5.0038342] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2020] [Accepted: 01/20/2021] [Indexed: 11/14/2022] Open
Abstract
Modeling the optical spectra of molecules in solution presents a challenge, so it is important to understand which of the solvation effects (i.e., electrostatics, mutual polarization, and hydrogen bonding interactions between solute and solvent molecules) are crucial in reproducing the various features of the absorption and fluorescence spectra and to identify a sufficient theoretical model that accurately captures these effects with minimal computational cost. In this study, we use various implicit and explicit solvation models, such as molecular dynamics coupled with non-polarizable and polarizable force fields, as well as Car-Parrinello molecular dynamics, to model the absorption and fluorescence spectra of indole in aqueous solution. The excited states are computed using the equation of motion coupled cluster with single and double excitations combined with the effective fragment potential to represent water molecules, which we found to be a computationally efficient approach for modeling large solute-solvent clusters at a high level of quantum theory. We find that modeling mutual polarization, compared to other solvation effects, is a dominating factor for accurately reproducing the position of the peaks and spectral line shape of the absorption spectrum of indole in solution. We present an in-depth analysis of the influence that different solvation models have on the electronic excited states responsible for the features of the absorption spectra. Modeling fluorescence is more challenging since it is hard to reproduce even the correct emitting state, and force field parameters need to be re-evaluated.
Collapse
Affiliation(s)
- Salsabil Abou-Hatab
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Vincenzo Carnevale
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
6
|
Gu B, Nenov A, Segatta F, Garavelli M, Mukamel S. Manipulating Core Excitations in Molecules by X-Ray Cavities. PHYSICAL REVIEW LETTERS 2021; 126:053201. [PMID: 33605757 DOI: 10.1103/physrevlett.126.053201] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Revised: 12/02/2020] [Accepted: 12/22/2020] [Indexed: 06/12/2023]
Abstract
Core excitations on different atoms are highly localized and therefore decoupled. By placing molecules in an x-ray cavity the core transitions become coupled via the exchange of cavity photons and form delocalized hybrid light-matter excitations known as core polaritons. We demonstrate these effects for the two inequivalent carbon atoms in 1,1-difluoroethylene. Polariton signatures in the x-ray absorption, two-photon absorption, and multidimensional four-wave mixing signals are predicted.
Collapse
Affiliation(s)
- Bing Gu
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari," Università degli studi di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari," Università degli studi di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari," Università degli studi di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
| |
Collapse
|
7
|
Borrego-Varillas R, Nenov A, Ganzer L, Oriana A, Manzoni C, Tolomelli A, Rivalta I, Mukamel S, Garavelli M, Cerullo G. Two-dimensional UV spectroscopy: a new insight into the structure and dynamics of biomolecules. Chem Sci 2019. [DOI: 10.1039/c9sc03871j] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Two-dimensional ultraviolet spectroscopy has the potential to deliver rich structural and dynamical information on biomolecules such as DNA and proteins.
Collapse
Affiliation(s)
| | - A. Nenov
- Dipartimento di Chimica Industriale
- Universitá degli Studi di Bologna
- I-40136 Bologna
- Italy
| | - L. Ganzer
- IFN-CNR
- Dipartimento di Fisica
- Politecnico di Milano
- I-20133 Milano
- Italy
| | - A. Oriana
- IFN-CNR
- Dipartimento di Fisica
- Politecnico di Milano
- I-20133 Milano
- Italy
| | - C. Manzoni
- IFN-CNR
- Dipartimento di Fisica
- Politecnico di Milano
- I-20133 Milano
- Italy
| | - A. Tolomelli
- Dipartimento di Chimica
- Universitá degli Studi di Bologna
- I-40126 Bologna
- Italy
| | - I. Rivalta
- Dipartimento di Chimica Industriale
- Universitá degli Studi di Bologna
- I-40136 Bologna
- Italy
| | - S. Mukamel
- Department of Chemistry
- Department of Physics and Astronomy
- University of California
- Irvine
- USA
| | - M. Garavelli
- Dipartimento di Chimica Industriale
- Universitá degli Studi di Bologna
- I-40136 Bologna
- Italy
| | - G. Cerullo
- IFN-CNR
- Dipartimento di Fisica
- Politecnico di Milano
- I-20133 Milano
- Italy
| |
Collapse
|
8
|
Nenov A, Conti I, Borrego-Varillas R, Cerullo G, Garavelli M. Linear absorption spectra of solvated thiouracils resolved at the hybrid RASPT2/MM level. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.07.025] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|
9
|
COBRAMM 2.0 — A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations. J Mol Model 2018; 24:271. [DOI: 10.1007/s00894-018-3769-6] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2018] [Accepted: 07/27/2018] [Indexed: 01/04/2023]
|
10
|
Towards Accurate Simulation of Two-Dimensional Electronic Spectroscopy. Top Curr Chem (Cham) 2018; 376:24. [DOI: 10.1007/s41061-018-0201-8] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Accepted: 04/24/2018] [Indexed: 10/14/2022]
|
11
|
Nenov A, Borrego-Varillas R, Oriana A, Ganzer L, Segatta F, Conti I, Segarra-Marti J, Omachi J, Dapor M, Taioli S, Manzoni C, Mukamel S, Cerullo G, Garavelli M. UV-Light-Induced Vibrational Coherences: The Key to Understand Kasha Rule Violation in trans-Azobenzene. J Phys Chem Lett 2018; 9:1534-1541. [PMID: 29504764 DOI: 10.1021/acs.jpclett.8b00152] [Citation(s) in RCA: 79] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
We combine sub-20 fs transient absorption spectroscopy with state-of-the-art computations to study the ultrafast photoinduced dynamics of trans-azobenzene (AB). We are able to resolve the lifetime of the ππ* state, whose decay within ca. 50 fs is correlated to the buildup of the nπ* population and to the emergence of coherences in the dynamics, to date unobserved. Nonlinear spectroscopy simulations call for the CNN in-plane bendings as the active modes in the subps photoinduced coherent dynamics out of the ππ* state. Radiative to kinetic energy transfer into these modes drives the system to a high-energy planar nπ*/ground state conical intersection, inaccessible upon direct excitation of the nπ* state, that triggers an ultrafast (0.45 ps) nonproductive decay of the nπ* state and is thus responsible for the observed Kasha rule violation in UV excited trans-AB. On the other hand, cis-AB is built only after intramolecular vibrational energy redistribution and population of the NN torsional mode.
Collapse
Affiliation(s)
- Artur Nenov
- Dipartimento di Chimica Industriale , Università degli Studi di Bologna , Viale del Risorgimento 4 , I-40136 Bologna , Italy
| | - Rocio Borrego-Varillas
- IFN-CNR, Dipartimento di Fisica , Politecnico di Milano , Piazza Leonardo da Vinci 32 , I-20133 Milano , Italy
| | - Aurelio Oriana
- IFN-CNR, Dipartimento di Fisica , Politecnico di Milano , Piazza Leonardo da Vinci 32 , I-20133 Milano , Italy
| | - Lucia Ganzer
- IFN-CNR, Dipartimento di Fisica , Politecnico di Milano , Piazza Leonardo da Vinci 32 , I-20133 Milano , Italy
| | - Francesco Segatta
- Dipartimento di Chimica Industriale , Università degli Studi di Bologna , Viale del Risorgimento 4 , I-40136 Bologna , Italy
- European Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT*-FBK) and Trento Institute for Fundamental Physics and Applications (TIFPA-INFN) , 38123 Trento , Italy
| | - Irene Conti
- Dipartimento di Chimica Industriale , Università degli Studi di Bologna , Viale del Risorgimento 4 , I-40136 Bologna , Italy
| | - Javier Segarra-Marti
- Laboratoire de Chimie UMR 5182 , Université Lyon, ENS de Lyon, CNRS, Université Lyon 1 , 46 Allée d'Italie , FR-69342 Lyon , France
| | - Junko Omachi
- Institute for Photon Science and Technology , University of Tokyo , 7-3-1 Hongo , Bunkyo-ku, Tokyo 113-0033 , Japan
| | - Maurizio Dapor
- European Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT*-FBK) and Trento Institute for Fundamental Physics and Applications (TIFPA-INFN) , 38123 Trento , Italy
| | - Simone Taioli
- European Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT*-FBK) and Trento Institute for Fundamental Physics and Applications (TIFPA-INFN) , 38123 Trento , Italy
- Faculty of Mathematics and Physics , Charles University , Praha 8 , 180 00 Prague , Czech Republic
| | - Cristian Manzoni
- IFN-CNR, Dipartimento di Fisica , Politecnico di Milano , Piazza Leonardo da Vinci 32 , I-20133 Milano , Italy
| | - Shaul Mukamel
- Department of Chemistry , University of California , Irvine , California 92697-2025 , United States
| | - Giulio Cerullo
- IFN-CNR, Dipartimento di Fisica , Politecnico di Milano , Piazza Leonardo da Vinci 32 , I-20133 Milano , Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale , Università degli Studi di Bologna , Viale del Risorgimento 4 , I-40136 Bologna , Italy
| |
Collapse
|
12
|
The highly excited-state manifold of guanine: calibration for nonlinear electronic spectroscopy simulations. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2225-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
13
|
Padnya PL, Khripunova IA, Mostovaya OA, Mukhametzyanov TA, Evtugyn VG, Vorobev VV, Osin YN, Stoikov II. Self-assembly of chiral fluorescent nanoparticles based on water-soluble L-tryptophan derivatives of p-tert-butylthiacalix[4]arene. BEILSTEIN JOURNAL OF NANOTECHNOLOGY 2017; 8:1825-1835. [PMID: 29046831 PMCID: PMC5629409 DOI: 10.3762/bjnano.8.184] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/02/2017] [Accepted: 08/07/2017] [Indexed: 05/27/2023]
Abstract
New water-soluble tetra-substituted derivatives of p-tert-butylthiacalix[4]arene containing fragments of L-tryptophan in cone and 1,3-alternate conformations were obtained. It was shown that the resulting compounds form stable, positively charged aggregates of 86-134 nm in diameter in water at a concentration of 1 × 10-4 M as confirmed by dynamic light scattering, scanning electron microscopy and transmission electron microscopy. It was established that these aggregates are fluorescently active and chiral. A distinctive feature of the compounds is the pronounced dependence of their spectral (emission and chiroptical) properties on the polarity of the solvent and the length of the linker between the macrocyclic and fluorophore parts of the molecule.
Collapse
Affiliation(s)
- Pavel L Padnya
- Kazan Federal University, 420008 Kremlevskaya, 18, Kazan, Russian Federation
- Peoples Friendship University of Russia (RUDN University), 117198 Miklukho-Maklaya St., 6, Moscow, Russian Federation
| | - Irina A Khripunova
- Kazan Federal University, 420008 Kremlevskaya, 18, Kazan, Russian Federation
| | - Olga A Mostovaya
- Kazan Federal University, 420008 Kremlevskaya, 18, Kazan, Russian Federation
| | | | - Vladimir G Evtugyn
- Kazan Federal University, 420008 Kremlevskaya, 18, Kazan, Russian Federation
| | | | - Yuri N Osin
- Kazan Federal University, 420008 Kremlevskaya, 18, Kazan, Russian Federation
| | - Ivan I Stoikov
- Kazan Federal University, 420008 Kremlevskaya, 18, Kazan, Russian Federation
| |
Collapse
|
14
|
Giussani A, Marcheselli J, Mukamel S, Garavelli M, Nenov A. On the Simulation of Two-dimensional Electronic Spectroscopy of Indole-containing Peptides. Photochem Photobiol 2017; 93:1368-1380. [PMID: 28380692 DOI: 10.1111/php.12770] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2017] [Accepted: 02/27/2017] [Indexed: 01/27/2023]
Abstract
A benchmark study of low-cost multiconfigurational CASSCF/CASPT2 schemes for computing the electronic structure of indole is presented. This facilitates the simulation of near-ultraviolet (UV) pump visible (VIS) probe (i.e. two-color) two-dimensional electronic spectra (2DES) of homo- and hetero-aggregates as well as for processing of multiple snapshots from molecular dynamics simulations. Fingerprint excited-state absorption signatures of indole are identified in a broad spectral window between 10 and 25 k cm-1 . The 18-24 k cm-1 spectral window which has no absorption of the monomer and noninteracting aggregates is ideally suited to embed charge-transfer signatures in stacked aggregates. The small peptide Trp-cage, containing a tryptophan and a tyrosine amino acids, having indole and phenol as side chains, respectively, serves to prove the concept. Clear charge-transfer signatures are found in the proposed spectral window for an interchromophore distance of 5 Å making near-UV pump VIS probe 2DES a suitable technique for resolving closely packed aggregates. We demonstrate that 2DES utilizing ultra-short pulses has the potential to resolve the nature of the spectroscopically resolved electronic states and that the line shapes of the excited-state absorption signals can be correlated to the polarity of the relevant states.
Collapse
Affiliation(s)
- Angelo Giussani
- Dipartimento di Chimica "G. Ciamician", Università degli Studi di Bologna, Bologna, Italy
| | | | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine, CA
| | - Marco Garavelli
- Dipartimento di Chimica "G. Ciamician", Università degli Studi di Bologna, Bologna, Italy.,Dipartimento di Chimica Industriale "Toso Montanari", Universita degli Studi di Bologna, Bologna, Italy
| | - Artur Nenov
- Dipartimento di Chimica "G. Ciamician", Università degli Studi di Bologna, Bologna, Italy.,Dipartimento di Chimica Industriale "Toso Montanari", Universita degli Studi di Bologna, Bologna, Italy
| |
Collapse
|
15
|
Nenov A, Mukamel S, Garavelli M, Rivalta I. Two-Dimensional Electronic Spectroscopy of Benzene, Phenol, and Their Dimer: An Efficient First-Principles Simulation Protocol. J Chem Theory Comput 2016; 11:3755-71. [PMID: 26574458 DOI: 10.1021/acs.jctc.5b00443] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
First-principles simulations of two-dimensional electronic spectroscopy in the ultraviolet region (2DUV) require computationally demanding multiconfigurational approaches that can resolve doubly excited and charge transfer states, the spectroscopic fingerprints of coupled UV-active chromophores. Here, we propose an efficient approach to reduce the computational cost of accurate simulations of 2DUV spectra of benzene, phenol, and their dimer (i.e., the minimal models for studying electronic coupling of UV-chromophores in proteins). We first establish the multiconfigurational recipe with the highest accuracy by comparison with experimental data, providing reference gas-phase transition energies and dipole moments that can be used to construct exciton Hamiltonians involving high-lying excited states. We show that by reducing the active spaces and the number of configuration state functions within restricted active space schemes, the computational cost can be significantly decreased without loss of accuracy in predicting 2DUV spectra. The proposed recipe has been successfully tested on a realistic model proteic system in water. Accounting for line broadening due to thermal and solvent-induced fluctuations allows for direct comparison with experiments.
Collapse
Affiliation(s)
- Artur Nenov
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via F. Selmi 2, 40126 Bologna, Italy
| | - Shaul Mukamel
- Department of Chemistry, University of California , Irvine, California 92697-2025, United States
| | - Marco Garavelli
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via F. Selmi 2, 40126 Bologna, Italy.,Université de Lyon , CNRS, Laboratoire de Chimie, École Normale Supérieure de Lyon, UMR 5182, 46 Allée d'Italie, 69364 Lyon, Cedex 07, France
| | - Ivan Rivalta
- Université de Lyon , CNRS, Laboratoire de Chimie, École Normale Supérieure de Lyon, UMR 5182, 46 Allée d'Italie, 69364 Lyon, Cedex 07, France
| |
Collapse
|
16
|
Giussani A, Segarra-Martí J, Nenov A, Rivalta I, Tolomelli A, Mukamel S, Garavelli M. Spectroscopic fingerprints of DNA/RNA pyrimidine nucleobases in third-order nonlinear electronic spectra. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1867-z] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
17
|
Nenov A, Giussani A, Fingerhut BP, Rivalta I, Dumont E, Mukamel S, Garavelli M. Spectral lineshapes in nonlinear electronic spectroscopy. Phys Chem Chem Phys 2016; 17:30925-36. [PMID: 26084213 DOI: 10.1039/c5cp01167a] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We outline a computational approach for nonlinear electronic spectra, which accounts for the electronic energy fluctuations due to nuclear degrees of freedom and explicitly incorporates the fluctuations of higher excited states, induced by the dynamics in the photoactive state(s). This approach is based on mixed quantum-classical dynamics simulations. Tedious averaging over multiple trajectories is avoided by employing the linearly displaced Brownian harmonic oscillator to model the correlation functions. The present strategy couples accurate computations of the high-lying excited state manifold with dynamics simulations. The application is made to the two-dimensional electronic spectra of pyrene, a polycyclic aromatic hydrocarbon characterized by an ultrafast (few tens of femtoseconds) decay from the bright S2 state to the dark S1 state. The spectra for waiting times t2 = 0 and t2 = 1 ps demonstrate the ability of this approach to model electronic state fluctuations and realistic lineshapes. Comparison with experimental spectra [Krebs et al., New Journal of Physics, 2013, 15, 085016] shows excellent agreement and allows us to unambiguously assign the excited state absorption features.
Collapse
Affiliation(s)
- Artur Nenov
- Dipartimento di Chimica G. Ciamician, Università di Bologna, V. F. Selmi 2, 40126 Bologna, Italy.
| | - Angelo Giussani
- Dipartimento di Chimica G. Ciamician, Università di Bologna, V. F. Selmi 2, 40126 Bologna, Italy.
| | - Benjamin P Fingerhut
- Max Born Institute for Nonlinear Optics and Short Pulse Spectroscopy, Max-Born-Strasse 2A, Berlin, 12489, Germany
| | - Ivan Rivalta
- Laboratoire de Chimie, Ecole Normale Suprieure de Lyon, 46, allée d'Italie, 69364 Lyon, France
| | - Elise Dumont
- Laboratoire de Chimie, Ecole Normale Suprieure de Lyon, 46, allée d'Italie, 69364 Lyon, France
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine, California 92697-2025, USA
| | - Marco Garavelli
- Dipartimento di Chimica G. Ciamician, Università di Bologna, V. F. Selmi 2, 40126 Bologna, Italy. and Laboratoire de Chimie, Ecole Normale Suprieure de Lyon, 46, allée d'Italie, 69364 Lyon, France
| |
Collapse
|
18
|
Segarra-Martí J, Garavelli M, Aquilante F. Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems. J Chem Theory Comput 2016; 11:3772-84. [PMID: 26574459 DOI: 10.1021/acs.jctc.5b00479] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A new flavor of the frozen natural orbital complete active space second-order perturbation theory method (FNO-CASPT2, Aquilante et al., J. Chem. Phys. 131, 034113) is proposed herein. In this new implementation, the virtual space in Cholesky decomposition-based CASPT2 computations (CD-CASPT2) is truncated by excluding those orbitals that contribute the least toward preserving a predefined value of the trace of an approximate density matrix, as that represents a measure of the amount of dynamic correlation retained in the model. In this way, the amount of correlation included is practically constant at all nuclear arrangements, thus allowing for the computation of smooth electronic states surfaces and energy gradients-essential requirements for theoretical studies in photochemistry. The method has been benchmarked for a series of relevant biochromophores for which large speed-ups have been recorded while retaining the accuracy achieved in the corresponding CD-CASPT2 calculations. Both vertical excitation energies and gradient calculations have been carried out to establish general guidelines as to how much correlation needs to be retained in the calculation for the results to be consistent with the CD-CASPT2 findings. Our results feature errors within a tenth of an eV for the most difficult cases and have been validated to be used for gradient computations where an up to 3-fold speed-up is observed depending on the size of the system and the basis set employed.
Collapse
Affiliation(s)
- Javier Segarra-Martí
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via Selmi 2, IT-40126 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via Selmi 2, IT-40126 Bologna, Italy.,Université de Lyon, CNRS , Institut de Chimie de Lyon, École Normale Supérieure de Lyon, 46 Allée d'Italie, F-69364 Lyon Cedex 07, France
| | - Francesco Aquilante
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via Selmi 2, IT-40126 Bologna, Italy
| |
Collapse
|
19
|
Nenov A, Segarra-Martí J, Giussani A, Conti I, Rivalta I, Dumont E, Jaiswal VK, Altavilla SF, Mukamel S, Garavelli M. Probing deactivation pathways of DNA nucleobases by two-dimensional electronic spectroscopy: first principles simulations. Faraday Discuss 2015; 177:345-62. [PMID: 25607949 DOI: 10.1039/c4fd00175c] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The SOS//QM/MM [Rivalta et al., Int. J. Quant. Chem., 2014, 114, 85] method consists of an arsenal of computational tools allowing accurate simulation of one-dimensional (1D) and bi-dimensional (2D) electronic spectra of monomeric and dimeric systems with unprecedented details and accuracy. Prominent features like doubly excited local and excimer states, accessible in multi-photon processes, as well as charge-transfer states arise naturally through the fully quantum-mechanical description of the aggregates. In this contribution the SOS//QM/MM approach is extended to simulate time-resolved 2D spectra that can be used to characterize ultrafast excited state relaxation dynamics with atomistic details. We demonstrate how critical structures on the excited state potential energy surface, obtained through state-of-the-art quantum chemical computations, can be used as snapshots of the excited state relaxation dynamics to generate spectral fingerprints for different de-excitation channels. The approach is based on high-level multi-configurational wavefunction methods combined with non-linear response theory and incorporates the effects of the solvent/environment through hybrid quantum mechanics/molecular mechanics (QM/MM) techniques. Specifically, the protocol makes use of the second-order Perturbation Theory (CASPT2) on top of Complete Active Space Self Consistent Field (CASSCF) strategy to compute the high-lying excited states that can be accessed in different 2D experimental setups. As an example, the photophysics of the stacked adenine-adenine dimer in a double-stranded DNA is modeled through 2D near-ultraviolet (NUV) spectroscopy.
Collapse
Affiliation(s)
- Artur Nenov
- Dipartimento di Chimica G. Ciamician, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
20
|
Nenov A, Giussani A, Segarra-Martí J, Jaiswal VK, Rivalta I, Cerullo G, Mukamel S, Garavelli M. Modeling the high-energy electronic state manifold of adenine: Calibration for nonlinear electronic spectroscopy. J Chem Phys 2015; 142:212443. [PMID: 26049463 DOI: 10.1063/1.4921016] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Pump-probe electronic spectroscopy using femtosecond laser pulses has evolved into a standard tool for tracking ultrafast excited state dynamics. Its two-dimensional (2D) counterpart is becoming an increasingly available and promising technique for resolving many of the limitations of pump-probe caused by spectral congestion. The ability to simulate pump-probe and 2D spectra from ab initio computations would allow one to link mechanistic observables like molecular motions and the making/breaking of chemical bonds to experimental observables like excited state lifetimes and quantum yields. From a theoretical standpoint, the characterization of the electronic transitions in the visible (Vis)/ultraviolet (UV), which are excited via the interaction of a molecular system with the incoming pump/probe pulses, translates into the determination of a computationally challenging number of excited states (going over 100) even for small/medium sized systems. A protocol is therefore required to evaluate the fluctuations of spectral properties like transition energies and dipole moments as a function of the computational parameters and to estimate the effect of these fluctuations on the transient spectral appearance. In the present contribution such a protocol is presented within the framework of complete and restricted active space self-consistent field theory and its second-order perturbation theory extensions. The electronic excited states of adenine have been carefully characterized through a previously presented computational recipe [Nenov et al., Comput. Theor. Chem. 1040-1041, 295-303 (2014)]. A wise reduction of the level of theory has then been performed in order to obtain a computationally less demanding approach that is still able to reproduce the characteristic features of the reference data. Foreseeing the potentiality of 2D electronic spectroscopy to track polynucleotide ground and excited state dynamics, and in particular its expected ability to provide conformational dependent fingerprints in dimeric systems, the performances of the selected reduced level of calculations have been tested in the construction of 2D electronic spectra for the in vacuo adenine monomer and the unstacked adenine homodimer, thereby exciting the Lb/La transitions with the pump pulse pair and probing in the Vis to near ultraviolet spectral window.
Collapse
Affiliation(s)
- Artur Nenov
- Dipartimento di Chimica "G. Ciamician," Università di Bologna, Via Selmi 2, IT-40126 Bologna, Italy
| | - Angelo Giussani
- Dipartimento di Chimica "G. Ciamician," Università di Bologna, Via Selmi 2, IT-40126 Bologna, Italy
| | - Javier Segarra-Martí
- Dipartimento di Chimica "G. Ciamician," Università di Bologna, Via Selmi 2, IT-40126 Bologna, Italy
| | - Vishal K Jaiswal
- Dipartimento di Chimica "G. Ciamician," Università di Bologna, Via Selmi 2, IT-40126 Bologna, Italy
| | - Ivan Rivalta
- Université de Lyon, CNRS, Institut de Chimie de Lyon, École Normale Supérieure de Lyon, 46 Allée d'Italie, F-69364 Lyon Cedex 07, France
| | - Giulio Cerullo
- Dipartimento di Fisica, Politecnico di Milano, IFN-CNR, Piazza Leonardo Da Vinci 32, IT-20133 Milano, Italy
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine, California 92697-2025, USA
| | - Marco Garavelli
- Dipartimento di Chimica "G. Ciamician," Università di Bologna, Via Selmi 2, IT-40126 Bologna, Italy
| |
Collapse
|