Toward a general mixed quantum/classical method for the calculation of the vibronic ECD of a flexible dye molecule with different stable conformers: Revisiting the case of 2,2,2-trifluoro-anthrylethanol

Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study

We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.

Taming the third order cumulant approximation to linear optical spectroscopy

The second order cumulant method offers a promising pathway to predicting optical properties in condensed phase systems. It allows for the computation of linear absorption spectra from excitation energy fluctuations sampled along molecular dynamics (MD) trajectories, fully accounting for vibronic effects, direct solute-solvent interactions, and environmental polarization effects. However, the second order cumulant approximation only guarantees accurate line shapes for energy gap fluctuations obeying Gaussian statistics. A third order correction has recently been derived but often yields unphysical spectra or divergent line shapes for moderately non-Gaussian fluctuations due to the neglect of higher order terms in the cumulant expansion. In this work, we develop a corrected cumulant approach, where the collective effect of neglected higher order contributions is approximately accounted for through a dampening factor applied to the third order cumulant term. We show that this dampening factor can be expressed as a function of the skewness and kurtosis of energy gap fluctuations and can be parameterized from a large set of randomly sampled model Hamiltonians for which exact spectral line shapes are known. This approach is shown to systematically remove unphysical contributions in the form of negative absorbances from cumulant spectra in both model Hamiltonians and condensed phase systems sampled from MD and dramatically improves over the second order cumulant method in describing systems exhibiting Duschinsky mode mixing effects. We successfully apply the approach to the coumarin-153 dye in toluene, obtaining excellent agreement with experiment.

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

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.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

The role of chlorine position in the electronic circular dichroism of chlorophenyl-ethanol investigated by vibronic calculations

(R)-1-phenyl-ethanol (PhEtOH) and the different isomers of (R)-1-(chlorophenyl)ethanol (ClPhEtOH) exhibit very interesting electronic circular dichroism (ECD) in methanol. In all cases, the spectrum shows clear vibronic features, but it is monosignated and negative for PhEtOH and meta-ClPhEtOH, positive for the ortho isomer and bisignated for the para isomer. We used computational chemistry to rationalise this behaviour adopting CAM-B3LYP/def2-TZVP, describing the bulk solvent effects with polarizable continuum models and solute-solvent specific interactions with clusters comprising the solute and two solvent molecules. We adopted harmonic vibronic models to compute the ECD spectral shapes of all stable conformers, and we obtained the room-temperature spectra by Boltzmann average. Simulated spectra are in very good agreement with experiment and allow us to rationalise their difference in terms of the relevance of Franck-Condon (FC) and Herzberg-Teller (HT) intensity-borrowing contributions, modulated by the substituent effect. The bisignated shape of the spectrum of para-ClPhEtOH arises from the competition of opposite-sign FC and HT bands, promoted by different vibrational modes. Due to the challenges we document in computing its ECD spectrum, para-ClPhEtOH represents a good test case to help the development of novel methodologies for an improved description of weak vibronic ECD spectra of flexible systems in explicit solvents.

Adiabatic-Molecular Dynamics Generalized Vertical Hessian Approach: A Mixed Quantum Classical Method To Compute Electronic Spectra of Flexible Molecules in the Condensed Phase

We present a general mixed quantum classical method that couples classical molecular dynamics (MD) and vibronic models to compute the shape of electronic spectra of flexible molecules in the condensed phase without, in principle, any phenomenological broadening. It is based on a partition of the nuclear motions of the solute + solvent system in "soft" and "stiff" vibrational modes and an adiabatic hypothesis that assumes that stiff modes are much faster than soft ones. In this framework, the spectrum is rigorously expressed as a conformational integral of quantum vibronic spectra along the stiff coordinates only. Soft modes enter at the classical level through the conformational distribution that is sampled with classical MD runs. In each configuration, reduced-dimensionality quadratic Hamiltonians are built in the space of the stiff coordinates only, thanks to a generalization of the Vertical Hessian harmonic model and an iterative application of projectors in internal coordinates to remove soft modes. Quantum vibronic spectra, specific for each sampled configuration of the soft coordinates, are then computed at the desired temperature with efficient time-dependent techniques, and the global spectrum simply arises from their average. For consistency of the whole procedure, classical MD runs are performed with quantum-mechanically derived force fields, parameterized at the same level of theory selected for generating the quadratic Hamiltonians along the stiff coordinates. Application to N-methyl-6-oxyquinolinium betaine in water, dithiophene in ethanol, and cyanidine in water is presented to show the performance of the method.

A computational study of the vibrationally-resolved electronic circular dichroism spectra of single-chain transoid and cisoid oligothiophenes in chiral conformations

We simulate the vibronic profile of the electronic circular dichroism (ECD) spectra of oligothiophenes in cisoid and transoid chiral arrangements. We consider oligomers of different lengths, from two to fifteen units, and investigate extensively how the ECD spectral shapes depend on the inter-ring torsions. In general, the molecular structures we consider are not stationary points of the ground state potential energy surface. Therefore, in order to perform vibronic calculations, we present a new computational protocol able to define reduced-dimensionality models where the effect of the off-equilibrium modes is removed. This is done adopting a description of the vibrational motions in curvilinear internal coordinates, and vertical harmonic models coupled with an iterative application of projectors to define energy Hessians, and therefore effective normal modes, in the space complementary to the one of the off-equilibrium coordinates. Although we consider both Franck-Condon and Herzberg-Teller contributions, the results show that transoid twisted ribbons always give rise to monosignated ECD spectra, while bi-signated and multi-signated spectra are expected for cisoid helices. These findings are explained on the basis of the different transition strengths of the lowest electronic states imparted by the different spatial arrangement, that is almost linear for transoid structures and more globular for cisoid ones. We predicted the chiroptical response of a large number of possible molecular arrangements. These data are employed to critically discuss the experimental ECD of polythiophenes in different experimental conditions, forming either aggregates or host-guest complexes. The method here proposed to perform vibronic calculations in reduced-dimensionality models is of general applicability and its potential interest goes beyond the practical application presented here.

Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water

Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.