Hydrogen bond dynamics and vibrational spectral diffusion in aqueous solution of acetone: A first principles molecular dynamics study#

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.

Dynamics of Ionic Liquid through Intrinsic Vibrational Probes Using the Dispersion-Corrected DFT Functionals

First principles molecular dynamics simulations have been utilized to study the spectral properties of the protic ionic liquid, methylammonium formate (MAF). All simulations were performed using density functional theory (DFT) and various van der Waals-corrected exchange-correlation functionals. We calculated the vibrational stretch frequency distributions, determined the time-frequency correlations of the intrinsic vibrational probes, the N-H and C-O modes in MAF, and the frequency-structure correlations. We also estimated the average hydrogen-bond lifetimes and orientation dynamics to capture the ultrafast spectral response. The spectroscopic signature of the N-H stretching vibrations using the Becke-Lee-Yang-Parr (BLYP) and Perdew-Burke-Ernzerhof (PBE) functionals displays a spectral shift in the lower frequency side, suggesting stronger hydrogen-bonding interactions represented by the gradient approximation functionals than the van der Waals (vdW)-corrected simulations. The carboxylate frequency profiles with the dispersion-corrected representations are almost similar without a significant difference in the normalized distributions. Besides, the COO stretching frequencies at the peak maxima positions of the PBE functionals exhibit a lesser deviation from the experimental data. Spectral diffusion dynamics of the intrinsic vibrational probes on the cationic and anionic sites of the ionic liquid proceed through a short time relaxation of the intact hydrogen bonds followed by an intermediate time constant and a longer time decay indicating the switchover of hydrogen bonds. Dispersion-corrected atom-centered one-electron potential (DCACP) correction added to the BLYP system slows down the picosecond time scales of frequency correlation and the time constants of rotational motion, lengthening the overall system dynamics. The observed trends in the time-dependent decays of frequency fluctuations and the orientation autocorrelation functions correlate with the structural interactions in liquid MAF and hydrogen-bond dynamics. In this study, we examine the predictions made by different density functional treatments comparing the results of the uncorrected BLYP and PBE representations with the semiempirical vdW methods of Grimme and matching our calculated data with the experimental observations.

Conformation-induced vibrational spectral dynamics of hydrogen peroxide and vicinal water molecules

We studied the conformation-induced spectral response of water molecules due to site-specific structural alterations of solvated hydrogen peroxide (H₂O₂) employing DFT-based first principles molecular dynamics (FPMD) simulations. Wavelet transform was used to determine the time-dependent frequencies of the hydroxyls of water molecules and the O-H stretch modes of H₂O₂. Shifts in the vibrational frequency of the hydrogen-bonded hydroxyls inside the solvation shell of H₂O₂ support multiple distinctive hydrogen bonding environments. This paper classifies two distinct hydrogen bond types inside the O-O_W solvation shell of H₂O₂, and the dynamical calculations provide a quantitative estimation of the relative hydrogen bond strength. We ascertain the reason for not observing the escape of water molecules from the hydrogen peroxide hydration shell, unlike the solvation shell of ionic solutions and neutral solutes. Besides, we provide a comprehensive analysis of the spectral shifts in the normalized frequency distribution, the time-dependent decay of frequency-frequency correlation functions, and the hydrogen bond length scale fluctuations. We also quantify the relative contribution of the cisoid and transoid conformers affecting the vibrational spectral signature of the vicinal water molecules. While the transoid conformers promote the hydrogen bonding interactions through the oxygen site (OH_W), the cisoid conformers facilitate hydrogen peroxide-water hydrogen bond formation through the hydrogen site (HO_W). These non-identical hydrogen bond associations stabilize hydrogen peroxide in water.

Decelerated Liquid Dynamics Induced by Component-Dependent Supercooling in Hydrogen and Deuterium Quantum Mixtures

Isotopic mixtures of p-H₂ and o-D₂ molecules have been an attractive binary system because they include two kinds of purely isotopic molecules which possess the same electronic potential but the twice different mass inducing differently pronounced nuclear quantum effects (NQEs). Accessing details of structures and dynamics in such quantum mixtures combining complex molecular dynamics with NQEs of different strengths remains a challenging problem. Taking advantage of the nonempirical molecular dynamics method which describes p-H₂ and o-D₂ molecules, we found that the liquid dynamics slows down at a specific mixing ratio, which can be connected to the observed anomalous slowdown of crystallization in the quantum mixtures. We attributed the decelerated dynamics to the component-dependent supercooling of p-H₂ taking place in the mixtures, demonstrating that there is an optimal mixing ratio to hinder crystallization. The obtained physical insights will help in experimentally controlling and achieving unknown quantum mixtures including superfluid.

Distinctive behavior and two-dimensional vibrational dynamics of water molecules inside glycine solvation shell

We present a first principles molecular dynamics study of a deuterated aqueous solution of a single glycine moiety to explore the structure, dynamics, and two-dimensional infrared spectra of water molecules found in the solvation shell of glycine.

Time-dependent vibrational spectral analysis of first principles trajectory of methylamine with wavelet transform

The fluctuation dynamics of amine stretching frequencies, hydrogen bonds, dangling N-D bonds, and the orientation profile of the amine group of methylamine (MA) were investigated under ambient conditions by means of dispersion-corrected density functional theory-based first principles molecular dynamics (FPMD) simulations. Along with the dynamical properties, various equilibrium properties such as radial distribution function, spatial distribution function, combined radial and angular distribution functions and hydrogen bonding were also calculated. The instantaneous stretching frequencies of amine groups were obtained by wavelet transform of the trajectory obtained from FPMD simulations. The frequency-structure correlation reveals that the amine stretching frequency is weakly correlated with the nearest nitrogen-deuterium distance. The frequency-frequency correlation function has a short time scale of around 110 fs and a longer time scale of about 1.15 ps. It was found that the short time scale originates from the underdamped motion of intact hydrogen bonds of MA pairs. However, the long time scale of the vibrational spectral diffusion of N-D modes is determined by the overall dynamics of hydrogen bonds as well as the dangling ND groups and the inertial rotation of the amine group of the molecule.