Six different members of the (BaxSr1-x)HAsO4·H2O solid solution were prepared and characterized, and then dissolution of the synthetic solids was studied at 25°C and pH 2 in a series of batch experiments for 4320 h. With the increase in the Ba/(Ba+Sr) mole fraction of the solids, the aqueous pH decreased and the aqueous barium concentration increased. The aqueous strontium and arsenic concentrations had the highest values at Ba/(Ba+Sr)=0.21. The solubility products for BaHAsO4·H2O and SrHAsO4·H2O were calculated to be 10-5.52 and 10-4.62, respectively. The corresponding free energies of formation were determined to be −1543.99±0.18 kJ/mol and −1537.94±0.02 kJ/mol. The solid solution had a minimum solubility product of 10-5.61 at Ba/(Ba+Sr)=0.82. The Guggenheim coefficients were determined to be ao=1.55 and a1=-4.35. The Lippmann diagram was a typical Lippmann diagram for a nonideal solid solution with a negative enthalpy of mixing. The system shows an “alyotropic” minimum at the aqueous Ba/(Ba+Sr) activity fraction of 0.87 where the solutus and solidus curves meet. At the end of the dissolution experiment, the dissolution followed the saturation curve for the pure endmember BaHAsO4·H2O and approached the intersection with the minimum stoichiometric saturation curve on the Lippmann diagram.
1. Introduction
Arsenic in natural waters is a worldwide problem. Arsenic has been known from antiquity to be highly toxic for animals and the majority of plants [1, 2]. Although arsenic has been classified at the top of the priority list of the most hazardous substances [3, 4], the crystal structures and the solubility as well as the thermodynamic properties of numerous arsenates remain poorly determined. This is an important handicap because an in-depth study of the arsenate behaviour in soils, sediments, and natural waters that have been subjected to pollution requires a precise knowledge of the possible precipitating phases, their crystal chemistry, and their thermodynamic properties. Furthermore, in multicomponent aqueous systems, the precipitation of solid solutions is always a possibility such that the study of arsenate solid solutions involving substitution between atoms of similar size and character is worthwhile [5, 6]. Previous studies indicate that, at high concentrations, compounds in this series could be a limiting phase for arsenic in natural aqueous environments [7, 8]. In solution, Ba tends to associate with As at pH 7.47–7.66, forming BaHAsO4·H2O and Ba3(AsO4)2 [9]. This chemical interaction is used to remove As from aqueous solutions, but its efficiency depends on the solubility of both elements, which is altered by the physicochemical conditions of water [10, 11]. Tiruta-Barna et al. [12] also found that the leaching behaviour of arsenic from a compacted coal fly ash was controlled by the weak soluble phase BaHAsO4·H2O. The only available data on this series correspond to the strontium endmember, whose space group (Pbca), cell parameters, and atomic positions of Sr and As were determined by Binas and Boll-Dornberger [13]. From powder diffraction data, Martin et al. [14] determine the same space group for the barium endmember and suggest that SrHAsO4·H2O and BaHAsO4·H2O are isomorphous and therefore good candidates to form solid solutions. The previous data indicate that SrHAsO4·H2O and BaHAsO4·H2O form a complete solid solution where Sr2+ ions substitute for Ba2+. Both endmembers crystallize in the same orthorhombic Pbca space group, with lattice parameters that vary significantly with composition [4].
The determination of the equilibrium behaviour in SS-AS systems requires knowledge of both the endmember solubility products and the degree of nonideality of the solid solution [5]. Unfortunately, the lack of thermodynamic data for arsenates is even greater than the lack of crystal-chemical data. The Ksp value for BaHAsO4·H2O was determined to be 10−4.70 by Robins [10], 10−24.64 by Essington [15], 10−5.31 by Itoh and Tozawa [16], 10−0.8 by Orellana et al. [17], 10−5.51 by Davis [18], and 10−5.60 by Zhu et al. [9], which showed a great inconsistency of the currently accepted solubility data. Moreover, the solubility of the strontium endmember as well as the thermodynamic properties of the solid solution is unknown [4]. The fact that there is a significant preferential partitioning of barium toward the solid phase seems to indicate a lesser solubility (≈ one order of magnitude) for the strontium endmember. A thermodynamic character of these compounds needs to be confirmed [4].
In the present study, a series of the (BaxSr1-x)HAsO4·H2O solid solution with different Ba/(Ba + Sr) atomic ratios was prepared by a precipitation method. The resulting solid solution particles were characterized by various techniques. This paper reports the results of a study that monitors the dissolution and release of constituent elements from synthetic (BaxSr1-x)HAsO4·H2O solid solutions using batch dissolution experiments. The solid solution aqueous solution reaction paths are also discussed using the Lippmann diagram to evaluate the potential impact of such solid solutions on the mobility of arsenic in the environment.
2. Experimental Methods2.1. Solid Preparation and Characterization
The experimental details for the preparation of the samples by precipitation were based on the following equation: X2+ + HAO42- + H2O = XHAsO4·H2O, where X = Ba2+ or Sr2+. The (BaxSr1-x)HAsO4·H2O solid solutions were synthesized by controlled mixing of a solution of 100 mL 0.5 M Na3AsO4 and a solution of 100 mL 0.5 M Ba(ClO4)2 and Sr(ClO4)2, so that the (Ba + Sr)/As molar ratio in the mixed solution was 1.00. The amounts of Ba(ClO4)2 and Sr(ClO4)2 were varied in individual syntheses to obtain synthetic solids with different mole fractions of Ba/(Ba + Sr) (Table 1). Reagent grade chemicals and ultrapure water were used for the synthesis and all experiments. The initial solutions were slowly mixed in a covered beaker in a course of 10 minutes at room temperature (25±1°C). The resulting solutions were kept at 70°C and stirred at a moderate rate (100 rpm) using a stir bar. After one week to crystallize, the precipitates were allowed to settle. The resultant precipitates were then washed thoroughly with ultrapure water and dried for 24 h at <110°C to avoid decomposition of the solid samples obtained.
Summary of synthesis and composition of the (BaxSr1-x)HAsO4·H2O solid solution.
Volumes of the precursors (mL)
Solid solution
Ba(ClO4)2
Sr(ClO4)2
Na3AsO4
Molecular formula
XBa
0.5 mol/L
0.5 mol/L
0.5 mol/L
0
100
100
SrHAsO4·H2O
0.00
20
80
100
(Ba0.21Sr0.79)HAsO4·H2O
0.21
40
60
100
(Ba0.42Sr0.58)HAsO4·H2O
0.42
60
40
100
(Ba0.61Sr0.39)HAsO4·H2O
0.61
80
20
100
(Ba0.82Sr0.18)HAsO4·H2O
0.82
100
0
100
BaHAsO4·H2O
1.00
The composition of the solid sample was determined. 50 mg of sample was digested in 10 mL of 2 M HNO3 solution and then diluted to 50 mL with 2% HNO3 solution. It was analysed for Ba, Sr, and As using an inductively coupled plasma atomic emission spectrometer (ICP-OES, Perkin Elmer Optima 7000DV). The synthetic solids were also characterised by powder X-ray diffraction (XRD) with an X’Pert PRO diffractometer using Cu Kα radiation (40 kV and 40 mA). Crystallographic identification of the synthesised solids was accomplished by comparing the experimental XRD pattern to standard compiled by the International Centre for Diffraction Data (ICDD), which were card 00-023-0823 for BaHAsO4·H2O and card 01-074-1622 for SrHAsO4·H2O. The morphology was analysed by scanning electron microscopy (SEM, Joel JSM-6380LV). Infrared transmission spectra (KBr) were recorded over the range of 4000–400 cm−1 using a Fourier transformed infrared spectrophotometer (FT-IR, Nicolet Nexus 470 FT-IR).
2.2. Dissolution Experiments
1.5 g of the synthetic (BaxSr1-x)HAsO4·H2O solid was placed in 250 mL polypropylene bottle. 150 mL of 0.01 M HNO3 (initial pH 2) was added to each bottle. The bottles were capped and placed in a temperature-controlled water bath (25±1°C). Water samples (3 mL) were taken from each bottle on 18 occasions (1 h, 3 h, 6 h, 12 h, 1 d, 2 d, 3 d, 5 d, 10 d, 20 d, 30 d, 40 d, 50 d, 60 d, 75 d, 90 d, 120 d, and 180 d). After each sampling, the sample volume was replaced with an equivalent amount of ultrapure water. The samples were filtered using 0.20 μm pore diameter membrane filters and stabilised with 0.2% HNO3 in 25 mL volumetric flask. Ba, Sr, and As were analysed by using an ICP-OES. After 180 d of dissolution, the solid samples were taken from each bottle, washed, dried, and characterised using XRD, FT-IR, and SEM in the same manner as described above.
2.3. Thermodynamic Calculations
Associated with each dissolution is an assemblage of solid phases, a solution phase containing dissolved calcium, phosphate, arsenate, and a pH value. Assuming equilibrium has been reached, the thermodynamic data can be calculated using established theoretical principles. In this study, the simulations were performed using PHREEQC (Version 3.1.1) together with the most complete literature database minteq.v4.dat, which bases on the ion dissociation theory. The input is free format and uses order-independent keyword data blocks that facilitate the building of models that can simulate a wide variety of aqueous-based scenarios [19]. The aqueous species considered in the calculations were Ba2+, BaOH+, BaAsO4-, BaHAsO4, BaH2AsO4+, Sr2+, SrOH+, SrAsO4-, SrHAsO4, SrH2AsO4+, H3AsO4, H2AsO4-, HAsO42-, and AsO43-. The activities of Ba2+(aq), Sr2+(aq), and HAsO42-(aq) were firstly calculated by using PHREEQC, and then the ion activity products (IAPs) for (BaxSr1-x)HAsO4·H2O, BaHAsO4·H2O, and SrHAsO4·H2O were determined according to the mass-action expressions by using Microsoft Excel.
3. Results and Discussion3.1. Solid Characterizations
The composition of the synthetic solid depends on the initial Ba : Sr : As mole ratio in the starting solution. To ensure that the (BaxSr1-x)HAsO4·H2O solid solution was formed, the precipitation was conducted by mixing barium solution, strontium solution, and arsenate solution at low rate. Results suggest that the crystal was the intended composition of (BaxSr1-x)HAsO4·H2O. The atomic (Ba + Sr)/As ratios were 1.00 which is a stoichiometric ratio of (BaxSr1-x)HAsO4·H2O. The atomic Ba/(Ba + Sr) ratios (XBa=0.00,0.21,0.42,0.61,0.82,1.00) were almost the same as those of the precursor solutions (x=0.00,0.20,0.40,0.60,0.80,1.00). No Na+ and NO3- were detected in the prepared solid (Table 1).
XRD, FT-IR, and SEM analyses were performed on the solid samples of (BaxSr1-x)HAsO4·H2O before and after dissolution (Figures 1, 2, and 3). As illustrated in the figures, the results of the analyses on materials before the dissolution were almost indistinguishable from the following reaction. No evidence of secondary mineral precipitation was observed in the dissolution experiment.
XRD patterns of the (BaxSr1 − x )HAsO4·H2O solid solution before (A) and after (B) dissolution.
FT-IR patterns of the (BaxSr1 − x )HAsO4·H2O solid solution before (A) and after (B) dissolution.
SEM images of the (BaxSr1 − x )HAsO4·H2O solid solution before (A) and after (B) dissolution.
The XRD patterns of the obtained solids indicated the formation of the (BaxSr1-x)HAsO4·H2O solid solution, which has the same type structure as BaHAsO4·H2O (ICDD PDF 00-023-0823) and SrHAsO4·H2O (ICDD PDF 01-074-1622) (Figure 1). The patterns correspond exactly with the database patterns and no impurities are observed. The solid solution is complete, with the space group Pbca (orthorhombic) being retained throughout. BaHAsO4·H2O and SrHAsO4·H2O are the two endmembers of a structural family series. When subjected to XRD, they produce the same reflections; but the reflections exist at different two-theta values; that is, the reflective planes are the same but “d” spacings are different. All the compounds have indicated the formation of a solid phase differing only in reflection location, reflection width, and absolute intensity of the diffraction patterns. The reflection peaks of BaHAsO4·H2O and SrHAsO4·H2O were slightly different from each other. The reflections of the (BaxSr1-x)HAsO4·H2O solid solutions shifted gradually to a lower-angle direction when the mole fraction XBa of the solids increased (Figure 1).
The FT-IR spectra of the (BaxSr1-x)HAsO4·H2O solid solutions are shown in Figure 2. The spectra may be divided into three sections: (a) hydroxyl-stretching region, (b) water HOH bending region, and (c) arsenate As–O stretching and OAsO bending region [9]. The free arsenate ion, AsO43-, belongs to the point group Td. The normal modes of the tetrahedral arsenate ion are ν1, symmetric As–O stretching; ν2, OAsO bending; ν3, As–O stretching; and ν4, OAsO bending. In the undistorted state, only the absorptions corresponding to ν3 and ν4 vibrations are observed. The two remaining fundamentals ν1 and ν2 become infrared active when the configuration of the AsO43- ions is reduced to some lower symmetry [9]. The degenerate modes are split by distortion of the arsenate groups through lack of symmetry in the lattice sites. As shown in Figure 2, the bands of AsO43- appeared around 815.74, 814.92, and 693.82 cm−1 (ν3) and 1747.66, 1662.27, 1620.03, and 1442.49 cm−1 (ν1) for the endmember BaHAsO4·H2O and around 854.23, 814.92, 704.10 cm−1 (ν3), and 1465.20 cm−1 (ν1) for the endmember SrHAsO4·H2O. For the (BaxSr1-x)HAsO4·H2O solid solutions, the ν1 and ν3 bands shifted slightly to a smaller wavenumber when the mole fraction XBa of the solids increased.
Well-crystallized solids formed (Figure 3). The cell parameters a and b and volume increased and c decreased in a nonlinear way with XBa. The flaky appearance increased; that is, the morphology of the (BaxSr1-x)HAsO4·H2O solid solutions changed from short prismatic or granular crystals (XBa=0) to platy or blades crystals (XBa=1).
3.2. Evolution of Aqueous Composition
The solution pH and element concentrations during the dissolution experiments at 25°C and initial pH 2 as a function of time are shown in Figure 4 for the (BaxSr1-x)HAsO4·H2O solid solution. The experimental results indicated that the dissolution could be stoichiometrical only at the very beginning of the process, and then dissolution became nonstoichiometrical and the system underwent a dissolution-recrystallization process that affects the ratio of the substituting ions in both the solid and the aqueous solution.
Change of the aqueous pH and component concentrations with dissolution time.
When dissolution progressed at the initial pH 2, the aqueous pHs increased rapidly from 2.00 to 5.51–7.71 within the first hour of the experiment. For the dissolution of the (BaxSr1-x)HAsO4·H2O solid solution with XBa = 0.00, 0.21, 0.42, and 0.61, the solution pHs increased gradually from 1 h to 120 h and after that decreased gradually until they reached the steady state after 2160 h. For the dissolution of the solids with XBa = 0.82 and 1.00, the aqueous pHs varied only slightly after 1 h and reached the steady state after 2160 h. The aqueous pHs decreased with the increasing XBa of the (BaxSr1-x)HAsO4·H2O solid solution.
For the dissolution of the solids with XBa = 0.21, 0.42, and 0.61, the aqueous Ba concentrations increased rapidly and reached the peak values within the first hour and then decreased gradually and reached the lowest peak values in 240 h. After that, they increased gradually and reached the steady state after 2160 h. For the dissolution of the solids with XBa = 0.82 and 1.00, the aqueous Ba concentrations increased rapidly within the first six hours and then increased gradually in 6–480 h. After that, they decreased gradually and reached the steady state after 2160 h. The aqueous Ba concentrations increased with the increasing XBa of the (BaxSr1-x)HAsO4·H2O solid solution.
When dissolution progressed at the initial pH 2, the aqueous Sr concentrations increased rapidly and reached the peak values within 120–480 h and then decreased gradually until they reached the steady state after 2880 h. The aqueous Sr concentrations increased with the increasing XBa of the (BaxSr1-x)HAsO4·H2O solid solution with 0<XBa<0.21 (or 1<XSr<0.79) and decreased with the increasing XBa of the solids with 0.21<XBa<1 (or 0.79<XSr<0); that is, the aqueous Sr concentrations had the highest value for the dissolution of (Ba0.21Sr0.79)HAsO4·H2O.
The aqueous As(V) concentrations increased rapidly and reached the peak values within the first hour and then decreased gradually and reached the lowest peak values in 48 h. After that, they increased gradually and reached the second peak values in 120–720 h. And then they decreased gradually and reached the steady state after 2880 h. The aqueous As(V) concentrations had the highest value for the dissolution of (Ba0.21Sr0.79)HAsO4·H2O.
3.3. Determination of Solubility and Free Energies of Formation
For the stoichiometric dissolution of the (BaxSr1-x)HAsO4·H2O solid solution according to
(1)(BaxSr1-x)HAsO4·H2O=xBa2+(aq)+(1-x)Sr2+(aq)+HAsO4-(aq)
a stoichiometric ion activity product, IAP, can be written as
(2)IAP={Ba2+}x{Sr2+}(1-x){HAsO4-}.
The activities of Ba2+(aq), Sr2+(aq), and HAsO42-(aq) were firstly calculated by using PHREEQC [19], and then the IAP values for (BaxSr1-x)HAsO4·H2O, BaHAsO4·H2O, and SrHAsO4·H2O were determined according to the mass-action expression (2). The aqueous pH, Ba, Sr, and As concentrations had reached stable values after 2880 h dissolution and the IAP values of 2880 h, 3600 h, and 4320 h were considered as Ksp of the solids (Figure 4 and Table 2).
Analytical data and solubility determination of the (BaxSr1-x)HAsO4·H2O solid solution.
Sample
Dissolution time (h)
pH
Concentration (mmol/L)
logKsp
Average logKsp
ΔGfo (kJ/mol)
Ba
Sr
As
SrHAsO4·H2O
2880
7.96
0.00
9.08
12.88
−4.62
−4.62
−1537.94
3600
7.90
0.00
9.12
12.83
−4.62
4320
7.95
0.00
9.16
12.53
−4.62
(Ba0.21Sr0.79)HAsO4·H2O
2880
7.77
0.91
13.34
15.38
−4.70
−4.69
−1538.54
3600
7.77
0.98
13.32
15.54
−4.69
4320
7.74
0.94
13.32
15.53
−4.69
(Ba0.42Sr0.58)HAsO4·H2O
2880
6.86
1.63
8.73
11.68
−5.11
−5.13
−1541.25
3600
6.85
1.56
8.82
11.07
−5.13
4320
6.83
1.54
8.69
10.81
−5.15
(Ba0.61Sr0.39)HAsO4·H2O
2880
6.18
3.92
5.18
10.15
−5.45
−5.49
−1543.50
3600
6.10
3.91
5.17
10.17
−5.51
4320
6.06
3.93
5.20
10.24
−5.53
(Ba0.82Sr0.18)HAsO4·H2O
2880
5.72
9.45
1.55
11.58
−5.57
−5.61
−1544.24
3600
5.67
9.51
1.56
11.16
−5.63
4320
5.68
9.57
1.57
11.07
−5.62
BaHAsO4·H2O
2880
5.61
10.53
0.00
11.42
−5.49
−5.52
−1543.99
3600
5.56
10.56
0.00
11.25
−5.53
4320
5.57
10.69
0.00
10.94
−5.53
For (1),
(3)ΔGro=xΔGfo[Ba2+]+(1-x)ΔGfo[Sr2+]+ΔGfo[HAsO42-]+ΔGfo[H2O]-ΔGfo[(BaxSr1-x)HAsO4·H2O].
The standard free energy of reaction (ΔGro), in kJ/mol, is related to Ksp at standard temperature (298.15 K) and pressure (0.101 MPa) by
(5)ΔGro=-5.708logKsp.
The solution chemistry representing equilibrium involving the solution phase and (BaxSr1-x)HAsO4·H2O, along with the calculated log Ksp using PHREEQC, is shown in Table 2. Based on the obtained literature data, ΔGfo[Ba2+]=-560.77 kJ/mol, ΔGfo[Sr2+]=-559.84 kJ/mol, ΔGfo[HAsO42-]=-714.59 kJ/mol, ΔGfo[H2O]=-237.141 kJ/mol, and the free energies of formation, ΔGfo[(BaxSr1-x)HAsO4·H2O], were also calculated (Table 2).
The mean Ksp values were calculated for BaHAsO4·H2O of 10−5.52 (10−5.49–10−5.53) at 25°C, for SrHAsO4·H2O of 10−4.62 at 25°C. The corresponding free energies of formation (ΔGfo) were determined to be −1543.99±0.18 kJ/mol and -1537.94±0.02 kJ/mol. The Ksp value of 10−5.52 for BaHAsO4·H2O is approximately 19.12log units higher than 10−24.64 reported by Essington [15] and approximately 4.72log units lower than 10−0.8 reported by Orellana et al. [17], but in accordance with those of Robins [10], Davis [18], and Zhu et al. [9] (10−4.70, 10−5.51, and 10−5.60, resp.). Essington [15] determined the solubility of BaHAsO4·H2O(c) based not on his own experimental measurement, but on the data of Chukhlantsev [20] for Ba3(AsO4)2(c). He took the Ba3(AsO4)2 solid used by Chukhlantsev [20] for BaHAsO4·H2O(c) and recalculated the original analytical data. Orellana et al. [17] determined the solubility product only from precipitation experiments. The aqueous solution in his experiment might not reach equilibrium and was still supersaturated with regard to the solid. Based on our experimental results and those available in the literature, the solubility product for BaHAsO4·H2O should be around 10−5.50.
BaHAsO4·H2O is less soluble than SrHAsO4·H2O. For the (BaxSr1-x)HAsO4·H2O solid solution, the solubility decreased as XBa increased when 0.00<XBa<0.82 and increased as XBa increased when 0.82<XBa<1.00. The (BaxSr1 − x )HAsO4·H2O solid solution had a minimum solubility product of 10−5.61 at XBa=0.82. This variation tendency is in accordance with the change of the unit cell parameters [4].
3.4. Saturation Index for BaHAsO4·H2O and SrHAsO4·H2O
The Ksp values of 10−5.52 for BaHAsO4·H2O and 10−4.62 for SrHAsO4·H2O were used in the calculation using the program PHREEQC in the present study. The calculated saturation indices for BaHAsO4·H2O show a trend of increasing values as the composition of the solid phases approaches that of the pure-phase endmember, BaHAsO4·H2O (Figure 5). At the beginning of the dissolution of the (BaxSr1-x)HAsO4·H2O solid solution, the BaHAsO4·H2O saturated index (SI) values increased with time until the aqueous solution was oversaturated with respect to BaHAsO4·H2O, and then the SI values decreased slowly. At the end of the dissolution experiment (2880–4320 h), the aqueous solution was saturated with respect to BaHAsO4·H2O for the (BaxSr1-x)HAsO4·H2O solid solution with SI = −0.06~0.10.
Saturation indices for BaHAsO4·H2O and SrHAsO4·H2O.
The calculated saturation indices for SrHAsO4·H2O show a distinctly different trend than those for BaHAsO4·H2O (Figure 5). For the (BaxSr1-x)HAsO4·H2O solid solution with x>0.2, the SrHAsO4·H2O saturated index (SI) values decreased as the SrHAsO4·H2O mole fraction decreased or the BaHAsO4·H2O mole fraction increased for the (BaxSr1-x)HAsO4·H2O solid solution with x>0.2. The SrHAsO4·H2O saturated index (SI) values increased with time for the (BaxSr1-x)HAsO4·H2O solid solution with x>0.2. For the (BaxSr1-x)HAsO4·H2O solid solution with x≤0.2, the SrHAsO4·H2O saturated index (SI) values increased as the SrHAsO4·H2O mole fraction increased or the BaHAsO4·H2O mole fraction decreased for the (BaxSr1-x)HAsO4·H2O solid solution with x≤0.2. At the beginning of the dissolution, the SrHAsO4·H2O saturated index (SI) values increased with time until the aqueous solution was oversaturated with respect to SrHAsO4·H2O, and then the SI values decreased slowly. At the end of the dissolution experiment (4320 h), the aqueous solution was saturated with respect to SrHAsO4·H2O for the (BaxSr1-x)HAsO4·H2O solid solution with x≤0.2, while the aqueous solution was undersaturated with respect to SrHAsO4·H2O for the (BaxSr1-x)HAsO4·H2O solid solution with x>0.2.
3.5. Lippmann Diagram3.5.1. Construction of Lippmann Diagram
Understanding solid solution aqueous solution (SSAS) processes is of fundamental importance. However, in spite of the numerous studies, the availability of thermodynamic data for SSAS systems is still scarce [6]. Lippmann extended the solubility product concept to solid solutions by developing the concept of “total solubility product ΣΠSS,” which is defined as the sum of the partial activity products contributed by the individual endmembers of the solid solution [5, 6]. At thermodynamic equilibrium, the total activity product ΣΠSS, expressed as a function of the solid composition, yields Lippmann’s “solidus” relationship. In the same way, the “solutus” relationship expresses ΣΠSS as a function of the aqueous solution composition. The graphical representation of “solidus” and “solutus” yields a phase diagram, usually known as a Lippmann diagram [5, 6]. A comprehensive methodology for describing reaction paths and equilibrium end points in solid solution aqueous solution systems had been presented and discussed in literatures [5, 6, 21–23].
In the case of the (BaxSr1-x)HAsO4·H2O solid solution, the term “total solubility product ΣΠSS” is defined as the ({Ba2+} + {Sr2+}){HAsO42-} at equilibrium and can be expressed by
(6)ΣΠSS=KBaXBaγBa+KSrXSrγSr,
where {·} designate aqueous activity. KBa and KSr, XBa and XSr, and γBa and γSr are the thermodynamic solubility products, the mole fractions (x,1-x), and the activity coefficients of the BaHAsO4·H2O and SrHAsO4·H2O components in the (BaxSr1-x)HAsO4·H2O solid solution, respectively. This relationship, called the solidus, defines all possible thermodynamic saturation states for the two-component solid solution series in terms of the solid phase composition.
The term “total solubility product ΣΠSS” can also be expressed by
(7)ΣΠSS=1XBa2+,aq/KBaγBa+XSr2+,aq/KSrγSr,
where XBa2+,aq and XSr2+,aq are the activity fractions of the aqueous Ba and Sr components. This relationship, called the solutus, defines all possible thermodynamic saturation states for the two-component solid solution series in terms of the aqueous phase composition.
For a solid solution with fixed composition XBa=1-XSr, a series of minimum stoichiometric saturation scenarios as a function of the aqueous activity fraction XBa2+,aq of the substituting ions in the aqueous solution can be described by
(8)ΣΠSS=IAP(XBa2+,aq)XBa(XSr2+,aq)XSr.
For the endmembers BaHAsO4·H2O and BaHAsO4·H2O and XBa=1 and XBa=0, the endmember saturation equations can be written as
(9)ΣΠBaHAsO4·H2O={Ba2+}{HAsO42-}(XBa2+,aq)XBa=KBa(XBa2+,aq)XBa,ΣΠSrHAsO4·H2O={Sr2+}{HAsO42-}(XSr2+,aq)XSr=KSr(XSr2+,aq)XSr.
A Lippmann phase diagram for the (BaxSr1-x)HAsO4·H2O solid solution is a plot of the solidus and solutus as logΣΠSS on the ordinate versus two superimposed aqueous and solid phase mole fraction scales on the abscissa.
3.5.2. Lippmann Diagram for the Ideal Solid Solution
In the (BaxSr1-x)HAsO4·H2O system, the endmember solubility products are −5.52 for BaHAsO4·H2O and −4.62 for SrHAsO4·H2O; that is, they differ by about 0.9 log units. The (BaxSr1 − x )HAsO4·H2O solid solutions are assumed ideal. In this case, the pure-phase solubility curves for BaHAsO4·H2O and SrHAsO4·H2O are quite distinct from the solutus, and the solution compositions along the solutus are clearly undersaturated with respect to both pure BaHAsO4·H2O and SrHAsO4·H2O solids. The dotted horizontal tie lines indicate the relation between the solid mole fractions and the aqueous activity fractions at thermodynamic equilibrium (points T2 and T1) or at primary saturation (points P2 and P1) with respect to a (Ba0.61Sr0.39)HAsO4·H2O solid. The dashed curve gives the series of possible (XBa2+,aq,ΣΠSS) aqueous compositions that satisfy the condition of stoichiometric saturation with respect to a (Ba0.61Sr0.39)HAsO4·H2O solid. Point M1 is the “minimum stoichiometric saturation” point for a (Ba0.61Sr0.39)HAsO4·H2O solid.
As with all the “minimum stoichiometric saturation” curves, the “pure endmember saturation” curves plot above the solutus for all aqueous activity fractions, concurring with the solutus only for XBa2+,aq or XSr2+,aq equal to one. This can be seen in Figure 6, in which the curves ΣΠBaHAsO4·H2O and ΣΠSrHAsO4·H2O have been represented for the same solid solution as in Figure 6. Note that the two curves intersect at the point “S.” This point represents an aqueous solution at simultaneous stoichiometric saturation with respect to both endmembers.
Lippmann diagram for the (BaxSr1− x )HAsO4·H2O ideal solid solution.
The hypothetical reaction path is also shown in Figure 6, in relation to Lippmann solutus and solidus curves for the (Ba0.61Sr0.39)HAsO4·H2O solid solution. The reaction path of a stoichiometrically dissolving solid solution moves vertically from the abscissa of a Lippmann diagram, originating at the mole fraction corresponding to the initial solid solution composition. The pathway shows initial stoichiometric dissolution up to the solutus curve, followed by nonstoichiometric dissolution along the solutus, towards the more soluble endmember.
The difference between the solubility products of the endmembers involves a strong preferential partitioning of the less soluble endmember towards the solid phase. Our dissolution data indicate a final enrichment in the BaHAsO4·H2O component in the solid phase and a persistent enrichment in the Sr2+ component in the aqueous phase. The possibility of formation of a phase close in composition to pure BaHAsO4·H2O seems unavoidable, given the lower solubility of BaHAsO4·H2O and the oversaturation with respect to BaHAsO4·H2O and the relatively high solubility of SrHAsO4·H2O and the undersaturation with respect to SrHAsO4·H2O (Figure 5). From the point of view of the crystallization behavior, all data seem to indicate that there is a preferential partitioning of barium into the solid phase [4].
3.5.3. Lippmann Diagram for the Nonideal Solid Solution
Complementary powder XRD measurements indicated that the cell parameters increased in a nonlinear way with XBa indicating the solid solution is complete but could be nonideal [4]. The solid phase activity coefficients of BaHAsO4·H2O (γBa) and SrHAsO4·H2O (γSr) as components of the solid solution can be calculated as a function of composition using the Redlich and Kister equations [5, 6], expressed in the form
(10)lnγBa=(1-x)2[a0-a1(3x-(1-x))+⋯],lnγSr=x2[a0-a1(3(1-x)-x)+⋯],
where x and 1-x are the mole fractions (XBa and XSr) of the BaHAsO4·H2O and SrHAsO4·H2O components in the (BaxSr1 − x )HAsO4·H2O solid solution; the Guggenheim coefficients a0 and a1 can be determined from an expansion of the excess Gibbs free energy of mixing GE [5, 6].
The excess Gibbs free energy of mixing GE has not been measured for the (BaxSr1 − x )HAsO4·H2O solid solution. Assuming that a stoichiometric saturation state was attained in the dissolution experiments, the Guggenheim coefficients a0 and a1 can be estimated from
(11)lnKsp=x(1-x)a0+x(1-x)(x-(1-x))a1+(1-x)ln[(1-x)KBa]+xln[xKSr].
Log IAP (logKsp) values of the samples taken after 2880 h dissolution are shown as a function of XBa in Figure 7. A plot of these log Ksp values versus solid mole fraction of BaHAsO4·H2O shows that the logKsp values are close to and slightly lower than what would be expected for an ideal solid solution. Fitting the Ksp values as a function of solid composition to (11) yields a best fit with a two-parameter Guggenheim model of a0=1.55 and a1=-4.35 (R2>0.95) (Figure 7).
Estimation of the Guggenheim coefficients for the (BaxSr1−x )HAsO4·H2O nonideal solid solution.
The diagram of Figure 8 is a typical Lippmann diagram for a solid solution with a negative enthalpy of mixing. The negative excess free energy of mixing produces a fall of the solutus curve with respect to the position of an equivalent ideal solutus. This means that the solubility of intermediate compositions is significantly smaller than that of an equivalent ideal solid solution. Indeed, for a certain range of solid compositions, the solutus values plot below the solubility product of the less soluble endmember. The system shows an “alyotropic” minimum at XBa2+,aq=0.87 where the solutus and solidus curves meet. Such a point represents a thermodynamic equilibrium state in which the mole fraction of the substituting ions in the solid phase equals its activity fraction in the aqueous solution. Obviously, at the alyotropic point, the equilibrium distribution coefficient is equal to unity [5, 6].
Lippmann diagram for the (BaxSr1−x )HAsO4·H2O nonideal solid solution.
A Lippmann diagram for the (BaxSr1 − x )HAsO4·H2O solid solution for the ideal case when a0=0.0 is shown in Figure 9. In addition to the solutus and the solidus, the diagrams contain the total solubility product curve at stoichiometric saturation for the (BaxSr1 − x )HAsO4·H2O solid solution at x = 0.21, 0.42, 0.61, and 0.82. The saturation curves for pure endmembers BaHAsO4·H2O (x=1.00) and SrHAsO4·H2O (x=0.00) have also been plotted in the chart. Also included are the data from our study, plotted as ({Ba2+} + {Sr2+}){HAsO42-} versus XBa2+,aq.
Plot of the experimental data on the Lippmann diagram.
In general, the location of data points on a Lippmann diagram will depend on the aqueous speciation, degree to which secondary phases are formed, and the relative rates of dissolution and precipitation. As (BaxSr1 − x )HAsO4·H2O dissolves in solution, aqueous Ba2+ is converted into Ba2+, BaOH+, BaAsO4-, BaHAsO4, and BaH2AsO4+, and aqueous Sr2+ is converted into Sr2+, SrOH+, SrAsO4-, SrHAsO4, and SrH2AsO4+; aqueous HAsO42- is converted primarily into AsO43-, H2AsO4-, and H3AsO4 and only a fraction remains as Ba2+, Sr2+, and HAsO42-. The smaller values of the activity fractions are a consequence of the Ba2+ and Sr2+ speciation. For the plot of the experimental data on the Lippmann diagram, the effect of the aqueous Ba2+ and Sr2+ speciation was considered by calculating the activity of Ba2+ and Sr2+ with the program PHREEQC.
The experimental data plotted on Lippmann phase diagrams show that the (Ba0.61Sr0.39)HAsO4·H2O precipitate dissolved stoichiometrically at the beginning and approached the Lippmann solutus curve and then overshot the Lippmann solutus curve, the saturation curves for pure endmembers BaHAsO4·H2O, and the stoichiometric saturation curve for x=0.61. After about 1 h dissolution, the aqueous solution was oversaturated with respect to BaHAsO4·H2O and the (BaxSr1 − x )HAsO4·H2O solid solution. And then, the XBa2+,aq decreased from 0.61 to 0.21 in 240 h with no obvious change in the logΣΠSS value, indicating the dissolution path for this precipitate may involve stoichiometric dissolution to the Lippmann solutus curve and overshot the Lippmann solutus curve followed by a possible exchange and recrystallization reaction. From 240 to 4320 h, the logΣΠSS value decreased further, but the XBa2+,aq increased from 0.21 to about 0.43, and the dissolution followed the saturation curve for pure endmembers BaHAsO4·H2O and approached the intersection with the minimum stoichiometric saturation curve for x=0.61.
4. Conclusions
With the increasing XBa, the morphology of the (BaxSr1 − x )HAsO4·H2O solid solutions changed from short prismatic or granular crystals (XBa=0) to platy or blades crystals (XBa=1). The synthetic solids used in the experiments were found to have no obvious variation after dissolution. During the dissolution of the solid solution, the aqueous pH values and component concentrations increased rapidly at the beginning and then varied slowly with time and finally exhibited stable state after 2880 h dissolution. With the increasing XBa, the aqueous pH value decreased and the aqueous Ba concentration increased. At XBa=0.21, the aqueous Sr and As concentrations had the highest values. The solubility products (Ksp) for BaHAsO4·H2O and SrHAsO4·H2O were calculated to be 10−5.52 and 10−4.62, respectively. The corresponding free energies of formation (ΔGfo) were -1543.99±0.18 kJ/mol and -1537.94±0.02 kJ/mol. BaHAsO4·H2O is less soluble than SrHAsO4·H2O. For the (BaxSr1 − x )HAsO4·H2O solid solution, the solubility decreased as XBa increased when 0.00<XBa<0.82 and increased as XBa increased when 0.82<XBa<1.00. The (BaxSr1 − x )HAsO4·H2O solid solution had a minimum solubility product of 10−5.61 at XBa=0.82. The Guggenheim coefficients were determined to be ao=1.55 and a1=-4.35 for the (BaxSr1 − x )HAsO4·H2O nonideal solid solution. The constructed Lippmann diagram was a typical Lippmann diagram for a nonideal solid solution with a negative enthalpy of mixing, which produced a fall of the solutus curve with respect to the position of an equivalent ideal solutus, which indicated that the solubility of intermediate compositions is significantly smaller than that of an equivalent ideal solid solution. The system shows an “alyotropic” minimum at XBa2+,aq=0.87. At the end of the experiment, the dissolution followed the saturation curve for the pure endmember BaHAsO4·H2O and approached the intersection with the minimum stoichiometric saturation curve on the Lippmann diagram.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors thank the financial supports from the National Natural Science Foundation of China (41263009, 40773059) and the Provincial Natural Science Foundation of Guangxi (2012GXNSFDA053022, 2011GXNSFF018003).
FulladosaE.MuratJ. C.MartínezM.VillaescusaI.Effect of pH on arsenate and arsenite toxicity to luminescent bacteria (vibrio fischeri)20044621761822-s2.0-1042301327MohanD.PittmanC. U.Jr.Arsenic removal from water/wastewater using adsorbents—a critical review20071421-21532-s2.0-3394717751210.1016/j.jhazmat.2007.01.006SmedleyP. L.KinniburghD. G.A review of the source, behaviour and distribution of arsenic in natural waters20021755175682-s2.0-003621960110.1016/S0883-2927(02)00018-5JiménezA.PrietoM.Salvadó MÁGarcía-GrandaS.Structure and crystallization behavior of the (Ba,Sr)HAsO4·H2O solid-solution in aqueous environments200489601609GlynnP. D.ReardonE. J.Solid-solution aqueous-solution equilibria: thermodynamic theory and representation199029021642012-s2.0-0025225521PrietoM.Thermodynamics of solid solution aqueous solution systems20097047852-s2.0-7514916767010.2138/rmg.2009.70.2LeeJ. S.NriaguJ. O.Stability constants for metal arsenates2007421231332-s2.0-3424720211610.1071/EN06070Planer-FriedrichB.ArmientaM. A.MerkelB. J.Origin of arsenic in the groundwater of the Rioverde basin, Mexico20014010129012982-s2.0-003481461610.1007/s002540100315ZhuY.ZhangX.XieQ.ChenY.WangD.LiangY.LuJ.Solubility and stability of barium arsenate and barium hydrogen arsenate at 25°C20051201–337442-s2.0-1634437469610.1016/j.jhazmat.2004.12.025RobinsR. G.The solubility of barium arsenate: sherritt's barium arsenate process19851624044062-s2.0-002207757310.1007/BF02679734Méndez-RodríguezL.Zenteno-SavínT.Acosta-VargasB.WurlJ.Imaz-LamadridM.Differences in arsenic, molybdenum, barium, and other physicochemical relationships in groundwater between sites with and without mining activities2013523824310.4236/ns.2013.52A035Tiruta-BarnaL.RakotoarisoaZ.MéhuJ.Assessment of the multi-scale leaching behaviour of compacted coal fly ash20061373146614782-s2.0-3374891050010.1016/j.jhazmat.2006.04.039BinasH.Boll-DornbergerK.The structure of CaHAsO4·H2O (haidingerite) and SrHAsO4·H2O196221450MartinC.DurifA.Adverbuch-PouchotM. T.Sur le groupe spatial de l'haidingerite. Données cristallographiques sur BaHAsO4·H2O197093397EssingtonM. E.Solubility of barium arsenate1988526156615702-s2.0-0024110970ItohC. T.TozawaK.Equilibria of the barium(II)-arsenic(III,V) water system at 25°C198945105OrellanaF.AhumadaE.SuarezC.CoteG.LizamaH.Thermodynamics studies of parameters involved in the formation of arsenic(V) precipitates with barium(II)200045415422DavisJ.Stability of metal-arsenic solids in drinking water systems20004131352-s2.0-003401066510.1061/(ASCE)1090-025X(2000)4:1(31)ParkhurstD. L.AppeloC. A. J.Description of input and examples for PHREEQC version 3—a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations2013chapter A431497ChukhlantsevV. G.Solubility-products of arsenates1956119751982GlynnP. D.ReardonE. J.PlummerL. N.BusenbergE.Reaction paths and equilibrium end-points in solid-solution aqueous-solution systems19905422672822-s2.0-0025244055GamsjägerH.KönigsbergerE.PreisW.Lippmann diagrams: theory and application to carbonate systems2000611913210.1023/A:1009690502299BaronD.PalmerC. D.Solid-solution aqueous-solution reactions between jarosite (KFe3(SO4)2(OH)6) and its chromate analog20026616284128532-s2.0-003668946510.1016/S0016-7037(02)00880-3