Ab Initio Molecular Dynamics Study of Hydrochloric Acid in Water

Acid Dissociation in (HX)n (H2 O)n Clusters (X=F, Cl, Br, I; n=2, 3)

In this work, we analyze the interactions between two or three hydrogen halide molecules and the same number of water moieties through a systematic exploration of their potential energy surfaces. Our results indicate that the most stable HF and HCl aggregates do not experience dissociation of any of the acid fragments, even with three water molecules. In contrast, in the HBr and HI clusters, one of the acid fragments does dissociate. While the global minimum of (HBr)₃ (H₂ O)₃ is a hydrogen-bridged bihalide anion (BrHBr^- ), which is persistent at temperatures up to 203 K, the lowest energy structure of (HI)₃ (H₂ O)₃ has a separated ion pair, but the motif with a bihalide anion (IHI^- ) is only 0.2 kcal mol^-1 above the global minimum. Among the more stable structures is a broad spectrum of contacts, including water⋯water, HX⋯water, and HX⋯HX hydrogen bonds, halogen bonds, ionic and long-range X⋯H contacts.

Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective

As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.

Catalytic Oxidation of Chlorinated Organics over Lanthanide Perovskites: Effects of Phosphoric Acid Etching and Water Vapor on Chlorine Desorption Behavior

In this article, the underlying effect of phosphoric acid etching and additional water vapor on chlorine desorption behavior over a model catalyst La₃Mn₂O₇ was explored. Acid treatment led to the formation of LaPO₄ and enhanced the mobility of lattice oxygen of La₃Mn₂O₇ evidenced by a range of characterization (i.e., X-ray diffraction, temperature-programmed analyses, NH₃-IR, etc.). The former introduced thermally stable Brönsted acidic sites that enhanced dichloromethane (DCM) hydrolysis while the latter facilitated desorption of accumulated chlorine at elevated temperatures. The acid-modified catalyst displayed a superior catalytic activity in DCM oxidation compared to the untreated sample, which was ascribed to the abundance of proton donors and Mn(IV) species. The addition of water vapor to the reaction favored the formation and desorption of HCl and avoided surface chlorination at low temperatures. This resulted in a further reduction in reaction temperature under humid conditions ( T₉₀ of 380 °C for the modified catalyst). These results provide an in-depth interpretation of chlorine desorption behavior for DCM oxidation, which should aid the future design of industrial catalysts for the durable catalytic combustion of chlorinated organics.

Insights into acid dissociation of HCl and HBr with internal electric fields

A critical electric field exerted by the solvent on the ionizable group leads to acid dissociation.

How Many Water Molecules Does it Take to Dissociate HCl?

The potential energy surfaces of the HCl(H2O)n (n is the number of water molecules) clusters are systematically explored using density functional theory and high-level ab initio computations. On the basis of electronic energies, the number of water molecules needed for HCl dissociation is four as reported by some experimental groups. However, this number is five owing to the inclusion of entropic factors. Wiberg bond indices are calculated and analyzed, and the results provide a quadratic correlation and classification of clusters according to the nondissociated, partially dissociated, and fully dissociated character of the H-Cl bond. Our computations show that if temperature is not controlled during the experiment, the values obtained for the dipole moment (or for any measurable property) are susceptible to change, providing a different picture of the number of water molecules needed for HCl dissociation in a nanoscopic droplet.

Aqueous solutions: state of the art in ab initio molecular dynamics

The simulation of liquids by ab initio molecular dynamics (AIMD) has been a subject of intense activity over the last two decades. The significant increase in computational resources as well as the development of new and efficient algorithms has elevated this method to the status of a standard quantum mechanical tool that is used by both experimentalists and theoreticians. As AIMD computes the electronic structure from first principles, it is free of ad hoc parametrizations and has thus been applied to a large variety of physical and chemical problems. In particular, AIMD has provided microscopic insight into the structural and dynamical properties of aqueous solutions which are often challenging to probe experimentally. In this review, after a brief theoretical description of the Born-Oppenheimer and Car-Parrinello molecular dynamics formalisms, we show how AIMD has enhanced our understanding of the properties of liquid water and its constituent ions: the proton and the hydroxide ion. Thereafter, a broad overview of the application of AIMD to other aqueous systems, such as solvated organic molecules and inorganic ions, is presented. We also briefly describe the latest theoretical developments made in AIMD, such as methods for enhanced sampling and the inclusion of nuclear quantum effects.

Hydrated alizarin complexes: hydrogen bonding and proton transfer

We investigated the hydrogen bonding structures and proton transfer for the hydration complexes of alizarin (Az) produced in a supersonic jet using fluorescence excitation (FE), dispersed laser induced fluorescence (LIF), visible-visible hole burning (HB), and fluorescence detected infrared (FDIR) spectroscopy. The FDIR spectrum of bare Az with two O-H groups exhibits two vibrational bands at 3092 and 3579 cm(-1), which, respectively, correspond to the stretching vibration of O1-H1 that forms a strong intramolecular hydrogen bond with the C9=O9 carbonyl group and the stretching vibration of O2-H2 that is weakly hydrogen-bonded to O1-H1. For the 1:1 hydration complex Az(H(2)O)(1), we identified three conformers. In the most stable conformer, the water molecule forms hydrogen bonds with the O1-H1 and O2-H2 groups of Az as a proton donor and proton acceptor, respectively. In the other conformers, the water binds to the C10=O10 group in two nearly isoenergetic configurations. In contrast to the sharp vibronic peaks in the FE spectra of Az and Az(H(2)O)(1), only broad, structureless absorption was observed for Az(H(2)O)(n) (n≥ 2), indicating a facile decay process, possibly due to proton transfer in the electronic excited state. The FDIR spectrum with the wavelength of the probe laser fixed at the broad band exhibited a broad vibrational band near the O2-H2 stretching vibration frequency of the most stable conformer of Az(H(2)O)(1). With the help of theoretical calculations, we suggest that the broad vibrational band may represent the occurrence of proton transfer by tunnelling in the electronic ground state of Az(H(2)O)(n) (n≥ 2) upon excitation of the O2-H2 vibration.

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.

Real-world predictions from ab initio molecular dynamics simulations

In this review we present the techniques of ab initio molecular dynamics simulation improved to its current stage where the analysis of existing processes and the prediction of further chemical features and real-world processes are feasible. For this reason we describe the relevant developments in ab initio molecular dynamics leading to this stage. Among them, parallel implementations, different basis set functions, density functionals, and van der Waals corrections are reported. The chemical features accessible through AIMD are discussed. These are IR, NMR, as well as EXAFS spectra, sampling methods like metadynamics and others, Wannier functions, dipole moments of molecules in condensed phase, and many other properties. Electrochemical reactions investigated by ab initio molecular dynamics methods in solution, on surfaces as well as complex interfaces, are also presented.

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.

Dynamical properties of liquid water from ab initio molecular dynamics performed in the complete basis set limit

Dynamical properties of liquid water were studied using Car-Parrinello [Phys. Rev. Lett. 55, 2471 (1985)] ab initio molecular dynamics (AIMD) simulations within the Kohn-Sham (KS) density functional theory employing the Becke-Lee-Yang-Parr exchange-correlation functional for the electronic structure. The KS orbitals were expanded in a discrete variable representation basis set, wherein the complete basis set limit can be easily reached and which, therefore, provides complete convergence of ionic forces. In order to minimize possible nonergodic behavior of the simulated water system in a constant energy (NVE) ensemble, a long equilibration run (30 ps) preceded a 60 ps long production run. The temperature drift during the entire 60 ps trajectory was found to be minimal. The diffusion coefficient [0.055 A2/ps] obtained from the present work for 32 D2O molecules is a factor of 4 smaller than the most up to date experimental value, but significantly larger than those of other recent AIMD studies. Adjusting the experimental result so as to match the finite-sized system used in the present study brings the comparison between theory and experiment to within a factor of 3. More importantly, the system is not observed to become "glassy" as has been reported in previous AIMD studies. The computed infrared spectrum is in good agreement with experimental data, especially in the low frequency regime where the translational and librational motions of water are manifested. The long simulation length also made it possible to perform detailed studies of hydrogen bond dynamics. The relaxation dynamics of hydrogen bonds observed in the present AIMD simulation is slower than those of popular force fields, such as the TIP4P potential, but comparable to that of the TIP5P potential.

Structure of liquid water at ambient temperature fromab initiomolecular dynamics performed in the complete basis set limit

Structural properties of liquid water at ambient temperature were studied using Car-Parrinello [Phys. Rev. Lett. 55, 2471 (1985)] ab initio molecular dynamics (CPAIMD) simulations combined with the Kohn-Sham (KS) density functional theory and the BLYP exchange-correlation functional for the electronic structure. Unlike other recent work on the same subject, where plane-wave (PW) or hybrid Gaussian/plane-wave basis sets were employed, in the present paper, a discrete variable representation (DVR) basis set is used to expand the KS orbitals, so that with the real-space grid adapted in the present work, the properties of liquid water could be obtained very near the complete basis set limit. Structural properties of liquid water were extracted from a 30 ps CPAIMD-BLYP/DVR trajectory at 300 K. The radial distribution functions (RDFs), spatial distribution functions, and hydrogen bond geometry obtained from the CPAIMD-BLYP/DVR simulation are generally in good agreement with the most up to date experimental measurements. Compared to recent ab initio MD simulations based on PW basis sets, less significant overstructuring was found in the RDFs and the distributions of hydrogen bond angles, suggesting that previous plane-wave and Gaussian basis set calculations have exaggerated the tendency toward overstructuring.

Structure and Dynamics of Br- Ion in Liquid Methanol

A Car-Parrinello molecular dynamics simulation has been performed on a solution of Br- in liquid methanol analyzing with particular attention charge transfer and polarization effects. The first solvation shell has been characterized in terms of H-bonds, and it has been found that the high polarization of the bromide gives rise to a stable solvent cage. The differences in the coordination number with the chloride can be ascribed to the ionic radius and to the stronger perturbations brought by the solvent to the bromide ion.

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.

Solvation Dynamics of Li+ and Cl- Ions in Liquid Methanol

Car-Parrinello molecular dynamics simulations have been performed on Li(+) and Cl(-) in fully deuterated liquid methanol. The results have been compared with available experimental and theoretical data. It has been found that the lithium cation has a stable tetrahedral coordination, whereas the chloride anion presents an average coordination number of 3.56. The polarization effects induced by the ion on the solvent have been analyzed in terms of Wannier function centers. Particular attention has been devoted to the charge transfer, which is particularly important in these types of systems. Evidence for the stability of the lithium cation solvent cage also has been found in the vibrational spectra.

Ions in water: The microscopic structure of a concentrated HCl solution

A neutron diffraction experiment with isotopic H/D substitution on a concentrated HCl/H2O solution is presented. The full set of partial structure factors is extracted, by combining the diffraction data with a Monte Carlo simulation. This allows us to investigate both the changes of the water structure in the presence of ions and their solvation shell, overcoming the limitations of standard diffraction experiments. It is found that the interaction with the solutes affects the tetrahedral network of hydrogen bonded water molecules, in a manner similar to the application of an external pressure to pure water, although HCl seems less effective than other solutes, such as NaOH, at the same concentration. Consistent with experimental and theoretical data, the number of water molecules in the solution is not sufficient to completely dissociate the acid molecule. As a consequence, both dissociated H+ and Cl- ions and undissociated HCl molecules coexist in the sample, and this mixture is correctly reproduced in the simulation box. In particular, the hydrated H+ ions, forming a H3O+ complex, participate in three strong and short hydrogen bonds, while a well-defined hydration shell is found around the chlorine ion. These results are not consistent with the findings of early diffraction experiments on the same system and could only be obtained by combining high quality experimental data with a proper computer simulation.