Liquid-Phase Activity Coefficients for Saturated HF/H2O Mixtures with Vapor-Phase Nonidealities Described by Molecular Simulation

New Molecular-Mechanics Model for Simulations of Hydrogen Fluoride in Chemistry and Biology

Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF)_n clusters, for n = 3-7; various clusters of HF and H₂O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H₂O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.

Association patterns in (HF)(m)(H2O)(n) (m + n = 2-8) clusters

In an attempt to understand the phase behavior of aqueous hydrogen fluoride, the clustering in the mixture is investigated at the molecular level. The study is performed at the mPW1B95/6-31+G(d,p) level of theory. Several previous studies attempted to describe the dissociation of HF in water, but in this investigation, the focus is only on the association patterns that are present in this binary mixture. A total of 214 optimized geometries of (HF)n(H2O)m clusters, with m + n as high as 8, were investigated. For each cluster combination, several different conformations are investigated, and the preferred conformations are presented. Using multiple linear regressions, the average strengths of the four possible H-bonding interactions are obtained. The strongest H-bond interaction is reported to be the H2O...H-F interaction. The most probable distributions of mixed clusters as a function of composition are also deduced. It is found that the larger (HF)n(H2O)m clusters are favored both energetically and entropically compared to the ones that are of size m + n < or = 3. Also, the clusters with equimolar contributions of HF and H2O are found to have the strongest interactions.