Quantum Anharmonic Calculations of Vibrational Spectra for Water Adsorbed on Titania Anatase(101) Surface: Dissociative versus Molecular Adsorption.

The interaction of water molecules and hydroxyl groups with titanium dioxide (TiO2) surfaces is ubiquitous and very important in anatase nanoparticle photocatalytic processes. Infrared spectroscopy, assisted by ab initio calculations of vibrational frequencies, can be a powerful tool to elucidate the mechanisms behind water adsorption. However, a straightforward comparison between measurements and calculations remains a challenging task because of the complexity of the physical phenomena occurring on nanoparticle surfaces. Consequently, severe computational approximations, such as harmonic vibrational ones, are usually employed. In the present work we partially address this complexity issue by overcoming some of the standard approximations used in theoretical simulations and employ the Divide and Conquer Semiclassical Initial Value Representation (DC-SCIVR) molecular dynamics. This method allows to perform simulations of vibrational spectra of large dimensional systems accounting not only for anharmonicities, but also for nuclear quantum effects. We apply this computational method to water and deuterated water adsorbed on the ideal TiO2 anatase(101) surface, contemplating both the molecular and the dissociated adsorption processes. The results highlight not only the presence of an anharmonic shift of the frequencies in agreement with the experiments, but also complex quantum mechanical spectral signatures induced by the coupling of molecular vibrational modes with the surface ones, which are different in the hydrogenated case from the deuterated one. These couplings are further analyzed by exploiting the mode subdivision performed during the divide and conquer procedure.

Representing molecular ground and excited vibrational eigenstates with nuclear densities obtained from semiclassical initial value representation molecular dynamics.

We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in the protonated glycine molecule [Aieta et al., Nat Commun 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in the correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, also to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond lengths, angles, and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching, the typical harmonic nodal pattern is absent in the anharmonic distribution.

Anharmonic quantum nuclear densities from full dimensional vibrational eigenfunctions with application to protonated glycine.

The interpretation of molecular vibrational spectroscopic signals in terms of atomic motion is essential to understand molecular mechanisms and for chemical characterization. The signals are usually assigned after harmonic normal mode analysis, even if molecular vibrations are known to be anharmonic. Here we obtain the quantum anharmonic vibrational eigenfunctions of the 11-atom protonated glycine molecule and we calculate the density distribution of its nuclei and its geometry parameters, for both the ground and the O-H stretch excited states, using our semiclassical method based on ab initio molecular dynamics trajectories. Our quantum mechanical results describe a molecule elongated and more flexible with respect to what previously thought. More importantly, our method is able to assign each spectral peak in vibrational spectroscopy by showing quantitatively how normal modes involving different functional groups cooperate to originate that spectroscopic signal. The method will possibly allow for a better rationalization of experimental spectroscopy.

Anharmonic calculations of vibrational spectra for molecular adsorbates: A divide-and-conquer semiclassical molecular dynamics approach.

The vibrational spectroscopy of adsorbates is becoming an important investigation tool for catalysis and material science. This paper presents a semiclassical molecular dynamics method able to reproduce the vibrational energy levels of systems composed by molecules adsorbed on solid surfaces. Specifically, we extend our divide-and-conquer semiclassical method for power spectra calculations to gas-surface systems and interface it with plane-wave electronic structure codes. The Born-Oppenheimer classical dynamics underlying the semiclassical calculation is full dimensional, and our method includes not only the motion of the adsorbate but also those of the surface and the bulk. The vibrational spectroscopic peaks related to the adsorbate are accounted together with the most coupled phonon modes to obtain spectra amenable to physical interpretations. We apply the method to the adsorption of CO, NO, and H2O on the anatase-TiO2 (101) surface. We compare our semiclassical results with the single-point harmonic estimates and the classical power spectra obtained from the same trajectory employed in the semiclassical calculation. We find that CO and NO anharmonic effects of fundamental vibrations are similarly reproduced by the classical and semiclassical dynamics and that H2O adsorption is fully and properly described in its overtone and combination band relevant components only by the semiclassical approach.

An effective semiclassical approach to IR spectroscopy.

We present a novel approach to calculate molecular infrared (IR) spectra based on semiclassical (SC) molecular dynamics. The main advance from a previous SC method [M. Micciarelli et al. J. Chem. Phys. 149, 064115 (2018)] consists of the possibility to avoid state-to-state calculations making applications to systems characterized by sizable densities of vibrational states feasible. Furthermore, this new method accounts not only for positions and intensities of the several absorption bands which make up the IR spectrum but also for their shapes. We show that accurate SC IR spectra including quantum effects and anharmonicities for both frequencies and intensities can be obtained starting from SC power spectra. The approach is first tested against the water molecule and then applied to the 10-atom glycine amino acid.

Unraveling the molecular mechanisms of color expression in anthocyanins.

Anthocyanins are a broad family of natural dyes, increasingly finding application as substitutes for artificial colorants in the food industry. In spite of their importance and ubiquity, the molecular principles responsible for their extreme color variability are poorly known. We address these mechanisms by computer simulations and photoabsorption experiments of cyanidin-3-O-glucoside in water solution, as a proxy for more complex members of the family. Experimental results are presented in the range of pH 1-9, accompanied by a comprehensive systematic computational study across relevant charge states and tautomers. The computed spectra are in excellent agreement with the experiments, providing unprecedented insight into the complex behavior underlying color expression in these molecules. Besides confirming the importance of the molecule's charge state, we also unveil the hitherto unrecognized role of internal distortions in the chromophore, which affect its degree of conjugation, modulating the optical gap and in turn the color. This entanglement of structural and electronic traits is also shared by other members of the anthocyanin family (e.g. pelargonidin and delphinidin) highlighting a common mechanism for color expression across this important family of natural dyes.

Anharmonic vibrational eigenfunctions and infrared spectra from semiclassical molecular dynamics.

We describe a new approach based on semiclassical molecular dynamics that allows simulating infrared absorption or emission spectra of molecular systems with inclusion of anharmonic intensities. This is achieved from semiclassical power spectra by computing first the vibrational eigenfunctions as a linear combination of harmonic states, and then the oscillator strengths associated with the vibrational transitions. We test the approach against a 1D Morse potential and apply it to the water molecule with results in excellent agreement with discrete variable representation quantum benchmarks. The method does not require any grid calculations, and it is directly extendable to high dimensional systems. The usual exponential scaling of the basis set size with the dimensionality of the system can be avoided by means of an appropriate truncation scheme. Furthermore, the approach has the advantage to provide IR spectra beyond the harmonic approximation without losing the possibility of an intuitive assignment of absorption peaks in terms of normal modes of vibration.

Sampling Molecular Conformers in Solution with Quantum Mechanical Accuracy at a Nearly Molecular-Mechanics Cost.

We introduce a method to evaluate the relative populations of different conformers of molecular species in solution, aiming at quantum mechanical accuracy, while keeping the computational cost at a nearly molecular-mechanics level. This goal is achieved by combining long classical molecular-dynamics simulations to sample the free-energy landscape of the system, advanced clustering techniques to identify the most relevant conformers, and thermodynamic perturbation theory to correct the resulting populations, using quantum-mechanical energies from density functional theory. A quantitative criterion for assessing the accuracy thus achieved is proposed. The resulting methodology is demonstrated in the specific case of cyanin (cyanidin-3-glucoside) in water solution.

Multimodel Approach to the Optical Properties of Molecular Dyes in Solution.

We introduce a multimodel approach to the simulation of the optical properties of molecular dyes in solution, whereby the effects of thermal fluctuations and of dielectric screening on the absorption spectra are accounted for by explicit and implicit solvation models, respectively. Thermal effects are treated by averaging the spectra of molecular configurations generated by an ab initio molecular-dynamics simulation where solvent molecules are treated explicitly. Dielectric effects are then dealt with implicitly by computing the spectra upon removal of the solvent molecules and their replacement with an effective medium, in the spirit of a continuum solvation model. Our multimodel approach is validated by comparing its predictions with those of a fully explicit-solvation simulation for cyanidin-3-glucoside (cyanin) chromophore in water. While multimodel and fully explicit-solvent spectra may differ considerably for individual configurations along the trajectory, their time averages are remarkably similar, thus providing a solid benchmark of the former and allowing us to save considerably on the computer resources needed to predict accurate absorption spectra. The power of the proposed methodology is finally demonstrated by the excellent agreement between its predictions and the absorption spectra of cyanin measured at strong and intermediate acidity conditions.

Photophysics and photochemistry of a DNA-protein cross-linking model: a synergistic approach combining experiments and theory.

The photophysical and photochemical properties of 5-benzyluracil and 5,6-benzyluracil, the latter produced by photocyclization of the former through irradiation with femtosecond UV laser pulses, are investigated both experimentally and theoretically. The absorption spectra of the two molecules are compared, and the principal electronic transitions involved are discussed, with particular emphasis on the perturbation induced on the two chromophore species (uracil and benzene) by their proximity. The photoproduct formation for different irradiation times was verified with high-performance liquid chromatography and nuclear magnetic resonance measurements. The steady-state fluorescence demonstrates that the fluorescence is a distinctive physical observable for detection and selective identification of 5- and 5,6-benzyluracil. The principal electronic decay paths of the two molecules, obtained through TDDFT calculations, explain the features observed in the emission spectra and the photoreactivity of 5-benzyluracil. The order of magnitude of the lifetime of the excited states is derived with steady-state fluorescence anisotropy measurements and rationalized with the help of the computational findings. Finally, the spectroscopic data collected are used to derive the photocyclization and fluorescence quantum yields. In obtaining a global picture of the photophysical and photochemical properties of the two molecules, our findings demonstrates that the use of 5-benzyluracil as a model system to study the proximity relations of a DNA base with a close-lying aromatic amino acid is valid at a local level since the main characteristics of the decay processes from the excited states of the uracil/thymine molecules remain almost unchanged in 5-benzyluracil, the main perturbation arising from the presence of the close-lying aromatic group.

Low-lying excited-states of 5-benzyluracil.

A numerical study is reported concerning the first and second singlet excited-states of 5-benzyluracil using the multireference self-consistent field (state-averaged CASSCF) method. The vertical excitation energies of low-lying excited-states were characterized using the SA-CASSCF method, as well as using higher-level methods, such as CASPT2, MRCI and EOM-CCSD. The local minima and conical intersections found on the potential energy surfaces (PESs) were characterized in terms of molecular geometry and natural population analysis. Different relaxation pathways on the PESs are identified and discussed by comparing with the similar pathways found for the individual monomers of uracil and benzene. The molecule can be thought of as a model system for the study of crosslink reaction between DNA and proteins induced by UV light.

Pushing the frontiers of first-principles based computer simulations of chemical and biological systems.

The Laboratory of Computational Chemistry and Biochemistry is active in the development and application of first-principles based simulations of complex chemical and biochemical phenomena. Here, we review some of our recent efforts in extending these methods to larger systems, longer time scales and increased accuracies. Their versatility is illustrated with a diverse range of applications, ranging from the determination of the gas phase structure of the cyclic decapeptide gramicidin S, to the study of G protein coupled receptors, the interaction of transition metal based anti-cancer agents with protein targets, the mechanism of action of DNA repair enzymes, the role of metal ions in neurodegenerative diseases and the computational design of dye-sensitized solar cells. Many of these projects are done in collaboration with experimental groups from the Institute of Chemical Sciences and Engineering (ISIC) at the EPFL.