Abstract
The solubilities and the relevant physicochemical properties of the ternary systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O at 308.15 K were investigated using an isothermal dissolution method. It was found that there is one invariant point, two univariant curves, and two crystallization regions of the systems. The systems belong to a simple co-saturated type, and neither double salts nor solid solutions were found. Based on the extended HW model and its temperature-dependent equations, the single-salt Pitzer parameters β(0), β(1), β(2) and CØ for MgCl2, MgSO4, and Mg(B6O7)(OH)6, the mixed ion-interaction parameters θCl,B6O10, θSO4,B6O10, ΨMg,Cl,B6O10, ΨMg,SO4,B6O10 of the systems at 308.15 K were fitted, In addition, the average equilibrium constants of the stable equilibrium solids at 308.15 K were obtained by a method using the activity product constant. Then the solubilities of the ternary systems are calculated. The calculated solubilities agree well with the experimental values.
Phase equilibrium; Phase diagram; Solubility; Pitzer modeling; Mcallisterite
THERMODYNAMICS
Solid-liquid stable phase equilibria of the ternary systems MgCl2 + MgB6O10+ H2O AND MgSO4 + MgB6O10 + H2O at 308.15 K
Lingzong MengI,II,* * To whom correspondence should be addressed ; Dan LiIII
ISchool of Chemistry and Chemical Engineering, Phone/Fax: + 86-539-8766600, Linyi University, Linyi, 276005, China.E-mail: menglingzong@lyu.edu.cn
IITianjin Key Laboratory of Marine Resources and Chemistry, Tianjin University of Science and Technology, Tianjin, 300457, China.
IIISchool of Media Communications, Linyi University, Linyi 276005, China
ABSTRACT
The solubilities and the relevant physicochemical properties of the ternary systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O at 308.15 K were investigated using an isothermal dissolution method. It was found that there is one invariant point, two univariant curves, and two crystallization regions of the systems. The systems belong to a simple co-saturated type, and neither double salts nor solid solutions were found. Based on the extended HW model and its temperature-dependent equations, the single-salt Pitzer parameters β(0), β(1), β(2) and CØ for MgCl2, MgSO4, and Mg(B6O7)(OH)6, the mixed ion-interaction parameters θCl,B6O10, θSO4,B6O10, ΨMg,Cl,B6O10, ΨMg,SO4,B6O10 of the systems at 308.15 K were fitted, In addition, the average equilibrium constants of the stable equilibrium solids at 308.15 K were obtained by a method using the activity product constant. Then the solubilities of the ternary systems are calculated. The calculated solubilities agree well with the experimental values.
Keywords: Phase equilibrium; Phase diagram; Solubility; Pitzer modeling; Mcallisterite.
INTRODUCTION
There are more than one thousand salt lakes in China, and a number of salt lakes with an abundance of lithium and boron resources are widely distributed in the Qaidam Basin of Qinghai-Tibet Plateau, China (Zheng et al., 1988). These brines belong mostly to the complex system Li + Na + K + Mg + Cl +SO4 + borate + H2O. It is well known that the thermodynamic phase equilibria and phase diagrams play an important role in exploiting the brine resources and describing the geochemical behavior of the brine and mineral system. Therefore, simulative experimental studies on stable phase equilibria containing magnesium and boron are essential for predicting the path of mineral crystallization for the effective separation and purification of the boron-containing mixture salts effectively. Although it has been pointed out that the structure of boron ion of mcallisterite is B6O7(OH)62-, B6O102- is used in this paper for convenience (Li and Gao, 1993).
The ternary systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O are subsystems of the seven-component system. Although the stable phase equilibrium of the system MgCl2 + MgB6O10 + H2O at 298.15 K and 323.15 K, the system MgSO4 + MgB6O10 + H2O at 323.15 K, and the system MgCl2 + MgSO4 +MgB6O10 + H2O at 323.15 K have been reported (Bi et al., 1997; Meng et al., 2012; Meng et al., 2011), the stable phase equilibria of the ternary systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O at 308.15 K are not reported in the literature. In this paper, the stable solubilities and the relevant physicochemical properties (densities and pH) of the ternary systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O at 308.15 K, as well as the predicted solubilities using the Harvie-Weare (HW) solubility modeling approach, which incorporates the Pitzer specific interaction formalism are presented.
EXPERIMENTAL SECTION
Apparatus and Reagents
A thermostatic shaker (model THZ-82) whose temperature can be controlled to ±0.1 K was used for the measurement of phase equilibria. When the solution temperature in the container was under (308.15 ± 0.1) K, the apparatus for controlling the temperature formed a circuit, and the heating apparatus began to heat with temperature at 308.15 K. Conversely, the circuit was broken, and the heating apparatus stopped working. Therefore, the temperature in the thermostatic shaker was always maintained at (308.15 ± 0.1) K. The solid phase minerals were identified using Schreinemaker's method (Song, 1991) and an X-ray diffractometer (X'pert PRO, Spectris. Pte. Ltd., The Netherlands). The chemicals used were of analytical grade and obtained from Tianjin Kermel Chemical Reagent Ltd. [epsomite (MgSO4·7H2O, 0.99 by mass fraction), bischofite (MgCl2·6H2O, 0.99 by mass fraction)], and were re-crystallized before use. Mcallisterite (MgB6O10·7.5H2O, 0.99 by mass fraction) was synthesized in our laboratory with the method described in the literature (Meng et al., 2010). Doubly deionized water (DDW) with conductivity less than 1.2×10-4 S·m-1 at 298.15 K was used to prepare the series of artificial synthesized brines and for chemical analysis.
Experimental Methods
The isothermal dissolution method was used in this study. The series of complexes of the ternary systems were loaded into clean polyethylene bottles and capped tightly. The bottles were placed in the thermostatic rotary shaker, whose temperature was controlled to (308.15 ± 0.1) K, and rotated at 120 rpm to accelerate the equilibration of those complexes. The rotary system rested for 2 h before sampling, and a 1.0 mL sample of the clarified solution was taken from each polyethylene bottle with a pipette at regular intervals and then used for a refractive index measurement. If the difference between two refractive index measurements was within ± 0.0002, then the equilibrium state was assumed. Otherwise, the solution was rotated continually until the equilibrium state was achieved. Generally, it took approximately 90 days to reach the equilibrium state. Twice, at different times, a sample of approximately 5.0 mL was taken, weighed and diluted to a final volume of 250 mL using a volumetric flask filled with DDW. Then, a quantitative analysis was performed. If the difference between the compositions of the two liquid phases was within ± 0.3% in mass fraction, then the solubility of this equilibrium point, which was the average of the experimental data, was obtained. Some other solution samples were taken to determine the physicochemical properties. Meanwhile, some of the wet residue was separated from the solution, and one portion of the wet residue was weighed and dissolved to a final volume of 250 mL using a volumetric flask filled with DDW for chemical analysis by Schreinemaker's wet residue method (i.e., the solid phase mineral point lies on the extension of the line on the phase diagram through the composition points of the liquid phase and the wet residue) (Song, 1991). Another portion was dried, ground into powder and analysed by X-ray diffraction. The solid phase could then be identified.
Analytical Methods
The solid phases were analyzed by Schreinemaker's method and X-ray diffraction; magnesium borate in the two systems is mcallisterite (MgB6O10·7.5H2O). It shows that mcallisterite is not converted, so B6O102- can be used to express the boron ions in the aqueous and solid phases. The Mg2+ ion concentration was determined by titration with EDTA standard solution in the presence of the indicator Eriochrome Black-T. The average deviation of the determination was less than ± 0.3% (by mass). For the ternary system MgCl2 + MgB6O10 + H2O, the Cl- ion concentration was determined by titration with a standard solution of Hg(NO3)2 (uncertainty ± 0.3%) and the concentration of B6O102- was evaluated using an ion balance (Zhai et al., 1988). For the ternary system MgSO4 + MgB6O10 + H2O, the concentrations of B6O102- were determined by basic titration in the presence of mannitol (uncertainty ± 0.3%) and SO42- was calculated by subtraction via charge balance (Zhai et al., 1988). The densities (ρ) were measured by means of a pyknometer with a precision of ± 0.2 × 10-6 kg·m-3. A PHS-3C precision pH meter (Shanghai Precision & Scientific Instrument Co. Ltd.) was used to measure the pH of the equilibrium aqueous solutions (precision of ± 0.01). The pH meter was calibrated with standard buffer solutions of a mixed phosphate of potassium dihydrogen phosphate and sodium dihydrogen phosphate (pH = 6.84), as well as potassium acid phthalate (pH = 4.02). All the measurements were maintained at (308.15 ± 0.1) K through control of the thermostat.
RESULTS AND DISCUSSION
The solubility data and the relevant physicochemical property data of the aqueous systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O at 308.15 K are presented in Tables 1 and 2. The ionic composition in the equilibrium solution of the isothermal dissolution is expressed as mass fraction. Using the experimental results in Tables 1 and 2, the phase diagrams of these systems at 308.15 K were plotted, as shown in Figures 1 and 2. Both of the phase diagrams of the ternary systems in Figures 1 and 2 consist of one invariant point, two univariant curves and two crystallization fields. In Figure 1, the point labeled E is invariant for bischofite (MgCl2·6H2O) and mcallisterite (MgB6O10·7.5H2O), and the compositions of MgCl2, MgB6O10 in the liquid phase with mass fraction (wb × 102) are 36.00 and 0.77, respectively. In Figure 2, the point labeled E is invariant for epsomite (MgSO4·7H2O) and mcallisterite (MgB6O10·7.5H2O), and the compositions of MgSO4, MgB6O10 in the liquid phase with mass fraction (wb × 102) are 29.29 and 1.33, respectively. Both of the ternary systems belong to a simple eutectic type, and neither double salts nor solid solutions were found.
Using the experimental results in Tables 1 and 2, the relationships of the solution physicochemical properties (densities and pH) with mass fraction of magnesium chloride or magnesium sulfate are plotted in Figures 3 and 4. It was found that the solution densities of the ternary systems were changed regularly with the increase of mass fraction of magnesium chloride or magnesium sulfate and reached the maximum value at the invariant, while the pH values of the stable equilibrium aqueous solution generally decreased gradually with increasing magnesium chloride or magnesium sulfate concentration, because the solution with magnesium chloride or magnesium sulfate is acidic and the pH values will decrease with increasing concentration. Using the diagrams (Figures 3 and 4), the physicochemical properties (density and refractive index) could be used to roughly estimate the magnesium concentration in the ternary system.
SOLUBILITY PREDICTION
Ion-Interaction Model
Pitzer has developed an ion-interaction model and published a series of papers (Pitzer 1973; Pitzer 1992), which gave a set of expressions for the osmotic coefficients of the solution and mean activity coefficient of electrolytes in the solution. Single ion activity coefficients, which are more convenient to use in solubility calculations, are given in the Harvie and Weare modeling approach (Harvie and Weare, 1980; Harvie et al., 1984). Because these equations are based on the excess free energy, all the activity expressions are consistent and ready for application to different types of data (e.g., solubility data) for parameter evaluations and calculations of other thermodynamic functions. Using the activity coefficients and the solubility products of the equilibrium solid phases allowed us to identify the coexisting solid phases and their compositions at equilibrium. Additional work has centered on developing variable temperature models, which will increase the applicability to a number of diverse geochemical systems (Møller, 1988; Greenberg and Møller, 1989; Pabalan and Pitzer, 1987).The following equations are the main expressions related to the solubility calculation of this model.
In expressions (1) to (3), M, c and c′ represent cations, and X, a and a′ represent anions, respectively. Additionally, mi and γi represent the molality (mol·kg−1) and the activity coefficient of the ions, respectively, zi refers to the valence of the ions, and ø refers to the permeability coefficient. The other symbols in expressions (1) to (3), for example, F, C, Z, AØ, Ψ, Φ, BØ, and B, are all described in the references (Harvie and Weare, 1980; Harvie et al., 1984).
Solubility equilibrium constants are used to calculate solubility, i.e., a salt solution is saturated at a given temperature and pressure when the ion activity product is equal to the solubility equilibrium constant. For a hydrated salt , the solubility equilibrium constant (Ksp) at a definite temperature and pressure for the dissolution reaction is expressed as follows:
The activity of water is related to the osmotic coefficient ø, by the equation
where, Mw is the molar mass of water and the sum covers all solute species. The activity coefficients, lnγi, and the osmotic coefficients ø, can be calculated using the the Pitzer ion-interaction model and the extended HW model.
According to the Pitzer and extended HW models, the activity and osmotic coefficients are parametric functions of β(0), β(1), β(2), CØ, θcc′ and ψcc′a, β(0), β(1), β(2), CØ are the parameters of a single salt, θcc′ represents the interaction of the two ions with the same sign and ψcc′a represents the interactions among the three ions, in which the sign of the third one is different from the first two ions. Solubility of the systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O at 308.15 K could be calculated from Eqs. (1)-(6) when the Pitzer parameters and Ksp of the system are all known.
Model Parameterization
A high-temperature thermo-chemical model for the system (Na + K + Mg + Cl + SO4 + OH + H2O) was presented by Pabalan and Pitzer (1987), including a range of temperature for the model. Take MgCl2 for example, the temperature-dependent equations for MgCl2 are as follows:
where P(T) in Eq. (7) refers to the binary parameters β(0), β(1),β(2) and CØ of MgCl2 at 308.15 K of the solution model and the values of the fitting constant terms of Q1, Q2,and Q3 were introduced (Pabalan and Pitzer, 1987). The DebyeHückel parameter AØ and the binary parameters of MgSO4 at 308.15 K can be calculated based on the temperature-dependent equations presented in the literature (Møller, 1988; Pabalan and Pitzer, 1987). The binary parameters β(0), β(1),β(2)and CØ of Mg(B6O7)(OH)6 and the mixing ion-interaction parameters ,
,
,
were available by a multiple linear regression procedure applied to the solubility data of the corresponding ternary system in this paper. All the parameters used in the prediction are presented in Table 3. According to the experimental solubility data in Tables 1 and 2, and Pitzer parameters in Table 3, the average equilibrium constants (lnKaver) for the solid phases were fitted in this study. These average equilibrium constants (lnKaver) of equilibrium solid salts for MgCl2·6H2O, MgSO4·7H2O, and Mg(B6O7)(OH)6·4.5H2O correspond to 10.7075, -4.114139 and -9.0057 at (308.15 ± 0.1) K. The uncertainties of the average equilibrium constants for MgCl2·6H2O, MgSO4·7H2O are ± 0.0278 and ± 0.0821, whereas the uncertainty for Mg(B6O7)(OH)6·4.5H2O is larger than ± 1.0000. The boron ions in the solution are very complicated, there are many different kinds of boron ions (Li and Gao, 1993), so the equilibrium constant of mcallisterite has a comparatively large uncertainty; a similar situation has been described for Li2(B4O5)(OH)4·H2O (Meng, 2012).
Calculated Solubility
To determine the accuracy of the parameters and the solubility equilibrium constants of the solids, the solubilities of the two ternary systems at 308.15 K have been calculated based on the Pitzer ion-interaction model and the extended HW model for aqueous electrolyte solutions. It should be noted that solutes below their saturated solution molalities were fixed at their experimental values and the saturated solutes' solubilities were calculated. The phase diagrams according to the experimental and calculated data were constructed and are shown in Figures 5 and 6, respectively.
It is shown that the calculated and the experimental data for the ternary system MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O are in good agreement. These results indicate that the model obtained in this work that combines the Pitzer parameters and the equilibrium constants is reliable and is capable of predicting equilibria in the system studied.
CONCLUSIONS
The experimental solubility data and the relevant physicochemical property data of the aqueous systems MgCl2 + MgB6O10 + H2O and MgSO4 + MgB6O10 + H2O at 308.15 K were determined. The experimental stable phase diagrams and the variation of the physicochemical properties with composition were constructed for the first time. The results show that the systems belong to a simple co-saturated type and neither solid solutions nor double salts were formed in the two ternary systems. The single-salt Pitzer parameters, the mixing ion interaction parameters, the equilibrium constant of the minerals were fitted, and the DebyeHückel parameter AØ in the systems were calculated based on the temperature-dependent equations presented in the literature (Møller, 1988; Pabalan and Pitzer, 1987). Solubility calculations for the ternary systems at 308.15 K have been made using the Pitzer parameters and average equilibrium constants of the solids. The calculated solubilities agree well with experimental values. These results indicate that the parameters and the equilibrium constants obtained in this work are reliable.
ACKNOWLEDGEMENTS
Financial support from the Foundation of Tianjin Key Laboratory of Marine Resources and Chemistry at Tianjin University of Science & Technology (Grant 201101) and the Science Foundation of Linyi University (Grant BS201005) are greatly acknowledged.
NOMENCLATURE
AØ
DebyeHückel parameter
BØ
parametric functions of Pitzer parameter
B
parametric functions of Pitzer parameter
B'
parametric functions of Pitzer parameter
CØ
single-salt Pitzer parameter
C
parametric functions of Pitzer parameter
fγ
parametric functions of Pitzer parameter
g(x), g'(x)
parametric functions of Pitzer parameter
I
ionic strength
mol/kg
Ksp
solubility equilibrium constant
m
molality
mol/kg
M,C
cation
n
neutral molecules
X, A
anion
z
valence of ions
Z
parametric functions of the molality of species i
Greek Letters
αw
water activity
β(0)
single-salt Pitzer parameter
β(1)
single-salt Pitzer parameter
β(2)
single-salt Pitzer parameter
υ
number of ions in the molecule
γ
activity coefficient of the ions
ø
osmotic coefficients
θ
mixing ion interaction Pitzer parameter
Ψ
mixing ion interaction Pitzer parameter
Submitted: January 15, 2013
Revised: March 5, 2013
Accepted: April 26, 2013
- Bi, W. B., Sun, B., Song, P. S. and Zhu, G. Q., Study on phase diagram of Mg2+/Cl-, borate-H2O ternary system at 25 °C. J. Salt Lake Res., 5, 42-46 (1997).
- Deng, T. L., Phase equilibrium for the aqueous system containing lithium, sodium, potassium, chloride, and borate ions at 298.15 K. J. Chem. Eng. Data. 49, 1295-1299 (2004).
- Greenberg, J. P. and Møller, N., The prediction of mineral solubilities in natural waters: A chemical equilibrium model for the Na-K-Ca-Cl-SO4-H2O system to high concentration from 0 to 250 °C. Geochim. Cosmochim. Acta, 53, 2503-2518 (1989).
- Harvie, C. E., Møller, N. and Weare, J. H., The prediction of mineral solubilities in nature waters: The Na-K-Mg- Ca-H-SO4-OH-HCO3-CO3-H2O system to high ionic strengths at 25 °C. Geochim. Cosmochim. Acta, 48, 723-751 (1984).
- Harvie, C. E. and Weare, J. H., The prediction of mineral solubilties in nature waters: The Na-K-Mg- Cl-SO4-H2O system from zero to high concentration at 25 °C. Geochim. Cosmochim. Acta, 44, 981-997 (1980).
- Li, J. and Gao, S. Y., Chemistry of borates. J. Salt Lake Res., 1, 62-67 (1993).
- Meng, L. Z., Solubility prediction for LiClLi2SO4Li2B4O7H2O system at 298.15 K using ion-interaction model. J. Chem. Eng. Jpn., 45, 563-567 (2012).
- Meng, L. Z., Deng, T. L., Duan, C. W., Li, D. and Guo, Y. F., A synthesis method for mcallisterite. World Sci-tech R&D (China), 32, 825-826 (2010).
- Meng, L. Z., Deng, T. L., Guo, Y. F. and Li, D., The phase diagrams of the systems (MgCl2 + MgB6O10 + H2O) and (MgSO4 + MgB6O10 + H2O) at 323.15 K and Pitzer model. Russ. J. Phys. Chem. A, 86, 1526-1532 (2012).
- Meng, L. Z., Li, D., Guo, Y. F. and Deng, T. L., Stable phase equilibrium of the aqueous quaternary system (MgCl2 + MgSO4 + MgB6O10 + H2O) at 323.15 K. J. Chem. Eng. Data, 56, 5060-5065 (2011).
- Møller, N., The prediction of mineral solubilities in natural waters: A chemical equilibrium model for the Na-Ca-Cl-SO4-H2O system to high temperature and concentration. Geochim. Cosmochim. Acta, 52, 821-837 (1988).
- Pabalan, R. T. and Pitzer, K. S., Thermodynamics of concentrated electrolyte mixtures and the prediction of mineral solubilities to high temperatures for mixtures in the system Na-K-Mg-Cl-SO4-OH-H2O. Geochim. Cosmochim. Acta, 51, 2429-2443 (1987).
- Pitzer, K. S., Thermodynamics of electrolytes I: Theoretical general equation. J. Phys. Chem., 77, 268-277 (1973).
- Pitzer, K. S., Activity coefficients in electrolyte solutions. 2nd Ed., CRC Press, London (1992).
- Song, P. S., Studies on the application of the wet residues method in phase Equilibrium for the salt-water systems. J. Salt Lake Res., (China), 1, 42-46 (1991).
- Zhai, Z. X., Cheng, Z. B., Hu, F., Wu, G. L., Xiao, Y. K. and Liu, F. M., Analytical Methods of Brines and Salts. 2nd Ed., Qinghai Institute of Salt Lakes, Chinese Academy of Science. Chin. Sci. Press, Beijing (1988).
- Zheng, X. Y., Tang, Y. and Xu, C., Tibet Saline Lake. Science Press, Beijing, China (1988).
Publication Dates
-
Publication in this collection
20 Mar 2014 -
Date of issue
Mar 2014
History
-
Received
15 Jan 2013 -
Accepted
26 Apr 2013 -
Reviewed
05 Mar 2013