Journal Information
Vol. 8. Issue 3.
Pages 2741-2746 (May - June 2019)
Download PDF
More article options
Vol. 8. Issue 3.
Pages 2741-2746 (May - June 2019)
Original article
DOI: 10.1016/j.jmrt.2019.04.010
Open Access
Theoretical confirmation of an established experimental phase diagram of ternary system (Cu(SO4)–Zn(SO4)–H2O) at 25°C: Pitzer parameters determination
Bouchra Mekkaoui
Corresponding author

Corresponding author.
, Mohamed Jouaiti, Abderrahim Barroug, Ahmed Jouaiti, Latifa Laallam
Sustainable Development Laboratory, Faculty of Sciences and Techniques, Sultan Moulay Slimane University, Beni Mellal, Morocco
This item has received

Under a Creative Commons license
Article information
Full Text
Download PDF
Figures (4)
Show moreShow less
Tables (5)
Table 1. The measured solubility data of the ternary system (Cu(SO4)–Zn(SO4)–H2O) at 25°C.
Table 2. Weight composition coordinate of eutectic points for the studied ternary system.
Table 3. Single-salts calculated and existing Pitzer ion interaction parameters.
Table 4. Calculated ternary Pitzer ion interaction parameters.
Table 5. Experimental and theoretical points of the studied system isotherm.
Show moreShow less

The isotherm phase diagram of the ternary system (Cu(SO4)–Zn(SO4)–H2O) was firstly established at 25°C using a synthetic method of conductivity measurement. Then, we made use of volumetric dosing, in order to identify solid phases. The obtained isotherm has revealed the presence of two solid phases that having a stable equilibrium with the liquid phase, which are Cu(SO4)·5H2O, and the Zn(SO4)·7H2O. This study was confirmed numerically while going through the calculation of the binary and ternary Pitzer parameters for the studied salts through a Nelder Mead simplex algorithm. Indeed, the comparison of these with the bibliographic data has confirmed the accuracy of our experimental results.

Zinc sulfate
Copper sulfate
Gibbs energy
Full Text

The need for the optoelectronic device with a low cost and high performance had driven by, to the development of thin films of transparent conductive oxides with improved optical and electrical properties. Therefore, electrical materials are currently at the heart of scientific and technological development and their uses are expanded across the most innovative fields of innovation. Such as electrical energy, computing, renewable energy, hydrometallurgy, geochemistry and oceanography. Undoubtedly, photovoltaic cells based on metal oxides are the object of several studies [1–5]. The rational development of such oxides from the corresponding metal salts needs a good knowledge of phase equilibrium in the systems thereby formed.

The first factor involved in the development of thin films of conductive oxide is the temperature. Absolutely, matter changes its physical state, according to this one, which makes diagrams development very useful in this domain because it promotes the explanation of mechanical and thermodynamic properties changes, of metals depending on temperature. In this context comes the present work, where the main objective is to establish the isothermal phase diagram of a ternary system based on metal salts. Furthermore, a theoretical study was approached based on the Pitzer model.

The Cu(SO4)–Zn(SO4)–H2O ternary system investigated in the present work has been the objective of many researches, because of its applicability in different fields. This is particularly relevant to catalyst oxides that are imposed in many applications in several domains. Indeed, the catalysts are chemical species used in most industrial synthesis processes, by allowing them to increase the reaction kinetics while minimizing energy costs [6], as in the case of methane steam reforming [7]. The present paper is, in fact, a start of a larger work about the establishment of the phase diagram for the ZnSO4–CuSO4–NiSO4–H2O quaternary system, whose knowledge is important in the synthesis of mixed oxides based on the three mentioned metallic ions.

As we know, models have been developed allowing an adequate description of the behavior of concentrated aqueous solutions [8–10]. Pitzer is the most known one proposed for the description of the properties of dense electrolyte solutions in a wide range of temperatures [11], pressures, and compositions. Thus, this one requires the knowledge of interaction parameters between particles. As a matter of fact, Pitzer (1973) has developed his model for an electrolyte, while introducing a set of interaction parameters whether binary or ternary ones.

Hence, comes the second part of our paper in which we have treated the numerical aspect of the established isotherm using the algorithm previously elaborated. We could be able to determine Pitzer parameters for the studied system. The method is based on the comparison between the experimental phase diagram and the one obtained using Pitzer parameters, this comparison has allowed us to confirm that the experimentally developed isotherm is effective.

On the other hand, several methods have been used for the Pitzer parameters determination, except that some are considered expensive in terms of time and money, because of the enormous need for materials, energy, and equipment [12,13].

Clearly, the most known and used way is the iterative calculation based on the determination of the osmotic coefficient experimentally through the isopiestic method. This is an effective method if we could manage to undertake the experimental protocol correctly with low uncertainty.

However, it is considered expensive since the process must be repeated for a large number of modalities’ values.


The analysis of the solution in equilibrium with the different solid phases has been performed with a synthetic method based on conductivity measurements [14–17]. The water is progressively added using a micro-burette Schott Gerate T82/50 under stirring, to a saturated mixture of a given initial composition. After the stirring was stopped and the conductivity of the decanted solution is measured.

The experiment was carried out using a thermostated cell, Fig. 1, in which flows the distilled water maintained at the desired room temperature of 25°C.

Fig. 1.

Illustration of the experimental device.


We followed for several poly-phases (liquid+solid) of appropriate composition, the variation of the liquid phase conductivity according to the initial mixture dilution. The electric conductivity–composition curve of the solution has a discontinuity at each phase change. When the identification of the solid phases was too difficult only with the preceding method, the solution and the solid phases of a given mixture composition were submitted to dosages. The compositions of double saturation points are obtained by direct determination of concentration using the wet residues [18] and the “ensembles” methods [19].

The salts used in this work are the Cu(SO4)·5H2O (Blue vitriol), the Zn(SO4)·H2O. Water was bi-distilled before use. While the conductivity measurements were performed using a Knick 703 device. And the system temperature is maintained constant at 25°C by using a Julabo FP50 thermostat.

The solid phase compositions were determined by a dosage of the Cu and Zn ions in the several sulfates. It is a chemical analysis method whose objective is to confirm the solid phases identified from the conductivity curve.

The Cu ions were measured by iodometry, it is an indirect dosage using potassium iodide (KI) and sodium thiosulfate (Na2S2O3). The Zn ions are titrated by volummetry using the EDTA method and finally the sulfate ions by barium sulfate (BaSO4) precipitation.

3Pitzer model

The other object of the present work is the Pitzer parameters determination of the ternary system (Cu(SO4)·Zn(SO4)·H2O). Hence, the algorithm introduced below was the main tool for this step. In fact, the simplex method has been used to determine both binary and ternary Pitzer parameters.

The Pitzer model is one of the most known models that has known a remarkable success. It is particularly popular with geochemists, waste chemists and engineers, for the prediction of mineral solubilities and phase equilibrium. Pitzer and his colleagues have developed the ionic interaction approach to predict the mean activity coefficients of salts in aqueous solutions of electrolytes as a function of molality as well as the theoretical solubilities of isotherms [8–10].

This ion-interaction model gives an equation for the activity and osmotic coefficient of electrolyte solution depending on four binary parameters, β(0), β(1), β(2), Cϕ and two ternary ones ψ and θ. These parameters allow us the calculation of the equilibrium constant of each component in the systems at different temperature. Indeed, the present work aims to predict Pitzer parameters for the studied system while using our approach that has been developed in a previous work [20].

The method is based on several steps whose most of them move the point of the simplex where the function is maximum to the point where it is minimal. First, we initialize the programs by a set of vectors chosen such that they are not collinear. Then, vectors are ordered while respecting the ascending order of the objective function values. We compute the centroid, these steps are repeated for each iteration. Thereafter, comes the fourth iterative operations that aim to find the minima of the objective function, Reflect, Expand, Contract and shrink; the method is wide explicit in previous work.

Fig. 2 illustrates the flowchart that generally summarizes all mentioned operations while respecting the conditions and sequences of the algorithm, since the initialization of starting vectors set until we accurate the correct Pitzer parameter values.

  • Initialization: we initialize the programs by a vectors set, chosen such that they are not collinear.

    where Ei=e1,e2,…ej……en is a vector of size n1 With …j=1?if?i=j?otherwise?ej=0n1

Fig. 2.

Flowchart of employed algorithm.


With ej=1 if i=j otherwise ej=0, si is an added step to each parameter of the vector, in order to create n+1 initial vectors, the si value is chosen such that each parameter Xi is always between two reasonable values. Previously set by the author.

  • Centroid: M=1n∑i=0n−1Pi

  • Reflect: PnR=2*M+Pn

  • Expand: PnE=2*PnR−M

  • Contract: Pnc=32(M−Pn)Pncc=12(M+Pn)

  • Shrink: Vi=12(P0−Pi)i=1,….n

The objective function f is in terms of the standard deviation between the fitting of experimental points and the theoretical function.

The calculated (Yexp) is obtained by fitting the experimental points, the theoretical function (Ythe) is established after the adjustment of the theoretical points found using the Pitzer model at 25°C. Indeed, the equilibrium constant (Kb) for the binary system, is calculated for an initial guess of Pitzer parameters. Then, the ternary equilibrium constant (Kt) is calculated for the same initial guess, in order to define the theoretical pairs of solubility limits of the ternary system. The method consists essentially of equalizing the equilibrium constant of the binary (Kb) and the ternary system (Kt) corresponding to the same solid phase [21], while varying the molality of cation and anion as well as the water quantity.

where Npts: number of data points in the solubility isotherm; mc,mc: molality of cation; ma: molality of anion; nw: water quantity.

4Results and discussion4.1The experimental isotherm

The weight composition coordinates of considered mixture points were calculated using the following expressions:

The solids phases encountered in this work are:

  • -

    Cu(SO4)·5H2O (Cu5)

  • -

    Zn(SO4)·H2O (Zn1)

The experimental results of the solubility isotherm obtained at 25°C are given in Table 1. The results that have been found are exploited to find the isothermal curves of the ternary systems. In fact, we have introduced a triangular representation based on the conversion of rectangular coordinates to triangular ones (Tables 2–4).

Table 1.

The measured solubility data of the ternary system (Cu(SO4)–Zn(SO4)–H2O) at 25°C.

W(Cu(SO4)), %  W(Zn(SO4)), %  Limit domains 
35.51  Zn7+Liq 
5.03  28.33  – 
7.61  24.99  Cu5+Zn7+Liq 
8.35  21.82  Cu5+Liq 
9.13  19.27  – 
11.02  15.51  – 
12.42  11.65  – 
13.81  8.33  – 
15.41  5.41  – 
16.53  2.58  – 
18.48  – 
Table 2.

Weight composition coordinate of eutectic points for the studied ternary system.

W(Cu(SO4)), %  W(Zn(SO4)), %  Domain limits 
7.61  24.99  Cu5+Zn7+Liq 
Table 3.

Single-salts calculated and existing Pitzer ion interaction parameters.

Salts    β0  β1  Cφ  β2  Reference 
Cu(SO4Calculated  0.195  2.4  −0.002  −50.12   
  Literature  0.2358  2.485  −0.012  −47.35  [24] 
    0.2340  2.527  0.0044  −48.33  [25] 
    0.21757  2.62597  0.013756  −56.2413  [26] 
Zn(SO4Calculated  0.14  2.88  0.05  −32.70   
  Literature  0.1949  2.883  0.02896  −32.81  [27] 
Table 4.

Calculated ternary Pitzer ion interaction parameters.

Ternary parameters  θCu,Zn  ψCu,Zn,SO4 
Calculated  0.314  −0.025 

The conversion of the orthogonal reference to a triangular one on a 45° angle is as follows:

Fig. 3 illustrates the isotherm of Cu(SO4)–Zn(SO4)–H2O system, according to the triangular representation. The existing domains in the curve are explicit as follows:

  • -

    (I): the liquidfield.

  • -

    (II): the saturated solution is in equilibrium with a stable Zn(SO4)·7H2O solid phase.

  • -

    (III): the solid phases in equilibrium with the saturated solution is Cu(SO4)·5H2O.

  • -

    (IV): the solid phases in equilibrium with their saturated solution are Zn(SO4)·7H2O and Cu(SO4)·5H2O.

Fig. 3.

The 25°C solubility isotherm of (Cu(SO4)–Zn(SO4)–H2O) ternary system.


Otherwise, only one invariant isotherm and isobar is found at 25°C with the following weight composition coordinates:

4.2Theoretical isotherm

In this part, our approach is to determine the Pitzer parameters for the studied solid phases using the algorithm previously confirmed while exploiting the isothermal curve of the stake ternary system, a comparison was made between the results of our program and existing data in literature as will be shown in tables.

The ternary Pitzer parameters of the studied system are not available in literature. Therefore, we have been able to find their values thanks to our used approach.

The main object of this part is to confirm our experimental work. As we have previously seen in this paper, the algorithm is based essentially on the isothermal curve, especially the experimental points of solubility limits. We made use of these points to find the Pitzer parameters numerically. Then, the comparison was made to discover how much we are near to the literature data.

Fig. 4 shows the solubility limits calculated from the developed method, it is a triangular representation of the mass fractions. As it is clearly shown in the curve, the results that we found from the Simplex method based on the Pitzer calculation model are identical to those we found experimentally, which confirms the certainty of our experimental curve.

Fig. 4.

The theoretical isotherm of the ternary system CuSO4–ZnSO4–H2O.


As it will be shown in Table 5, a comparison was made between experimental data and theoretical points established from the used method, it has proved the efficiency of the approach based on Nelder Mead Simplex.

Table 5.

Experimental and theoretical points of the studied system isotherm.

Experimental pointsTheoretical points
W(Cu(SO4)), %  W(Zn(SO4)), %  W(Cu(SO4)), %  W(Zn(SO4)), % 
35.51  35.5 
2.51  31.87  2.97  31.62 
5.03  28.33  5.5  28.3 
7.61  24.99  7.7  25 
8.54  22.32  8.3  22.5 
9.33  19.70  9.2  19.5 
10.87  15.30  11  14.2 
12.33  11.56  12.4  11.6 
13.81  8.33  13.8 
15.41  5.41  15.4 
16.82  2.63  16.5  3.7 
18.48  18.28 

The most interesting issue about this method is that it can also be used inversely. Thus, this work has been dedicated to determine the binary and ternary Pitzer parameters of the metallic salt system, while basing it on experimental results that has been established before. In a reversible way, the experimental points that define the isotherm can be elaborated from known Pitzer parameters while making these last the fixed starting values of the simplex method, that has been used in a previous work [22] and with the same protocol we will search for the function that represents the isotherm, by minimizing the Gibbs energy [23]. The discovery of this option about our program represents a major gain in the field of experimental and theoretical chemistry, because being able to simulate the experiment without having to invest in terms of time, energy as well as chemicals, is already a remarkable progress in the present discipline, which comes to reduce and optimize significant expenses.


With this investment, we have proved that the used algorithm is efficient for the studied system. In our case, we have opted for two different aspects: In the first place the isotherm of the ternary system (Cu(SO4)–Zn(SO4)–H2O) was determined using the experimental protocol seen above, then a numerical study was conducted to confirm our experimental result while comparing the obtained PITZER parameters with literature.

In fact, this identification allowed us to conclude that the obtained isotherm is 100 percent correct. Furthermore, we have justified that the algorithm previously established can be useful, to elaborate a phase diagram from Pitzer parameters, which means less investment in terms of time, money and energy.

Conflicts of interest

The authors declare no conflicts of interest.


The authors would like to thank all people that have contributed to achieving this work.

F. Le Peltier, P. Chaumette, J. Saussey, M.M. Bettahar, J.C. Lavalley.
In situ FT-IR and kinetic study of methanol synthesis from CO2/H2 over ZnAl2O4 and Cu–ZnAl2O4 catalysts.
J Mol Catal A Chem, 132 (1998), pp. 91-100
A. Monzón, E. Romeo, C. Royo, R. Trujillano, F.M. Labajos, V. Rives.
Use of hydrotalcites as catalytic precursors of multimetallic mixed oxides. Application in the hydrogenation of acetylene.
Appl Catal A: Gen, 185 (1999), pp. 53-63
Hala G. El-Shobaky, M. Mokhtar, Gamil A. El-Shobaky.
Physicochemical surface and catalytic properties of CuO–ZnO/Al2O3 system.
Appl Catal A: Gen, 180 (1999), pp. 335-344
Velinov, Nikolay, et al.
Mechanochemical synthesis and characterization of nanocrystalline copper–cobalt ferrites.
Acta Metal Sin (Engl Lett), 28 (2015), pp. 367-372
G.A. Gomes, G.L. da Costa, A.B.-H. da Silva Figueiredo.
Synthesis of ferrite nanoparticles Cu1−xAg×Fe2O4 and evaluation of potential antibacterial activity.
J Mater Res Technol, 7 (2018), pp. 381-386
Q. Zhang, S. Wen, D. Wu, Q. Feng, S. Li.
Dissolution kinetics of hemimorphite in trichloroacetic acid solutions.
J Mater Res Technol, (2019),
D. Hammoud.
Synthèses et caractérisations d’oxydes mixtes à base de cuivre, zinc et aluminium issus de précurseurs de type hydrotalcite: Application pour la réaction de vaporeformage du biométhanol.
Université du Littoral Côte d’Opale, (2015),
Kenneth S. Pitzer.
Thermodynamics of electrolytes. I. Theoretical basis and general equations.
J Phys Chem, 77 (1973), pp. 268-277
Kenneth S. Pitzer, G. Mayorga.
Thermodynamics of electrolytes. II. Activity and osmotic coefficients for strong electrolytes with one or both ions univalent.
J Phys Chem, 77 (1973), pp. 2300-2308
T. Vielma, J. Salminen, U. Lassi.
Thermodynamics of the ZnSO4–H2SO4–H2O system.
Calphad, 60 (2018), pp. 126-133
R.G. da Silva, M. Seckler, S.D.F. Rocha, D. Saturnino, É.D. de Oliveira.
Thermodynamic modeling of phases equilibrium in aqueous systems to recover potassium chloride from natural brines.
J Mater Res Technol, 6 (2017), pp. 57-64
H. Park, P. Englezos.
Osmotic coefficient data for NaOH–NaCl–NaAl(OH)–H2O system measured by an isopiestic method and modeled using Pitzer's model at 298.15K.
Fluid Phase Equilibria, 155 (1999), pp. 251-259
M.El Guendouzi, A. Dinane, A. Mounir.
Hygrometric determination of water activities and osmotic and activity coefficients of NH4Cl–LiCl–H2O at 25°C.
J Sol Chem, 31 (2002), pp. 119-129
R. Tenu, J.J. Counioux.
Conductimetrically determined 45-degrees-C isothermal section of the system H2O–CaCl2–SrCl2–BaCl2.
Bull Soc Chim France I – Physic Syst Liqu Electrochim Catal Genie Chim, 5–6 (1979), pp. I155-I162
L. Laallam, M. Kaddami, S. El Allali, L. Missali, R. Tenu, J. Berthet, et al.
Diagramme de phases du système ternaire H2O–Al(NO3)3–Zn(NO3)2 entre-25 et 40°C.
Annal Chim, 29 (2004),
L. Laallam, A. Jouaiti, S. El Allali, M. Ferhat.
Phase diagram of the ternary Cu (NO3)2–Zn(NO3)2–H2O system.
Mater Lett, 57 (2003), pp. 2707-2709
L. Lâallam, A. Jouaiti, A. Barroug, S. El Allali, S. Berthet, R. Tenu, et al.
Phase diagram of the quaternary Al(NO3)3–Cu(NO3)2–Zn(NO3)2–H2O system.
J Phase Equil Diffusion, 32 (2011), pp. 543-550
A. Chrétien.
Etude du système quaternaire eau, nitrate de sodium, chlorure de sodium, sulfate de sodium, par André Chrétien.
Masson, (1929),
F.A.H. Schreinemakers.
Die Alkalichromate.
Zeitsch Phys Chem, 55 (1906), pp. 71-98
M. Bouchra, J. Mohamed, L. Latifa, L. Ahmed.
Determination of Pitzer parameters from experimental isotherm of a ternary system using Nelder mead method.
Der Pharma Chem, 10 (2018),
A. Jouaiti, B. Mekkaoui, et al.
Estimation of Pitzer parameters for aqueous salts equilibrium. Ternary system with a common ion. Al(NO3)3–Zn(NO3)2–H2O.
Int J Appl Eng Res (IJAER), 10 (2015), pp. 41326-41328
J.C. Lagarias, J.A. Reeds, M.H. Wright, P.E. Wright.
Convergence properties of the Nelder–Mead simplex method in low dimensions.
SIAM J Optim, 9.1 (1998), pp. 112-147
K.K. Alaneme, M.O. Bodunrin, S.R. Oke.
Processing, alloy composition and phase transition effect on the mechanical and corrosion properties of high entropy alloys: a review.
J Mater Res Technol, 5 (2016), pp. 384-393
Kenneth S. Pitzer, G. Mayorga.
Thermodynamics of electrolytes: III. Activity and osmotic coefficients for 2–2 electrolytes.
Mol Struct Stat Thermodyn: Selected Papers of Kenneth S Pitzer, (1993), pp. 405-412
Colin J. Downes, Kenneth S. Pitzer.
Thermodynamics of electrolytes. Binary mixtures formed from aqueous NaCl, Na2SO4, CuCl2, and CuSO4, at 25°C.
J Sol Chem, 5 (1976), pp. 389-398
D.G. Miller, J.A. Rard, L.B. Eppstein, R.A. Robinson.
Mutual diffusion coefficients, electrical conductances, osmotic coefficients, and ionic transport coefficientsl ij for aqueous CuSO4 at 25°C.
J Sol Chem, 9 (1980), pp. 467-496
Kenneth S. Pitzer, G. Mayorga.
Thermodynamics of electrolytes: III. Activity and osmotic coefficients for 2-2 electrolytes.
Mol Struct Stat Thermodyn: Selected Papers of Kenneth S Pitzer, (1993), pp. 405-412
Copyright © 2019. Brazilian Metallurgical, Materials and Mining Association
Journal of Materials Research and Technology

Subscribe to our newsletter

Article options
Cookies policy
To improve our services and products, we use cookies (own or third parties authorized) to show advertising related to client preferences through the analyses of navigation customer behavior. Continuing navigation will be considered as acceptance of this use. You can change the settings or obtain more information by clicking here.