On the Quantum Nature of an Excess Proton in Liquid Hydrogen Fluoride

Proton transfer 200 years after von Grotthuss: insights from ab initio simulations

In the last decade, ab initio simulations and especially Car-Parrinello molecular dynamics have significantly contributed to the improvement of our understanding of both the physical and chemical properties of water, ice, and hydrogen-bonded systems in general. At the heart of this family of in silico techniques lies the crucial idea of computing the many-body interactions by solving the electronic structure problem "on the fly" as the simulation proceeds, which circumvents the need for pre-parameterized potential models. In particular, the field of proton transfer in hydrogen-bonded networks greatly benefits from these technical advances. Here, several systems of seemingly quite different nature and of increasing complexity, such as Grotthuss diffusion in water, excited-state proton-transfer in solution, phase transitions in ice, and protonated water networks in the membrane protein bacteriorhodopsin, are discussed in the realms of a unifying viewpoint.

Insights into the Structure and Dynamics of a Room-Temperature Ionic Liquid: Ab Initio Molecular Dynamics Simulation Studies of 1-n-Butyl-3-methylimidazolium Hexafluorophosphate ([bmim][PF6]) and the [bmim][PF6]−CO2 Mixture

Ab initio molecular dynamics (AIMD) studies have been carried out on liquid 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6]) and its mixture with CO2 using the Car-Parrinello molecular dynamics (CPMD) method. Results from AIMD and empirical potential molecular dynamics (MD) have been compared and were found to differ in some respects. With a strong resemblance to the crystal, the AIMD simulated neat liquid exhibits many cation-anion hydrogen bonds, a feature that is almost absent in the MD results. The anions were observed to be strongly polarized in the condensed phase. The addition of CO2 increased the probability of this hydrogen bond formation. CO2 molecules in the vicinity of the ions of [bmim][PF6] exhibit larger deviations from linearity in their instantaneous configurations. The polar environment of the liquid induces a dipole moment in CO2, lifting the degeneracy of its bending mode. The calculated splitting in the vibrational mode compares well with infrared spectroscopic data. The solvation of CO2 in [bmim][PF6] is primarily facilitated by the anion, as seen from the radial and spatial distribution functions. CO2 molecules were found to be aligned tangential to the PF6 spheres with their most probable location being the octahedral voids of the anion. The structural data obtained from AIMD simulations can serve as a benchmark to refine interaction potentials for this important room-temperature ionic liquid.

First Principles and Experimental 1H NMR Signatures of Solvated Ions: The Case of HCl(aq)

A combined experimental and ab initio study is presented of the 1H NMR chemical shift distribution of aqueous hydrogen chloride solution as a function of acid concentration, based on Car-Parrinello molecular dynamics simulations and fully periodic NMR chemical-shift calculations. The agreement of computed and experimental spectra is very good. From first-principles calculations, we can show that the individual contributions of Eigen and Zundel ions, regular water molecules, and the chlorine solvation shell to the NMR line are very distinct and almost independent of the acid concentration. From the computed instantaneous NMR distributions, it is further possible to characterize the average variation in hydrogen-bond strength of the different complexes.

Ab Initio Molecular Dynamics Study of a Mixture of HF(aq) and HCl(aq)†

We have studied a mixture of HF and HCl molecules in water using Car-Parrinello ab initio molecular dynamics (CPMD). We have done simulations with 1 HF and 3 HCl molecules, 3 HF and 4 HCl, 6 HF and 8 HCl (6/8 simulation), and 14 HF molecules. All simulations consist of 32 molecules, and they were 10-96 ps long. The HF dissociation probability was around 30%, and HCl's was more than 90%. The solvation of the HF molecule was much better than the solvation of HCl. The solvation environment of F, both the F- ion and the F in HF, did not depend much on the acids concentration, whereas the Cl coordination numbers were rather sensitive to the concentration. In the 6/8 simulation, all XH-Y (X, Y = F, Cl) type molecules were observed and the FH-F was the most probable. In general, the molecular structures in mixed aqueous acid systems were similar to the pure HF(aq) and HCl(aq) systems.

Ab Initio Molecular Dynamics Simulation of a Water-Hydrogen Fluoride Equimolar Mixture

Hydrogen fluoride and water can be mixed in any proportion. The resulting solutions have unique acidic properties. In particular, hydrogen fluoride undergoes a weak-to-strong acidity transition with increasing concentration of HF To supplement the knowledge already obtained on dilute or moderately concentrated solutions and gas-phase aggregates, an equimolar mixture is studied here by Car-Parrinello molecular dynamics. The natures of the ions and of the complexes formed in the equimolar liquid were determined. Specifically, H3O+, H5O2+, FH-OH2, and HF2- were spontaneously obtained while only hydronium and fluoride ions pre-exist in the equimolar crystal. The behaviour of the proton in the equimolar liquid was compared with mixtures of other proportions simulated previously in an attempt to relate proton dynamics to acidity. In the same way, the behaviour of HF2- was also examined. In this case, proton localization and transfer appeared to be driven by the fluctuating environment of the solvated ion.