Hydrogen fluoride phase behavior and molecular structure: Ab initio derived potential models

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.

Bond cooperativity and ring formation in hydrogen fluoride thermodynamic properties: A two-density formalism framework

In this work, we develop a thermodynamic perturbation theory using a two-density formalism framework to model the bond cooperativity effect for associating hard sphere and Lennard-Jones fluids. The theory predictions are compared with Monte Carlo simulation results and they are in excellent agreement. We incorporate bond angle dependent ring formation into the theory to calculate hydrogen fluoride thermodynamic properties. The liquid density and vapor pressure obtained by the theory are in good agreement with the experimental data. Comparing the thermo-physical properties of hydrogen fluoride calculated by this theory with previous studies reveals the importance of bond angle dependent ring formation and cooperative hydrogen bonding to capture its anomalous behavior especially in the vapor phase. The cooperativity ratio obtained in our model is close to the values reported by previous quantum studies.

Correlated atomic motions in liquid deuterium fluoride studied by coherent quasielastic neutron scattering

The collective dynamics of liquid deuterium fluoride are studied by means of high-resolution quasielastic and inelastic neutron scattering over a range of four decades in energy transfer. The spectra show a low-energy coherent quasielastic component which arises from correlated stochastic motions as well as a broad inelastic feature originating from overdamped density oscillations. While these results are at variance with previous works which report on the presence of propagating collective modes, they are fully consistent with neutron diffraction, nuclear magnetic resonance, and infrared/Raman experiments on this prototypical hydrogen-bonded fluid.

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.

On the variation of the structure of liquid deuterium fluoride with temperature

The structure of liquid deuterium fluoride has been measured using pulsed neutron diffraction and high energy x-ray diffraction techniques as a function of temperature. The neutron experiments were performed at T=296+/-2 K, 246+/-2 K, and 193+/-2 K and the x-ray measurements carried out at 296+/-2 K and 195+/-2 K. The x-ray pair correlation functions, which are dominated by fluorine-fluorine interactions, show the first peak at approximately 2.53+/-0.05 A remains very nearly invariant with decreasing temperature. Peaks around 4.5 and 5.0 A also appear at both temperatures in the x-ray data. In contrast, the intermolecular peaks in the total neutron pair correlation function show that significant systematic local structural changes occur as the temperature is lowered. The first intermolecular peak position shortens from 1.64+/-0.05 A at 296 K to 1.56+/-0.05 A at 195 K. Although there are overlapping contributions from the intermolecular hydrogen-fluorine and hydrogen-hydrogen correlations, it is clear that the temperature dependent structural changes are largely due to a rearrangement of the deuterium atom positions in the fluid. By comparison with partial structure factor data the hydrogen bonds appear to become more linear at lower temperatures.

Computational quantum mechanics: An underutilized tool in thermodynamics

In this paper, we highlight the various ways computational quantum mechanics (QM) can be used in applied thermodynamics. We start with the most rigorous procedures of calculating the interactions between molecules that can then be used in simulation and progress, in steps, to less rigorous but easily used methods, including the very successful continuum solvation models.

Evidence of the presence of opticlike collective modes in a liquid from neutron scattering experiments

Inelastic neutron scattering data from liquid DF close to the melting point show, in addition to spectra comprising quasielastic and heavily damped acoustic motions, an intense, nondispersive band centered at about 27 meV along with a broader higher energy feature. Observation of the former band provides the first direct verification of the existence within the liquid state of collective opticlike excitations as predicted by molecular dynamics simulations. The latter corresponds to mainly reorientational motions assigned from mode eigenvector analysis carried out by computer simulations.

Development of a new polarizable potential model of hydrogen fluoride and comparison with other effective models in liquid and supercritical states

Development of a new polarizable potential of hydrogen fluoride through the reparametrization of the JV-P model is presented: The length of the H-F bond has been shortened and the other parameters of the model have been readjusted accordingly. The structural, thermodynamic, and liquid-vapor equilibrium properties of the new model are compared with those of other effective potential models of HF as well as with experimental data in a broad range of thermodynamic states, from near-freezing to supercritical conditions. It is found that although the reparametrization does not change the structural properties of the HF model noticeably at the level of the pair correlations, it improves the reproduction of the thermodynamic properties of hydrogen fluoride over the entire range of existence of a thermodynamically stable liquid phase and also that of the vapor-liquid coexistence curve. However, the new model, which still overestimates the close-contact separation of the HF molecules, underestimates the density of the coexisting liquid phase and overestimates the saturation pressure, probably due to the too steep repulsion of the potential function.

Anomalous corresponding-states surface tension of hydrogen fluoride and of the Onsager model

In a corresponding-states analysis of the liquid-vapor surface tension originally suggested by Guggenheim, we study the behavior of different simple (i.e., nonpolar), polar and ionic fluids. The results are compared to the corresponding ones for model fluids of each of the three types. For simple and weakly polar fluids (both real and model), the data map onto a master curve, as demonstrated by Guggenheim. For strongly dipolar, associating fluids, which also exhibit hydrogen bonding, one finds deviations from the master curve at low temperatures and, thus, observes the characteristic sigmoid behavior of the reduced surface tension as a function of temperature. The same is obtained for the model ionic fluid, the restricted primitive model. Truly exceptionally low values of the reduced surface tension are found for hydrogen fluoride and for the Onsager model of dipolar fluids, the surface tension of which we evaluate using an approximate hypernetted chain relation to obtain the square-gradient term in a modified van der Waals theory. Remarkably, in the corresponding-states plot, the surface tensions of HF and of the Onsager model agree very closely, while being well separated from the values for the other fluids. We also study the gradual transition of a model fluid from a simple fluid to a strongly dipolar one by varying the relative strength of dipolar and dispersion forces.