Thermodynamic evaluation and optimization of the (MgCl2 + CaCl2 + MnCl2 + FeCl2 + CoCl2 + NiCl2) system

A complete critical evaluation of all available phase diagram and thermodynamic data has been performed for all condensed phases of the (MgCl₂ + CaCl₂ + MnCl₂ + FeCl₂ + CoCl₂ + NiCl₂) system, and optimized model parameters have been found. The model parameters obtained for the binary subsystems can be used to predict thermodynamic properties and phase equilibria for the multicomponent system. The Modified Quasichemical Model was used for the molten salt phase, and the (MgCl₂ + MnCl₂ + FeCl₂ + CoCl₂ + NiCl₂) solid solution was modeled using a cationic substitutional model with an ideal entropy and an excess Gibbs free energy expressed as a polynomial in the component mole fractions. This is the first of two articles on the optimization of the (NaCl + KCl + MgCl₂ + CaCl₂ + MnCl₂ + FeCl₂ + CoCl₂ + NiCl₂) system.     

Thermodynamic evaluation and optimization of the (NaCl + KCl + MgCl2 + CaCl2 + ZnCl2) system

A complete critical evaluation of all available phase diagram and thermodynamic data has been performed for all condensed phases and relevant gaseous species of the (NaCl + KCl + MgCl₂ + CaCl₂ + ZnCl₂) system, and optimized model parameters have been found. The (NaCl + KCl + MgCl₂ + CaCl₂) subsystem has been critically evaluated in a previous article. The model parameters obtained for the binary and ternary subsystems can be used to predict thermodynamic properties and phase equilibria for the multicomponent system. The Modified Quasichemical Model for short-range ordering was used for the molten salt phase.     

Thermodynamic evaluation and optimization of the (LiF + NaF + KF + MgF2 + CaF2 + SrF2) system

A complete critical evaluation of all available phase diagram and thermodynamic data has been performed for all condensed phases of the (LiF + NaF + KF + MgF₂ + CaF₂ + SrF₂) system, and optimized model parameters have been found. The (LiF + NaF + KF + MgF₂ + CaF₂) subsystem has been critically evaluated in a previous article. The model parameters obtained for the binary and ternary subsystems can be used to predict thermodynamic properties and phase equilibria for the multicomponent system. The Modified Quasichemical Model for short-range ordering was used for the molten salt phase, and the low-temperature and high-temperature (CaF₂ + SrF₂) solid solutions were modelled using a cationic substitutional model with an ideal entropy and an excess Gibbs free energy expressed as a polynomial in the component mole fractions. Finally, the (Li, Na, K)(Mg, Ca, Sr)F₃ perovskite phase was modelled using the Compound Energy Formalism.     

Pitzer ion-interaction parameters for Fe(II) and Fe(III) in the quinary {Na + K + Mg + Cl + SO4 + H2O} system at T=298.15 K

This paper describes a chemical model that calculates (solid + liquid) equilibria in the {m₁FeCl₂ + m₂FeCl₃}(aq), {m₁FeSO₄ + m₂Fe₂(SO₄)₃}(aq), {m₁NaCl + m₂FeCl₃}(aq), {m₁Na₂SO₄ + m₂FeSO₄}(aq), {m₁NaCl + m₂FeCl₂}(aq), {m₁KCl + m₂FeCl₃}(aq), {m₁K₂SO₄ + m₂Fe₂(SO₄)₃}(aq), {m₁KCl + m₂FeCl₂}(aq), {m₁K₂SO₄ + m₂FeSO₄}(aq), and {m₁MgCl₂ + m₂FeCl₂}(aq) systems, where m denotes molality at T=298.15 K. The Pitzer ion-interaction model has been used for thermodynamic analysis of the experimental activity data in binary FeCl₂(aq) and FeCl₃(aq) solutions, and ternary solubility data, presented in the literature. The thermodynamic functions needed (binary and ternary parameters of ionic interaction, thermodynamic solubility products) have been calculated and the theoretical solubility isotherms have been plotted. The mixed solution model parameters {θ(MN) and ψ(MNX)} have been chosen on the basis of the compositions of saturated ternary solutions and data on the pure water solubility of the K₂SO₄ · FeSO₄ · 6H₂O double salt. The standard chemical potentials of four ferrous {FeCl₂ · 4H₂O, Na₂SO₄ · FeSO₄ · 4H₂O, K₂SO₄ · FeSO₄ · 6H₂O, and MgCl₂ · FeCl₂ · 8H₂O} and three ferric {FeCl₃ · 6H₂O, 2KCl · FeCl₃ · H₂O, and 2K₂SO₄ · Fe₂(SO₄)₃ · 14H₂O} solid phases have been determined. Comparison of solubility predictions with experimental data not used in model parameterization is given. The component activities of the saturated {m₁MgSO₄ + m₂FeSO₄}(aq) and in the mixed crystalline phase were determined and the change of the molar Gibbs free energy of mixing Δ_mixG^∘_m(s) of crystals was determined as a function of the solid phase composition. It is established that at T=298.15 K the mixed (Mg,Fe)SO₄ · 7H₂O and (Fe,Mg)SO₄ · 7H₂O crystals show small positive deviations from the ideal mixed crystals. Limitations of the {Fe(II) + Fe(III)} model due to data insufficiencies are discussed.     

Thermodynamic study of the system (LiCl + CaCl2 + H2O)

Solubility isotherms of the ternary system (LiCl + CaCl₂ + H₂O) were elaborately determined at T = (283.15 and 323.15) K. Several thermodynamic models were applied to represent the thermodynamic properties of this system. By comparing the predicted and experimental water activities in the ternary system, an empirical modified BET model was selected to represent the thermodynamic properties of this system. The solubility data determined in this work at T = (283.15 and 323.15) K, as well as those from the literature at other temperatures, were used for the model parameterization. A complete phase diagram of the ternary system was predicted over the temperature range from (273.15 to 323.15) K. Subsequently, the Gibbs free energy of formation of the solid phases CaCl₂ · 4 H₂O_(s), CaCl₂ · 2 H₂O_(s), LiCl · 2H₂O_(s), and LiCl · CaCl₂ · 5H₂O_(s) was estimated and compared with the literature data.     

Thermodynamic evaluation and optimization of the (NaF + AlF3 + CaF2 + BeF2 + Al2O3 + BeO) system

A complete critical evaluation of all available phase diagram and thermodynamic data has been performed for all condensed phases and relevant gaseous species of the (NaF + AlF₃ + CaF₂ + BeF₂ + Al₂O₃ + BeO) system, and optimized model parameters have been found. The (NaF + AlF₃ + CaF₂ + Al₂O₃) subsystem, which is the base electrolyte used for the electro-reduction of alumina in Hall–Héroult cells, has been critically evaluated in a previous article. The Modified Quasichemical Model in the Quadruplet Approximation for short-range ordering was used for the molten salt phase. The thermodynamic database developed is a first step towards a quantitative study of the beryllium mass balance in an electrolysis cell. In particular, the predominant Be-containing species in the gas phase evolved at the anode were identified; and, for a given beryllium content of the alumina, the beryllium content of the electrolytic bath at steady state was assessed under several approximations.     

Isopiestic measurement and solubility evaluation of the ternary system (CaCl2 + SrCl2 + H2O) at T = 298.15 K

Water activities in the ternary system (CaCl₂ + SrCl₂ + H₂O) and its sub-binary system (CaCl₂ + H₂O) at T = 298.15 K have been elaborately measured by an isopiestic method. The data of the measured water activity were used to justify the reliability of solubility isotherms reported in the literature by correlating them with a thermodynamic Pitzer–Simonson–Clegg (PSC) model. The model parameters for representing the thermodynamic properties of the (CaCl₂ + H₂O) system from (0 to 11) mol ⋅ kg⁻¹ at T = 298.15 K were determined, and the experimental water activity data in the ternary system were compared with those predicted by the parameters determined in the binary systems. Their agreement indicates that the PSC model parameters can reliably represent the properties of the ternary system. Under the assumption that the equilibrium solid phases are the pure solid phases (SrCl₂ ⋅ 6H₂O and CaCl₂ ⋅ 6H₂O)_(s) or the ideal solid solution consisting of CaCl₂ ⋅ 6H₂O_(s) and SrCl₂ ⋅ 6H₂O_(s), the solubility isotherms were predicted and compared with experimental data from the literature. It was found that the predicted solubility isotherm agrees with experimental data over the entire concentration range at T = 298.15 K under the second assumption described above; however, it does not under the first assumption. The modeling results reveal that the solid phase in equilibrium with the aqueous solution in the ternary system is an ideal solid solution consisting of SrCl₂ ⋅ 6H₂O_(s) and CaCl₂ ⋅ 6H₂O_(s). Based on the theoretical calculation, the possibility of the co-saturated points between SrCl₂ ⋅ 6H₂O_(s) and the solid solution (CaCl₂ ⋅ 6H₂O + SrCl₂ ⋅ 6H₂O)_(s) and between CaCl₂ ⋅ 6H₂O_(s) and the solid solution (CaCl₂ ⋅ 6H₂O + SrCl₂ ⋅ 6H₂O)_(s), which were reported by experimental researchers, has been discussed, and the Lippann diagram of this system has been presented.     

Thermodynamic investigation of the (LiF + NaF + CaF2 + LaF3) system

In this work a thermodynamic assessment of the (LiF + NaF + CaF₂ + LaF₃) system is reported. For the thermodynamic modeling of the liquid phase, the classical polynomial model, and the modified quasi-chemical model were used in parallel and compared. The extrapolation to higher order systems was done according to the Toop mathematical formalism. Furthermore, differential-scanning calorimetry data of the ternary (LiF + CaF₂ + LaF₃), (NaF + CaF₂ + LaF₃), and the quaternary (LiF + NaF + CaF₂ + LaF₃) mixtures are presented. Good agreement between the experimental data and the thermodynamic assessment was obtained.     

Thermodynamic evaluation and optimization of the (NaCl + Na2SO4 + Na2CO3 + KCl + K2SO4 + K2CO3) system

A complete, critical evaluation of all phase diagrams and thermodynamic data was performed for all condensed phases of the (NaCl + Na₂SO₄ + Na₂CO₃ + KCl + K₂SO₄ + K₂CO₃) system, and optimized parameters for the thermodynamic solution models were obtained. The Modified Quasichemical Model in the Quadruplet Approximation was used for modelling the liquid phase. The model evaluates first- and second-nearest-neighbour short-range order, where the cations (Na⁺ and K⁺) were assumed to mix on a cationic sublattice, while anions $({\text{CO}}_3^{2-}\text{,}{\text{SO}}_4^{2-}\text{,}\text{and}{\text{Cl}}^{-})$ were assumed to mix on an anionic sublattice. The thermodynamic properties of the solid solutions of (Na,K)₂(SO₄,CO₃) were modelled using the Compound Energy Formalism, and (Na,K)Cl was modelled using a substitutional model in previous studies. Phase transitions in the common-cation ternary systems (NaCl + Na₂SO₄ + Na₂CO₃) and (KCl + K₂SO₄ + K₂CO₃) were studied experimentally using d.s.c./t.g.a. The experimental results were used as input for evaluating the phase equilibrium in the common-cation ternary systems. The models can be used to predict the thermodynamic properties and phase equilibria in multicomponent heterogeneous systems. The experimental data from the literature are reproduced within experimental error limits.     

Thermodynamic evaluation and optimization of the (NaNO3 + KNO3 + Na2SO4 + K2SO4) system

A complete critical evaluation of all available phase diagram and thermodynamic data has been performed for all condensed phases of the (NaNO₃ + KNO₃ + Na₂SO₄ + K₂SO₄) ternary reciprocal system, and optimised model parameters have been found. The model parameters obtained for the four binary common-ion subsystems (i.e. (NaNO₃ + Na₂SO₄), (KNO₃ + K₂SO₄), (NaNO₃ + KNO₃) and (Na₂SO₄ + K₂SO₄)) are used to predict thermodynamic properties and phase equilibria for the entire system. The Modified Quasichemical Model in the Quadruplet Approximation for short-range ordering was used for the molten salt phase, and the Compound Energy Formalism was used for the various solid solutions.     

Thermodynamic study of the mixed system (CsCl + CaCl2 + H2O) by EMF Measurements at T = 298.15 K

This work reports the results of a thermodynamic investigation of the ternary mixed-electrolyte system (CsCl + CaCl₂ + H₂O). The activity coefficients of this mixed aqueous electrolyte system have been studied with the electromotive force measurement (EMF) of the cell: Cs ion-selective electrode (ISE)|CsCl(m_A), CaCl₂(m_B), H₂O|Ag/AgCl at T = 298.15 K and over total ionic strengths from (0.01 to 1.50) mol · kg^?1 for different ionic strength fractions y_B of CaCl₂ with y_B = (0, 0.2, 0.4, 0.6, and 0.8). The cesium ion-selective electrode (Cs-ISE) and the Ag/AgCl electrode used in this work were made in our laboratory and had a good Nernst response. The experimental results obey the Harned rule, and the Pitzer model can be used to describe this ternary system satisfactorily. The osmotic coefficients, excess Gibbs free energies and activities of water of the mixtures were also calculated.     

Phase equilibria for H2 + CO2 + H2O system containing gas hydrates

Isothermal phase equilibrium (pressure–composition in the gas phase) for the ternary system of H₂ + CO₂ + H₂O has been investigated in the presence of gas hydrate phase. Three-phase equilibrium pressure increases with the H₂ composition of gas phase. The Raman spectra suggest that H₂ is not enclathrated in the hydrate-cages and behaves only like the diluent gas toward the formation of CO₂ hydrate. This fact is also supported by the thermodynamic analysis using Soave–Redlich–Kwong equation of state.     

Salt effect on (liquid + liquid) equilibrium of (water + tert-butanol + 1-butanol) system: Experimental data and correlation

(Liquid + liquid) equilibrium data for the quaternary systems (water + tert-butanol + 1-butanol + KBr) and (water + tert-butanol + 1-butanol + MgCl₂) were experimentally determined at T = 293.15 K and T = 313.15 K. For mixtures with KBr, the overall salt concentrations were 5 and 10 mass percent; for mixtures with MgCl₂, the overall salt concentrations were 2 and 5 mass percent. The experimental results were used to estimate molecular interaction parameters for the NRTL activity coefficient model, using the Simplex minimization method and a concentration-based objective function. The correlation results are extremely satisfactory, with deviations in phase compositions below 1.7%.     

Thermodynamic optimization of the (Na2O + SiO2 + NaF + SiF4) reciprocal system using the Modified Quasichemical Model in the Quadruplet Approximation

All available thermodynamic and phase diagram data for the condensed phases of the ternary reciprocal system (NaF + SiF₄ + Na₂O + SiO₂) have been critically assessed. Model parameters for the unary (SiF₄), the binary systems and the ternary reciprocal system have been found, which permit to reproduce the most reliable experimental data. The Modified Quasichemical Model in the Quadruplet Approximation was used for the oxyfluoride liquid solution, which exhibits strong first-nearest-neighbor and second-nearest-neighbor short-range ordering. This thermodynamic model takes into account both types of short-range ordering as well as the coupling between them. Model parameters have been estimated for the hypothetical high-temperature liquid SiF₄.     

Phase equilibrium data and thermodynamic modeling of the system (CO2 + biodiesel + methanol) at high pressures

The main objective of this work was to investigate the high pressure phase behavior of the binary systems {CO₂(1) + methanol(2)} and {CO₂(1) + soybean methyl esters (biodiesel)(2)} and the ternary system {CO₂(1) + biodiesel(2) + methanol(3)} were determined. Biodiesel was produced from soybean oil, purified, characterized and used in this work. The static synthetic method, using a variable-volume view cell, was employed to obtain the experimental data in the temperature range of (303.15 to 343.15) K and pressures up to 21 MPa. The mole fractions of carbon dioxide were varied according to the systems as follows: (0.2383 to 0.8666) for the binary system {CO₂(1) + methanol(2)}; (0.4201 to 0.9931) for the binary system {CO₂(1) + biodiesel(2)}; (0.4864 to 0.9767) for the ternary system {CO₂(1) + biodiesel(2) + methanol(3)} with a biodiesel to methanol molar ratio of (1:3); and (0.3732 to 0.9630) for the system {CO₂ + biodiesel + methanol} with a biodiesel to methanol molar ratio of (8:1). For these systems, (vapor + liquid), (liquid + liquid), (vapor + liquid + liquid) transitions were observed. The phase equilibrium data obtained for the systems were modeled using the Peng–Robinson equation of state with the classical van der Waals (PR-vdW2) and Wong-Sandler (PR–WS) mixing rules. Both thermodynamic models were able to satisfactorily correlate the phase behavior of the systems investigated and the PR–WS presented the best performance.     

Phase behavior of (CO2 + methanol + lauric acid) system

In this study the phase equilibrium behaviors of the binary system (CO₂ + lauric acid) and the ternary system (CO₂ + methanol + lauric acid) were determined. The static synthetic method, using a variable-volume view cell, was employed to obtain the experimental data in the temperature range of (293 to 343) K and pressures up to 24 MPa. The mole fractions of carbon dioxide were varied according to the systems as follows: (0.7524 to 0.9955) for the binary system (CO₂ + lauric acid); (0.4616 to 0.9895) for the ternary system (CO₂ + methanol + lauric acid) with a methanol to lauric acid molar ratio of (2:1); and (0.3414 to 0.9182) for the system (CO₂ + methanol + lauric acid) with a methanol to lauric acid molar ratio of (6:1). For these systems (vapor + liquid), (liquid + liquid), (vapor + liquid + liquid), and (solid + fluid) transitions were observed. The phase equilibrium data obtained for the systems were modeled using the Peng–Robinson equation of state with the classical van der Waals mixing rule with a satisfactory correlation between experimental and calculated values.     

Thermodynamic evaluation and optimization of the (NaCl + KCl + AlCl3) system

All available phase equilibrium and thermodynamic data for the (NaCl + KCl + AlCl₃) system were collected and critically evaluated. An optimization was performed to obtain the parameters of one set of model equations for each phase (solids, liquid, gas) in order to best reproduce all the data simultaneously. In this way the data are rendered self-consistent, discrepancies among the data are identified, and extrapolations and interpolations can be performed. For the molten phase the Modified Quasichemical Model for short-range ordering was used, with monomeric Al³⁺ ions (corresponding to AlCl₄⁻ complexes in earlier models) predominating in alkali-rich melts, and dimeric aluminum species (corresponding to Al₂Cl₇⁻ complexes in previous models) predominating in AlCl₃-rich melts. No ternary model parameters were required for the liquid phase; the binary parameters suffice. The models can be used with Gibbs free energy minimization software to calculate phase diagram sections, vapor pressures, and all thermodynamic properties at all compositions and over extended ranges of temperature and pressure.     

Thermodynamic evaluation and optimization of the (Na2CO3 + Na2SO4 + Na2S + K2CO3 + K2SO4 + K2S) system

A critical evaluation of all phase diagram and thermodynamic data were performed for the solid and liquid phases of the (Na₂CO₃ + Na₂SO₄ + Na₂S + K₂CO₃ + K₂SO₄ + K₂S) system and optimized model parameters were obtained. The Modified Quasichemical Model in the Quadruplet Approximation was used for modelling the liquid phase. The model evaluates first- and second-nearest-neighbour short-range ordering, where the cations (Na⁺ and K⁺) are assumed to mix on a cationic sublattice, while anions $({\text{CO}}_3^{2-}\text{,}{\text{SO}}_4^{2-}\text{,}\text{and}{\text{S}}^{2-})$ are assumed to mix on an anionic sublattice. The Compound Energy Formalism was used for modelling the solid solutions of (Na, K)₂(CO₃, SO₄, S). The models can be used to predict the thermodynamic properties and phase equilibria in multicomponent heterogeneous systems. The experimental data from the literature were reproduced within experimental error limits.     

(Liquid + solid) metastable equilibria in quinary system Li2CO3 + Na2CO3 + K2CO3 + Li2B4O7 + Na2B4O7 + K2B4O7 + H2O at T=288 K for Zhabuye salt lake

An experimental study on metastable equilibria at T=288 K in the quinary system Li₂CO₃ + Na₂CO₃ + K₂CO₃ + Li₂B₄O₇ + Na₂B₄O₇ + K₂B₄O₇ + H₂O was done by isothermal evaporation method. Metastable equilibrium solubilities and densities of the solution were determined experimentally. According to the experimental data, the metastable equilibrium phase diagram under the condition saturated with Li₂CO₃ was plotted, in which there are four invariant points; nine univariant curves; six fields of crystallization: K₂CO₃ · 3/2H₂O, K₂B₄O₇ · 5H₂O, Li₂B₂O₄ · 16H₂O, Na₂B₂O₄ · 8H₂O, Na₂CO₃ · 10H₂O, NaKCO₃ · 6H₂O. Some differences were found between the stable phase diagram at T=298 K and the metastable one at T=288 K.