Electronic spectroscopy of a solvatochromic dye in water: comparison of static cluster/implicit and dynamical/explicit solvent models on structures and energies

Unexpected longer T1 lifetime of 6-sulfur guanine than 6-selenium guanine: the solvent effect of hydrogen bonds to brake the triplet decay

The decay of the T1 state to the ground state is an essential property of photosensitizers because it decides the lifetime of excited states and, thus, the time window for sensitization. The sulfur/selenium substitution of carbonyl groups can red-shift absorption spectra and enhance the triplet yield because of the large spin-orbit coupling, modifying nucleobases to potential photosensitizers for various applications. However, replacing sulfur with selenium will also cause a much shorter T1 lifetime. Experimental studies found that the triplet decay rate of 6-seleno guanine (6SeGua) is 835 times faster than that of 6-thio guanine (6tGua) in aqueous solution. In this work, we reveal the mechanism of the T1 decay difference between 6SeGua and 6tGua by computing the activation energy and spin-orbit coupling for rate calculation. The solvent effect of water is treated with explicit microsolvation and implicit solvent models. We find that the hydrogen bond between the sulfur atom of 6tGua and the water molecule can brake the triplet decay, which is weaker in 6SeGua. This difference is crucial to explain the relatively long T1 lifetime of 6tGua in an aqueous solution. This insight emphasizes the role of solvents in modulating the excited state dynamics and the efficiency of photosensitizers, particularly in aqueous environments.

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.

An Expedited Route to Optical and Electronic Properties at Finite Temperature via Unsupervised Learning

Electronic properties and absorption spectra are the grounds to investigate molecular electronic states and their interactions with the environment. Modeling and computations are required for the molecular understanding and design strategies of photo-active materials and sensors. However, the interpretation of such properties demands expensive computations and dealing with the interplay of electronic excited states with the conformational freedom of the chromophores in complex matrices (i.e., solvents, biomolecules, crystals) at finite temperature. Computational protocols combining time dependent density functional theory and ab initio molecular dynamics (MD) have become very powerful in this field, although they require still a large number of computations for a detailed reproduction of electronic properties, such as band shapes. Besides the ongoing research in more traditional computational chemistry fields, data analysis and machine learning methods have been increasingly employed as complementary approaches for efficient data exploration, prediction and model development, starting from the data resulting from MD simulations and electronic structure calculations. In this work, dataset reduction capabilities by unsupervised clustering techniques applied to MD trajectories are proposed and tested for the ab initio modeling of electronic absorption spectra of two challenging case studies: a non-covalent charge-transfer dimer and a ruthenium complex in solution at room temperature. The K-medoids clustering technique is applied and is proven to be able to reduce by ∼100 times the total cost of excited state calculations on an MD sampling with no loss in the accuracy and it also provides an easier understanding of the representative structures (medoids) to be analyzed on the molecular scale.

Water-Mediated Excited State Proton Transfer of Pyranine-Acetate in Aqueous Solution: Vibrational Fingerprints from Ab Initio Molecular Dynamics

In this work, we simulate the excited state proton transfer (ESPT) reaction involving the pyranine photoacid and an acetate molecule as proton acceptor, connected by a bridge water molecule. We employ ab initio molecular dynamics combined with an hybrid quantum/molecular mechanics (QM/MM) framework. Furthermore, a time-resolved vibrational analysis based on the wavelet-transform allows one to identify two low frequency vibrational modes that are fingerprints of the ESPT event: a ring wagging and ring breathing. Their composition suggests their key role in optimizing the structure of the proton donor–acceptor couple and promoting the ESPT event. We find that the choice of the QM/MM partition dramatically affects the photoinduced reactivity of the system. The QM subspace was gradually extended including the water molecules directly interacting with the pyranine–water–acetate system. Indeed, the ESPT reaction takes place when the hydrogen bond network around the reactive system is taken into account at full QM level.

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

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.

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.