Search

We study the absorption and emission electronic spectra in an aqueous solution of N-methyl-6-oxyquinolinium betaine (MQ), an interesting dye characterized by a large change of polarity and H-bond ability between the ground (S0) and the excited (S1) states. To that end we compare alternative approaches based either on explicit solvent models and density functional theory (DFT)/molecular-mechanics (MM) calculations or on DFT calculations on clusters models embedded in a polarizable continuum (PCM). In the first approach (ClMD), the spectrum is computed according to the classical Franck-Condon principle, from the dispersion of the time-dependent (TD)-DFT vertical transitions at selected snapshots of molecular dynamics (MD) on the initial state. In the cluster model (Qst) the spectrum is simulated by computing the quantum vibronic structure, estimating the inhomogeneous broadening from state-specific TD-DFT/PCM solvent reorganization energies. While both approaches provide absorption and emission spectral shapes in nice agreement with experiment, the Stokes shift is perfectly reproduced by Qst calculations if S0 and S1 clusters are selected on the grounds of the MD trajectory. Furthermore, Qst spectra better fit the experimental line shape, mostly in absorption. Comparison of the predictions of the two approaches is very instructive: the positions of Qst and ClMD spectra are shifted due to the different solvent models and the ClMD spectra are narrower than the Qst ones, because MD underestimates the width of the vibrational density of states of the high-frequency modes coupled to the electronic transition. On the other hand, both Qst and ClMD approaches highlight that the solvent has multiple and potentially opposite effects on the spectral width, so that the broadening due to solute-solvent vibrations and electrostatic interaction with bulk solvent is (partially) counterbalanced by a narrowing of the contribution due to the solute vibrational modes. Qst analysis evidences a pure quantum broadening effect of the spectra in water due to vibronic progressions along the solute/solvent H-bonds.

We simulated the intrinsic reaction path of the Green Fluorescent Protein (GFP) proton shuttle in both the ground state (S0) and first singlet excited state (S1), accounting for the main energetic and steric effects of the protein in a convenient model including the chromophore, the crystallographic water, and the residues directly involved in the proton transfer event. We adopted density functional theory (DFT) and time-dependent density functional theory (TD-DFT) levels to define the potential energy surfaces of the two electronic states, and we compared results obtained by the Damped Velocity Verlet and the Hessian-based Predictor-Corrector integrators of the intrinsic reaction coordinate, which gave a comparable and consistent picture of the mechanism. We show that, at S1, the GFP proton transfer becomes favored, with respect to S0, as suggested by the experimental evidence. As an important finding, this change is strictly related to the rearrangement of the hydrogen bond network composing the reaction path, which, in S1, relaxes to a tighter and planar configuration, as a consequence of the photoinduced relaxation in the GFP chromophore structure, thus prompting more effectively for the proton shuttle. Therefore, we give an unprecedented direct proof of the key role played by the photoinduced structural relaxation of the GFP on the chromophore photoacidity, validating, in particular, the hypothesis of Fang and co-workers [Nature 2009, 462, 200].

In-chain donor/acceptor block copolymers comprised of alternating electron rich/poor moieties are emerging as promising semiconducting chromophores for use in organic photovoltaic devices. The mobilities of charge carriers in these materials are experimentally probed using gated organic field-effect transistors to quantify electron and hole mobilities, but a mechanistic understanding of the relevant charge diffusion pathways is lacking. To elucidate the mechanisms of electron and hole transport following excitation to optically accessible low-lying valence states, we utilize mean-field quantum electronic dynamics in the TDDFT formalism to explicitly track the evolution of these photo-accessible states. From the orbital pathway traversed in the dynamics, p- and n-type conductivities can be distinguished. The electronic dynamics of the studied polymers show the time-resolved transitions between the initial photoexcited state, a tightly-bound excitonic state that is dark to the ground state, and a partially charge separated state indicated by long-lived, out-of-phase charge oscillations along the polymer backbone. The frequency of these charge oscillations yields an insight into the characteristic mobilities of charge carriers in these materials. When the barycenters of the electron and hole densities are followed during the dynamics, a pseudo-classical picture for the translation of charge carrier densities along the polymer backbone emerges that clarifies a crucial aspect in the design of efficient organic photovoltaic materials.

In this study, we analyze the photophysics of the N-methyl-6-oxyquinolinium betaine (MQ) solvatochromic dye in aqueous solution, focussing on the important structural rearrangement of its first solvation shell following the electronic excitation, which is characterized by a strong charge transfer. To this aim, we compare the results provided by ground- and excited-state ab-initio molecular dynamics with that of full QM calculations on clusters including a small number of solvent molecules (from 1 to 4), taking into account bulk solvent effects by the Polarizable Continuum Model. The two methods agree in predicting that while in the ground electronic state between three and four water molecules are strongly coordinated to the oxygen atom of MQ, in the excited state two water molecules are present in the first solvation layer of the MQ oxygen. Vertical excitation and emission energies computed on the structures provided by the two approaches allow for estimates of the Stokes shift consistent with the experimental results. On the ground of the present results, some general considerations on the advantages and limitations of the dynamical and static approaches in describing a photoactivated process in solution are proposed.

The detailed interpretation of time-resolved spectroscopic signals in terms of the molecular rearrangement during a photoreaction or a photophysical event is one of the most important challenges of both experimental and theoretical chemistry. Here we simulate a time-resolved fluorescence spectrum of a dye in aqueous solution, the N-methyl-6-oxyquinolinium betaine, and analyze it in terms of far IR and THz frequency contributions, providing a direct connection to specific molecular motions. To obtain this result, we build up an innovative and general approach based on excited state ab-initio molecular dynamics and a wavelet-based time-dependent frequency analysis of nonstationary signals. We obtain a nice agreement with key parameters of the solvent dynamics, such as the total Stokes shift and the Stokes shift relaxation times. As an important finding, we observe a strong change of specific solute-solvent interactions upon the electronic excitation, with the migration of about 1.5 water molecules from the first solvation shell toward the bulk. In spite of this event, the Stokes shift dynamics is ruled by collective solvent motions in the THz and far IR, which guide and modulate the strong rearrangement of the dye microsolvation. By the relaxation of THz and IR contributions to the emission signal, we can follow and understand in detail the molecularity of the process. The protocol presented here is, in principle, transferable to other time-resolved spectroscopic techniques.

In this work, we consider a new combination of vibrational analysis and normal-like mode decomposition of Debye-Waller factors of solvated ions entirely based on molecular dynamics data. Such a novel time-dependent analysis procedure provides a direct link between x-ray absorption fine structure parameters and normal mode contributions for an ion-solvent system. The potentialities of such a methodology rely on two fundamental aspects which distinguish it from already available tools. First, a general vibrational analysis that does not require any Gaussian or harmonic model for describing atomic fluctuations in liquids. Second, a very accurate sampling of the short range motions around the structural probe via the recently developed atom centered density matrix propagation/general liquid optimized boundary method. This novel molecular dynamics methodology is based on an integrated ab initio/classical potential using localized basis functions and nonperiodic boundary conditions. As a case study we have chosen the Zn(II) ion in aqueous solution. The consistency of our results and the observed good agreement with experiments show how the key support to advanced structural techniques from molecular dynamics can be further expanded and investigated.