Abstract

A thermodynamic method of the global Gibbs energy variation calculation for describing heterogeneous equilibria of transformations of calcium, strontium, and barium polyvanadates, that occur in systems MeO-V2O5-H2O, where Me is the alkaline earth element, has been developed and used. Its quintessence consists in the thermodynamic analysis of the real conditions of various processes on the basis of their total thermodynamic characteristics. On the basis of the selected thermodynamic data for involved species, the thermodynamic stability areas of solid polyvanadates towards the solution pH and vanadium and alkaline earth metal ion concentrations in heterogeneous mixtures have been established, taking into account the complex formation reactions in multicomponent heterogeneous systems. The existing experimental data confirm the results on the thermodynamic stability of polyvanadates obtained in this paper.

1. Introduction

Polyvanadates obtained by solid-phase synthesis and subsequent study of their chemical properties in water-salt systems have a huge practical significance, because regardless of the type of source of raw materials, their processing includes stages of the synthesis of final products with desired properties (high or low solubility, etc.), with the transfer of vanadium in the solution and its subsequent precipitation. The abovementioned processes are associated with extremely complex chemical equilibria, the study of which would allow the development of technological basis of a number of combined productions: obtaining technical and pure vanadium pentoxide, vanadium catalysts, dyes, etc. [1, 2]. Up to 85–87% of the total amount of polyvanadates is used in ferrous metallurgy as highly efficient and habitually an indispensable alloying additive in the production of diverse steels. About 10–12% of polyvanadates is used in nonferrous metallurgy, mainly in the form of aluminum-vanadium alloys for alloying structural materials based on titanium, applied in space technology. Inorganic materials containing polyvanadates are used for production of optical quantum generators, ionic conductors, ferroelectrics, dielectrics, ceramics, etc. [1]. The demand for polyvanadates is projected to grow by seven percent a year till 2025 [2], thanks to the usage of vanadium-containing steels in both traditional areas and the implementation of new technologies for the manufacture of batteries [2]. The concentrations of vanadium in surface waters caused by industrial wastes are small and, in general, are within the natural content (up to 65 μg/L) [3]. Literature sources provide information on vanadium concentrations in surface waters of industrial wastes up to 2 mg/L [3]. The chemistry of aqueous solutions of vanadium is complicated by the existence in solutions of numerous ions and polyions, the transition of which into an equilibrium state depends on the pH value and oxidation-reduction potentials in the systems, as well as its concentration and the content of other elements. Thus, the composition and stability of solid phases formed in the MeO–V2O5–H2O system, where Me is the alkaline earth element, exhibit complex functions of the chemical composition of solution, pH, etc. Experimentally, it has been demonstrated [46] that depending on the conditions in the MeO–V2O5–H2O system, different polyvanadates MeV12O31, MeV6O16, Me3V10O28, Me(VO4)2, Me2V2O7, and Me3(VO4)2 can precipitate. For the compounds of Ca2+, Sr2+, and Ba2+ the same authors determined the solubility products [4, 5], taking into account the mutual transformation reactions of the vanadium (V) ionic species. The coexistence areas of these polyvanadates were also calculated in the cited papers as a function of the solution pH and alkaline earth metal concentrations. At the same time, the reciprocal relationship between the solution composition and solid phases was established in the form of diagrams and . This procedure is not always efficient. In particular, the solubility diagrams provide valuable information only if the solubility of solid phase is low.

In this paper, based on calculations of the global Gibbs energy variation of the dissolution-precipitation process of the polyvanadates of alkaline earth metals as a function of pH and total concentrations of vanadium and alkaline earth metals, their thermodynamic stability areas are determined. It has been proved that the value of global Gibbs energy variation under real conditions constitutes a more objective criterion for estimation of the stability areas of polyvanadates of alkaline earth metals than the solubility diagrams.

2. Theoretical Part

In a series of papers [68], it has been shown that as a strict criterion of solid-phase stability serves the value of the Gibbs energy variation of the solid-phase formation-dissolution process. This criterion will be applied for the determination of the thermodynamic stability areas of the alkaline earth metal polyvanadates.

The essence of the calculation method will be explained by a concrete example of equilibrium of calcium dodecavanadate with the saturated aqueous solution:

The Gibbs energy variation under standard conditions is equal:where is the standard Gibbs energy of formation of the i species. The values used (kJ/mol) are recalculated from the selected equilibrium constants [4, 5] and shown in Table 1. The value cannot, however, serve as a characteristic of the thermodynamic stability of calcium dodecavanadate under real conditions. In the latter case, the equilibrium is described by the equation of reaction isotherm, which for Reaction (1) takes the following form:

Note that the Gibbs energy variation in Reaction (1) strongly depends on the pH value. At the same time, Equation (3) cannot also serve as a characteristic of the thermodynamic stability of calcium dodecavanadate, because the vanadyl ion, in function of pH and CV, is subjected to complex chemical transformations in the solution (Table 2).

Obviously, in calculating ΔGS, all these equilibria must be taken into account. A rigorous thermodynamic analysis, developed in a series of papers [7, 8], shows the global Gibbs energy variation of Reaction (1) under real conditions, considering the equilibria (Table 2), as described by the following equation:where is the coefficient that takes into account the contribution of the equilibria (Table 2), which is determined in the following way:where is the equilibrium concentration of vanadyl ion , computed for certain values of pH and CV from the mass balance conditions:

Analogous relations are valid for other vanadates. Under this approach, for −ΔGS> 0, the solid phase is thermodynamically unstable to dissolution according to Reaction (1) and vice versa; if −ΔGS< 0, the dissolution of the solid phase takes place.

3. Results and Discussion

Within the pH range where the solid phase (oxide) is thermodynamically unstable (dissolves), it is necessary to determine the areas of predominance of the soluble species in solution. In the presence of polynuclear species, it is necessary to calculate the partial molar fractions of the species as a function of pH and CV. The last ones for vanadium compounds (V), for example, are determined in the following way [7]:

The results of calculating these functions for vanadium (V) compounds for CV = 1 × 10−3 mol·L−1 as a function of pH are shown in Figure 1. Based on these results, it is easy to determine the thermodynamic stability of the species in the solution as a function of pH and CV. As one can see, in acidic media, within the pH range between 2.4 and 5.5, the species H2V10O284− and HV10O285− are dominant. In alkaline solutions, after pH > 8, the HVO42− species prevail.

The obtained data (Figure 1) correlate well with existing experimental evidences. It is well known that high-field NMR gives high-quality structural data and quantitative information about speciation and concentration of a large variety of species in aqueous solvents [9]. Consequently, NMR 51V spectroscopy measurements [10] confirm the speciation of vanadium (V) in aqueous vanadate solutions as a function of pH and vanadate concentrations. As seen in Figure 1, in alkaline solutions, tetrahedrally coordinated vanadates, metavanadate, and pyrovanadate are abundant. Decavanadate only forms when metavanadate is added to solutions of pH 3 or less. Acidification of metavanadate solutions to pH 4 or lower polymerizes all the oligomers to form V10. Readjusting the pH near to 8 produces a depolymerization of V10 to form V2, V4, and V5, as it has been stated in [1113]. As the pH of metavanadate solutions increases, V1 and V2 become the predominant species [14]. As reported Ralston et al. [15], in the basic pH (between 8 and 9) solution, and (or ) coexist. Orthovanadates dominate at a solution pH greater than 13 [16]. As highly basic solutions are acidified to pH values between 9 and 12, orthovanadate tetrahedral units can combine to form pyrovanadates , , and , which are colorless [16]. Continued acidification to pH values between 6 and 9 leads to the formation of colorless or yellow metavanadates , , , and [16, 17]. Additional acidification to pH between 2 and 6 causes the coordination of vanadium to change from tetrahedral to octahedral. A combination of ten octahedral units forms the decavanadate ions , , and , leading to orange or red solutions [16, 17]. Pervanadyl is yellow and dominates at dilute concentrations in strongly acidic solutions, at pH 1.0–2.5 (Figure 1).

The calculated dependence of ΔGS on pH for CV = CCa = 1 mol·L−1 for different calcium polyvanadates are shown in Figure 2.

The intersection points of the curves correspond to the boundaries of the stability areas of the respective salts. For two salts, e.g., 1 and 2, the Gibbs energy values in the given conditions are ΔGS(1) and ΔGS(2); in the case −ΔGS(1) < −ΔGS(2), salt 1 is thermodynamically more stable, and vice versa; for −ΔGS(1) > −ΔGS(2), salt 1 must transform into salt 2. In the case ΔGS(1) = ΔGS(2), equilibrium is established. In Figures 35, the ΔGS(pH) dependence for different calcium, strontium, and barium polyvanadates are displayed for several CV and CMe values. For simplicity, only the lower boundaries of the ΔGS values are shown.

As it is evident from Figures 25, the thermodynamic stability of polyvanadates towards each other strongly depends on pH and total alkaline earth metal and vanadium concentrations. The most insoluble compounds are barium polyvanadates in comparison with calcium and strontium ones. With decrease in the CV and CMe concentrations, the polyvanadates become thermodynamically unstable to dissolution. From Figures 35, one can see that alkaline earth metal polyvanadates are thermodynamically instable when CV = CMe ≤ 1 × 10−3 mol·L−1 for all the examined pH values, excepting barium polyvanadates in very acidic solutions (pH < 3.3). The last fact is not reflected in the and diagrams, that is, the main priority of the criterion proposed in this paper for assessing the thermodynamic stability of alkaline earth metal polyvanadates. The existing experimental data [4, 5, 1820] confirm the results on the thermodynamic stability of polyvanadates obtained in this paper.

The effect of temperature on the thermodynamics stability of hydrocomplexes and solubility of polyvanadates for different temperatures are estimated on the basis of van’t Hoff equation [7]:

The appropriate values of the standard enthalpies ΔH can be found in handbooks available everywhere [2127]. Usually, it is supposed that and the temperature insignificantly influences the ΔH values within the investigated temperature interval.

Vanadium in inorganic media (a mixture of an inorganic acids and their salts) may form complexes with such anions such as , etc. Note that the sulphate medium is important due to its industrial applications [2832]. Schematically, if L- denotes the ligand, the reaction of complex formation may be written as follows:

The thermodynamic equilibrium formation constant of at 25°C is equal to [28]. Consequently, the complexes may form in the presence of a large excess of ligand, . The contribution of anions of aqueous electrolytes on the thermodynamics stability of soluble and insoluble species can be taken into account through the alpha coefficient, in the form described in Equation (6) [7, 8].

4. Conclusions

(i)A thermodynamic method for calculating the global Gibbs energy variation for describing heterogeneous equilibria of transformations of polyvanadates of calcium, strontium, and barium, that occur in systems MeO-V2O5-H2O, where Me is the alkaline earth element, has been developed and applied.(ii)Its quintessence consists in the thermodynamic analysis of the real conditions of various processes on the basis of their total thermodynamic characteristics. On the basis of the selected thermodynamic data for involved soluble and insoluble species, the thermodynamic stability areas of solid polyvanadates towards the solution pH and vanadium and alkaline earth metal ion concentrations in heterogeneous mixtures have been established, taking into account the complex formation reactions in multicomponent heterogeneous systems.(iii)The calcium, strontium, and barium polyvanadates are dissolved for CV = CMe ≤ 1 × 10−3 mol·L−1 for all the examined pH values.(iv)The thermodynamic stability of soluble species as function of pH and CV has been determined.(v)The obtained results, based on the thermodynamic analysis and graphical design of the calculated data, are in good agreement with the available experimental data.

Data Availability

The standard Gibbs energies of formation of chemical species used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partially supported by the Academy of Sciences of Moldova (grant number 15.817.02.15A).