Infrared spectroscopy in the gas and liquid phase from first principle molecular dynamics simulations: application to small peptides

After 50 Years of Vibrational Circular Dichroism Spectroscopy: Challenges and Opportunities of Increasingly Accurate and Complex Experiments and Computations

VCD research continues to thrive, driven by ongoing experimental and theoretical advances. Modern studies deal with increasingly complex samples featuring weak intermolecular interactions and shallow potential energy surfaces. Likewise, the combination of VCD measurements with, for instance, cryo-spectroscopic techniques has significantly increased their sensitivity. The extent to which such modern measurements enhance the informative value of VCD depends significantly on the quality of the theoretical models, which must adequately account for anharmonicity, solvation and molecular dynamics. We herein discuss how experimental advancements engage in a stimulating interplay with recent theoretical developments, pursuing either the static or the dynamic computational route. Both paths have their own strengths and limitations, each addressing fundamentally different problems. We give an outlook on future challenges of VCD research, including the possibility to combine static and dynamic approaches to obtain a full picture of the sample.

Boosting the Modeling of Infrared and Raman Spectra of Bulk Phase Chromophores with Machine Learning

In the field of vibrational spectroscopy simulation, hybrid approximations to Kohn-Sham density-functional theory (KS-DFT) are often considered computationally prohibitive due to the significant effort required to evaluate the exchange-correlation potential in planewave codes. In this Letter, we show that by leveraging the porting of KS-DFT on GPU and incorporating machine-learning techniques, simulating IR and Raman spectra of real-life chromophores in bulk aqueous solution becomes a routine application at this level of theory.

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.

A Critical Assessment on Calculating Vibrational Spectra in Nanostructured Materials

Vibrational spectroscopy is an omnipresent spectroscopic technique to characterize functional nanostructured materials such as zeolites, metal-organic frameworks (MOFs), and metal-halide perovskites (MHPs). The resulting experimental spectra are usually complex, with both low-frequency framework modes and high-frequency functional group vibrations. Therefore, theoretically calculated spectra are often an essential element to elucidate the vibrational fingerprint. In principle, there are two possible approaches to calculate vibrational spectra: (i) a static approach that approximates the potential energy surface (PES) as a set of independent harmonic oscillators and (ii) a dynamic approach that explicitly samples the PES around equilibrium by integrating Newton's equations of motions. The dynamic approach considers anharmonic and temperature effects and provides a more genuine representation of materials at true operating conditions; however, such simulations come at a substantially increased computational cost. This is certainly true when forces and energy evaluations are performed at the quantum mechanical level. Molecular dynamics (MD) techniques have become more established within the field of computational chemistry. Yet, for the prediction of infrared (IR) and Raman spectra of nanostructured materials, their usage has been less explored and remain restricted to some isolated successes. Therefore, it is currently not a priori clear which methodology should be used to accurately predict vibrational spectra for a given system. A comprehensive comparative study between various theoretical methods and experimental spectra for a broad set of nanostructured materials is so far lacking. To fill this gap, we herein present a concise overview on which methodology is suited to accurately predict vibrational spectra for a broad range of nanostructured materials and formulate a series of theoretical guidelines to this purpose. To this end, four different case studies are considered, each treating a particular material aspect, namely breathing in flexible MOFs, characterization of defects in the rigid MOF UiO-66, anharmonic vibrations in the metal-halide perovskite CsPbBr3, and guest adsorption on the pores of the zeolite H-SSZ-13. For all four materials, in their guest- and defect-free state and at sufficiently low temperatures, both the static and dynamic approach yield qualitatively similar spectra in agreement with experimental results. When the temperature is increased, the harmonic approximation starts to fail for CsPbBr3 due to the presence of anharmonic phonon modes. Also, the spectroscopic fingerprints of defects and guest species are insufficiently well predicted by a simple harmonic model. Both phenomena flatten the potential energy surface (PES), which facilitates the transitions between metastable states, necessitating dynamic sampling. On the basis of the four case studies treated in this Review, we can propose the following theoretical guidelines to simulate accurate vibrational spectra of functional solid-state materials: (i) For nanostructured crystalline framework materials at low temperature, insights into the lattice dynamics can be obtained using a static approach relying on a few points on the PES and an independent set of harmonic oscillators. (ii) When the material is evaluated at higher temperatures or when additional complexity enters the system, e.g., strong anharmonicity, defects, or guest species, the harmonic regime breaks down and dynamic sampling is required for a correct prediction of the phonon spectrum. These guidelines and their illustrations for prototype material classes can help experimental and theoretical researchers to enhance the knowledge obtained from a lattice dynamics study.

Conformational preferences of the flexible galactofuranose sugar in gas-phase

In contrast with the predominant pyranose form of galactose, galactofuranose is known to be highly flexible. Such flexibility poses a remarkable challenge in terms of structural studies, thus hindering the in depth understanding of the structure/function relationship in this rare sugar. A thorough computational study based on molecular dynamics and density functional theory supported by vibrational spectroscopy in the gas phase was carried out to provide a better understanding of the instrinsic conformational preferences of galactofuranose. Based on energetic and spectroscopic criteria, we report a subtantially reduced conformational landscape: methyl α-D-galactofuranose adopts E₂/¹E conformations and methyl β-D-galactofuranose adopts ¹T₂/¹E conformations.

Influence of the environment on the infrared spectrum of alanine: An effective mode analysis

The vibrational spectrum of the alanine amino acid was computationally determined in the infrared range 1000-2000 cm^-1, under various environments encompassing the gas, hydrated, and crystalline phases, by means of classical molecular dynamics trajectories, carried out with the Atomic Multipole Optimized Energetics for Biomolecular Simulation polarizable force field. An effective mode analysis was performed, in which the spectra are optimally decomposed into different absorption bands arising from well-defined internal modes. In the gas phase, this analysis allows us to unravel the significant differences between the spectra obtained for the neutral and zwitterionic forms of alanine. In condensed phases, the method provides invaluable insight into the molecular origins of the vibrational bands and further shows that peaks with similar positions can be traced to rather different molecular motions.

How Crystal Symmetry Dictates Non-Local Vibrational Circular Dichroism in the Solid State

Solid-State Vibrational Circular Dichroism (VCD) can be used to determine the absolute structure of chiral crystals, but its interpretation remains a challenge in modern spectroscopy. In this work, we investigate the effect of a twofold screw axis on the solid-state VCD spectrum in a combined experimental and theoretical analysis of P2₁ crystals of (S)-(+)-1-indanol. Even though the space group is achiral, a single proper symmetry operation has an important impact on the VCD spectrum, which reflects the supramolecular chirality of the crystal. Distinguishing between contributions originating from molecular chirality and from chiral crystal packing, we find that while IR absorption hardly depends on the symmetry of the space group, the situation is different for VCD, where completely new non-local patterns emerge. Understanding the two underlying mechanisms, namely gauge transport and direct coupling, will help to use VCD to distinguish polymorphic forms.

Molecular Dynamics of Linear Molecules in Strong Magnetic Fields

Molecular rotations and vibrations have been extensively studied by chemists for decades, both experimentally using spectroscopic methods and theoretically with the help of quantum chemistry. However, the theoretical investigation of molecular rotations and vibrations in strong magnetic fields requires computationally more demanding tools. As such, proper calculations of rotational and vibrational spectra were not feasible up until very recently. In this work, we present rotational and vibrational spectra for two small linear molecules, H2 and LiH, in strong magnetic fields. By treating the nuclei as classical particles, trajectories for rotations and vibrations are simulated from ab initio molecular dynamics. Born-Oppenheimer potential energy surfaces are calculated at the Hartree-Fock and MP2 levels of theory, using London atomic orbitals to ensure gauge origin invariance. For the calculation of nuclear trajectories, a highly efficient Tajima propagator is introduced, incorporating the Berry curvature tensor accounting for the screening of nuclear charges.

Some opinions on MD-based vibrational spectroscopy of gas phase molecules and their assembly: An overview of what has been achieved and where to go

We hereby review molecular dynamics simulations for anharmonic gas phase spectroscopy and provide some of our opinions of where the field is heading. With these new directions, the theoretical IR/Raman spectroscopy of large (bio)-molecular systems will be more easily achievable over longer time-scale MD trajectories for an increase in accuracy of the MD-IR and MD-Raman calculated spectra. With the new directions presented here, the high throughput 'decoding' of experimental IR/Raman spectra into 3D-structures should thus be possible, hence advancing e.g. the field of MS-IR for structural characterization by spectroscopy. We also review the assignment of vibrational spectra in terms of anharmonic molecular modes from the MD trajectories, and especially introduce our recent developments based on Graph Theory algorithms. Graph Theory algorithmic is also introduced in this review for the identification of the molecular 3D-structures sampled over MD trajectories.

Direct observation of the solvent organization and nuclear vibrations of [Ru(dcbpy)2(NCS)2]4-, [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], via ab initio molecular dynamics

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)₂(NCS)₂]^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₂ 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.

The important role of non-covalent interactions for the vibrational circular dichroism of lactic acid in aqueous solution

We present a computational study of vibrational circular dichroism (VCD) in solutions of (S)-lactic acid, relying on ab initio molecular dynamics (AIMD) and full solvation with bulk water. We discuss the effect of the hydrogen bond network on the aggregation behaviour of the acid: while aggregates of the solute represent conditions encountered in a weakly interacting solvent, the presence of water drastically interferes with the clusters - more strongly than originally anticipated. For both scenarios we computed the VCD spectra by means of nuclear velocity perturbation theory (NVPT). The comparison with experimental data allows us to establish a VCD-structure relationship that includes the solvent network around the chiral solute. We suggest that fundamental modes with strong polarisation such as the carbonyl stretching vibration can borrow VCD from the chirally restructured solvent cage, which extends the common explanatory models of VCD generation in aqueous solution.

Computation of Solid-State Vibrational Circular Dichroism in the Periodic Gauge

We introduce a new theoretical formalism to compute solid-state vibrational circular dichroism (VCD) spectra from molecular dynamics simulations. Having solved the origin-dependence problem of the periodic magnetic gauge, we present IR and VCD spectra of (1S,2S)-trans-1,2-cyclohexanediol obtained from first-principles molecular dynamics calculations and nuclear velocity perturbation theory, along with the experimental results. Because the structure model imposes periodic boundary conditions, the common origin of the rotational strength has hitherto been ill-defined and was approximated by means of averaging multiple origins. The new formalism reconnects the periodic model with the finite physical system and restores gauge freedom. It nevertheless fully accounts for nonlocal spatial couplings from the gauge transport term. We show that even for small simulation cells the rich nature of solid-state VCD spectra found in experiments can be reproduced to a very satisfactory level.

Ab Initio molecular dynamics with screened Lorentz forces. II. Efficient propagators and rovibrational spectra in strong magnetic fields

Strong magnetic fields have a large impact on the dynamics of molecules. In addition to the changes in the electronic structure, the nuclei are exposed to the Lorentz force with the magnetic field being screened by the electrons. In this work, we explore these effects using ab initio molecular dynamics simulations based on an effective Hamiltonian calculated at the Hartree-Fock level of theory. To correctly include these non-conservative forces in the dynamics, we have designed a series of novel propagators that show both good efficiency and stability in test cases. As a first application, we analyze simulations of He and H₂ at two field strengths characteristic of magnetic white dwarfs (0.1 B₀ = 2.35 × 10⁴ T and B₀ = 2.35 × 10⁵ T). While the He simulations clearly demonstrate the importance of electron screening of the Lorentz force in the dynamics, the extracted rovibrational spectra of H₂ reveal a number of fascinating features not observed in the field-free case: couplings of rotations/vibrations with the cyclotron rotation, overtones with unusual selection rules, and hindered rotations that transmute into librations with increasing field strength. We conclude that our presented framework is a powerful tool to investigate molecules in these extreme environments.

Probing anharmonic phonons by quantum correlators: A path integral approach

We devise an efficient scheme to determine vibrational properties from Path Integral Molecular Dynamics (PIMD) simulations. The method is based on zero-time Kubo-transformed correlation functions and captures the anharmonicity of the potential due to both temperature and quantum effects. Using analytical derivations and numerical calculations on toy-model potentials, we show that two different estimators built upon PIMD correlation functions fully characterize the phonon spectra and the anharmonicity strength. The first estimator is associated with the force-force quantum correlators and, in the weak anharmonic regime, yields reliable zero-point motion frequencies and thermodynamic properties of the quantum system. The second one is instead connected to displacement-displacement correlators and accurately probes the lowest-energy phonon excitations, regardless of the anharmonicity strength of the system. We also prove that the use of generalized eigenvalue equations, in place of the standard normal mode equations, leads to a significant speed-up in the PIMD phonon calculations, both in terms of faster convergence rate and smaller time step bias. Within this framework, using ab initio PIMD simulations, we compute phonon dispersions of diamond and of the high-pressure I4₁/amd phase of atomic hydrogen. We find that in the latter case, the anharmonicity is stronger than previously estimated and yields a sizeable red-shift in the vibrational spectrum of atomic hydrogen.

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole-dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

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.

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.

Vibrational Spectra of Zeolite Y as a Function of Ion Exchange

Zeolite Y is one of the earliest known and most widely used synthetic zeolites. Many experimental investigations verify the valuable ion exchange capability of this zeolite. In this study, we assessed the effects of ion exchange on its vibrational spectra. We applied classical lattice dynamics methods for IR and Raman intensity calculations. Computed spectra of optimized zeolite Y structures with different cations were compared with experimental data. The spectra obtained in this study are in agreement with previous experimental and computational studies on zeolites from the faujasite group.

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.

Time-Resolved Vibrational Analysis of Excited State Ab Initio Molecular Dynamics to Understand Photorelaxation: The Case of the Pyranine Photoacid in Aqueous Solution

We present a novel time-resolved vibrational analysis for studying photoinduced nuclear relaxation. Generalized modes velocities are defined from ab initio molecular dynamics and wavelet transformed, providing the time localization of vibrational signals in the electronic excited state. The photoexcited pyranine in aqueous solution is presented as a case study. The transient and sequential activation of the simulated vibrational signals is in good agreement with vibrational dynamics obtained from femtosecond stimulated Raman spectroscopy data.

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.

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.

Improving Vibrational Mode Interpretation Using Bayesian Regression

To streamline the interpretation of vibrational spectra, this work introduces the use of Bayesian linear regression with automatic relevance determination as a viable approach to decompose the atomic motions along any vibrational mode as a weighted combination of displacements along chemically meaningful internal coordinates. This novel approach denominated vibrational mode automatic relevance determination (VMARD) is presented and compared with the well-established potential energy decomposition (PED) scheme. Good agreement is generally attained between the two methods. VMARD returns a decomposition of the atomic displacement using only a small number of internal coordinates, thus aiding the interpretation of the vibrational spectra. Moreover, the results show that the VMARD descriptions are resilient toward the addition of additional internal coordinates, achieving a concise description of the vibrational modes despite the use of redundant internal coordinates. Potential applications of VMARD involving the gathering of physical insights on the atomic motions along the reaction coordinate at transition state structures, as well as the improvement of theoretically predicted vibrational frequencies, are also presented under a proof-of-concept perspective.

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.

Structure and dynamics of propylammonium nitrate-acetonitrile mixtures: An intricate multi-scale system probed with experimental and theoretical techniques

In this article, we report the study of structural and dynamical properties for a series of acetonitrile/propylammonium nitrate mixtures as a function of their composition. These systems display an unusual increase in intensity in their X-ray diffraction patterns in the low-q regime, and their ¹H-NMR diffusion-ordered NMR spectroscopy (DOSY) spectra display unusual diffusivities. However, the magnitude of both phenomena for mixtures of propylammonium nitrate is smaller than those observed for ethylammonium nitrate mixtures with the same cosolvent, suggesting that the cation alkyl tail plays an important role in these observations. The experimental X-ray scattering data are compared with the results of molecular dynamics simulations, including both ab initio studies used to interpret short-range interactions and classical simulations to describe longer range interactions. The higher level calculations highlight the presence of a strong hydrogen bond network within the ionic liquid, only slightly perturbed even at high acetonitrile concentration. These strong interactions lead to the symmetry breaking of the NO₃^- vibrations, with a splitting of about 88 cm^-1 in the ν₃ antisymmetric stretch. The classical force field simulations use a greater number of ion pairs, but are not capable of fully describing the longest range interactions, although they do successfully account for the observed concentration trend, and the analysis of the models confirms the nano-inhomogeneity of these kinds of samples.

Application of far-infrared spectroscopy to the structural identification of protein materials

Far-IR spectroscopy was applied to monitor the structure of two types of silk fibroins and the results indicate that they both show several sharp characteristic peaks, which are totally different from those of globular proteins.

Efficient and Accurate Born-Oppenheimer Molecular Dynamics for Large Molecular Systems

An efficient scheme for the calculation of Born-Oppenheimer molecular dynamics (BOMD) simulations is introduced. It combines the corrected small basis set Hartree-Fock (HF-3c) method by Sure and Grimme [J. Comput. Chem. 2013, 43, 1672], extended Lagrangian BOMD (XL-BOMD) by Niklasson et al. [J. Chem. Phys. 2009, 130, 214109], and the calculation of the two electron integrals on graphics processing units (GPUs) [J. Chem. Phys. 2013, 138, 134114; J. Chem. Theory Comput. 2015, 11, 918]. To explore the parallel performance of our strong scaling implementation of the method, we present timings and extract, as its validation and first illustrative application, high-quality vibrational spectra from simulated trajectories of β-carotene, paclitaxel, and liquid water (up to 500 atoms). We conclude that the presented BOMD scheme may be used as a cost-efficient and reliable tool for computing vibrational spectra and thermodynamics of large molecular systems including explicit solvent molecules containing 500 atoms and more. Simulating 50 ps of maitotoxin (nearly 500 atoms) employing time steps of 0.5 fs requires ∼3 weeks on 12 CPUs (Intel Xeon E5 2620 v3) with 24 GPUs (AMD FirePro 3D W8100).

Kinetic theory for DNA melting with vibrational entropy

By treating DNA as a vibrating nonlinear lattice, an activated kinetic theory for DNA melting is developed to capture the breakage of the hydrogen bonds and subsequent softening of torsional and bending vibration modes. With a coarse-grained lattice model, we identify a key bending mode with GHz frequency that replaces the hydrogen vibration modes as the dominant out-of-phase phonon vibration at the transition state. By associating its bending modulus to a universal in-phase bending vibration modulus at equilibrium, we can hence estimate the entropic change in the out-of-phase vibration from near-equilibrium all-atom simulations. This and estimates of torsional and bending entropy changes lead to the first predictive and sequence-dependent theory with good quantitative agreement with experimental data for the activation energy of melting of short DNA molecules without intermediate hairpin structures.

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.

Ab Initio Transient Vibrational Spectral Analysis

Pump probe spectroscopy techniques have enabled the direct observation of a variety of transient molecular species in both ground and excited electronic states. Time-resolved vibrational spectroscopy is becoming an indispensable tool for investigating photoinduced nuclear dynamics of chemical systems of all kinds. On the other hand, a complete picture of the chemical dynamics encoded in these spectra cannot be achieved without a full temporal description of the structural relaxation, including the explicit time-dependence of vibrational coordinates that are substantially displaced from equilibrium by electronic excitation. Here we present a transient vibrational analysis protocol combining ab initio direct molecular dynamics and time-integrated normal modes introduced in this work, relying on the recent development of analytic time-dependent density functional theory (TDDFT) second derivatives for excited states. Prototypical molecules will be used as test cases, showing the evolution of the vibrational signatures that follow electronic excitation. This protocol provides a direct route to assigning the vibrations implicated in the (photo)dynamics of several (photoactive) systems.