Infrared Spectroscopy of N-Methylacetamide Revisited by ab Initio Molecular Dynamics Simulations

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Unravelling Mn4Ca cluster vibrations in the S1, S2 and S3 states of the Kok-Joliot cycle of photosystem II

Vibrational spectroscopy serves as a powerful tool for characterizing intermediate states within the Kok-Joliot cycle. In this study, we employ a QM/MM molecular dynamics framework to calculate the room temperature infrared absorption spectra of the S1, S2, and S3 states via the Fourier transform of the dipole time auto-correlation function. To better analyze the computational data and assign spectral peaks, we introduce an approach based on dipole-dipole correlation function of cluster moieties of the reaction center. Our analysis reveals variation in the infrared signature of the Mn4Ca cluster along the Kok-Joliot cycle, attributed to its increasing symmetry and rigidity resulting from the rising oxidation state of the Mn ions. Furthermore, we successfully assign the debated contributions in the frequency range around 600 cm-1. This computational methodology provides valuable insights for deciphering experimental infrared spectra and understanding the water oxidation process in both biological and artificial systems.

Vibrational Spectra Simulations in Amino Acid-Based Imidazolium Ionic Liquids

We present maximally localized Wannier functions and Voronoi tessellation to obtain dipole moment distributions for vibrational spectra in several important ionic liquids calculated by using ab initio molecular dynamics simulations. IR and Raman spectra of various imidazolium-based ionic liquids (ILs) paired with six amino acid anions are shown herein. For IR spectra, two approaches (Wannier and Voronoi) are in agreement with respect to the relative intensities and the overall shapes for the main peaks. Under Raman spectra, the polarizability of the covalent bonds is shown to affect the strength of the Raman scattering signal. The advantage of the Voronoi tessellation method, being that it does not have strong spikes in its time development, is demonstrated by the comparison of two theoretical methods (Wannier and Voronoi) with experimental data. We analyze the errors between theoretical and experimental spectroscopic data, with the Voronoi method shown to accurately reproduce experimental values. In addition, theoretical spectroscopy shows the ability to accurately separate components of a mixture. The combination of theoretical and experimental methods is utilized to understand the spectroscopic properties of amino acid-based imidazolium ILs.

Evidence of Excited-State Vibrational Mode Governing the Photorelaxation of a Charge-Transfer Complex

Modern, nonlinear, time-resolved spectroscopic techniques have opened new doors for investigating the intriguing but complex world of photoinduced ultrafast out-of-equilibrium phenomena and charge dynamics. The interaction between light and matter introduces an additional dimension, where the complex interplay between electronic and vibrational dynamics needs the most advanced theoretical-computational protocols to be fully understood on the molecular scale. In this study, we showcase the capabilities of ab initio molecular dynamics simulation integrated with a multiresolution wavelet protocol to carefully investigate the excited-state relaxation dynamics in a noncovalent complex involving tetramethylbenzene (TMB) and tetracyanoquinodimethane (TCNQ) undergoing charge transfer (CT) upon photoexcitation. Our protocol provides an accurate description that facilitates a direct comparison between transient vibrational analysis and time-resolved spectroscopic signals. This molecular level perspective enhances our understanding of photorelaxation processes confined in the adiabatic regime and offers an improved interpretation of vibrational spectra. Furthermore, it enables the quantification of anharmonic vibrational couplings between high- and low-frequency modes, specifically the TCNQ "rocking" and "bending" modes. Additionally, it identifies the primary vibrational mode that governs the adiabaticity between the ground state and the CT state. This comprehensive understanding of photorelaxation processes holds significant importance in the rational design and precise control of more efficient photovoltaic and sensor devices.

Polarisation effects on the H-bond acceptor properties of secondary amides

H-bonding interactions in networks are stabilised by cooperativity, but the relationship between the chemical structures of the interacting functional groups and the thermodynamic consequences is not well-understood. We have used compounds with an intramolecular H-bond between a pyridine H-bond acceptor and an amide NH group to quantify cooperative effects on the H-bond acceptor properties of the amide carbonyl group. 1H NMR experiments in n-octane confirm the presence of the intramolecular H-bond and show that this interaction is intact in the 1 : 1 complex formed with perfluoro-tert-butanol (PFTB). UV-vis absorption titrations were used to measure the relationship between the association constant for formation of this complex and the H-bond acceptor properties of the pyridine involved in the intramolecular H-bond. Electron-donating substituents on the pyridine increase the strength of the intermolecular H-bond between PFTB and the amide. There is a linear relationship between the H-bond acceptor parameter β measured for the amide carbonyl group and the H-bond acceptor parameter for the pyridine. The cooperativity parameter κ determined from this relationship is 0.2, i.e. β for an amide carbonyl group is increased by one fifth of the value of β of an acceptor that interacts with the NH group. This result is reproduced by DFT calculations of H-bond parameters for the individual molecules in the gas phase, which implies that the observed cooperativity can be understood as polarisation of the electron density in the amide π-system in response to formation of a H-bond. The cooperativity parameter κ measured for the secondary amide H-bond donor and H-bond acceptor is identical, which implies that polarisation of an amide mediates the interaction between an external donor or acceptor in a reciprocal manner.

DFT and ab initio molecular dynamics simulation study of the infrared spectrum of the protic ionic liquid 2-hydroxyethylammonium formate

Protic ionic liquids (PILs) typically show a complex band shape in their infrared (IR) spectra in the high-frequency range due to the hydrogen stretching vibrations of functional groups forming rather strong hydrogen bonds (H-bonds). In the low-frequency range, the intermolecular stretching mode of the H-bond leaves a mark in the far-IR spectrum of PILs. In this study, the IR spectrum of the PIL 2-hydroxyethylammonium formate, [HOCH2CH2NH3][HCOO], is investigated in order to identify the different modes that contribute to the high-frequency band shape, i.e. the cation ν(NH), ν(OH), and ν(CH) modes, and the anion ν(CH) mode, as well as the intermolecular mode of the strongest H-bond in the far-IR spectrum. The assignment is validated by quantum chemistry calculations of clusters at the density functional theory (DFT) level for four ionic pairs and by ab initio molecular dynamics (AIMD) simulations of ten ionic pairs. There is good agreement between the vibrational frequencies obtained from DFT and AIMD simulations for both the high- and low-frequency ranges. Based on the calculations, the strong H-bond interaction between the cation -NH3 group and [HCOO]- gives a broad band envelope associated with the ν(NH) mode in the high-frequency range of the IR spectrum on which there are narrower peaks corresponding to the ν(OH) and ν(CH) modes. In the far-IR (FIR) spectrum, the anions' rattling motion gives a broad feature with a maximum at 160 cm-1, while the H-bond's intermolecular NH⋯O stretching mode appears as a peak at 255 cm-1.

A Molecular Dynamic Study of the Effects of Surface Partitioning on the OH Radical Interactions with Solutes in Multicomponent Aqueous Aerosols

The surface-bulk partitioning of small saccharide and amide molecules in aqueous droplets was investigated using molecular dynamics. The air-particle interface was modeled using a 80 Å cubic water box containing a series of organic molecules and surrounded by gaseous OH radicals. The properties of the organic solutes within the interface and the water bulk were examined at a molecular level using density profiles and radial pair distribution functions. Molecules containing only polar functional groups such as urea and glucose are found predominantly in the water bulk, forming an exclusion layer near the water surface. Substitution of a single polar group by an alkyl group in sugars and amides leads to the migration of the molecule toward the interface. Within the first 2 nm from the water surface, surface-active solutes lose their rotational freedom and adopt a preferred orientation with the alkyl group pointing toward the surface. The different packing within the interface leads to different solvation shell structures and enhanced interaction between the organic molecules and absorbed OH radicals. The simulations provide quantitative information about the dimension, composition, and organization of the air-water interface as well as about the nonreactive interaction of the OH radicals with the organic solutes. It suggests that increased concentrations, preferred orientations, and decreased solvation near the air-water surface may lead to differences in reactivities between surface-active and surface-inactive molecules. The results are important to explain how heterogeneous oxidation mechanisms and kinetics within interfaces may differ from those of the bulk.

Experimental Detection and Simulation of Terahertz Spectra of Aqueous L-Arginine

Terahertz (THz) wave is a good candidate for biological sample detection, because vibration and rotation energy levels of biomolecule are in THz band. However, the strong absorption of THz wave by water in biological samples hinders its development. In this paper, a method for direct detection of THz absorption spectra of L-arginine suspension was proposed by using a strong field THz radiation source combined with a polyethylene cell with micrometer thickness in a THz time-domain spectroscopy system. And the THz absorption spectrum of L-arginine solution was simulated by the density functional theory and the simulation result is in good agreement with the experimental results. Finally, the types of chemical bond interaction that cause the absorption peak are identified based on the experimental and simulation results. This work paves a way to investigate the THz absorption spectra and intramolecular interactions of aqueous biological samples.

H-bond cooperativity: polarisation effects on secondary amides

Formation of a H-bond with an amide carbonyl oxygen atom increases the strength of subsequent H-bonds formed by the amide NH, due to polarisation of the bond. The magnitude of this effect has been quantified by measuring association constants for the formation of 1 : 1 complexes of 2-hydroxylbenzamides with tri-n-butyl phosphine oxide. In 2-hydroxybenzamides, there is an intramolecular H-bond between the phenol OH group and the carbonyl oxygen atom. Comparison of the association constants measured for compounds with and without the 2-hydroxy group allows direct quantification of the effect of the intramolecular H-bond on the H-bond donor properties of the amide NH group. Substituents were used to modulate the strength of the intramolecular and intermolecular H-bonds. The presence of an intramolecular H-bond increases the strength of the intermolecular H-bond by more than one order of magnitude in n-octane solution. The increase in the H-bond donor parameter used to describe the amide NH group is directly proportional to the H-bond donor parameter of the phenol OH group that makes the intramolecular H-bond. These polarisation effects will lead to substantial cooperativity in complex systems that feature networks of non-covalent interactions, and the measurements described here provide a quantitative basis for understanding such phenomena.

Formation of an intramolecular phenol-amide H-bond leads to a dramatic increase in the H-bond donor strength of the amide NH group. Polarisation of the amide group is directly proportional to the polarity of the phenol H-bond donor.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

Multipolar Force Fields for Amide-I Spectroscopy from Conformational Dynamics of the Alanine Trimer

The dynamics and spectroscopy of N-methyl-acetamide (NMA) and trialanine in solution are characterized from molecular dynamics simulations using different energy functions, including a conventional point charge (PC)-based force field, one based on a multipolar (MTP) representation of the electrostatics, and a semiempirical DFT method. For the 1D infrared spectra, the frequency splitting between the two amide-I groups is 10 cm^-1 from the PC, 13 cm^-1 from the MTP, and 47 cm^-1 from self-consistent charge density functional tight-binding (SCC-DFTB) simulations, compared with 25 cm^-1 from experiment. The frequency trajectory required for the frequency fluctuation correlation function (FFCF) is determined from individual normal mode (INM) and full normal mode (FNM) analyses of the amide-I vibrations. The spectroscopy, time-zero magnitude of the FFCF C(t = 0), and the static component Δ₀² from simulations using MTP and analysis based on FNM are all consistent with experiments for (Ala)₃. Contrary to this, for the analysis excluding mode-mode coupling (INM), the FFCF decays to zero too rapidly and for simulations with a PC-based force field, the Δ₀² is too small by a factor of two compared with experiments. Simulations with SCC-DFTB agree better with experiment for these observables than those from PC-based simulations. The conformational ensemble sampled from simulations using PCs is consistent with the literature (including P_II, β, α_R, and α_L), whereas that covered by the MTP-based simulations is dominated by P_II with some contributions from β and α_R. This agrees with and confirms recently reported Bayesian-refined populations based on 1D infrared experiments. FNM analysis together with a MTP representation provides a meaningful model to correctly describe the dynamics of hydrated trialanine.

Effect of the Hydration Shell on the Carbonyl Vibration in the Ala-Leu-Ala-Leu Peptide

The vibrational spectrum of the Ala-Leu-Ala-Leu peptide in solution, computed from first-principles simulations, shows a prominent band in the amide I region that is assigned to stretching of carbonyl groups. Close inspection reveals combined but slightly different contributions by the three carbonyl groups of the peptide. The shift in their exact vibrational signature is in agreement with the different probabilities of these groups to form hydrogen bonds with the solvent. The central carbonyl group has a hydrogen bond probability intermediate to the other two groups due to interchanges between different hydrogen-bonded states. Analysis of the interaction energies of individual water molecules with that group shows that shifts in its frequency are directly related to the interactions with the water molecules in the first hydration shell. The interaction strength is well correlated with the hydrogen bond distance and hydrogen bond angle, though there is no perfect match, allowing geometrical criteria for hydrogen bonds to be used as long as the sampling is sufficient to consider averages. The hydrogen bond state of a carbonyl group can therefore serve as an indicator of the solvent's effect on the vibrational frequency.

TAO-DFT-Based Ab Initio Molecular Dynamics

Recently, AIMD (ab initio molecular dynamics) has been extensively employed to explore the dynamical information of electronic systems. However, it remains extremely challenging to reliably predict the properties of nanosystems with a radical nature using conventional electronic structure methods (e.g., Kohn-Sham density functional theory) due to the presence of static correlation. To address this challenge, we combine the recently formulated TAO-DFT (thermally-assisted-occupation density functional theory) with AIMD. The resulting TAO-AIMD method is employed to investigate the instantaneous/average radical nature and infrared spectra of n-acenes containing n linearly fused benzene rings (n = 2-8) at 300 K. According to the TAO-AIMD simulations, on average, the smaller n-acenes (up to n = 5) possess a nonradical nature, and the larger n-acenes (n = 6-8) possess an increasing radical nature, showing remarkable similarities to the ground-state counterparts at 0 K. Besides, the infrared spectra of n-acenes obtained with the TAO-AIMD simulations are in qualitative agreement with the existing experimental data.

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.

The Amide I Spectrum of Proteins-Optimization of Transition Dipole Coupling Parameters Using Density Functional Theory Calculations

The amide I region of the infrared spectrum is related to the protein backbone conformation and can provide important structural information. However, the interpretation of the experimental results is hampered because the theoretical description of the amide I spectrum is still under development. Quantum mechanical calculations, for example, using density functional theory (DFT), can be used to study the amide I spectrum of small systems, but the high computational cost makes them inapplicable to proteins. Other approaches that solve the eigenvalues of the coupled amide I oscillator system are used instead. An important interaction to be considered is transition dipole coupling (TDC). Its calculation depends on the parameters of the transition dipole moment. This work aims to find the optimal parameters for TDC in three major secondary structures: α-helices, antiparallel β-sheets, and parallel β-sheets. The parameters were suggested through a comparison between DFT and TDC calculations. The comparison showed a good agreement for the spectral shape and for the wavenumbers of the normal modes for all secondary structures. The matching between the two methods improved when hydrogen bonding to the amide oxygen was considered. Optimal parameters for individual secondary structures were also suggested.

Gas phase dynamics, conformational transitions and spectroscopy of charged saccharides: the oxocarbenium ion, protonated anhydrogalactose and protonated methyl galactopyranoside

Protonated intermediates are postulated to be involved in the rate determining step of many sugar reactions. This paper presents a study of protonated sugar species, isolated in the gas phase, using a combination of infrared multiple photon dissociation (IRMPD) spectroscopy, classical ab initio molecular dynamics (AIMD) and quantum mechanical vibrational self-consistent field (VSCF) calculations. It provides a likely identification of the reactive intermediate oxocarbenium ion structure in a d-galactosyl system as well as the saccharide pyrolysis product anhydrogalactose (that suggests oxocarbenium ion stabilization), along with the spectrum of the protonated parent species: methyl d-galactopyranoside-H+. Its vibrational fingerprint indicates intramolecular proton sharing. Classical AIMD simulations for galactosyl oxocarbenium ions, conducted in the temperature range ∼300-350 K (using B3LYP potentials on-the-fly) reveal efficient transitions on the picosecond timescale. Multiple conformers are likely to exist under the experimental conditions and along with static VSCF calculations, they have facilitated the identification of the individual structural motifs of the galactosyl oxocarbenium ion and protonated anhydrogalactose ion conformers that contribute to the observed experimental spectra. These results demonstrate the power of experimental IRMPD spectroscopy combined with dynamics simulations and with computational spectroscopy at the anharmonic level to unravel conformer structures of protonated saccharides, and to provide information on their lifetimes.

Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain

Gas-phase, double resonance IR spectroscopy has proven to be an excellent approach to obtain structural information on peptides ranging from single amino acids to large peptides and peptide clusters. In this review, we discuss the state-of-the-art of infrared action spectroscopy of peptides in the far-IR and THz regime. An introduction to the field of far-IR spectroscopy is given, thereby highlighting the opportunities that are provided for gas-phase research on neutral peptides. Current experimental methods, including spectroscopic schemes, have been reviewed. Structural information from the experimental far-IR spectra can be obtained with the help of suitable theoretical approaches such as dynamical DFT techniques and the recently developed Graph Theory. The aim of this review is to underline how the synergy between far-IR spectroscopy and theory can provide an unprecedented picture of the structure of neutral biomolecules in the gas phase. The far-IR signatures of the discussed studies are summarized in a far-IR map, in order to gain insight into the origin of the far-IR localized and delocalized motions present in peptides and where they can be found in the electromagnetic spectrum.

Computing Bulk Phase Resonance Raman Spectra from ab Initio Molecular Dynamics and Real-Time TDDFT

We present our novel approach for computing resonance Raman (RR) spectra of periodic bulk phase systems from ab initio molecular dynamics, including solvent influence and some anharmonic effects. Based on real-time time-dependent density functional theory, we obtain the RR spectra for all laser wavelengths in one pass. We compute the RR spectrum of uracil in aqueous solution, which is in good agreement with experiment. This is the first simulation of a bulk phase RR spectrum.

A fragmented, permutationally invariant polynomial approach for potential energy surfaces of large molecules: Application to N-methyl acetamide

We describe and apply a method to extend permutationally invariant polynomial (PIP) potential energy surface (PES) fitting to molecules with more than 10 atoms. The method creates a compact basis of PIPs as the union of PIPs obtained from fragments of the molecule. An application is reported for trans-N-methyl acetamide, where B3LYP/cc-pVDZ electronic energies and gradients are used to develop a full-dimensional potential for this prototype peptide molecule. The performance of several fragmented bases is verified against a benchmark PES using all (66) Morse variables. The method appears feasible for much larger molecules.

Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water

The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.

Conformational assignment of gas phase peptides and their H-bonded complexes using far-IR/THz: IR-UV ion dip experiment, DFT-MD spectroscopy, and graph theory for mode assignment

Graph theory based vibrational modes as new entities for vibrational THz spectroscopy.

Aqueous solvation of Mg(ii) and Ca(ii): A Born-Oppenheimer molecular dynamics study of microhydrated gas phase clusters

The hydration features of [Mg(H₂O)_n]²⁺ and [Ca(H₂O)_n]²⁺ clusters with n = 3-6, 8, 18, and 27 were studied by means of Born-Oppenheimer molecular dynamics simulations at the B3LYP/6-31+G** level of theory. For both ions, it is energetically more favorable to have all water molecules in the first hydration shell when n ≤ 6, but stable lower coordination average structures with one water molecule not directly interacting with the ion were found for Mg²⁺ at room temperature, showing signatures of proton transfer events for the smaller cation but not for the larger one. A more rigid octahedral-type structure for Mg²⁺ than for Ca²⁺ was observed in all simulations, with no exchange of water molecules to the second hydration shell. Significant thermal effects on the average structure of clusters were found: while static optimizations lead to compact, spherically symmetric hydration geometries, the effects introduced by finite-temperature dynamics yield more prolate configurations. The calculated vibrational spectra are in agreement with infrared spectroscopy results. Previous studies proposed an increase in the coordination number (CN) from six to eight water molecules for [Ca(H₂O)_n]²⁺ clusters when n ≥ 12; however, in agreement with recent measurements of binding energies, no transition to a larger CN was found when n > 8. Moreover, the excellent agreement found between the calculated extended X-ray absorption fine structure spectroscopy spectra for the larger cluster and the experimental data of the aqueous solution supports a CN of six for Ca²⁺.

Structure and dynamics of the peptide strand KRFK from the thrombospondin TSP-1 in water

Theoretical investigations of a solute in liquid water at normal temperature and pressure can be performed at different levels of theory. Static quantum calculations as well as classical and ab initio molecular dynamics are used to completely explore the conformational space for large solvated molecular systems. In the classical approach, it is essential to describe all of the interactions of the solute and the solvent in detail. Water molecules are very often described as rigid bodies when the most commonly used interaction potentials, such as the SPCE and the TIP4P models, are employed. Recently, a physical model based upon a cluster of rigid water molecules with a tetrahedral architecture (AB₄) was proposed that describes liquid water as a mixture of both TIP4P and SPCE molecular species that occur in the proportions implied by the tetrahedral architecture (one central molecule versus four outer molecules; i.e., 20% TIP4P versus 80% SPCE molecules). In this work, theoretical spectroscopic data for a peptide strand were correlated with the structural properties of the peptide strand solvated in water, based on data calculated using different theoretical approaches and physical models. We focused on a particular peptide strand, KRFK (lysine-arginine-phenylalanine-lysine), found in the thrombospondin TSP-1, due to its interesting properties. As the activity and electronic structure of this system is strongly linked to its structure, we correlated its structure with charge-density maps obtained using different semi-empirical charge Q_eq equations. The structural and thermodynamic properties obtained from classical simulations were correlated with ab initio molecular dynamics (AIMD) data. Structural changes in the peptide strand were rationalized in terms of the motions of atoms and groups of atoms. To achieve this, conformational changes were investigated using calculated infrared spectra for the peptide in the gas phase and in water solvent. The calculated AIMD infrared spectrum for the peptide was correlated with static quantum calculations of the molecular system based on a harmonic approach as well as the VDOS (vibrational density of states) spectra obtained using various classical solvent models (SPCE, TIP4P, and AB₄) and charge maps.

The dynamics of peptide-water interactions in dialanine: An ultrafast amide I 2D IR and computational spectroscopy study

We present a joint experimental and computational study of the dynamic interactions of dialanine (Ala-Ala) with water, comparing the results of ultrafast 2D IR and infrared transient absorption spectroscopy of its amide I vibration with spectra modeled from molecular dynamics (MD) simulations. The experimental data are analyzed to describe vibrational frequency fluctuations, vibrational energy relaxation, and chemical exchange processes. The origin of these processes in the same underlying fluctuating forces allows a common description in terms of the fluctuations and conformational dynamics of the peptide and associated solvent. By comparing computational spectroscopy from MD simulations with multiple force fields and water models, we describe how the dynamics of water hydrogen bond fluctuations and switching processes act as a source of friction that governs the dephasing and vibrational relaxation, and provide a description of coupled water and peptide motions that give rise to spectroscopic exchange processes.

Static vs dynamic DFT prediction of IR spectra of flexible molecules in the condensed phase: The (ClCF2CF(CF3)OCF2CH3) liquid as a test case

First-principles molecular dynamics (FPMD) simulations in the framework of Density Functional Theory (DFT) are carried out for the prediction of the infrared spectrum of the fluorinated molecule ClCF₂CF(CF₃)OCF₂CH₃ in liquid and gas phase. This molecule is characterized by a flexible structure, allowing the co-existence of several stable conformers, that differ by values of the torsional angles. FPMD computed spectra are compared to the experimental ones, and to Boltzmann weighted IR spectra based on gas phase calculations.

Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations

We present our novel methodology for computing Raman optical activity (ROA) spectra of liquid systems from ab initio molecular dynamics (AIMD) simulations. The method is built upon the recent developments to obtain magnetic dipole moments from AIMD and to integrate molecular properties by using radical Voronoi tessellation. These techniques are used to calculate optical activity tensors for large and complex periodic bulk phase systems. Only AIMD simulations are required as input, and no time-consuming perturbation theory is involved. The approach relies only on the total electron density in each time step and can readily be combined with a wide range of electronic structure methods. To the best of our knowledge, these are the first computed ROA spectra for a periodic bulk phase system. As an example, the experimental ROA spectrum of liquid (R)-propylene oxide is reproduced very well.

Correlations in the Solute-Solvent Dynamics Reach Beyond the First Hydration Shell of Ions

While the real-space structure of solvation shells has been explored for decades, a dynamical perspective that directly relies on changes in the H-bond network became accessible more recently mainly via far-infrared (THz) spectroscopies. A remaining key question is how many hydration shells are affected by ion-induced network perturbations. We disclose that theoretical THz difference spectra of aqueous salt solutions can be deciphered in terms of only a handful of dipolar auto- and cross-correlations, including the second solvation shell. This emphasizes the importance of cross-correlations being often neglected in multicomponent models. Analogously, experimental THz responses of simple ions can be deciphered in a similar way. Dramatic intensity cancellations due to large positive and negative contributions are found to effectively shift intensity maxima. Thus, THz spectroscopy provides an unprecedented view on the details of hydration dynamics, which can be understood by a combination of experiment and theory.

On-the-Fly ab Initio Semiclassical Calculation of Glycine Vibrational Spectrum

We present an on-the-fly ab initio semiclassical study of vibrational energy levels of glycine, calculated by Fourier transform of the wavepacket correlation function. It is based on a multiple coherent states approach integrated with monodromy matrix regularization for chaotic dynamics. All four lowest-energy glycine conformers are investigated by means of single-trajectory semiclassical spectra obtained upon classical evolution of on-the-fly trajectories with harmonic zero-point energy. For the most stable conformer I, direct dynamics trajectories are also run for each vibrational mode with energy equal to the first harmonic excitation. An analysis of trajectories evolved up to 50 000 atomic time units demonstrates that, in this time span, conformers II and III can be considered as isolated species, while conformers I and IV show a pretty facile interconversion. Therefore, previous perturbative studies based on the assumption of isolated conformers are often reliable but might be not completely appropriate in the case of conformer IV and conformer I for which interconversion occurs promptly.

Local Mode Analysis: Decoding IR Spectra by Visualizing Molecular Details

Integration of experimental and computational approaches to investigate chemical reactions in proteins has proven to be very successful. Experimentally, time-resolved FTIR difference-spectroscopy monitors chemical reactions at atomic detail. To decode detailed structural information encoded in IR spectra, QM/MM calculations are performed. Here, we present a novel method which we call local mode analysis (LMA) for calculating IR spectra and assigning spectral IR-bands on the basis of movements of nuclei and partial charges from just a single QM/MM trajectory. Through LMA the decoding of IR spectra no longer requires several simulations or optimizations. The novel approach correlates the motions of atoms of a single simulation with the corresponding IR bands and provides direct access to the structural information encoded in IR spectra. Either the contributions of a particular atom or atom group to the complete IR spectrum of the molecule are visualized, or an IR-band is selected to visualize the corresponding structural motions. Thus, LMA decodes the detailed information contained in IR spectra and provides an intuitive approach for structural biologists and biochemists. The unique feature of LMA is the bidirectional analysis connecting structural details to spectral features and vice versa spectral features to molecular motions.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Correcting the record: the dimers and trimers of trans-N-methylacetamide

Raman jet spectroscopy reveals three N-methylacetamide molecules organizing into a ring structure, previously overlooked in computations.

Infrared Action Spectroscopy of Low-Temperature Neutral Gas-Phase Molecules of Arbitrary Structure

We demonstrate a technique for IR action spectroscopy that enables measuring IR spectra in a background-free fashion for low-temperature neutral gas-phase molecules of arbitrary structure. The method is exemplified experimentally for N-methylacetamide molecules in the mid-IR spectral range of 1000-1800  cm^{-1}, utilizing the free electron laser FELIX. The technique involves the resonant absorption of multiple mid-IR photons, which induces molecular dissociation. The dissociation products are probed with 10.49 eV vacuum ultraviolet photons and analyzed with a mass spectrometer. We also demonstrate the capability of this method to record, with unprecedented ease, mid-IR spectra for the molecular associates, such as clusters and oligomers, present in a molecular beam. In this way the mass-selected spectra of low-temperature gas-phase dimers and trimers of N-methylacetamide are measured in the full amide I-III range.

A polarizable QM/MM approach to the molecular dynamics of amide groups solvated in water

The infrared (IR) spectra of polypeptides are dominated by the so-called amide bands. Because they originate from the strongly polar and polarizable amide groups (AGs) making up the backbone, their spectral positions sensitively depend on the local electric fields. Aiming at accurate computations of these IR spectra by molecular dynamics (MD) simulations, which derive atomic forces from a hybrid quantum and molecular mechanics (QM/MM) Hamiltonian, here we consider the effects of solvation in bulk liquid water on the amide bands of the AG model compound N-methyl-acetamide (NMA). As QM approach to NMA we choose grid-based density functional theory (DFT). For the surrounding MM water, we develop, largely based on computations, a polarizable molecular mechanics (PMM) model potential called GP6P, which features six Gaussian electrostatic sources (one induced dipole, five static partial charge distributions) and, therefore, avoids spurious distortions of the DFT electron density in hybrid DFT/PMM simulations. Bulk liquid GP6P is shown to have favorable properties at the thermodynamic conditions of the parameterization and beyond. Lennard-Jones (LJ) parameters of the DFT fragment NMA are optimized by comparing radial distribution functions in the surrounding GP6P liquid with reference data obtained from a "first-principles" DFT-MD simulation. Finally, IR spectra of NMA in GP6P water are calculated from extended DFT/PMM-MD trajectories, in which the NMA is treated by three different DFT functionals (BP, BLYP, B3LYP). Method-specific frequency scaling factors are derived from DFT-MD simulations of isolated NMA. The DFT/PMM-MD simulations with GP6P and with the optimized LJ parameters then excellently predict the effects of aqueous solvation and deuteration observed in the IR spectra of NMA. As a result, the methods required to accurately compute such spectra by DFT/PMM-MD also for larger peptides in aqueous solution are now at hand.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

Fundamental frequency from classical molecular dynamics

We give a theoretical validation for calculating fundamental frequencies of a molecule from classical molecular dynamics (MD) when its anharmonicity is small enough to be treated by perturbation theory. We specifically give concrete answers to the following questions: (1) What is the appropriate initial condition of classical MD to calculate the fundamental frequency? (2) From that condition, how accurately can we extract fundamental frequencies of a molecule? (3) What is the benefit of using ab initio MD for frequency calculations? Our analytical approaches to those questions are classical and quantum normal form theories. As numerical examples we perform two types of MD to calculate fundamental frequencies of H2O with MP2/aug-cc-pVTZ: one is based on the quartic force field and the other one is direct ab initio MD, where the potential energies and the gradients are calculated on the fly. From those calculations, we show comparisons of the frequencies from MD with the post vibrational self-consistent field calculations, second- and fourth-order perturbation theories, and experiments. We also apply direct ab initio MD to frequency calculations of C-H vibrational modes of tetracene and naphthalene. We conclude that MD can give the same accuracy in fundamental frequency calculation as second-order perturbation theory but the computational cost is lower for large molecules.