Strategy for Modeling the Infrared Spectra of Ion-Containing Water Drops

Tensorial Properties via the Neuroevolution Potential Framework: Fast Simulation of Infrared and Raman Spectra

Infrared and Raman spectroscopy are widely used for the characterization of gases, liquids, and solids, as the spectra contain a wealth of information concerning, in particular, the dynamics of these systems. Atomic scale simulations can be used to predict such spectra but are often severely limited due to high computational cost or the need for strong approximations that limit the application range and reliability. Here, we introduce a machine learning (ML) accelerated approach that addresses these shortcomings and provides a significant performance boost in terms of data and computational efficiency compared with earlier ML schemes. To this end, we generalize the neuroevolution potential approach to enable the prediction of rank one and two tensors to obtain the tensorial neuroevolution potential (TNEP) scheme. We apply the resulting framework to construct models for the dipole moment, polarizability, and susceptibility of molecules, liquids, and solids and show that our approach compares favorably with several ML models from the literature with respect to accuracy and computational efficiency. Finally, we demonstrate the application of the TNEP approach to the prediction of infrared and Raman spectra of liquid water, a molecule (PTAF-), and a prototypical perovskite with strong anharmonicity (BaZrO3). The TNEP approach is implemented in the free and open source software package gpumd, which makes this methodology readily available to the scientific community.

Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases

Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.

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.

Theoretical Sum Frequency Generation Spectra of Protein Amide with Surface-Specific Velocity-Velocity Correlation Functions

Vibrational sum frequency generation (vSFG) spectroscopy is widely used to probe the protein structure at interfaces. Because protein vSFG spectra are complex, they can only provide detailed structural information if combined with computer simulations of protein molecular dynamics and spectra calculations. We show how vSFG spectra can be accurately modeled using a surface-specific velocity-velocity scheme based on ab initio normal modes. Our calculated vSFG spectra show excellent agreement with the experimental sum frequency spectrum of LTα14 peptide and provide insight into the origin of the characteristic α-helical amide I peak. Analysis indicates that the peak shape can be explained largely by two effects: (1) the uncoupled response of amide groups located on opposite sides of the α-helix will have different orientations with respect to the interface and therefore different local environments affecting the local mode vibrations and (2) vibrational splitting from nearest neighbor coupling evaluated as inter-residue vibrational correlation. The conclusion is consistent with frequency mapping techniques with an empirically based ensemble of peptide structures, thus showing how time correlation approaches and frequency mapping techniques can give independent yet complementary molecular descriptions of protein vSFG. These models reveal the sensitive relationship between protein structure and their amide I response, allowing exploitation of the complicated molecular vibrations and their interference to derive the structures of proteins under native conditions at interfaces.

Manifolds of low energy structures for a magic number of hydrated sulfate: SO42-(H2O)24

Low energy structures of SO₄^2-(H₂O)₂₄ have been obtained using a combination of classical molecular dynamics simulations and refinement of structures and energies by quantum chemical calculations. Extensive exploration of the potential energy surface led to a number of low-energy structures, confirmed by accurate calibration calculations. An overall analysis of this large set was made after devising appropriate structural descriptors such as the numbers of cycles and their combinations. Low energy structures bear common motifs, the most prominent being fused cycles involving alternatively four and six water molecules. The latter adopt specific conformations which ensure the appropriate surface curvature to form a closed cage without dangling O-H bonds and at the same time provide 12-coordination of the sulfate ion. A prominent feature to take into account is isomerism via inversion of hydrogen bond orientations along cycles. This generates large families of ca. 100 isomers for this cluster size, spanning energy windows of 10-30 kJ mol^-1. This relatively ignored isomerism must be taken into account to identify reliably the lowest energy minima. The overall picture is that the magic number cluster SO₄^2-(H₂O)₂₄ does not correspond to formation of a single, remarkable structure, but rather to a manifold of structural families with similar stabilities. Extensive calculations on isomerization mechanisms within a family indicate that large barriers are associated to direct inversion of hydrogen bond networks. Possible implications of these results for magic number clusters of other anions are discussed.

Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches

Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.

Infrared Spectra of Deprotonated Dicarboxylic Acids: IRMPD Spectroscopy and Empirical Valence‐Bond Modeling

Experimental infrared multiple-photon dissociation (IRMPD) spectra recorded for a series of deprotonated dicarboxylic acids, HO₂ (CH₂ )_n CO 2 - (n=2-4), are interpreted using a variety of computational methods. The broad bands centered near 1600 cm^-1 can be reproduced neither by static vibrational calculations based on quantum chemistry nor by a dynamical description of individual structures using the many-body polarizable AMOEBA force field, strongly suggesting that these molecules experience dynamical proton sharing between the two carboxylic ends. To confirm this assumption, AMOEBA was combined with a two-state empirical valence-bond (EVB) model to allow for proton transfer in classical molecular dynamics simulations. Upon suitable parametrization based on ab initio reference data, the EVB-AMOEBA model satisfactorily reproduces the experimental infrared spectra, and the finite temperature dynamics reveals a significant amount of proton sharing in such systems.