While the exact concentrations of these species remain undetermined from spectroscopic data, it's evident that the first equilibrium is shifted significantly to the right, with cryoscopic data suggesting a dissociation of the complex around 15%.
The unique interplay of molecular interactions, particularly involving HF₂⁻ and H₃O⁺ ions, enhances proton transport in HF solutions. The observed increase in conductivity suggests a Grotthuss-like hopping mechanism facilitated by hydrogen-bonded HF₂⁻ networks, consistent with prior spectroscopic studies.5 While direct experimental confirmation of proton hopping in HF solutions remains limited, previous spectroscopic studies have identified structural features consistent with this transport mechanism. This behavior, starkly contrasting with other salt/solvent electrolyte solutions, emphasizes the profound impact of these molecular dynamics on the macroscopic property of conductivity. Interestingly, the absence of a conductivity maximum, as typically expected in Abu-Lebdeh’s model, suggests that the solution's structure remains unaltered and is not heterogeneous which contradicts the presence of the eutectic in the phase diagram. The answer might lie in the formation of unique microdomains with enhanced stability that is likely due to hydrofluoric acid's small size, which seamlessly integrates into the hydrogen-bonded chains of water, forming (HF)x(H2O)y clusters akin and intertwined to the (H2O)x clusters. These extended structures or clusters have been theoretically predicted where the larger clusters with equimolar amounts of HF and H2O are the most stable due the very strong H2O --- HF interactions.27 Such a structural formation facilitates a Grotthus-type hopping mechanism that prevails across varying concentrations, further elucidating hydrofluoric acid's distinctive conductivity attributes in contrast to conventional solvent/electrolyte solutions as depicted in the phase diagram of Fig. 2b.
Modeling the behavior of HF in water at high concentrations is challenging due to several limiting factors. First, the ionic strength of concentrated HF solutions lies outside the range where simple electrolyte solution models (e.g., the Debye–Hückel theory) are accurate. Early researchers addressed this issue by introducing a concentration-dependent dissociation degree or “stoichiometric” equilibrium constant for HF to fit experimental data.28 This approach acknowledged that the apparent extent of HF dissociation decreases as total HF concentration increases. Modern thermodynamic models (such as those based on Pitzer ion-interaction equations or specific ion-pairing schemes) improve upon simple theories and have been applied to HF–H2O systems.29 However, these models still require reliable experimental data for calibration, and beyond a certain concentration (above roughly 6–10 M HF) such data become scarce or uncertain. Another limitation is that conventional speciation models may not include all relevant fluoride species present at high acid concentrations. In reality, at high [HF] the formation of HF2− (and potentially higher-order fluoride complexes) significantly affects the solution chemistry.30 If a model neglects such associations and assumes that HF dissociates only into H+ and F−, it will underestimate the true acidity and mis predict the species distribution. In practice, treating HF with a fixed dissociation constant or as fully dissociated leads to errors when predicting thermodynamic properties of concentrated solutions.31 This is because effective equilibrium constants themselves change with concentration due to medium effects and ion pairing that are not captured by simplistic models.
Modeling Ionic Transport: Speciation, Mobility, and κ(T,x) Analysis
Ionic conductivity is a product of the number of charge carriers, by their charge, and charge mobility as represented in Eq.
1
A
where
ni is number of free charged species,
qi denotes charge number and 'µ
i' is mobility of free charged species. To perform the theoretical calculation of ionic conductivity based on Eq.
1, one must first understand speciation of ions present in the aqueous solutions of both acids, and obtain ionic mobility and activity data for each. We purposefully ignored the contribution of OH
−, and H
+ from water auto dissociation, as it is few orders of magnitude lower than the contributions of the ions dissociated from most acids, even weak ones.
ni which represents the effective concentration of ions in molarity,
i.e. activity,
ai. For acetic acid, data for activity is only available up to 1 M from Cohn et al.
32 While that for hydrofluoric acid was obtained to higher concentrations from the National Bureau of standards.
28 On the other hand, 25°C values of mobility at infinite dilution, µ
io, were sourced from the literature to be 349.8 Scm
2 for H
+, 40.9 Scm
2 for CH
3COO
−, 55.4 Scm
2 for F
− ions,
33 and 84.6 Scm
2 for HF
2−.
30 pH values versus concentration for acetic acid and hydrofluoric acid are tabulated in table S.1 and S.2, respectively. To the best of our knowledge, mobility values at elevated temperatures and concentrations are scarce in literature for most ions. We have used these extracted species activity and mobility data to theoretically calculate conductivity and the results of this exercise are presented graphically in Fig S.1. For acetic acid, one can see clearly the great discrepancy between the theoretically calculated line versus the experimental reported conductivity curves. This is mainly due to underestimating the reduced ionic mobility as concentration increases. Several studies have focused solely on the this effect and introduced empirical equations to correct for mobilities at higher concentration (< 0.1 M) (Eq. 2).
34,35
Where d is a parameter with a value > 0.5, z is the charge number and I is ionic strength. But at higher concentration it can be directly related to x (empirically to simplify the equation due to lack of literature values above 0.1 M).
36
Ion mobility refers to the average speed of ions per unit electric field strength, determined by the interplay between the external electric field and resistance to ion movement. When considering the ion and its hydration layer as a single entity, the resistance encountered during migration includes ion–ion, ion–solvent, and solvent–solvent interactions. Ion–ion forces are long-range electrostatic interactions, while ion–solvent and solvent–solvent forces are short-range interactions. At low concentrations, long-range interactions dominate, allowing short-range interactions to be neglected. However, as concentration of electrolyte solution increases, the molecular spacing decreases, making short-range interactions significant and causing a rapid increase in movement resistance. Consequently, ion mobility typically decreases with higher electrolyte solution concentrations.36 This helps explain the presence of a maximum in conductivity isotherms of most aqueous electrolyte solutions. For hydrofluoric acid, one expects that this is not the case. In fact, we found evidence that the mobility of the HF2− triple ion actually increases with concentration!30
Our research group has introduced a semi-empirical equation (Eq. 4) that combines equations
1 and 3 and introduces temperature effect (that causes an increase in the number of free charges and their mobility) hence takes into account mobility and activity changes with temperature and concentration.
5
where A, B, C, and D are parameters and the pre-exponential term corrects activity while the exponential term is related to mobility. Fitting both acids’ κ vs x and T data lead to perfect surface fits with both fits’ coefficient of determination (R2) being 100% as shown in Fig. 3 below.
Figure 4(a) elucidates more clearly in 2D the variation of ionic conductivity (k) with molar fraction (x) at different temperatures (25°C, 45°C, and 65°C) for acetic acid solutions. In the extremely dilute region, specifically below x = 0.03, a noteworthy observation is the minimal deviation in conductivity across the three temperatures, with the curves nearly overlapping. This subtlety in temperature influence can be attributed to the pronounced effects of low molar volume and viscosity in water-rich mixtures, optimizing ion mobility and conductivity. However, as the solution approaches xmax, the disparity in conductivity between the different temperatures becomes significantly pronounced. Interestingly, the increment in conductivity observed when transitioning from 25°C to 45°C is nearly double that of the shift from 45°C to 65°C. Beyond xmax, while the temperature effect persists, it is less pronounced, with the conductivity values converging at higher concentrations. Figure 4(b) presents the ionic conductivity (k) of hydrofluoric acid (HF) solutions across varying molar fractions (x) at temperatures of 0°C, 18°C, 37.8°C, and 93.3°C. Remarkably, the isotherms depicted are linear, with a near overlap in the extremely dilute region, similar to the behavior observed in acetic acid solutions (Fig. 4(a)). However, a distinct characteristic of hydrofluoric acid solutions is the continuous and increasing influence of temperature on conductivity, with no discernible xmax. The temperature effect strengthens as x increases, and the increment in conductivity is greater between 0°C and 18°C than between subsequent temperature intervals. This pattern mirrors the non-linear temperature effect observed in acetic acid solutions, albeit in the absence of a conductivity peak. The pronounced temperature sensitivity in hydrofluoric acid solutions especially in concentrated regions, was observed previously for other solutions such as nitrates and acids, and can be attributed to the presence of a higher number of ion pairs and clusters (IPs and ICs) and the consequent increase in mobile "free" ions as temperature rises, a phenomenon previously discussed.4 The expansion of the liquid and the increase in free volume at higher temperatures further enhance ion mobility, contributing to the substantial increase in conductivity.
Figure 5 outlines the variation of the natural logarithm of ionic conductivity (ln(k)) against the reciprocal of temperature (1000/T) for both hydrofluoric acid and acetic acid solutions. The linear trend exhibited by both acids in this plot signifies their adherence to the classical Arrhenius behavior, described by Eq. (1)
(5) where E
a is the activation energy, R is the gas constant, T is the temperature in Kelvin, and “A” is the pre-exponential factor. While the original conductivity data for hydrofluoric acid displayed linearity with respect to x, this graph is valuable for extracting the variation of activation energy and the pre-exponential factor across different molar fractions.
Figure 6(a) and (b) presents the variation of activation energy (Ea) with molar fraction (x). In the case of acetic acid, a notable increase in Ea is observed in the dilute region, reaching a plateau around xmax. Beyond this point, Ea exhibits a further increase until x = 0.15, followed by a sudden drop to values close to the value near xmax. Subsequently, Ea experiences another increase, stabilizing between x = 0.27 and 0.33, before rising once more at higher molar fractions. A similar trend has been observed in another weak acid (H₃PO₄),5 though further research is needed to determine whether this behavior is common among weak acids. This trend in the Ea vs x graph is very similar to that of aqueous phosphoric acid solutions. For hydrofluoric acid, while data in the extremely dilute region is unavailable, we hypothesize that Ea initiates at a value comparable to most acids (approximately 8.5 kJ mol− 1) and undergoes a decrease, aligning with the first available data point at x = 0.08 where Ea is around 3.3 kJ mol − 1. Following this, a continuous decrease in Ea is observed until a minimum is reached at x = 0.14, after which a slight increment and subsequent decrement are noted at higher molar fractions. This Ea vs x trend is very similar to that observed in aqueous nitric acid solutions except that nitric acid’s Ea does not drop as much as that of hydrofluoric acid. The Ea drops at the minimum point and also it continues to increase beyond the minimum unlike that of hydrofluoric acid. This may be related to the unique structural features of hydrofluoric acid (e.g., HF₂⁻ formation and hydrogen bonding), though we emphasize that this remains a hypothesis due to the lack of direct structural validation. However, it is worth noting that the overall decrease in Ea is less than the thermal activation energy of 25°C at ~ 2 kJ mol − 1 indicating a facile transport of ions with little activation necessary.
Ionic conductivity in liquid electrolyte solutions shares some conceptual similarities with amorphous materials, such as glass-forming liquids or solid ionic glasses above the glass transition temperature (Tg). Regardless of the conduction mechanism, ionic conductivity can be expressed as the product of the concentration of charge carriers, their mobility, and their charge, as described in Eq. 1. While we reference the empirical Vogel-Tammann-Fulcher (VTF) equation as a useful tool to describe conductivity-temperature relationships in some materials, we clarify that we are not asserting that our data strictly follows the VTF equation. Instead, we included the VTF discussion as a conceptual framework to highlight the transition from Arrhenius to non-Arrhenius behavior in many ionic systems, especially in glassy or polymeric electrolyte. It is also important to note that direct application of the VTF equation to liquid aqueous electrolyte solutions is non-trivial because defining key parameters such as the ideal vitreous temperature (T₀) and free volume (Vf) in these systems remains challenging. While VTF theory has successfully described transport in polymeric or glassy systems, its role in liquid electrolyte solutions remains largely qualitative. Previous studies have drawn analogies between VTF-type behavior and ionic conduction in highly concentrated or supercooled liquids, but these analogies are limited by the lack of a well-defined T₀. In light of this, we treated VTF as a conceptual tool rather than an explicit fitting model, which is why we did not perform a VTF analysis on our data. To analyze these behaviours, we primarily apply the Arrhenius model, while referencing the VTF framework conceptually to discuss its potential relevance to ionic transport in certain cases focusing on the role of activation energy, ion mobility, and the potential for free volume influence in determining conductivity. In glass electrolytes, two different mechanisms govern temperature dependence: below Tg, conductivity follows the Arrhenius behavior, while above Tg, chain movement and free volume rearrangement become significant, resulting in Vogel-Tammann-Fulcher (VTF) dependence.
Nascimento et al.
37 summarized this behavior and provided an expression for the concentration of charge carriers (n⁺):
where
is the Gibbs energy of cationic vacancy and interstitial pair formation, calculated from the corresponding formation enthalpy (
and entropy
) Nascimento also described the equation to calculate the mobility (µ+) of charge carriers in the presence of an electric field at low temperatures (below T
g):
where λ is the mean distance between cationic sites,
is the jump attempt frequency and
is the migration Gibbs free energy calculated the migration enthalpy (
and entropy
).
Using Equations
1,
6, and
7, Nascimento derived the equation to describe the cationic conductivity (σT) in the low-temperature range for glass electrolytes below T
g:
This is essentially the Arrhenius equation, where the experimental value of activation energy (
) and the pre-exponential factor (
) are identified as:
In hydrofluoric acid, the Arrhenius-like behavior observed at low to moderate concentrations aligns with Eq. 8, suggesting that ion mobility remains relatively unaffected by free volume changes. Instead, conductivity relies on stable hydrogen-bonded clusters and proton transfer, allowing for a unique linear increase without transitioning to VTF behavior.
Above T
g, when free volume is activated, another conduction mechanism is observed which involve the transfer of defects to neighboring positions because of the local deformations of the macromolecular chains due to fluctuations in the free volume of the chain segments.
38 However, in aqueous solutions, the concept of free volume is not as well-defined as in polymers or glasses, making the direct application of VTF theory challenging. We have therefore avoided a formal VTF analysis, as it would require detailed viscosity and relaxation data not presently available. Moreover, the conceptual application of VTF theory to aqueous systems remains nontrivial due to the limited understanding of free volume behaviour in liquids. Instead, we focus on the conceptual aspects of ion transport and their implications for observed conductivity trends. To arrive to an equation equivalent to Eq.
8 but now accounting for this new conduction mechanism, one must first define a couple of parameters. First, this second mechanism only appears above the ideal vitreous temperature (T
0), where at T > T
0, the cationic conduction can occur by either a successful activated jump or by the entropic free volume mechanism. The probability of an activated jump is given by P
1 which is related the enthalpy of migration by the following equation:
38
and P
2 is the probability that minimum free volume is available for the second mechanism:
38
where Vf
* is the minimum required free volume for conduction and V
f is the chain segment mean free volume. V
f temperature dependence is given by the following equation:
38
where
is the free volume thermal expansion coefficient, P
2 can now be re-arranged to be:
38
This equation confirms that conduction due to defect transfer can happen only at T > T
0. At this high T range, cationic displacement can occur by either mechanisms and the total probability of a successful would be statistically
38
and we can now arrive to a new equation for mobility:
meaning that mobility at the high temperatures depends on the probabilities P
1 and P
2 and a free volume migration is only possible above T
0. However, it is only noticeable for T > T
g. Souquet et al.
38 suggest that at the following temperature range T
0 < T < T
g, P
2 remains very small compared to P
1 and the mobility equation remains the same as Eq.
7. At temperatures above T
g, P
2 is not negligible and the conductivity equation now becomes:
38
This equation is the product of two exponentials where the first represents the activated mechanism and the second is associated with the VTF behavior. Notably, the pre-exponential factor is not changed in this equation which means that the limit of the conductivity-temperature product at infinite T does not depend on any conduction mechanism. VTF behavior typically requires sufficient free volume, enabling large-scale ion migration facilitated by chain motion or defect movement. In acetic acid, we observe features that may suggest a transition beyond xmax where ion mobility is restricted—possibly due to viscosity and ion pairing—which is conceptually consistent with VTF-type behavior described in Eq. 17. However, hydrofluoric acid does not exhibit VTF-like behavior due to its unique proton hopping mechanism along hydrogen-bonded networks, a Grotthuss-like process that sustains high mobility without the need for free volume reorganization. Thus, while acetic acid’s conductivity behavior transitions as expected with increasing free volume, hydrofluoric acid’s proton-hopping mechanism maintains a stable structure that bypasses typical VTF dependencies, even at higher concentrations
Figure 6(c) and (d) shows the variation of the pre-exponential factor (A) with molar fraction (x). A can also be described as the conductivity when charge carriers encounter no energy barriers, meaning the activation energy (Ea) is non-existent.39,40 For acetic acid, (A) exhibits an initial increase with increasing x, reflecting the enhancement in ion jump frequency and available conductive sites due to higher ion concentrations. This trend continues until reaching a plateau around xmax, where the solution's structural constraints begin to counterbalance the benefits of increased ion concentration. Beyond x = 0.15, a notable sudden drop followed by a consistent decline in (A) is observed, that could be attributed to the increased formation of non-conductive ion pairs or clusters and a significant rise in viscosity, which collectively impede ion mobility. For hydrofluoric acid, (A) demonstrates a slow and gradual continuous increase, suggesting that the solution's structural changes facilitate better total ion mobility across all concentrations. We hypothesize that this behaviour may arise from the formation of complex ions that facilitate charge transport more efficiently, without strong hindrance from ion pairing or aggregation. However, this remains a proposed explanation based on indirect inference.
Lastly, it is important to note that we have based our calculations of Ea and A on the simple Arrhenius equation, assuming its validity in the κ vs T data range despite dealing with aqueous solutions that actually behave like glasses or polymers that are well above their Tg. However, one can say that the isotherms follow pseudo-Arrhenius behavior due to the narrowness of the temperature range we are discussing especially for acetic acid, where the difference between the highest temperature and the lowest is just 40°C. if the temperature range was to be larger, it would be wiser to have considered the VTF behavior when linearizing the data.
We emphasize that many of the mechanistic suggestions made in this section—particularly those involving ion pairing, complex formation, or structural transport modes—are working hypotheses motivated by literature and consistent with our data trends, but they are not experimentally verified. This reflects a broader challenge in the field, as probing the microstructure of highly concentrated or structurally unconventional liquid systems remains a major experimental and theoretical difficulty.