A theoretical study on the size-dependence of ground-state proton transfer in phenol-ammonia clusters

Geometries and infrared (IR) spectra in the mid-IR region of phenol-(ammonia)_n (PhOH-(NH₃)_n) (n = 0-10) clusters have been studied using density functional theory (DFT) to investigate the critical number of solvent molecules necessary to promote ground-state proton transfer (GSPT). For n ≤ 8 clusters, the most stable isomer is a non-proton-transferred (non-PT) structure, and all isomers found within 1.5 kcal mol^-1 from it are also non-PT structures. For n = 9, the most stable isomer is also a non-PT structure; however, the second stable isomer is a PT structure, whose relative energy is within the experimental criterion of population (0.7 kcal mol^-1). For n = 10, the PT structure is the most stable one. We can therefore estimate that the critical size of GSPT is n = 9. This is confirmed by the fact that these calculated IR spectra are in good accordance with our previous experimental results of mid-IR spectra. It is demonstrated that characteristic changes of the ν_9a and ν₁₂ bands in the skeletal vibrational region provide clear information that the GSPT reaction has occurred. It was also found that the shortest distance between the π-ring and the solvent moiety is a good indicator of the PT reaction.

Acid Dissociation in HCl-Water Clusters is Temperature Dependent and Cannot be Detected Based on Dipole Moments

The dissociation of acids in aqueous environments at low temperatures in the presence of a limited amount of water is underlying a wealth of processes from atmospheric to interstellar science. For the paradigmatic case of HCl(H_{2}O)_{n} clusters, our extensive ab initio path integral simulations quantify in terms of free energy differences and barriers that n=4 water molecules are indeed required to dissociate HCl at low temperatures. Increasing the temperature, however, reverses the process and thus counteracts dissociation by fluctuation-driven recombination. The size of the electric dipole moment is shown to not correlate with the acid being in its dissociated or molecular state, thus rendering its measurement as a function of n unable to detect the dissociation transition.

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.

Aggregation-induced chemical reactions: acid dissociation in growing water clusters

Understanding chemical reactivity at ultracold conditions, thus enabling molecular syntheses via interstellar and atmospheric processes, is a key issue in cryochemistry. In particular, acid dissociation and proton transfer reactions are ubiquitous in aqueous microsolvation environments. Here, the full dissociation of a HCl molecule upon stepwise solvation by a small number of water molecules at low temperatures, as relevant to helium nanodroplet isolation (HENDI) spectroscopy, is analyzed in mechanistic detail. It is found that upon successive aggregation of HCl with H(2)O molecules, a series of cyclic heteromolecular structures, up to and including HCl(H(2)O)(3), are initially obtained before a precursor state for dissociation, HCl(H(2)O)(3)···H(2)O, is observed upon addition of a fourth water molecule. The latter partially aggregated structure can be viewed as an "activated species", which readily leads to dissociation of HCl and to the formation of a solvent-shared ion pair, H(3)O(+)(H(2)O)(3)Cl(-). Overall, the process is mostly downhill in potential energy, and, in addition, small remaining barriers are overcome by using kinetic energy released as a result of forming hydrogen bonds due to aggregation. The associated barrier is not ruled by thermal equilibrium but is generated by athermal non-equilibrium dynamics. These "aggregation-induced chemical reactions" are expected to be of broad relevance to chemistry at ultralow temperature much beyond HENDI spectroscopy.