The dissociation mechanism and thermodynamic properties of HCl(aq) in hydrothermal fluids (to 700 °C, 60 kbar) by ab initio molecular dynamics simulations
Introduction
Hydrogen chloride (HCl) is one of the most significant volatiles in the Earth’s crust. In aqueous solutions, the neutral HCl(aq) ion pair dissociates to the hydrogen ion (H+) and chloride anion (Cl−) through a reversible chemical reaction:
Both H+ and Cl− are important in controlling mineral solubility and element mobility in natural and man-made hydrothermal fluids. The chloride anion (Cl−) is one of the most abundant anions in hydrothermal fluids, and plays an important role in promoting element transport in hydrothermal fluids by forming stable chloro-complexes with many metals (e.g., Cu+, Co2+, Zn2+, Pb2+, Pd2+, etc.), whereas pH is an essential control of mineral solubility and element speciation in aqueous solutions (as reviewed in Seward and Barnes, 1997, Brugger et al., 2016). A recent study also reveals the importance of the pH of subduction zone fluids in controlling metal and volatile cycles (Galvez et al., 2016). Hence, the knowledge of the dissociation mechanism and thermodynamic properties of HCl(aq) over a wide range of temperature and pressure is crucial for the understanding of element transport and deposition in hydrothermal fluids.
The logarithm of the dissociation constant of HCl(aq) (Eq. (1), logKd) can be expressed as;where a refers to the activities of the subscripted species. As summarized in Table 1, experiments up to 700 °C/4 kbar (one experiment to 4 kbar, and two to 2 kbar) and thermodynamic extrapolations have quantified the dissociation of HCl in aqueous fluids. At low temperature, significant discrepancies in the dissociation constant (up to 6 orders of magnitude) exist between different studies (see Table 1); as a result, some studies assume that HCl is significantly associated at room-temperature in concentrated HCl solutions, whereas other see no significant association. For instance, a few experimental studies obtained the dissociation constant of HCl(aq) by measuring the partial vapor pressure to determine the activity of HCl(aq) () (e.g., Robinson and Bates, 1971, Marsh and McElroy, 1985). In these studies, the equilibrium of HCl(aq) ⇌ HCl(g) was calculated based on Raoult’s Law, yielding a value of logKd (HCl(aq)) around 6.2 at 25 °C. Textbooks (e.g., p. 952 of Greenwood and Earnshaw, 1984) and the IUPAC dataset (p. 46 of Perrin, 1982) list a similar logKd value of 7, making HCl a strong acid with no association at room temperature. However, Clegg and Brimblecombe (1986) pointed out that the approach based on Raoult’s Law is inappropriate for strong acids like HCl, and the Henry’s Law constant is a better approach to describe equilibrium between the vapor pressure of HCl solution and dissociated H+ and Cl− ions.
In contrast to the above, several studies based on potentiometric and conductance data (Table 1) gave values as low as 0.9 (e.g., Johnson and Pytkowicz, 1978, Sretenskaya, 1992), implying significant ion pairing even at room temperature and suggesting that HCl is a relatively weak acid. The extrapolation of Ruaya and Seward (1987)’s AgCl(s) solubility measurements at 100–350 °C also indicated a logKd value of 0.7 at room temperature, and Sverjensky et al. (1991) recommended a logKd value at 0.86 based on a fit of the experimental solubility data for alkali minerals.
Simonson et al. (1990) showed that the discrepancies in logKd of HCl(aq) can be due to the different chemical models employed to calculate the dissociation constants based on the available experimental data. Based on Holmes et al. (1987)’s excess thermodynamic data, Simonson et al. (1990) calculated two sets of the dissociation constants of HCl(aq) using two types of data treatment (as listed in Table 1), and obtained a logKd at 25 °C of 7.6 based on the ion association-interaction model (Pitzer’s model), and 0.2 based on the activity expansion-chemical equilibrium model.
The logKd from different studies generally converge with increasing temperature to mostly overlap at temperatures ≥ 300 °C (Tagirov et al., 1997, Pokrovskii, 1999, Ho et al., 2001). This reflects the increased association of HCl making measurements more achievable. There are still some small discrepancies, in particular between mineral solubility and conductance measurements, and in low-density (<0.6 g cm−3) fluids. For example, there is up to 1 order of magnitude difference among the logKd values reported by Franck, 1956, Frantz and Marshall, 1984 and Sverjensky et al. (1991).
The Helgeson-Kirkham-Flowers (HKF) model has been widely used to fit logKd(HCl(aq)) at high temperature (e.g., Sverjensky et al., 1991, Tagirov et al., 1997, Pokrovskii, 1999). The logKd values extrapolated by the HKF model, based on mineral equilibrium data for K-feldspar-muscovite-quartz (Sverjensky et al., 1991) and AgCl(s) solubility experiments (Tagirov et al., 1997), are within a range of half a log unit at high temperature (300–600 °C, 0.5–2 kbar).
Over the past few decades, ab initio MD has been playing an increasingly significant role in understanding the molecular structures and thermodynamics of aqueous systems over a wide range of T, P, and solution compositions. In particular, the dissociation of acids at room temperature has been investigated quantitatively using such ab initio MD approaches. Trout and Parrinello (1998) first calculated the free energy of H2O dissociation using ab initio distance-constrained thermodynamic integration. Although their results are 16.7 kJ/mol less than the experimental values due to the choice of the density functional and finite size (32 water molecules) of the system due to limited computing power at the time, they first demonstrated the capability of quantitative prediction of thermodynamic properties for aqueous species. Later on, Sprik (2000) calculated the logKd of liquid water at ambient conditions using ab initio thermodynamic integration based on the coordination number of a given oxygen by hydrogen atoms and obtained values that are within one log unit of the experimental value.
In recent years, new techniques have been used to provide a molecular-level understanding of the association of HCl(aq). Chialvo et al., 2002, Chialvo et al., 2003, Chialvo and Simonson, 2007 conducted potential of mean force (PMF) calculations based on classical force-field molecular dynamics simulation to predict the dissociation constants as a function of solution density. They used Lennard-Jones potentials between the oxygen in rigid H2O and H3O+ molecules and the Cl− ion to model intermolecular interactions, and obtained reasonable representations of properties such as dielectric constant of water at near critical conditions and the logKd values as a function of density changes. However, those studies did not provide direct association/dissociation constants as a function of temperature and pressure. Murakhtina et al. (2006) used density functional theory (DFT) based ab initio MD to calculate the H NMR chemical shift of HCl(aq), and showed good agreement with experimental results. Sulpizi and Sprik, 2008, Sulpizi and Sprik, 2010 applied ab initio thermodynamic integration by constraining a dummy proton to calculate the free energy of dissociation of a series of acids; they obtained excellent agreement for logKd (57 °C) values to within 0.4 log units for several small acid molecules (HCl, H2S, formic acid), and within 2 log units for the wide range of acids investigated. They obtained logKd values of 6.7–7.1 for HCl(aq) at 57 °C (330 K), in excellent agreement with the IUPAC value of 7 at 25 °C (Perrin, 1982). Tummanapelli and Vasudevan (2014) calculated pKa values of weak organic acids at a temperature of 57 °C (330 K) using ab initio metadynamics simulations and achieved accurate predictions (within 0.2 log units).
Recently ab initio MD has also been widely used in molecular-level understanding of the geometry of metal complexes and the energetics of ligand exchange reactions relevant for metal transport in hydrothermal fluids (e.g., Cu(I)–Cl−, Sherman, 2007, Mei et al., 2014; U(VI)-Cl−, Bühl et al., 2008; Au(I)–HS−, Liu et al., 2011; Cu(I)–HS−–Cl−, Mei et al., 2013a; Ag(I)–Cl−, Liu et al., 2012, Pokrovski et al., 2013; Ag(I)-HS−/OH−, He et al., 2016; Au(I)–HS−/OH−/S3−, Mei et al., 2013b, Pokrovski et al., 2015; Au(I)–Cl−, Mei et al., 2014, Mei et al., 2017). Using distance-constrained thermodynamic integration, our previous studies successfully predicted the quantitative thermodynamic properties for the formation constants for Cu(I)–Cl−–HS− (Mei et al., 2013a), Zn(II)–Cl− (Mei et al., 2015b), Zn(II)–HS− (Mei et al., 2016) and Pd(II)–Cl−–HS− (Mei et al., 2015a). Two studies have successfully measured acid dissociation constants using ab initio MD simulations under hydrothermal conditions. Liu et al. (2013) obtained agreement to within 2 pKa units for molybdic acid to 300 °C, and Liu et al. (2015) obtained agreements of within 1 unit for pKa1 of As(OH)3(aq) to 300 °C, and for pKa1 to 3 of As(HS)3(aq) at 25 °C (no experimental values exist at T ≥ 25 °C).
Despite all of the experimental and theoretical studies dedicated to HCl association, there is no ab initio MD study of HCl(aq) dissociation at elevated temperature, and the mechanism of the increasing ion pairing/association of HCl(aq) with increasing temperatures has not been clarified. There is also no data for HCl dissociation constants at ultra-high P-T relevant to mantle metasomatism and ultra-high pressure subduction environments (up to 60 kbar, 1000 °C). Under such conditions the role of HCl in controlling pH and metal complexation is unknown.
The present study aims to address these knowledge gaps in order to improve our understanding of HCl(aq) dissociation/association, provide an independent check of the dissociation constants, and provide the first dataset for HCl dissociation at ultra-high P-T fluids. Hence, we conducted ab initio molecular dynamics (MD) simulations to investigate the association mechanism of HCl(aq) as a function of temperature and pressure (25–700 °C, 1 bar to 60 kbar), and then calculated dissociation constants of HCl(aq) based on distance-constraint thermodynamic integration. The dissociation constants obtained from ab initio MD simulations were then used to fit thermodynamic properties using the HKF model and DEW model (Sverjensky et al., 2014). These properties allow estimation of the role of HCl in controlling fluid pH and element mobility in deep earth fluids.
Section snippets
Ab initio molecular dynamics simulations
Ab initio MD simulations were performed using the Car–Parrinello (CP) molecular dynamics code CPMD (Car and Parrinello, 1985). This method implements density functional theory using plane-wave basis sets and pseudo-potentials for the core electrons and the nucleus. The BLYP exchange-correlation functional was employed with a cutoff for the gradient correction of 5 × 10−8 a.u. (Lee et al., 1988, Becke, 1988). A plane-wave cutoff of 80 Ry (1088.46 eV) was used together with Martins–Troullier
Dissociation as a function of temperature from ab initio MD simulations
Unconstrained ab initio MD simulations of HCl(aq) were conducted at 25–700 °C at pressures up to 60 kbar (Table 2) with an initial H–Cl distance of 1.27 Å (Fig. 1). At 25 °C, 1 bar, the HCl molecule dissociated very early in the simulation (in 0.3 ps), and chloride existed as the free Cl– ion during the remaining simulation time of 18.86 ps. Fig. 3a shows the time evolution of the H–Cl distance during the first 2 ps of the simulation. From data in the range 100–250 °C (100 bar), the HCl ion
Dissociation and hydration of HCl in aqueous fluids as a function of temperature and pressure
The ab initio MD simulations conducted in this study revealed the dissociation/association of HCl(aq) in aqueous fluids over a wide range of temperatures and pressures. At temperatures < 300 °C and P ≥ Psat, the HCl(aq) ion pair was unstable and dissociated quickly into H+ and Cl−. At temperatures ≥ 300 °C, the simulations reveal the strong pressure dependence of the association/dissociation of chloride. Upon an increase in pressure from Psat to 5 kbar, the predominant form of chloride changed
Conclusion and geological implications
The MD simulations in this study provide insights into the molecular structure and dissociation mechanism of the HCl(aq) molecule and the corresponding products. Despite the uncertainties and limitations of the MD simulation approach, our results provide a reasonable prediction of the dissociation constants for HCl(aq) compared with the existing experimental data. Furthermore, the MD simulations provide the first thermodynamic data for HCl speciation at high pressures relevant to conditions
Acknowledgements
Research funding was provided by the CSIRO OCE fellowship to Y.M., and by the Australian Research Council (ARC) to W.L. (FT130100510), J.B. (DP130100471) and J.G. (DP160100677). D.S. was supported by UK NERC grant NE/P002196/1. The MD calculations in this work were supported by resources provided by the Pawsey Supercomputing Centre with funding from the Australian Government and the Government of Western Australia and the high-performance computers in CSIRO. We are also grateful for the helpful
References (85)
- et al.
A review of the coordination chemistry of hydrothermal systems, or do coordination changes make ore deposits?
Chem. Geol.
(2016) - et al.
H3O+/Cl− ion pairing in hydrothermal solutions by simulation and electrical conductance. A review
J. Mol. Liquids
(2003) - et al.
The dissociation constant and henry's law constant of HCl in aqueous solution
Atmos. Environ.
(1986) - et al.
Molecular dynamics simulation study of ionic hydration and ion association in dilute and 1 molal aqueous sodium chloride solutions from ambient to supercritical conditions
Geochim. Cosmochim. Acta
(1998) - et al.
The relationship between subduction zone redox budget and arc magma fertility
Earth Planet. Sci. Lett.
(2011) - et al.
In situ Raman study and thermodynamic model of aqueous carbonate speciation in equilibrium with aragonite under subduction zone conditions
Geochim. Cosmochim. Acta
(2014) - et al.
The enthalpy of dilution of HCl(aq) to 648 K and 40 MPa thermodynamic properties
J. Chem. Thermodyn.
(1987) - et al.
Speciation of gold in hydrosulphide-rich ore-forming fluids: Insights from first-principles molecular dynamics simulations
Geochim. Cosmochim. Acta
(2011) - et al.
Silver speciation in chloride-containing hydrothermal solutions from first principles molecular dynamics simulations
Chem. Geol.
(2012) - et al.
Structures and acidity constants of arsenite and thioarsenite species in hydrothermal solutions
Chem. Geol.
(2015)
The chemistry of subduction-zone fluids
Earth Planet. Sci. Lett.
The dissociation constant and Henry's law constant of HCl in aqueous solution
Atmos. Environ. (1967)
Palladium complexation in chloride- and bisulfide-rich fluids: Insights from ab initio molecular dynamics simulations and X-ray absorption spectroscopy
Geochim. Cosmochim. Acta
Speciation and thermodynamic properties of zinc in sulfur-rich hydrothermal fluids: Insights from ab initio molecular dynamics simulations and X-ray absorption spectroscopy
Geochim. Cosmochim. Acta
Metal complexation and ion hydration in low density hydrothermal fluids: Ab initio molecular dynamics simulation of Cu(I) and Au(I) in chloride solutions (25–1000 °C, 1–5000 bar)
Geochim. Cosmochim. Acta
Ab initio molecular dynamics simulation and free energy exploration of copper(I) complexation by chloride and bisulfide in hydrothermal fluids
Geochim. Cosmochim. Acta
Complexation of gold in S3−-rich hydrothermal fluids: Evidence from ab-initio molecular dynamics simulations
Chem. Geol.
Zinc complexation in chloride-rich hydrothermal fluids (25–600 °C): A thermodynamic model derived from ab initio molecular dynamics
Geochim. Cosmochim. Acta
Calculation of the standard partial molal thermodynamic properties and dissociation constants of aqueous HCl0 and HBr 0 at temperatures to 1000°C and pressures to 5 kbar
Geochim. Cosmochim. Acta
Silver in geological fluids from in situ X-ray absorption spectroscopy and first- principles molecular dynamics
Geochim. Cosmochim. Acta
The ion-pair constant and other thermodynamic properties of HCl up to 350°C
Geochim. Cosmochim. Acta
Volumes and heat capacities of aqueous solutions of hydrochloric acid at temperatures from 298.15 K to 623 K and pressures to 28 MPa
J. Chem. Thermodyn.
Complexation of Cu+ in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics
Geochim. Cosmochim. Acta
Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Correlation algorithms for ionic species and equation of state predictions to 5 kb and 1000°C
Geochim. Cosmochim. Acta
Inorganic species in geologic fluids: correlations among standard molal thermodynamic properties of aqueous ions and hydroxide complexes
Geochimica et Cosmochimica Acta
A suite of programs, OptimA, OptimB, OptimC, and OptimS compatible with the Unitherm database, for deriving the thermodynamic properties of aqueous species from solubility, potentiometry and spectroscopy measurements
Appl. Geochem.
Computation of the pK of liquid water using coordination constraints
Chem. Phys.
Thermodynamic assessment of hydrothermal alkali feldspar-mica-aluminosilicate equilibria
Geochim. Cosmochim. Acta
Prediction of the thermodynamic properties of aqueous metal complexes to 1000 °C and 5 kb
Geochim. Cosmochim. Acta
Water in the deep Earth: The dielectric constant and the solubilities of quartz and corundum to 60 kb and 1200 °C
Geochim. Cosmochim. Acta
Experimental study of dissociation of HCl from 350 to 500°C and from 500 to 2500 bars: Thermodynamic properties of HCl°(aq)
Geochim. Cosmochim. Acta
The dissociation mechanism of H2O in water studied by first-principles molecular dynamics
Chem. Phys. Lett.
Prediction of the PVT properties of water over wide range of temperatures and pressures from molecular dynamics simulation
Phys. Earth Planet Inter.
Roles of water for chemical reactions in high-temperature water
Chem. Rev.
Computer simulation of liquids
Oxford University Press
Structure, dynamics, and spectral diffusion of water from first-principles molecular dynamics
J. Phys. Chem. C
A systematic study of chloride ion solvation in water using van der Waals inclusive hybrid density functional theory
Mol. Phys.
Density-functional exchange-energy approximation with correct asymptotic behavior
Phys. Rev. A, General Phys.
Zur quantentheorie der molekeln
Ann. Phys.
Binding modes of oxalate in UO2 (oxalate) in aqueous solution studied with first-principles molecular dynamics simulations. Implications for the chelate effect
Dalton Trans.
Density functional theory study of uranium(VI) aquo chloro complexes in aqueous solution
J. Phys. Chem. A
Unified approach for molecular dynamics and density-functional theory
Phys. Rev. Lett.
Cited by (30)
Molecular dynamics simulations of LiCl ion pairs in high temperature aqueous solutions by deep learning potential
2022, Journal of Molecular LiquidsCitation Excerpt :Molecular dynamics (MD, both classical MD and ab initio MD), an efficient method to study the association behavior of metal complexes at the micro-level, provides a detailed molecular picture of the ion pairs [4–7]. Ion association constants could be successfully calculated from the MD potential of mean force (PMF) [8–12]. The accuracy of classical molecular dynamics (CMD) simulations, is force fields sensitive, which counts on the precision of the applied empirical potential [13].
An experimental and thermodynamic study of sphalerite solubility in chloride-bearing fluids at 300–450 °C, 500 bar: implications for zinc transport in seafloor hydrothermal systems
2022, Geochimica et Cosmochimica ActaCitation Excerpt :However, the effect of chlorinity is more significant at higher fluid pH (Fig. 7). We suggest that this is caused by the fact that at elevated temperature, the effects of in-situ pH on Zn solubility decreases as a result of strong complexation of HCl(aq) (Tagirov et al., 1997; Ho et al., 2001; Mei et al., 2018), such that the formation of high-order Zn chloride complexes becomes relatively more important in controlling Zn solubility. The experimental data from Cygan et al. (1994) indicate similar trends (Fig. 7; Table S6), in that fluids with lower quench pH contain higher Zn concentrations, and the variation of salinity plays a negligible role on Zn solubility at low pH (Fig. 7a,b).