Car-Parrinello Molecular Dynamics Study of DCl Hydrate Crystals

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

HCl hydrates as model systems for protonated water

Ab initio molecular dynamics simulations are presented of vibrational dynamics and spectra of crystal HCl hydrates. Depending on the composition, the hydrates include distinct protonated water forms, which in their equilibrium structures approximate either the Eigen ion H3O+(H2O)3 (in the hexahydrate) or the Zundel H2O...H+...OH2 ion (in the di- and trihydrate). Thus, the hydrates offer the opportunity to study spectra and dynamics of distinct species of protonated water trapped in a semirigid solvating environment. The experimentally measured spectra are reproduced quite well by BLYP/DZVP-level calculations employing Fourier transform of the system dipole. The large overall width (800-1000 cm-1) of structured proton bands reflects a broad range of solvating environments generated by crystal vibrations. The aqueous HCl solution was also examined in search of an objective criterion for separating the contributions of "Zundel-like" and "Eigen-like" protonated forms. It is suggested that no such criterion exists since distributions of proton-related structural properties appear continuous and unimodal. Dipole derivatives with respect to OH and O...H+ stretches in water and protonated water were also investigated to advance the understanding of the corresponding IR intensities. The effects of H bonding and solvation on the intensities were analyzed with the help of the Wannier centers' representation of electron density.

Elusive structure of HCl monohydrate

The study addresses the structure of crystalline HCl monohydrate which is composed of H3O+ and Cl-. The published x-ray diffraction patterns indicate an element of disorder, the nature of which is debated in the literature and is addressed in the present study. The computational investigations include searches for alternative crystal structures employing an empirical potential, and on-the-fly simulations as implemented in the density functional code QUICKSTEP employing Gaussian basis sets. The experimental work focuses on Fourier-transform infrared (FTIR) spectra of crystal nanoparticles. Simulations of FTIR spectra and of the x-ray diffraction patterns are consistent with crystal monohydrate structure composed of ferroelectric domains, joined by "boundary tissue" of antiferroelectric structure.

A new glance at HCl-monohydrate spectroscopy, using on-the-fly dynamics

On-the-fly dynamics is used to analyze the remarkably anharmonic infrared spectroscopy of crystalline HCl monohydrate, an ionic solid composed of H3O+ and Cl-. The dominant intense infrared feature is shown to originate from specific sections of the hydronium trajectory, in which one of the H-atoms interacts strongly with a neighboring Cl-.