1
|
Bao S, Raymond N, Zeng T, Nooijen M. Vibrational Electronic-Thermofield Coupled Cluster (VE-TFCC) Theory for Quantum Simulations of Vibronic Coupling Systems at Thermal Equilibrium. J Chem Theory Comput 2024. [PMID: 38950345 DOI: 10.1021/acs.jctc.4c00338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/03/2024]
Abstract
A vibrational electronic-thermofield coupled cluster (VE-TFCC) approach is developed to calculate thermal properties of non-adiabatic vibronic coupling systems. The thermofield (TF) theory and a mixed linear exponential ansatz based on second-quantized Bosonic construction operators are introduced to propagate the thermal density operator as a "pure state" in the Bogoliubov representation. Through this compact representation of the thermal density operator, the approach is basis-set-free and scales classically (polynomial) as the number of degrees of freedoms (DoF) in the system increases. The VE-TFCC approach is benchmarked with small test models and a real molecular compound (CoF4- anion) against the conventional sum over states (SOS) method and applied to calculate thermochemistry properties of a gas-phase reaction: CoF3 + F- → CoF4-. Results shows that the VE-TFCC approach, in conjunction with vibronic models, provides an effective protocol for calculating thermodynamic properties of vibronic coupling systems.
Collapse
Affiliation(s)
- Songhao Bao
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Neil Raymond
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Tao Zeng
- Department of Chemistry, York University, Toronto, Ontario M3J 1P3, Canada
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
2
|
García-González F, Otero JC, Ávila Ferrer FJ, Santoro F, Aranda D. Linear Vibronic Coupling Approach for Surface-Enhanced Raman Scattering: Quantifying the Charge-Transfer Enhancement Mechanism. J Chem Theory Comput 2024; 20:3850-3863. [PMID: 38687961 PMCID: PMC11099975 DOI: 10.1021/acs.jctc.4c00061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Revised: 04/12/2024] [Accepted: 04/15/2024] [Indexed: 05/02/2024]
Abstract
The outstanding amplification observed in surface-enhanced Raman scattering (SERS) is due to several enhancement mechanisms, and standing out among them are the plasmonic (PL) and charge-transfer (CT) mechanisms. The theoretical estimation of the enhancement factors of the CT mechanism is challenging because the excited-state coupling between bright plasmons and dark CT states must be properly introduced into the model to obtain reliable intensities. In this work, we aim at simulating electrochemical SERS spectra, considering models of pyridine on silver clusters subjected to an external electric field E⃗ that represents the effect of an electrode potential Vel. The method adopts quantum dynamical propagations of nuclear wavepackets on the coupled PL and CT states described with linear vibronic coupling models parametrized for each E⃗ through a fragment-based maximum-overlap diabatization. By presenting results at different values of E⃗, we show that indeed there is a relation between the population transferred to the CT states and the total scattered intensity. The tuning and detuning processes of the CT states with the bright PLs as a function of the electric field are in good agreement with those observed in experiments. Finally, our estimations for the CT enhancement factors predict values in the order of 105 to 106, meaning that when the CT and PL states are both in resonance with the excitation wavelength, the CT and PL enhancements are comparable, and vibrational bands whose intensity is amplified by different mechanisms can be observed together, in agreement with what was measured by typical experiments on silver electrodes.
Collapse
Affiliation(s)
- Francisco García-González
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Juan Carlos Otero
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Francisco J. Ávila Ferrer
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Fabrizio Santoro
- Istituto
di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca
del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Daniel Aranda
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| |
Collapse
|
3
|
Cerezo J, Gierschner J, Santoro F, Prampolini G. Explicit Modelling of Spectral Bandshapes by a Mixed Quantum-Classical Approach: Solvent Order and Temperature Effects in the Optical Spectra of Distyrylbenzene. Chemphyschem 2024:e202400307. [PMID: 38728539 DOI: 10.1002/cphc.202400307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Revised: 04/29/2024] [Accepted: 05/08/2024] [Indexed: 05/12/2024]
Abstract
The absorption and emission spectral shapes of a flexible organic probe, the distyrylbenzene (DSB) dye, are simulated accounting for the effect of different environments of increasing complexity, ranging from a homogeneous, low-molecular- weight solvent, to a long-chain alkane, and, eventually, a channel-forming organic matrix. Each embedding is treated explicitly, adopting a mixed quantum-classical approach, the Adiabatic Molecular Dynamics - generalized vertical Hessian (Ad-MD|gVH) model, which allows a direct simulation of the environment-induced constraining effects on the vibronic spectral shapes. In such a theoretical framework, the stiff modes of the dye are described at a quantum level within the harmonic approximation, including Duschinsky mixing effects, while flexible degrees of freedom of the solute (e. g. torsions) and those of the solvent are treated classically by means of molecular dynamics sampling. Such a setup is shown to reproduce the distinct effects exerted by the different environments in varied thermodynamic conditions. Besides allowing for a first-principles rationale on the supramolecular mechanism leading to the experimental spectral features, this result represents the first successful application of the Ad-MD|gVH method to complex embeddings and supports its potential application to other heterogeneous environments, such as for instance, pigment-protein complexes or organic dyes adsorbed into metal-organic frameworks.
Collapse
Affiliation(s)
- Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049, Madrid, Spain
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM), Consiglio Nazionale delle Ricerche (CNR), 50019, Sesto Fiorentino, Italien
| | - Johannes Gierschner
- Madrid Institute for Advanced Studies, IMDEA Nanoscience, C/Faraday 9, Ciudad Universitaria de Cantoblanco, 28049, Madrid, Spain
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM), Consiglio Nazionale delle Ricerche (CNR), 50019, Sesto Fiorentino, Italien
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM), Consiglio Nazionale delle Ricerche (CNR), 50019, Sesto Fiorentino, Italien
| |
Collapse
|
4
|
Hunter KE, Mao Y, Chin AW, Zuehlsdorff TJ. Environmentally Driven Symmetry Breaking Quenches Dual Fluorescence in Proflavine. J Phys Chem Lett 2024; 15:4623-4632. [PMID: 38647005 DOI: 10.1021/acs.jpclett.4c00677] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
Nonadiabatic couplings between several electronic excited states are ubiquitous in many organic chromophores and can significantly influence optical properties. A recent experimental study demonstrated that the proflavine molecule exhibits surprising dual fluorescence in the gas phase, which is suppressed in polar solvent environments. Here, we uncover the origin of this phenomenon by parametrizing a linear-vibronic coupling Hamiltonian from spectral densities of system-bath coupling constructed along molecular dynamics trajectories, fully accounting for interactions with the condensed-phase environment. The finite-temperature absorption, steady-state emission, and time-resolved emission spectra are then computed using powerful, numerically exact tensor network approaches. We find that the dual fluorescence in vacuum is driven by a single well-defined coupling mode but is quenched in solution due to dynamic solvent-driven symmetry breaking that mixes the two low-lying electronic states. We expect the computational framework developed here to be widely applicable to the study of non-Condon effects in complex condensed-phase environments.
Collapse
Affiliation(s)
- Kye E Hunter
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, Paris 75005, France
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| |
Collapse
|
5
|
Bao S, Raymond N, Nooijen M. Time dependent vibrational electronic coupled cluster (VECC) theory for non-adiabatic nuclear dynamics. J Chem Phys 2024; 160:094105. [PMID: 38426527 DOI: 10.1063/5.0190034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2023] [Accepted: 01/28/2024] [Indexed: 03/02/2024] Open
Abstract
A time-dependent vibrational electronic coupled-cluster (VECC) approach is proposed to simulate photo-electron/UV-VIS absorption spectra as well as time-dependent properties for non-adiabatic vibronic models, going beyond the Born-Oppenheimer approximation. A detailed derivation of the equations of motion and a motivation for the ansatz are presented. The VECC method employs second-quantized bosonic construction operators and a mixed linear and exponential ansatz to form a compact representation of the time-dependent wave-function. Importantly, the method does not require a basis set, has only a few user-defined inputs, and has a classical (polynomial) scaling with respect to the number of degrees of freedom (of the vibronic model), resulting in a favorable computational cost. In benchmark applications to small models and molecules, the VECC method provides accurate results compared to multi-configurational time-dependent Hartree calculations when predicting short-time dynamical properties (i.e., photo-electron/UV-VIS absorption spectra) for non-adiabatic vibronic models. To illustrate the capabilities, the VECC method is also successfully applied to a large vibronic model for hexahelicene with 14 electronic states and 63 normal modes, developed in the group by Aranda and Santoro [J. Chem. Theory Comput. 17, 1691, (2021)].
Collapse
Affiliation(s)
- Songhao Bao
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Neil Raymond
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
6
|
Wiethorn ZR, Hunter KE, Zuehlsdorff TJ, Montoya-Castillo A. Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations. J Chem Phys 2023; 159:244114. [PMID: 38153146 DOI: 10.1063/5.0180405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 12/06/2023] [Indexed: 12/29/2023] Open
Abstract
While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute-solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.
Collapse
Affiliation(s)
- Zachary R Wiethorn
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Kye E Hunter
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | | |
Collapse
|
7
|
Bousquet MHE, Papineau TV, Veys K, Escudero D, Jacquemin D. Extensive Analysis of the Parameters Influencing Radiative Rates Obtained through Vibronic Calculations. J Chem Theory Comput 2023; 19:5525-5547. [PMID: 37494031 DOI: 10.1021/acs.jctc.3c00191] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2023]
Abstract
Defining a theoretical model systematically delivering accurate ab initio predictions of the fluorescence quantum yields of organic dyes is highly desirable for designing improved fluorophores in a systematic rather than trial-and-error way. To this end, the first required step is to obtain reliable radiative rates (kr), as low kr typically precludes effective emission. In the present contribution, using a series of 10 substituted phenyls with known experimental kr, we analyze the impact of the computational protocol on the kr determined through the thermal vibration correlation function (TVCF) approach on the basis of time-dependent density functional theory (TD-DFT) calculations of the energies, structures, and vibrational parameters. Both the electronic structure (selected exchange-correlation functional, application or not of the Tamm-Dancoff approximation) and the vibronic parameters (line-shape formalism, coordinate system, potential energy surface model, and dipole expansion) are tackled. Considering all possible combinations yields more than 3500 cases, allowing to extract statistically-relevant information regarding the impact of each computational parameter on the magnitude of the estimated kr. It turns out that the selected vibronic model can have a significant impact on the computed kr, especially the potential energy surface model. This effect is of the same order of magnitude as the difference noted between B3LYP and CAM-B3LYP estimates. For the treated compounds, all evaluated functionals do deliver reasonable trends, fitting the experimental values.
Collapse
Affiliation(s)
| | | | - Koen Veys
- Department of Chemistry, KU Leuven, B-3001 Leuven, Belgium
| | | | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| |
Collapse
|
8
|
Segatta F, Ruiz DA, Aleotti F, Yaghoubi M, Mukamel S, Garavelli M, Santoro F, Nenov A. Nonlinear Molecular Electronic Spectroscopy via MCTDH Quantum Dynamics: From Exact to Approximate Expressions. J Chem Theory Comput 2023; 19:2075-2091. [PMID: 36961952 PMCID: PMC10100531 DOI: 10.1021/acs.jctc.2c01059] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/26/2023]
Abstract
We present an accurate and efficient approach to computing the linear and nonlinear optical spectroscopy of a closed quantum system subject to impulsive interactions with an incident electromagnetic field. It incorporates the effect of ultrafast nonadiabatic dynamics by means of explicit numerical propagation of the nuclear wave packet. The fundamental expressions for the evaluation of first- and higher-order response functions are recast in a general form that can be used with any quantum dynamics code capable of computing the overlap of nuclear wave packets evolving in different states. Here we present the evaluation of these expressions with the multiconfiguration time-dependent Hartree (MCTDH) method. Application is made to pyrene, excited to its lowest bright excited state S2 which exhibits a sub-100-fs nonadiabatic decay to a dark state S1. The system is described by a linear vibronic coupling Hamiltonian, parametrized with multiconfiguration electronic structure methods. We show that the ultrafast nonadiabatic dynamics can have a remarkable effect on the spectral line shapes that goes beyond simple lifetime broadening. Furthermore, a widely employed approximate expression based on the time scale separation of dephasing and population relaxation is recast in the same theoretical framework. Application to pyrene shows the range of validity of such approximations.
Collapse
Affiliation(s)
- Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Daniel Aranda Ruiz
- ICMol, Universidad de Valencia, Catedrático José Beltrán Martínez, 2, 46980 Paterna, Spain
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Flavia Aleotti
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Martha Yaghoubi
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, United States
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", University of 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
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| |
Collapse
|
9
|
Penfold TJ, Eng J. Mind the GAP: quantifying the breakdown of the linear vibronic coupling Hamiltonian. Phys Chem Chem Phys 2023; 25:7195-7204. [PMID: 36820783 DOI: 10.1039/d2cp05576g] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
Excited state dynamics play a critical role across a broad range of scientific fields. Importantly, the highly non-equilibrium nature of the states generated by photoexcitation means that excited state simulations should usually include an accurate description of the coupled electronic-nuclear motion, which often requires solving the time-dependent Schrödinger equation (TDSE). One of the biggest challenges for these simulations is the requirement to calculate the PES over which the nuclei evolve. An effective approach for addressing this challenge is to use the approximate linear vibronic coupling (LVC) Hamiltonian, which enables a model potential to be parameterised using relatively few quantum chemistry calculations. However, this approach is only valid provided there are no large amplitude motions in the excited state dynamics. In this paper we introduce and deploy a metric, the global anharmonicity parameter (GAP), which can be used to assess the accuracy of an LVC potential. Following its derivation, we illustrate its utility by applying it to three molecules exhibiting different rigidity in their excited states.
Collapse
Affiliation(s)
- Thomas J Penfold
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK.
| | - Julien Eng
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK.
| |
Collapse
|
10
|
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
|
11
|
Cerezo J, García-Iriepa C, Santoro F, Navizet I, Prampolini G. Unraveling the contributions to the spectral shape of flexible dyes in solution: insights on the absorption spectrum of an oxyluciferin analogue. Phys Chem Chem Phys 2023; 25:5007-5020. [PMID: 36722876 DOI: 10.1039/d2cp05701h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
We present a computational investigation of the absorption spectrum in water of 5,5-spirocyclopropyl-oxyluciferin (5,5-CprOxyLH), an analogue of the emitter compound responsible for the bioluminescence in fireflies. Several factors participate in determining the 5,5-CprOxyLH's spectral shape: (i) the contribution of the four close-energy excited states, which show significant non-adiabatic couplings, (ii) the flexible molecular structure and (iii) the specific interactions established with the surrounding environment, which strongly couple the protic solvent dynamics with the dye's spectral response. To tackle the challenge to capture and dissect the role of all these effects we preliminarily investigate the role of non-adiabatic couplings with quantum dynamics simulations and a linear vibronic coupling model in the gas phase. Then, we account for both the molecular flexibility and solvent interactions by resorting to a mixed quantum classical protocol, named Adiabatic Molecular Dynamics generalized Vertical Gradient (Ad-MD|gVG), which is built on a method recently proposed by some of us. It is rooted in the partition between stiff degrees of freedom of the dye, accounted for at the vibronic level within the harmonic approximation, and flexible degrees of freedom of the solute (and of the solvent), described classically through a sampling based on Molecular Dynamics (MD). Ad-MD|gVG avoids spurious effects arising in the excited state Hessians due to non-adiabatic couplings, and can therefore be applied to account for the contributions of the first four excited states to the 5,5-CprOxyLH absorption spectrum. The final simulated spectrum is in very good agreement with the experiment, especially when the MD is driven by a refined quantum-mechanically derived force-field. More importantly, the origin of each separate contribution to the spectral shape is appropriately accounted for, paving the way to future applications of the method to more complex systems or alternative spectroscopies, as emission or circular dichroism.
Collapse
Affiliation(s)
- Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain. .,CNR - Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - Cristina García-Iriepa
- Universidad de Alcalá, Departamento de Química Analítica, Química Física e Ingeniería Química, Grupo de Reactividad y Estructura Molecular (RESMOL), 28806 Alcalá de Henares (Madrid), Spain. .,Universidad de Alcalá, Instituto de Investigación Química "Andrés M. del Río" (IQAR), 28806 Alcalá de Henares (Madrid), Spain
| | - Fabrizio Santoro
- CNR - Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - Isabelle Navizet
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Giacomo Prampolini
- CNR - Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy.
| |
Collapse
|
12
|
Green JA, Gómez S, Worth G, Santoro F, Improta R. Solvent Effects on Ultrafast Charge Transfer Population: Insights from the Quantum Dynamics of Guanine-Cytosine in Chloroform. Chemistry 2022; 28:e202201731. [PMID: 35950519 PMCID: PMC9828530 DOI: 10.1002/chem.202201731] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Indexed: 01/12/2023]
Abstract
We study the ultrafast photoactivated dynamics of the hydrogen bonded dimer Guanine-Cytosine in chloroform solution, focusing on the population of the Guanine→Cytosine charge transfer state (GC-CT), an important elementary process for the photophysics and photochemistry of nucleic acids. We integrate a quantum dynamics propagation scheme, based on a linear vibronic model parameterized through time dependent density functional theory calculations, with four different solvation models, either implicit or explicit. On average, after 50 fs, 30∼40 % of the bright excited state population has been transferred to GC-CT. This process is thus fast and effective, especially when transferring from the Guanine bright excited states, in line with the available experimental studies. Independent of the adopted solvation model, the population of GC-CT is however disfavoured in solution with respect to the gas phase. We show that dynamical solvation effects are responsible for this puzzling result and assess the different chemical-physical effects modulating the population of CT states on the ultrafast time-scale. We also propose some simple analyses to predict how solvent can affect the population transfer between bright and CT states, showing that the effect of the solute/solvent electrostatic interactions on the energy of the CT state can provide a rather reliable indication of its possible population.
Collapse
Affiliation(s)
- James A. Green
- Istituto di Biostrutture e Bioimmagini-CNRVia De Amicis 95I-80145Napoli
| | - Sandra Gómez
- Department of ChemistryUniversity College London20 Gordon StreetLondonWC1H 0AJUnited Kingdom,Departamento de Química FísicaUniversity of SalamancaSalamanca37008Spain
| | - Graham Worth
- Department of ChemistryUniversity College London20 Gordon StreetLondonWC1H 0AJUnited Kingdom
| | - Fabrizio Santoro
- Istituto di Chimica die Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNRVia Moruzzi 1I-56124Pisa
| | - Roberto Improta
- Istituto di Biostrutture e Bioimmagini-CNRVia De Amicis 95I-80145Napoli
| |
Collapse
|
13
|
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
|
14
|
Ozcelik A, Aranda D, Pereira-Cameselle R, Talavera M, Covelo B, Santoro F, Peña-Gallego Á, Alonso-Gómez JL. ON/OFF Spiroconjugation through Peripheral Functionalization: Impact on the Reactivity and Chiroptical Properties of Spirobifluorenes. Chempluschem 2022; 87:e202100554. [PMID: 35415974 DOI: 10.1002/cplu.202100554] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Revised: 03/16/2022] [Indexed: 02/03/2023]
Abstract
Spirobifluorenes are an important class of spiro compounds frequently used in the field of organic electronics. However, harnessing spiroconjugation to obtain high-performance in such structural motifs remains unexplored. We herein propose that peripheral functionalization may serve as a useful tool to control spiroconjugation in an ON/OFF manner on both chemical reactivity and photophysical properties. In particular, the ratio of mono- and di-functionalized spirobifluorenes found experimentally during their synthesis were found to be 3/2, 7/2, and 12/2 for phenyl, nitro-phenyl and amino-phenyl analogues, respectively. These remarkable reactivity differences correlate with the spiroconjugation character evaluated theoretically at the CAM-B3LYP/6-31G(d,p) level of theory. Additionally, comparison of experimental and predicted optical and chiroptical responses shows that spiroconjugated molecular orbitals have a significant or negligible involvement on the main electronic transitions depending on the peripheral functionality of the spirobifluorene.
Collapse
Affiliation(s)
- Ani Ozcelik
- Departamento de Química Orgánica, Universidad de Vigo, Campus Universitario, 36310, Vigo, Spain
| | - Daniel Aranda
- Istituto di Chimica dei Composti Organometallici, Consiglio Nazionale delle Ricerche, Via Madonna del Piano 10, 50019, Sesto Fiorentino, Pisa, Italy
| | | | - María Talavera
- Department of Chemistry, Humboldt Universität zu Berlin, Brook-Taylor-Strasse 2, 12489, Berlin, Germany
| | - Berta Covelo
- CACTI (Centro de Apoyo Científico-Tecnológico a la Investigación), Universidad de Vigo, Campus Universitario, 36310, Vigo, Spain
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organometallici, Consiglio Nazionale delle Ricerche, Via Madonna del Piano 10, 50019, Sesto Fiorentino, Pisa, Italy
| | - Ángeles Peña-Gallego
- Departamento de Química Física, Universidad de Vigo, Campus Universitario, 36310, Vigo, Spain
| | - J Lorenzo Alonso-Gómez
- Departamento de Química Orgánica, Universidad de Vigo, Campus Universitario, 36310, Vigo, Spain
| |
Collapse
|
15
|
Segalina A, Aranda D, Green JA, Cristino V, Caramori S, Prampolini G, Pastore M, Santoro F. How the Interplay among Conformational Disorder, Solvation, Local, and Charge-Transfer Excitations Affects the Absorption Spectrum and Photoinduced Dynamics of Perylene Diimide Dimers: A Molecular Dynamics/Quantum Vibronic Approach. J Chem Theory Comput 2022; 18:3718-3736. [PMID: 35377648 PMCID: PMC9202308 DOI: 10.1021/acs.jctc.2c00063] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
![]()
In this contribution
we present a mixed quantum-classical dynamical
approach for the computation of vibronic absorption spectra of molecular
aggregates and their nonadiabatic dynamics, taking into account the
coupling between local excitations (LE) and charge-transfer (CT) states.
The approach is based on an adiabatic (Ad) separation between the
soft degrees of freedom (DoFs) of the system and the stiff vibrations,
which are described by the quantum dynamics (QD) of wave packets (WPs)
moving on the coupled potential energy surfaces (PESs) of the LE and
CT states. These PESs are described with a linear vibronic coupling
(LVC) Hamiltonian, parameterized by an overlap-based diabatization
on the grounds of time-dependent density functional theory computations.
The WPs time evolution is computed with the multiconfiguration time-dependent
Hartree method, using effective modes defined through a hierarchical
representation of the LVC Hamiltonian. The soft DoFs are sampled with
classical molecular dynamics (MD), and the coupling between the slow
and fast DoFs is included by recomputing the key parameters of the
LVC Hamiltonians, specifically for each MD configuration. This method,
named Ad-MD|gLVC, is applied to a perylene diimide (PDI) dimer in
acetonitrile and water solutions, and it is shown to accurately reproduce
the change in the vibronic features of the absorption spectrum upon
aggregation. Moreover, the microscopic insight offered by the MD trajectories
allows for a detailed understanding of the role played by the fluctuation
of the aggregate structure on the shape of the vibronic spectrum and
on the population of LE and CT states. The nonadiabatic QD predicts
an extremely fast (∼50 fs) energy transfer between the two
LEs. CT states have only a moderate effect on the absorption spectrum,
despite the fact that after photoexcitation they are shown to acquire
a fast and non-negligible population, highlighting their relevance
in dictating the charge separation and transport in PDI-based optical
devices.
Collapse
Affiliation(s)
- Alekos Segalina
- Université de Lorraine and CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Daniel Aranda
- Instituto de Ciencia Molecular (ICMol), Universidad de Valencia, Catedrático J. Beltrán 2, 46980 Paterna, Valencia, Spain
| | - James A Green
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Vito Cristino
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie, Via Fossato di Mortara 17, 44121 Ferrara, Italy
| | - Stefano Caramori
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie, Via Fossato di Mortara 17, 44121 Ferrara, Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | | | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
16
|
Tissier RC, Rigaud B, Thureau P, Huix-Rotllant M, Jaber M, Ferré N. Stressing the differences in alizarin and purpurin dyes through UV-visible light absorption and 1H-NMR spectroscopies. Phys Chem Chem Phys 2022; 24:19452-19462. [DOI: 10.1039/d2cp00520d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Three anthraquinone-based chromophores (9,10-anthraquinone, alizarin, purpurin) are compared from the point of view of their experimental and computed NMR and UV-visible light absorption spectra. Using an hybrid (explicit/implicit) solvent model,...
Collapse
|
17
|
Vibrationally resolved absorption and emission spectral shapes of one 5-carbohelicene derivative: A theoretical study. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.139067] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
18
|
Dunnett AJ, Gowland D, Isborn CM, Chin AW, Zuehlsdorff TJ. Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue. J Chem Phys 2021; 155:144112. [PMID: 34654312 DOI: 10.1063/5.0062950] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute-solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute-solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.
Collapse
Affiliation(s)
- Angus J Dunnett
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Duncan Gowland
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
19
|
Jouybari MY, Green JA, Improta R, Santoro F. The Ultrafast Quantum Dynamics of Photoexcited Adenine-Thymine Basepair Investigated with a Fragment-based Diabatization and a Linear Vibronic Coupling Model. J Phys Chem A 2021; 125:8912-8924. [PMID: 34609880 PMCID: PMC9281421 DOI: 10.1021/acs.jpca.1c08132] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
In
this contribution we present a quantum dynamical study of the
photoexcited hydrogen bonded base pair adenine–thymine (AT)
in a Watson–Crick arrangement. To that end, we parametrize
Linear Vibronic Coupling (LVC) models with Time-Dependent Density
Functional Theory (TD-DFT) calculations, exploiting a fragment diabatization
scheme (FrD) we have developed to define diabatic states on the basis
of individual chromophores in a multichromophoric system. Wavepacket
propagations were run with the multilayer extension of the Multiconfiguration
Time-Dependent Hartree method. We considered excitations to the three
lowest bright states, a ππ* state of
thymine and two ππ* states (La and Lb) of adenine, and we found that on the 100 fs time
scale the main decay pathways involve intramonomer population transfers
toward nπ* states of the same nucleobase. In AT this transfer
is less effective than in the isolated nucleobases, because hydrogen
bonding destabilizes the nπ* states. The population transfer
to the A → T charge transfer state is negligible, making the
ultrafast (femtosecond) decay through the proton coupled electron
transfer mechanism unlikely, in line with experimental results in
apolar solvents. The excitation energy transfer is also very small.
We carefully compare the predictions of LVC Hamiltonians obtained
with different sets of diabatic states, defined so to match either
local states of the two separated monomers or the base pair adiabatic
states in the Franck–Condon region. To that end we also extend
the flexibility of the FrD-LVC approach, introducing a new strategy
to define fragments diabatic states that account for the effect of
the rest of the multichromohoric system through a Molecular Mechanics
potential.
Collapse
Affiliation(s)
- 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
| | - James A Green
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - 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
|
20
|
Green JA, Yaghoubi Jouybari M, Asha H, Santoro F, Improta R. Fragment Diabatization Linear Vibronic Coupling Model for Quantum Dynamics of Multichromophoric Systems: Population of the Charge-Transfer State in the Photoexcited Guanine-Cytosine Pair. J Chem Theory Comput 2021; 17:4660-4674. [PMID: 34270258 DOI: 10.1021/acs.jctc.1c00416] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
We introduce a method (FrD-LVC) based on a fragment diabatization (FrD) for the parametrization of a linear vibronic coupling (LVC) model suitable for studying the photophysics of multichromophore systems. In combination with effective quantum dynamics (QD) propagations with multilayer multiconfigurational time-dependent Hartree (ML-MCTDH), the FrD-LVC approach gives access to the study of the competition between intrachromophore decays, like those at conical intersections, and interchromophore processes, like exciton localization/delocalization and the involvement of charge-transfer (CT) states. We used FrD-LVC parametrized with time-dependent density functional theory (TD-DFT) calculations, adopting either CAM-B3LYP or ωB97X-D functionals, to study the ultrafast photoexcited QD of a guanine-cytosine (GC) hydrogen-bonded pair, within a Watson-Crick arrangement, considering up to 12 coupled diabatic electronic states and the effect of all of the 99 vibrational coordinates. The bright excited states localized on C and, especially, on G are predicted to be strongly coupled to the G → C CT state, which is efficiently and quickly populated after an excitation to any of the four lowest energy bright local excited states. Our QD simulations show that more than 80% of the excited population on G and ∼50% of that on C decay to this CT state in less than 50 fs. We investigate the role of vibronic effects in the population of the CT state and show that it depends mainly on its large reorganization energy so that it can occur even when it is significantly less stable than the bright states in the Franck-Condon region. At the same time, we document that the formation of the GC pair almost suppresses the involvement of dark nπ* excited states in the photoactivated dynamics.
Collapse
Affiliation(s)
- James A Green
- Istituto di Biostrutture e Bioimmagini (IBB-CNR), Consiglio Nazionale delle Ricerche, via Mezzocannone 16, I-80136 Napoli, Italy
| | - Martha Yaghoubi Jouybari
- Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), Consiglio Nazionale delle Ricerche, SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Haritha Asha
- Istituto di Biostrutture e Bioimmagini (IBB-CNR), Consiglio Nazionale delle Ricerche, via Mezzocannone 16, I-80136 Napoli, Italy
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), Consiglio Nazionale delle Ricerche, SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Roberto Improta
- Istituto di Biostrutture e Bioimmagini (IBB-CNR), Consiglio Nazionale delle Ricerche, via Mezzocannone 16, I-80136 Napoli, Italy
| |
Collapse
|
21
|
Green JA, Jouybari MY, Aranda D, Improta R, Santoro F. Nonadiabatic Absorption Spectra and Ultrafast Dynamics of DNA and RNA Photoexcited Nucleobases. Molecules 2021; 26:1743. [PMID: 33804640 PMCID: PMC8003674 DOI: 10.3390/molecules26061743] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2021] [Revised: 03/12/2021] [Accepted: 03/14/2021] [Indexed: 11/16/2022] Open
Abstract
We have recently proposed a protocol for Quantum Dynamics (QD) calculations, which is based on a parameterisation of Linear Vibronic Coupling (LVC) Hamiltonians with Time Dependent (TD) Density Functional Theory (TD-DFT), and exploits the latest developments in multiconfigurational TD-Hartree methods for an effective wave packet propagation. In this contribution we explore the potentialities of this approach to compute nonadiabatic vibronic spectra and ultrafast dynamics, by applying it to the five nucleobases present in DNA and RNA. For all of them we computed the absorption spectra and the dynamics of ultrafast internal conversion (100 fs timescale), fully coupling the first 2-3 bright states and all the close by dark states, for a total of 6-9 states, and including all the normal coordinates. We adopted two different functionals, CAM-B3LYP and PBE0, and tested the effect of the basis set. Computed spectra are in good agreement with the available experimental data, remarkably improving over pure electronic computations, but also with respect to vibronic spectra obtained neglecting inter-state couplings. Our QD simulations indicate an effective population transfer from the lowest energy bright excited states to the close-lying dark excited states for uracil, thymine and adenine. Dynamics from higher-energy states show an ultrafast depopulation toward the more stable ones. The proposed protocol is sufficiently general and automatic to promise to become useful for widespread applications.
Collapse
Affiliation(s)
- James A. Green
- CNR—Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via Mezzocannone 16, I-80136 Napoli, Italy;
| | - Martha Yaghoubi Jouybari
- CNR—Consiglio Nazionale Delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area Della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy; (M.Y.J.); (D.A.)
| | - Daniel Aranda
- CNR—Consiglio Nazionale Delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area Della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy; (M.Y.J.); (D.A.)
| | - Roberto Improta
- CNR—Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via Mezzocannone 16, I-80136 Napoli, Italy;
| | - Fabrizio Santoro
- CNR—Consiglio Nazionale Delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area Della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy; (M.Y.J.); (D.A.)
| |
Collapse
|
22
|
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
|
23
|
Liu Y, Aranda D, Santoro F. A computational study of the vibronic effects on the electronic spectra and the photophysics of aza[7]helicene. Phys Chem Chem Phys 2021; 23:16551-16563. [PMID: 34338704 DOI: 10.1039/d1cp00822f] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report a computational study on vibronic effects in the spectroscopy, photoinduced processes and decay back to the ground state of aza[7]helicene, a helicene with an unusually high fluorescence quantum yield (QY = 0.39). In a first step, we compute and assign the absorption and electronic circular dichroism (ECD) spectra in its full frequency range from 2.7 to 5.0 eV, accounting for nonadiabatic effects. Then we compute the quantum dynamics of the cascade of ultrafast internal conversions of the highly-excited singlet states to the lowest-energy one S1. Finally we adopt Fermi golden rule rates to compute the QY of the dye, taking into account the competition between the radiative decay and the nonradiative decays to the ground state and to the energy-accessible triplet states. We use time-dependent density functional theory (TD-DFT), including solvent (dichloromethane) effects within the polarizable continuum model, to parameterize a linear vibronic coupling (LVC) model involving the first lowest 12 singlet states and all the normal coordinates. Nonadiabatic spectra and internal conversions dynamics are then computed through wavepacket propagations with the Multilayer (ML) extension of the Multiconfigurational Time Dependent Hartree method (ML-MCTDH). We highlight the molecular vibrations playing a major role in determining the shape of the spectra and analyse the effect of inter-state couplings. At the same time we report a breakdown of perturbative Herzberg-Teller approach. The computed QY is in perfect agreement with experiment and allows us to ascertain that intersystem crossings are the processes limiting the fluorescence from S1. They involve the three lowest triplet states and are made effective by spin-orbit coupling and vibronic effects.
Collapse
Affiliation(s)
- Yanli Liu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, P. R. China.
| | | | | |
Collapse
|