Determination of the Ionic Association Constants of Na+ with CO32− and HCO3− Ions, in NaCl-NaHCO3-H2O Ternary Systems, at 25 °C

The ionic association constants of sodium with carbonate ion (K1C′) and acidic carbonate ions (K2C′) were measured in NaCl-NaHCO3-H2O ternary systems to determine the distribution of sodium among the chemical species present in the growth medium of Chlorella homosphaera 424 algae. The mean activity coefficients of sodium chloride (in pure sodium chloride and in a mixture of electrolytes) were determined experimentally using two electrochemical cells, namely Ag, AgCl| KCl (3 M)|| NH4NO3 (1 M)| NaCl (mNaCl)| Na+-ISE and Ag, AgCl|KCl (3 M)|| NH4NO3 (1 M)| NaCl (mNaCl)| Cl−-ISE. The studies carried out show that the values of the association constants of K1C′ and K2C′ do not depend on the composition of the medium, but only on the effective ionic strength. The experimentally obtained γNaCl0 values in the binary system are comparable to the mean activity coefficients values for NaCl, calculated using data from the literature, with −0.9 to 0.1% relative standard deviation. The obtained results show that the experimentally determined mean activity coefficient in the ternary system, γNaCl, is smaller than γNaCl0 in the binary system over the entire field of ionic strengths studied. The ternary system NaCl-NaHCO3-H2O obeys Harned’s rule.


Introduction
The studies carried out and presented in this paper are part of the general topic of interest of reducing the greenhouse effect that is generated globally by emission of gases such as carbon dioxide, methane, nitrous oxide, chlorofluorocarbons (CFCs), and sulfur dioxide.Currently, the most important greenhouse gas is carbon dioxide [1,2].
The accumulation of carbon dioxide in the atmosphere, with all its negative consequences for human society, is both a reality and a challenge for science and technology today.The warming of the atmosphere is due to carbon dioxide emissions in a proportion of 52%.The main cause for the increase in the level of carbon dioxide in the atmosphere is the use of fossil fuels for energy generation, transport, industry, and domestic utilities.
To reduce carbon dioxide emissions to the atmosphere, many technological solutions have been advanced, including geological sequestration, oceanic sequestration, biological sequestration, and mineral sequestration [1,[3][4][5][6][7].A relatively recent biological solution for reducing the harmful effects of carbon dioxide emissions in the atmosphere is its sequestration in the form of biomass, using microalgae.Microalgae are present in the entire ecosystem, not only in the aquatic environment, but also on land, presenting great variety.It is estimated that there are more than 50,000 species, but only around 30,000 have been studied and analyzed [8,9].
Microalgae have an interesting and completely unexploited potential in biotechnology.They synthesize products such as natural dyes, polyunsaturated fatty acids, polysaccharides, and vitamins.As one of the largest groups of oxygen producers worldwide, microalgae play an important role in ecosystem balance [3].After extraction, the spent material can be used as fertilizer, fish food, etc.
Biologically, the selection of microalgae species that can be used for the efficient sequestration of carbon dioxide from industrial emissions is not at all simple, because they are different both in morphology, and during the growth period [10].
The cultivation of microalgae requires real-time analytical control of the chemical species present in the specific culture media, which are complex synthetic solutions of electrolytes with a composition that is significantly different to that of natural waters/media, especially due to the high concentration of acidic carbonate ions; these can vary in the 0.10-0.50mol/L range, in the case of a 0.2-0.55mol/L total ionic strength [10], when the ionic species are closely associated.
The nutrient growth medium for Chlorella homosphaera 424 is medium Z (Zarrouk) [10], consisting of high levels of NaHCO 3 , accompanied by K 2 HPO 4 , NaNO 3 , K 2 SO 4 , NaCl, MgSO 4 ×7H 2 O, and CaCl 2 ×2H 2 O.In situ monitoring of the species distribution is a critical step for avoiding inactivation of microorganisms and offering optimal development.Potentiometric measurements might deliver real-time information that is useful for building a species distribution diagram, provided that all possible interactions and associations are taken into account.Thermodynamic constants describing such processes have been reported in the literature, generally at ionic strengths higher than 0.5 mol/L.Developing a procedure for anionic and cationic species determination, by direct or indirect potentiometric measurements, in the culture media of Chlorella homosphaera 424 algae could start with studying the media influence on the values of the association constants for ternary mixtures of the compounds present (NaCl-NaHCO 3 -H 2 O NaCl-Na 2 SO 4 -H 2 O, and NaCl-Na 2 HPO 4 -H 2 O), before approaching the multicomponent system (NaCl-NaHCO 3 -Na 2 HPO 4 -Na 2 SO 4 -H 2 O).
Software packages have often been used for chemical speciation calculation and graphical representation in the literature (HySS 2009 [11,12], Aq_solution 3.2 [13], KEV 3.0 [14]).Generally, input data for the specialized applications are the reactants' total concentrations, the association constants, and the solubility products.
There are no values reported in the literature for the association constants in the NaCl-NaHCO 3 -H 2 O ternary system at ionic strengths lower than 0.5 mol/L.Kester and Pytkowicz [15,16] used a method based on electromotive voltage measurement to determine the association constants in sodium chloride-and sodium sulfate-containing solutions at hydrogen ion concentrations where the HSO 4 − type ionic species are negligible, using a 0.638-0.690mol/L effective ionic strength.The activity of the free cation (Na + ) was measured in standard sodium chloride solutions and solutions containing NaCl-Na 2 SO 4 .The authors demonstrated that the association constants do not depend on the composition of the medium in which they were determined, but on the effective ionic strength of the medium.The reproducibility of their results decreases with increasing sodium chloride concentration.
Butler and Huston [17][18][19] determined the ionic association constants of NaHCO 3 and NaCO 3 − in the presence of sodium chloride in cells without liquid junction, in the 0.5-3 mol/L ionic strength range.
Pytkowicz and Hawley [20] measured the apparent association constants for the formation of acidic carbonate and carbonate ion pairs with Na + , Ca 2+ , and Mg 2+ at 25 • C and 0.72 mol/L ionic strength.The results were used, together with the sulfate ion association constants obtained by Kester and Pytkowicz [16], to determine the major ionic species in seawater at 25 • C and 34.8% salinity.Capewell, Hefter, and May [21] determined the association constant of NaCO 3 − in CsCl medium in the 0.5-7.0mol/L ionic strength range and at 25 • C, using sodium ion-selective electrodes.Blum [22] proposed a new method for determining ion association constants in KCl-Me 4 NCl and Me 4 NCl media, using ion-selective electrodes at various ionic strengths (Table 1).Chung Yan-Chun [23]  A synthesis of the experimental values of the association constants of sodium with carbonate and carbonate acid ions available in the literature is presented in Table 1.The purpose of this work is to evaluate the ionic association constants of sodium ion with carbonate and acidic carbonate ions in ternary systems of the NaCl-NaHCO 3 -H 2 O type at 0.1-0.5 mol/kg total ionic strength, to determine the distribution of sodium among the chemical species present in the growth medium of Chlorella homosphaera 424 algae, also known as "sodium speciation".This paper is part of a larger study which has as main objective to develop a procedure for the determination of species (anionic and cationic) by direct or indirect potentiometric measurements, in situ or on samples taken from the culture media of Chlorella homosphaera 424 algae.
Sodium speciation in the NaCl-NaHCO 3 -H 2 O ternary system has been approached through the ionic association model as modified by Pytkowicz and Kester [16].It is based on the hypothesis that the change in the mean activity coefficients of the salts present occurs only as a result of the formation of ion pairs.The verification of this hypothesis has been carried out by preliminarily studying the behavior of the NaCl-NaHCO 3 -H 2 O ternary systems in relation to Harned's rule.

Results and Discussion
The analytical control of a biological carbon dioxide sequestration process involves two stages.The first stage of the carbon dioxide sequestration process is its capture by absorption in an aqueous solution of sodium carbonate.The chemical species present in the absorption medium are CO 3 2− , HCO 3 − NaCO 3 − , NaHCO 3 , and H 2 CO 3 , and physically dissolved CO 2 , Na + , H 3 O + , and OH − .The control of the absorption process involves the monitoring of the chemical species in the solution, and sometimes in the gas phase.The experimental data obtained in the first stage were the subject of a previous study [24].
The second stage of the carbon dioxide sequestration process is storage of CO 2 in the form of microalgal biomass.In this step, the acidic carbonate solution obtained in the first step is diluted and fertilized to create a suitable culture medium for the growth of Chlorella homosphaera 424 algae.A typical Chorellahomosphaera 424 growth media composition [10] consisted of 16.80 g/L NaHCO 3 , 0.5 g/L K 2 HPO 4 , 2.5 g/L NaNO 3 , 1 g/L K 2 SO  The association constants of sodium ions with the anionic species present in the algal media growth were determined simultaneously with the concentrations of all species, by solving the mathematical model of the equilibrium state for the electrolytic system.The equilibrium state of the electrolyte solution is described by a system of equations including mass balance equations, neutrality conditions (electric charge balance), and chemical equilibrium relations for the present species.
The model is based on three simplifying hypotheses: (i) the anionic species with which sodium ions can associate in the solution are CO and the influence of the anions present in micro-concentrations is negligible; (ii) the salts of strong acids such as sodium chloride, sodium nitrate, and sodium sulfate are considered completely dissociated; and (iii) mean activity coefficients of the salts change only as a result of the formation of ion pairs, provided that the total ionic strength is kept constant.Standard sodium solutions, also used by J. N. Butler [18] and R. M. Pytkowicz [16], were considered as reference systems, because they present a number of advantages.The associations between sodium and chloride ions are negligible and there are commercial and reliable ion-selective electrodes, both for sodium and chloride.In addition, sodium chloride is an intensively studied electrolyte, so there is information in the literature on its mean activity coefficients, obtained by various procedures, with appropriate uncertainties.
In order to evaluate the magnitude of the associations between sodium ion and the CO 3 2− and HCO 3 − species, the influence of the composition on the activity coefficients of sodium chloride was studied at constant ionic strengths, and Harned's coefficients (α 12 ) [25] were determined in NaCl-NaHCO 3 -H 2 O ternary electrolyte systems.Both the mean activity coefficients, obtained experimentally, and database [26] values were used.
The association constants were determined as functions of the effective ionic strength in NaCl-NaHCO 3 -H 2 O ternary systems.In the pH range close to neutral, in which algae growth occurs, carbonic acid is completely dissociated and presents two differently protonated forms, CO 3 2− and HCO 3 − .An iterative procedure was used to calculate the association constants.In the first iteration, it was assumed that the effective ionic strength, I E , was equal to the total ionic strength, I T , and the values for the (n − 1) association constants were calculated, where n is the number of differently protonated species of the studied acid.The procedure was repeated until two successive I E values did not differ by more than 0.001 units.The experimental values of the association constants of sodium with carbonate ion (K 1C ) and acidic carbonate ions (K 2C ) for NaCl-NaHCO 3 -H 2 O ternary systems will be used later for sodium speciation in the NaCl-NaHCO 3 -Na 2 HPO 4 -Na 2 SO 4 -H 2 O multicomponent electrolyte system.Speciation software packages available in the specialized literature were used for the calculations and graphic representations of the compositions of multicomponent mixtures [11][12][13][14]26].

Determination of Mean Activity Coefficients
Table 2 presents the experimental values of the mean activity coefficients of sodium chloride, γ NaCl , in the NaCl-NaHCO 3 -H 2 O ternary system for different compositions and constant total ionic strengths, and the experimental values of the mean activity coefficients of sodium chloride, γ 0 NaCl , in the NaCl-H 2 O binary system.To facilitate the comparison of the obtained results with other data from the literature, concentrations and ionic strengths are expressed in mol/kg.The γ 0 NaCl values obtained experimentally are comparable to the mean activity coefficients of sodium chloride values calculated with the Pitzer-Millero equation and using the SIT model [27], with a −0.9 to 0.1% relative standard deviation.The results show that the activity coefficients of sodium chloride in the ternary system are smaller than the mean experimental activity coefficient in pure sodium chloride, γ NaCl < γ 0 NaCl , over the entire field of ionic states studied.This is due to the presence of sodium acidic carbonate in the system.

Determination of Harned's Coefficient
The association model is applicable with good results when the ionic associations are significant.An assessment of the magnitude of the ionic associations is made based on the criterion of compliance with Harned's rule.When this rule is respected, the extent of ionic associations is small [16].
Harned's rule is given in Equation ( 1): where: I T -total ionic strength; γ NaCl − -activity coefficient of sodium chloride in a mixture of electrolytes at I T ionic strength; γ 0 NaCl − -mean activity coefficient of sodium chloride in the NaCl-H 2 O binary system, at I T ionic strength; α 12 -Harned's coefficient; X NaCl -ionic strength fraction of sodium chloride (X NaCl = m NaCl /I T ).
Figure 1 shows that the logγ NaCl − logγ 0 NaCl term varies linearly with I T *(1 − X NaCl ) at a 0.1-0.5 mol/kg constant total ionic strength.The NaCl-NaHCO 3 -H 2 O ternary system obeys Harned's rule.The compliance of Harned's rule in the NaCl-NaHCO 3 -H 2 O ternary system (Table 3) indicates weak associations of the sodium ion with the acidic carbonate ions.

Determination of Harned's Coefficient
The association model is applicable with good results when the ionic associations are significant.An assessment of the magnitude of the ionic associations is made based on the criterion of compliance with Harned's rule.When this rule is respected, the extent of ionic associations is small [16].
Harned`s rule is given in Equation ( 1): where: IT -total ionic strength;  NaCl − -activity coefficient of sodium chloride in a mixture of electrolytes at ITionic strength; NaCl − -mean activity coefficient of sodium chloride in the NaCl-H2O binary system, at IT ionic strength;  -Harned's coefficient; XNaCl -ionic strength fraction of sodium chloride ( NaCl =  NaCl  ⁄ ).
Figure 1 shows that the log NaCl -log NaCl term varies linearly with IT*(1 − XNaCl) at a 0.1-0.5 mol/kg constant total ionic strength.The NaCl-NaHCO3-H2O ternary system obeys Harned's rule.The compliance of Harned's rule in the NaCl-NaHCO3-H2O ternary system (Table 3) indicates weak associations of the sodium ion with the acidic carbonate ions.The values of Harned's coefficient (α 12 ), obtained experimentally at a 0.5 mol/kg total ionic strength, are similar to those obtained by Butler and Huston [18].There are no values of Harned's coefficient available in the literature for total ionic strengths lower than 0.5 mol/kg.Harned's coefficient depends significantly on the total ionic strength at low ionic strength values, but it becomes almost invariable at ionic strength values larger than 0.4 mol/kg (Figure 2), a behavior also reported by Sirbu et al. [28] for the NaCl-Na 2 SO 4 -H 2 O ternary system.
The values of Harned's coefficient (α12), obtained experimentally at a 0.5 mol/kg total ionic strength, are similar to those obtained by Butler and Huston [18].There are no values of Harned's coefficient available in the literature for total ionic strengths lower than 0.5 mol/kg.
Harned's coefficient depends significantly on the total ionic strength at low ionic strength values, but it becomes almost invariable at ionic strength values larger than 0.4 mol/kg (Figure 2), a behavior also reported by Sirbu et al. [28] for the NaCl-Na2SO4-H2O ternary system.The influence of the medium composition on the sodium chloride activity coefficients,  NaCl , in the NaCl-NaHCO3-H2O ternary system at constant total ionic strength and different compositions is presented in Figure 3.The influence of the medium composition on the sodium chloride activity coefficients, γ NaCl , in the NaCl-NaHCO 3 -H 2 O ternary system at constant total ionic strength and different compositions is presented in Figure 3.
The values of Harned's coefficient (α12), obtained experimentally at a 0.5 mol/kg total ionic strength, are similar to those obtained by Butler and Huston [18].There are no values of Harned's coefficient available in the literature for total ionic strengths lower than 0.5 mol/kg.
Harned's coefficient depends significantly on the total ionic strength at low ionic strength values, but it becomes almost invariable at ionic strength values larger than 0.4 mol/kg (Figure 2), a behavior also reported by Sirbu et al. [28] for the NaCl-Na2SO4-H2O ternary system.The influence of the medium composition on the sodium chloride activity coefficients,  NaCl , in the NaCl-NaHCO3-H2O ternary system at constant total ionic strength and different compositions is presented in Figure 3.The mean activity coefficient values in the NaCl-NaHCO 3 -H 2 O ternary system increase linearly with the increase in the sodium chloride concentration and the decrease in the sodium acidic carbonate level in the medium.

Determination of Association Constants
The association constant of sodium ions with carbonate ions (K 1C ) and the association constant of sodium ions with acidic carbonate ions (K 2C ) for the NaCl-NaHCO 3 -H 2 O ternary system were determined as functions of the effective ionic strength.
The dependence of the association constants upon the effective ionic strength (Figures 4 and 5) is described by Equations ( 2) and (3).
The mean activity coefficient values in the NaCl-NaHCO3-H2O ternary system increase linearly with the increase in the sodium chloride concentration and the decrease in the sodium acidic carbonate level in the medium.

Determination of Association Constants
The association constant of sodium ions with carbonate ions ( ) and the association constant of sodium ions with acidic carbonate ions ( ) for the NaCl-NaHCO3-H2O ternary system were determined as functions of the effective ionic strength.
The dependence of the association constants upon the effective ionic strength (Fig-   The mean activity coefficient values in the NaCl-NaHCO3-H2O ternary system increase linearly with the increase in the sodium chloride concentration and the decrease in the sodium acidic carbonate level in the medium.

Determination of Association Constants
The association constant of sodium ions with carbonate ions ( ) and the association constant of sodium ions with acidic carbonate ions ( ) for the NaCl-NaHCO3-H2O ternary system were determined as functions of the effective ionic strength.
The dependence of the association constants upon the effective ionic strength (Fig-    The experimental values of the association constant of sodium ions with carbonate ions (K 1C ) and the association constant of sodium ions with acidic carbonate ions (K 2C ) are presented in Table 4. Extrapolated to infinite dilution, the values of the association constants are K 1C = 14.41 and K 2C = 1.71.The possibility of using the K 1C and K 2C association constant values measured in the NaCl-NaHCO 3 -H 2 O ternary system to determine the distribution of sodium among the chemical species present in the growth medium of Chlorella homosphaera 424 algae was verified.This was carried out by demonstrating the association constants' independence of the composition of the medium, at constant total ionic strength, as demonstrated in Figures 6 and 7.
Molecules 2023, 28, x FOR PEER REVIEW 9 of 16 The experimental values of the association constant of sodium ions with carbonate ions ( ) and the association constant of sodium ions with acidic carbonate ions ( ) are presented in Table 4. Extrapolated to infinite dilution, the values of the association constants are  = 14.41 and  = 1.71.The possibility of using the  and  association constant values measured in the NaCl-NaHCO3-H2O ternary system to determine the distribution of sodium among the chemical species present in the growth medium of Chlorella homosphaera 424 algae was verified.This was carried out by demonstrating the association constants' independence of the composition of the medium, at constant total ionic strength, as demonstrated in Figures 6 and 7.The studies carried out show that the values of the association constant of sodium ions with carbonate ions ( ) and the association constant of sodium ions with acidic carbonate ions ( ) do not depend on the composition of the medium, but only on the effective ionic strength.The obtained results prove that the association constants of sodium with carbonate ion ( ) and acidic carbonate ions ( ) for NaCl-NaHCO3-H2O ternary systems can be used for sodium speciation in the NaCl-NaHCO3-Na2HPO4-Na2SO4-H2O multicomponent electrolyte system, and similarly with algal media growth.
Chemical speciation can be visualized with the aid of Hyperquad simulation and speciation (HySS2009) [11,12], a utility program for the investigation of equilibria involving soluble and partially soluble species, available in electronic format at http://www.hyperquad.co.uk/hyss.htm(accessed on 20 September 2023).HySS, a computer program compatible with the Windows operating system, allows both the simulation of titration curves and the obtaining of speciation diagrams.Input data for speciation are the reactants' total concentrations and association constants.The distribution of free ions and ion pairs present in the NaCl-NaHCO3-H2O ternary system, at different pH values, is presented in  The studies carried out show that the values of the association constant of sodium ions with carbonate ions (K 1C ) and the association constant of sodium ions with acidic carbonate ions (K 2C ) do not depend on the composition of the medium, but only on the effective ionic strength.The obtained results prove that the association constants of sodium with carbonate ion (K 1C ) and acidic carbonate ions (K 2C ) for NaCl-NaHCO 3 -H 2 O ternary systems can be used for sodium speciation in the NaCl-NaHCO 3 -Na 2 HPO 4 -Na 2 SO 4 -H 2 O multicomponent electrolyte system, and similarly with algal media growth.
Chemical speciation can be visualized with the aid of Hyperquad simulation and speciation (HySS2009) [11,12], a utility program for the investigation of equilibria involving soluble and partially soluble species, available in electronic format at http://www.hyperquad.co.uk/hyss.htm(accessed on 20 September 2023).HySS, a computer program compatible with the Windows operating system, allows both the simulation of titration curves and the obtaining of speciation diagrams.Input data for speciation are the reactants' total concentrations and association constants.The distribution of free ions and ion pairs present in the NaCl-NaHCO 3 -H 2 O ternary system, at different pH values, is presented in Figures 8 and 9.
Molecules 2023, 28, x FOR PEER REVIEW 10 of 16 Figure 7. Variation in the association constant of sodium ions with acidic carbonate ions values with sodium chloride concentration at 0.1-0.5 mol/kg total ionic strength range.
The studies carried out show that the values of the association constant of sodium ions with carbonate ions ( ) and the association constant of sodium ions with acidic carbonate ions ( ) do not depend on the composition of the medium, but only on the effective ionic strength.The obtained results prove that the association constants of sodium with carbonate ion ( ) and acidic carbonate ions ( ) for NaCl-NaHCO3-H2O ternary systems can be used for sodium speciation in the NaCl-NaHCO3-Na2HPO4-Na2SO4-H2O multicomponent electrolyte system, and similarly with algal media growth.
Chemical speciation can be visualized with the aid of Hyperquad simulation and speciation (HySS2009) [11,12], a utility program for the investigation of equilibria involving soluble and partially soluble species, available in electronic format at http://www.hyperquad.co.uk/hyss.htm(accessed on 20 September 2023).HySS, a computer program compatible with the Windows operating system, allows both the simulation of titration curves and the obtaining of speciation diagrams.Input data for speciation are the reactants' total concentrations and association constants.The distribution of free ions and ion pairs present in the NaCl-NaHCO3-H2O ternary system, at different pH values, is presented in   In the 4-12 pH range (Figure 8), starting from H2CO3, the CO3 2− species are present in the form of HCO3 − , NaCO3 − , and free CO3 2− , meaning that at pH = 12, the carbonate is associated with sodium in a proportion of 54%, while 44% CO3 2− is free.It can be seen from Figure 9 that in the NaCl-NaHCO3-H2O system, as the pH increases, the percentage of associated sodium in the form of NaCO3 − and NaHCO3 increases, reaching up to 12% at pH 12, and the unassociated sodium decreases to 88%.At pH 8, it is observed that there is 90% free sodium, and 10% associated sodium in the form of NaHCO3.
Class A glassware was used for the preparation of all solutions.
The solutions pH values were measured with a P901 CONSORT portable pH-meter (Turnhout, Belgium), equipped with DG111 (Mettler Toledo, Schwerzenbach, Switzerland) combined glass electrode.A 2-point calibration routine with 7.00 and 10.00 pH buffers was employed to adjust the slope and the asymmetry potential of the measuring device.The hydrogen ion activity scale (pH = −log aH) was used.
Density determinations were performed with a Gay-Lussac pycnometer (Roth).The pycnometer was calibrated with double-distilled water.In the 4-12 pH range (Figure 8), starting from H 2 CO 3 , the CO 3 2− species are present in the form of HCO 3 − , NaCO 3 − , and free CO 3 2− , meaning that at pH = 12, the carbonate is associated with sodium in a proportion of 54%, while 44% CO 3 2− is free.It can be seen from Figure 9 that in the NaCl-NaHCO 3 -H 2 O system, as the pH increases, the percentage of associated sodium in the form of NaCO 3 − and NaHCO 3 increases, reaching up to 12% at pH 12, and the unassociated sodium decreases to 88%.At pH 8, it is observed that there is 90% free sodium, and 10% associated sodium in the form of NaHCO 3 .
Class A glassware was used for the preparation of all solutions.
The solutions pH values were measured with a P901 CONSORT portable pH-meter (Turnhout, Belgium), equipped with DG111 (Mettler Toledo, Schwerzenbach, Switzerland) combined glass electrode.A 2-point calibration routine with 7.00 and 10.00 pH buffers was employed to adjust the slope and the asymmetry potential of the measuring device.The hydrogen ion activity scale (pH = −log a H ) was used.
Density determinations were performed with a Gay-Lussac pycnometer (Roth).The pycnometer was calibrated with double-distilled water.
All determinations were carried out at 25 • C. The cell temperature was monitored with a calibrated thermometer with an expanded uncertainty of 0.5 • C for a coverage factor k = 2 and a probability P = 95%.The mean activity coefficients of sodium chloride, both in pure sodium chloride and in a mixture of electrolytes, were determined according to Bates and Dickson [29], through two consecutive measurements of the electromotive voltage in cells with a liquid junction consisting of the indicating electrode and the same reference electrode and immersed in the electrolytic medium.
For the correct measurement of the mean activity coefficients of sodium chloride, the sodium and chloride ion-selective electrodes were calibrated each time before use in sodium chloride standard solutions, using two electrochemical cells of the I and II type [29].
The electromotive force corresponding to each cell is described by Equations ( 4) and ( 5): where: E 0 Na + , E 0 Cl − -the standard potentials of cells I and II; S Na + , S Cl − -the slopes of the sodium and chloride electrodes, respectively; E j -the junction potential; E ref -the potential of the reference electrode; a Na + -sodium ion activity; a Cl − -chloride ion activity.The difference between the electromotive force of the I and II cells is equal to the electromotive force of cell III, without liquid junction: From a stock solution of 1 mol/L sodium chloride, a series of fourteen standard solutions in the 0.0001-0.5 mol/L concentration range was prepared.The electromotive force of the two cells I and II were measured consecutively, practically under the same conditions: the same solution, the same reference electrode, the same temperature, with a time difference between determinations of maximum 5 min.The temperature of the solutions was kept constant at 25 • C, with a precision of ±0.1 • C.
Calibration data were fitted to first-degree equations: a NaCl = log m NaCl × γ 0 NaCl (7) where: m NaCl − -sodium chloride molal concentration, γ 0 NaCl − -mean activity coefficient of sodium chloride in pure sodium chloride, S Na + , S Cl − -slopes of the sodium and chloride electrodes, respectively, E 0 Na + , E 0 Cl − -standard potentials of cells I and II, E 0 NaCl -standard potential of cell III.At concentrations higher than 0.01 mol/L, the mean activity coefficients of sodium chloride were calculated with the extended Debye-Hückel Equation (10).Electrochemical cells of the type IV and V were used to determine the mean activity coefficients of sodium chloride (γ NaCl ) in the NaCl-NaHCO 3 -H 2 O ternary system at 0.1-0.5 mol/kg total ionic strength.m NaHCO 3 -the sodium acidic carbonate molal concentration.The total ionic strength of the medium is described by the equation: The mean activity coefficients of sodium chloride in the presence of sodium acidic carbonate (γ NaCl ) were calculated with the equation: The association constants of sodium ions with carbonate ions (K 1C ) and acid carbonate ions (K 2C ), respectively, were determined in a series of NaCl-NaHCO 3 -H 2 O solutions of different compositions at 0.1-0.5 mol/kg total ionic in the 8.34-8.95pH range (Table 1).
An iterative procedure was used to calculate the association constants [17][18][19].In the first iteration, the effective ionic strength, I E , was assumed to be equal to the total ionic strength, I T , and the constant K 1C was considered equal to 3.162 [18].Iterations were developed using a computer and were repeated until two successive values of I E did not differ by more than 0.001 units [18].
Starting from the hypothesis that at constant total ionic strength, the modification of the mean activity coefficients occurs only as a result of ion pairing, it can be written that NaCl can be expressed as the mean activity coefficient of sodium chloride at an effective ionic strength, I E = I T .This ionic strength is expressed based on the equilibrium concentrations of the various ionic species: The equilibrium concentration values are obtained from the protonation and association equilibria: [NaHCO 3 ] = K 2C Na + ] F * HCO − 3 (21) and from the mass balance equations for the three components: HCO − 3 + CO 2− 3 where: K a 2 C -the second dissociation constant of H 2 CO 3 ; K 1C -the association constant of sodium ions with carbonate ions; K 2C -the association constant of sodium ions with acidic carbonate ions.
determined the stoichiometric association constant of NaCO 3 − in the 10-40 • C temperature range and the 0.20-0.62mol/L ionic strength range in NaCl-[(CH 3 ) 4 N] 2 BPDS-Me 4 NCl medium, and the stoichiometric association constant of NaHCO 3 in the 10-40 • C temperature range and 0.20-0.62mol/L ionic strength range in NaCl-Me 4 NCl medium.

Figure 1 .
Figure 1.Compliance with Harned's rule verification for ionic strengths in the 0.1-0.5 mol/kg range.

Figure 1 .
Figure 1.Compliance with Harned's rule verification for ionic strengths in the 0.1-0.5 mol/kg range.

Figure 2 .
Figure 2. Variation in Harned's coefficient with total ionic strength for the NaCl-NaHCO3-H2O ternary system at 25 °C.

Figure 3 .Figure 2 .
Figure 3. Variation in mean activity coefficients with composition in the NaCl-NaHCO3-H2O ternary system, at different constant total ionic strength values.

Figure 2 .
Figure 2. Variation in Harned's coefficient with total ionic strength for the NaCl-NaHCO3-H2O ternary system at 25 °C.

Figure 3 .Figure 3 .
Figure 3. Variation in mean activity coefficients with composition in the NaCl-NaHCO3-H2O ternary system, at different constant total ionic strength values.

Figure 6 .
Figure 6.Variation in the association constant of sodium ions with carbonate ion ( ) values with sodium chloride concentration at 0.1-0.5 mol/kg total ionic strength range.

Figure 6 .
Figure 6.Variation in the association constant of sodium ions with carbonate ion (K 1C ) values with sodium chloride concentration at 0.1-0.5 mol/kg total ionic strength range.

Figure 7 .
Figure 7. Variation in the association constant of sodium ions with acidic carbonate ions values with sodium chloride concentration at 0.1-0.5 mol/kg total ionic strength range.

Table 1 .
Existing data on the association constants of sodium with carbonate log K 1C and acid carbonate ions log K 2C , in different media, at 25 • C, available in the literature.

Table 2 .
The experimental values of mean activity coefficients of sodium chloride in the ternary system, γ NaCl , and NaCl-H 2 O system, γ 0 NaCl , against the ionic strength fraction of sodium chloride (X NaCl ), at constant total ionic strength, I T , and 25 • C.

Table 3 .
Values of Harned's coefficient in the NaCl-NaHCO3-H2O ternary system, at different total ionic strengths, and 25 °C.

Table 3 .
Values of Harned's coefficient in the NaCl-NaHCO 3 -H 2 O ternary system, at different total ionic strengths, and 25 • C.

Table 4 .
Experimental values of association constant of sodium ions with carbonate ions (K 1C ) and the association constant of sodium ions with acidic carbonate ions (K 2C ).

Table 4 .
Experimental values of association constant of sodium ions with carbonate ions ( ) and the association constant of sodium ions with acidic carbonate ions ( ).