Search

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.

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(II) complex, [Ru(dcbpy) 2 (NCS) 2 ] 4- [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], in water solution, since this system belongs to a broader class of transition-metal compounds undergoing upon photo-excitation rapid and complex charge transfer (CT) dynamics, which can be dictated by structural rearrangement and solvent environment. Ab initio molecular dynamics (AIMD) relying on a hybrid quantum/molecular mechanics scheme is used to probe the equilibrium microsolvation around the metal complex in terms of radial distribution functions of the main solvation sites and solvent effects on the overall equilibrium structure. Then, using our AIMD-based generalized normal mode approach, we investigate how the ligand vibrational spectroscopic features are affected by water solvation, also contributing to the interpretation of experimental Infra-Red spectra. Two solvation sites are found for the ligands: the sulfur and the oxygen sites can interact on average with ∼4 and ∼3 water molecules, respectively, where a stronger interaction of the oxygen sites is highlighted. On average an overall dynamic distortion of the C 2 symmetric gas-phase structure was found to be induced by water solvation. Vibrational analysis reproduced experimental values for ligand symmetric and asymmetric stretchings, linking the observed shifts with respect to the gas-phase to a complex solvent distribution around the system. This is the groundwork for future excited-state nuclear and electronic dynamics to monitor non-equilibrium processes of CT excitation in complex environments, such as exciton migration in photovoltaic technologies.

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

In this work, excited state ab initio molecular dynamics together with a time resolved vibrational analysis is employed to shed light on the vibrational photoinduced dynamics of a well-known diarylethene molecule experiencing a ring opening reaction upon electronic excitation. The photoreactivity of diarylethenes is recognized to be controlled by a non-adiabatic intersection point between the ground and the first excited state surfaces. The computation of an energy scan, along a suitable reaction coordinate, allows us to identify the region of potential energy surfaces in which the ground (S0) and the first excited (S1) state are well separated. The adiabatic sampling of that region in S1 shows that in the first 3 picoseconds, the central CC bond, which is subject to break, oscillates in an antiphase with respect to the energy gap ΔE(S1 - S0). A multiresolution analysis based on the wavelet transform was then applied to the structural parameters extracted from the excited state dynamics. The wavelet maps show characteristic oscillations of the frequencies, mainly CC stretching and CCC bending localized on the central 4-ring moiety. Moreover, we have identified the main frequency (methyl wagging motion) involved in the modulation of these oscillations. The anharmonic coupling within a group of vibrational modes was therefore highlighted, in good agreement with experimental evidence. For the first time, a quantitative analysis of time resolved signals from a wavelet transform/ab initio molecular dynamics approach was performed.

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.

We individuate a photoinduced electron transfer (PeT) as a quenching mechanism affecting rhodamine B photophysics in solvent. The PeT involves an electron transfer from the carboxylate group to the xanthene ring of rhodamine B. This is finely modulated by the subtle balance of coulombic and non-classical interactions between the carboxyphenyl and xanthene rings, also mediated by the solvent. We propose the use of an electronic density based index, the so called DCT index, as a new tool to assess and quantify the nature of the excited states involved in non-radiative decays near the region of their intersection. In the present case, this analysis allows us to gain insight on the interconversion process from the bright state to the dark state responsible for the quenching of rhodamine B fluorescence. Our findings encourage the use of density based indices to study the processes affecting excited state reactions that are characterized by a drastic change in the excitation nature, in order to rationalize the photophysical behavior of complex molecular systems.

The mechanism of base to base intermolecular proton shuttling occurring in the excited state proton transfer reaction between 7-hydroxy-4-(trifluoromethyl)coumarin (CouOH) and concentrated 1-methylimidazole base (1-MeId) in toluene solution is disclosed here by means of a computational approach based on Density Functional Theory (DFT) and Time Dependent DFT (TD-DFT). These methods allow us to characterize both the ground and excited state potential energy surfaces along the proton shuttling coordinate, and to assess the nature of the emitting species in the presence of an excess of 1-MeId. As a result, the tautomerism of CouOH is found to be photo-activated and, from a mechanistic point of view, the calculations clearly show that the overall driving force of the entire shuttling is the coumarin photoacidity, which is responsible for both the first proton transfer event and the strengthening of the following chain mechanism of base to base proton hopping.

In spite of the large number of experimental and theoretical studies, the optical absorption trend of the green fluorescent protein (GFP) chromophore in several environments has not been fully understood. We calculated at the same level of time dependent density functional theory the vertical excitation energy of the anionic GFP chromophore in the protein and in ethanol, dioxane, methanol and water solutions. As result, we reproduced for the first time the experimental trend of the absorption peaks with 0.015 eV as the standard deviation of the accuracy. This systematic error allowed us to analyze with confidence the relative weight of several solvation effects on the vertical excitation energy. Experimental trends not correlated with the solvent polarity were therefore explained with a fine balance of different steric and electronic effects on the photophysics of the chromophore. As an indirect and remarkable result, the present analysis confirms that the optical absorption of the chromophore in the gas phase is close to the value of 2.84 eV extrapolated by Dong et al. (J. Am. Chem. Soc., 2006, 128, 12038), and, as a consequence, that the protein environment induces a red shift of 0.23 eV.

ESR experimental data combined with quantum mechanical calculations and first-principle molecular dynamics simulations have allowed the unequivocal identification and characterization of the canonical form of the uracil anion radical as the prevalent tautomer in aqueous solution, contrary to both gas-phase results and recently reported results in solution based on thermodynamic analyses. The present study confirms the effectiveness of spectroscopic/theoretical investigations of short lived molecular species in different environments.