Next Article in Journal
Association Measure and Compact Prediction for Chemical Process Data from an Information-Theoretic Perspective
Next Article in Special Issue
Interaction between Iron Fluoride and Molten FLiBe
Previous Article in Journal
The Optimized Roadway Layouts and Surrounding Rock Control Technology of the Fully Mechanized Mining Surface with Large Mining Heights in High-Gas Mines
Previous Article in Special Issue
Synthesis, Hydration Processes and Ionic Conductivity of Novel Gadolinium-Doped Ceramic Materials Based on Layered Perovskite BaLa2In2O7 for Electrochemical Purposes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ionic Equilibria in Polytungstate Melts

by
Alexander V. Kosov
*,
Olga L. Semerikova
,
Sergey V. Vakarin
,
Olga V. Grishenkova
,
Alexey S. Vorob’ev
,
Anastasia O. Khudorozhkova
and
Yury P. Zaikov
Department of Electrolysis, Institute of High Temperature Electrochemistry, Ural Branch of the Russian Academy of Sciences, Yekaterinburg 620990, Russia
*
Author to whom correspondence should be addressed.
Processes 2022, 10(12), 2658; https://doi.org/10.3390/pr10122658
Submission received: 10 November 2022 / Revised: 2 December 2022 / Accepted: 6 December 2022 / Published: 9 December 2022

Abstract

:
Polytungstate melts are used for the electrodeposition of oxide tungsten bronzes (OTBs). The scarce information on the ionic composition and properties of these electrolytes hinders effective control of the electrochemical synthesis of OTBs with desired electrical and optical properties. In this work, a comprehensive study of Na2WO4–WO3 melts that contained up to 55 mol% of tungsten trioxide was performed in the temperature range from 983 to 1073 K. Melt densities were measured using the Archimedes method. DFT calculations were carried out for various tungsten-containing compounds, including W x O 3 x 1 2 + , W x O 3 x + 1 2 , NaW x O 3 x + 1 , and Na 2 W x O 3 x + 1 . The calculated values of the W–O bond energy indicate that the tested compounds are stable in the specified temperature range, and the   WO 2 2 + cation is the most stable. The experimental dependences of the redox potential on the mole fraction of tungsten trioxide in the Na2WO4–WO3 melt were obtained using the EMF method. A model that considers the processes of interaction between tungsten-containing ions and O2− ions was proposed for the quantitative interpretation of these dependences. The equilibrium constants were found through fitting according to the Levenberg–Marquardt algorithm. The effect of the WO3 mole fraction and temperature on the concentrations of WO 4 2 , W 2 O 7 2 , W 3 O 10 2 , W 4 O 13 2 , WO 2 2 + , and O2− ions was analyzed.

1. Introduction

Oxide tungsten bronzes (OTBs), i.e., nonstoichiometric compounds such as MxWO3 (here, M is an alkali metal and 0 < x < 1), have a wide range of composition-dependent properties [1,2]. The diverse electrophysical and optical properties of MxWO3 with different x and M [1,2,3,4,5,6] determine the broad use of OTBs in many fields of catalysis [7,8,9] and photocatalysis [10,11,12], including air and wastewater purification, anti-virus sterilization and hydrogen and oxygen production. OTB-based materials are also in demand for electrochromic [13], plasmonic [14], innovative biomedical [15] and field emission [16] applications, and NIR shielding [17].
Electrodeposition from polytungstate melts is one of the most promising methods for the synthesis of OTBs [18,19,20,21] and multilayered hybrid systems with nanocrystalline OTBs [9,22,23,24]. The main advantages of this method are the high process rate, the ability to change the composition and structure of deposits by varying the electrolysis parameters, and low capital and operating costs. Evidently, understanding the regularities of the processes that occur in the polytungstate electrolyte and at the electrolyte/electrode interface is necessary for efficient control of electrodeposition and obtaining OTBs with the desired physicochemical properties. However, the mechanism of OTB electrocrystallization is difficult to determine without the availability of sufficient information on the structure and ionic composition of polytungstate melts. The limited data lead to the emergence of unfounded hypotheses about the mechanisms of OTB formation in these melts [21,25,26,27,28].
The structures of molten Na2WO4 at 1033 K [29] and Na2W2O7 at 1087 K [30] were investigated using the radial distribution function based on X-ray scattered intensity data. The presence of tetrahedral (WO4)2− anions in molten Na2WO4 was proven [29]. In the Na2WO4–WO3 (1:1) system [30], it was found that the pattern of the [(W2O7)2−] anion chain in the crystalline state is practically preserved in the molten state within a small range (<0.45 nm), but the long-range order within one anion chain is almost lost in the melt and the neighboring chains are randomly oriented relative to each other. A Raman scattering study of molten alkali–metal tungstates [31,32] and transparent glasses obtained by rapid quenching of ditungstate melts [32] showed the presence of isolated [WO4]2− tetrahedral anions in Na2WO4 and K2WO4 melts and complexes that consisted of two corner-sharing tetrahedra (i.e., [W2O7]2−) in Na2W2O7 and K2W2O7 melts. The vibration modes of [W2O7]2− dimers were also found during the study of the microstructure of molten M2W2O7 (M = Li, Na, K) by in situ high-temperature Raman spectroscopy, in combination with density functional theory (DFT) analysis [33,34]. The cation influence is mainly limited by the change in the characteristics of the W–Onb bond (nb is non-bridging oxygen) in [W2O7]2− [31,33]. A joint analysis of the data of in situ high-temperature Raman spectroscopy and a quantum-chemical ab initio calculation for the evolution of stably existing structural units in Li2O–WO3 melts was carried out in [35]. It was concluded that complexes ([WO4]2−) and chains ([W2O7]2−, [W3O10]2−, and [W4O13]2−) composed of two, three, or four [WO4]2− are formed by sharing their corner oxygen atoms when the mole ratio of Li2O:WO3 in the melt is 1:1, 1:2, 1:3, and 1:4, respectively [35]. Raman spectroscopy or X-ray scattering data for the M2O–WO3 melts with any intermediate component ratio (for example, 1:1.5) are not available, which does not allow us to predict changes in the structure or structure-dependent properties of M2WO4–WO3 melts.
The data obtained by the EMF method can be used to estimate the ratio of anions in polytungstate melts [18,36,37,38]. The approach is based on the analysis of the experimental dependence of the potential difference between two platinum–oxygen electrodes semi-immersed in melts with different concentrations of tungsten trioxide (the melts are separated by a porous diaphragm); assumptions about ionic equilibria are required for analysis. In [18], this method was applied to determine the WO 4 2 / W 2 O 7 2 ratios in Na2WO4–M’2WO4–WO3 (M* = Li or K) melts that contained up to 20 mol% WO3. The results of [18] satisfactorily explained the change in the OTB composition with an increase in the tungsten trioxide content in the above concentration range. The M2WO4–WO3 melts with a WO3 mole fraction up to 0.5 were studied using the EMF method in [38]. The authors considered these melts as completely dissociated mutual solutions of Na2WO4 and Na2W2O7 that consisted of an ion mixture (Na+,   WO 4 2 , and W 2 O 7 2 ). A conclusion was made about the weak interionic interaction in this system. An acceptable agreement was achieved between the experimental and calculated results at a WO3 mole fraction up to 0.4 by introducing the activity coefficients of ditungstate ions [38]. According to [18,38], a model that considers the formation of W 3 O 10 2 and W 4 O 13 2 is required to more accurately describe the experimental results at high WO3 concentrations in melts.
The above literature data confirm that polytungstate melts contain various tungsten-containing anions and alkali metal cations. However, it is unlikely that the anions are directly involved in the cathodic process that leads to the electrodeposition of OTB and/or tungsten on the cathode. An assumption was made in [18] about the formation of W6+ cations as a result of the following processes:
4 W 2 O 7 2     7 W 2 O 7 2 + W 6 + ,
7 W 4 O 13 2     13 W 2 O 7 2 + 2 W 6 +
However, the authors provide no evidence to support this hypothesis. Thus, the problem of the ionic composition and its changes as the proportion of tungsten trioxide in the melt increases remains topical.
In this work, the stability of various tungsten-containing anions and cations is evaluated, and a model is proposed for calculating the ionic composition of Na2WO4–WO3 melts with a WO3 mole fraction up to 0.55.

2. Materials and Methods

2.1. Preparation of Electrolytes

Sodium tungstate and tungsten trioxide (purity 99.9 wt%, Vecton, RF) were used to prepare the melts. The dried (523 K, 2 h) reagents were weighed using a VK-600 balance (Massa-K, RF) with an accuracy of ±0.01 g, and then the electrolyte components were mixed in the required proportions in a porcelain container. The WO3 mole fraction in the mixtures, ν, ranged from 0 to 0.55. The prepared mixtures were placed in platinum or alumina crucibles (see Section 2.2 and Section 2.3) and heated to the experimental temperature (983, 1023 or 1073 K). In accordance with the phase diagram of the Na2WO4–WO3 system [18], the composition range was limited to 0 ≤ ν ≤ 0.40 at 983 K and 0 ≤ ν ≤ 0.55 at 1023 and 1073 K, respectively.

2.2. Electrochemical Measurements

All the electrochemical measurements were performed using Autolab PGSTAT302N (Metrohm, The Netherlands) with Nova 1.9 software. Figure 1 shows the electrochemical cell designs described in detail in Section 2.2.1 and Section 2.2.2.
All high-temperature experiments were carried out in a shaft furnace with a quartz protective vessel installed inside. The melt temperature was controlled using a Pt/Pt-Rh thermocouple and a Varta TP703 temperature controller (Varta, RF) with an accuracy of ±1 K. The operating temperature range corresponded to the usual conditions for OTB electrodeposition from Na2WO4–WO3 melts.

2.2.1. Measurement of the Oxygen Function of the Pt Electrode

The dependences of the potential difference between two Pt electrodes half-immersed in the melt on the partial pressure of oxygen over one of the electrodes were measured at 983 K. This was carried out to confirm that the reaction
1 2 O 2 + 2 e     O 2
is the potential-determining process of the platinum–oxygen electrode. Both electrodes were made of Pt foil; the surface area of each electrode was S = 0.9 cm2. Before measurements, one of the electrodes and a quartz test tube equipped with a vacuum rubber stopper, a spout for gas removal, and a bottom channel for electrolyte inflow were lowered into the 0.8Na2WO4–0.2WO3 melt in a platinum crucible (Figure 1a). The Pt electrode and a quartz tube for supplying a mixture of air with high-purity argon (99.999%, Uralcryogas, RF) were inserted into the stopper holes beforehand. The volume fraction of air in the gas flow was changed using a mass flow controller (LOW-ΔP-FLOW F-201DV) and a mass flow meter (LOW-ΔP-FLOW F-101D) for low pressure drops or corrosive gas (Bronkhorst High-Tech B.V., Ruurlo village, The Netherlands).

2.2.2. EMF Measurement

The EMF of the cell
Pt ,   O 2 air | 0 . 8 Na 2 WO 4 0 . 2 WO 3 ¦ 1 ν Na 2 WO 4 ν WO 3 | O 2 ( air ) ,   Pt
was measured at 983, 1023 and 1073 K in air ( p O 2 = 21.3 kPa). The cell design is presented in Figure 1b. A small thin-walled alumina crucible 1.8 cm in diameter with the (1 − ν)Na2WO4–νWO3 melt was placed on the bottom of a platinum crucible 3.5 cm in diameter, which served as a container for the standard melt (0.8Na2WO4–0.2WO3 in this work). The electrodes (Pt foil, S = 0.9 cm2) were half-immersed in the above melts before measurements. Alumina tubes protected the platinum current leads; air was blown through the second opening of these tubes. The recording of the potential difference was stopped when the ΔE value remained constant (within ±0.5 mV) for 30 min.

2.3. Density Measurement

The densities of the Na2WO4–WO3 melts with a mole fraction of tungsten trioxide up to 0.5 were measured using the Archimedean method, according to the technique described in [39,40]. A platinum sphere was suspended on a platinum wire 0.5 mm in diameter and 0.6 m long attached to a Mettler AT20 analytical balance (Mettler Toledo, Billerica, MA, USA). A weighed amount of a prepared oxide–salt mixture was placed in an alumina crucible. The measurements were performed within a temperature range from 983 K to 1073 K. The platinum sphere was placed into the melt at 1073 K, and temperature dependence of the sphere weight was registered. The density was calculated according to the equation
ρ = Δm/Vs,
where Δm (g) is the difference between the sphere weight in the melt and in air, and Vs (cm3) is the sphere volume. The cargo weight in the air was 14.5126 g. The temperature dependence of the sphere volume was determined by calibration using molten chlorides [40]. The experimental error for the density did not exceed ±0.5%.

2.4. Density Functional Theory (DFT) Calculations

Calculations of the stability of tungsten-containing compounds ( W x O 3 x 1 2 + , W x O 3 x + 1 2 , NaW x O 3 x + 1 , and Na 2 W x O 3 x + 1 ) were performed using the Siesta software package [41] on a cluster-type hybrid computer with 1864 CPUs and a peak performance of 216 Tflops. Geometric optimization was carried out using the general gradient approximation method in the PBE form [42] for all proposed systems. The dynamic relaxation of atoms continued until the change in the total energy of the system became less than 0.1 meV. The cutoff energy of the plane wave basis was 400 Ry. A cubic unit cell with a face translation vector of 20 Å was used. After geometric optimization, the resulting systems were tested for temperature stability using ab initio molecular dynamics with a Nose–Hoover thermostat [43] at 983, 1023 and 1073 K for 1000 time steps with a step of 1 fs.

3. Experimental Results

3.1. Electrochemical Measurements

The data obtained in the experiment described in Section 2.2.1 and their analysis are shown in Figure 2. Figure 2a demonstrates the experimental dependence of the potential difference, ΔE (V), between two Pt electrodes half-immersed in the 0.8Na2WO4–0.2WO3 melt with a stepwise change in the volume fraction of air in the gas mixture, Wair (vol%), above one of the electrodes. Evidently, an increase in Wair causes a decrease in the ΔE value. The dependence of ΔE on the logarithm of the air volume fraction is linear (Figure 2b) and the slope is close to theoretical one, as shown in the following equation:
Δ E = k T 4 e ln [ O 2 ] 1 [ O 2 ] 2 = k T 4 e ln p O 2   air W air p O 2   air = 21.54 10 3 T ln W air ,
where k (1.38 × 10−23 J K−1) is the Boltzmann constant, T (K) is the absolute temperature, e (C) is the elementary electric charge, [O2]1 and [O2]2 (mol L−1) are the concentrations of oxygen in the atmosphere above the first and second Pt electrodes, respectively (see Figure 1a), and p O 2   air is the partial pressure of oxygen in air. The error is small and equal to 0.39%. Formula (6) was derived using the Nernst equation for Reaction (3) and is as follows:
E = E 0 + k T 2 e ln [ O 2 ] 1 2 [ O 2 ] ,
where E0 (V) is the standard electrode potential and [O2] and [O2−] (mol L−1) are the concentrations of oxygen (above the melt) and O2− ions (in this melt), respectively. This means that the potential of the platinum electrode is indeed determined by Reaction (3).
In this case, the potential difference measured in Cell (4) is described by the following equation:
Δ E = k T 2 e ln [ O 2 ] melt   1 [ O 2 ] melt   2 ,
This allows us to find the ratio of oxide–ion concentrations in melts with different WO3 contents. Equation (8) can be rewritten as follows:
Δ E = k T 2 e ln g 1 V m 2 g 2 V m 1 ,
where g (mol) is the amount of O2− ions in the melt, Vm (cm3 mol−1) is the molar volume of the melt, and subscripts 1 and 2 refer to the (1 − ν)Na2WO4–νWO3 melt and the standard melt (0.8Na2WO4–0.2WO3 in this work), respectively.
The experimental concentration dependences of ΔE measured in Cell (4) at 983, 1023 and 1073 K are shown in Figure 3. Since the decrease in ΔE is associated with a change in the O2− concentration in the melt due to an increase in the WO3 mole fraction (see Equation (8)), these experimental dependences can be interpreted as follows. In pure molten sodium tungstate, the O2− concentration is significantly higher than in the standard melt (by six orders of magnitude at 983 K). A sharp decrease in the O2− concentration occurs when the WO3 mole fraction is only 0.01 (about 4 orders of magnitude at 983 K). In the range of ν from 0.05 to 0.55, the concentration of O2− ions decreases approximately according to the logarithmic law, since the ΔE(ν) dependences decrease almost linearly. At 1073 K, the range of the O2− concentration becomes noticeably narrower, as demonstrated by the ΔE values at ν = 0 and ν > 0.4 (see Figure 3).
The quantitative interpretation of these data is difficult because the dependences of g and V on the WO3 mole fraction are unknown. Dependences V(ν) can be found from the measured densities of melts with different contents of tungsten trioxide (see Section 3.2). The situation for g(ν) is much more complicated, since information is required on the processes in the melt and equilibria with the participation of O2− ions, as well as on the corresponding equilibrium constants. The absolute values of the equilibrium concentrations of ions in the melt can be found only in this case.
The problem is that there is still no unambiguous opinion regarding the ionic composition of polytungstate melts. Published data [29,30,31,32,33,34] show that various polytungstate anions with the general formula W x O 3 x 1 2 + can exist in the melt. At the same time, the results obtained by electrochemical methods (see, for example, [9,18,21]) indicate the presence of tungsten-containing cations in the melt, but do not allow one to determine their composition and charge. Nevertheless, the hypothesis of the presence of W6+ ions in the melt seems unlikely, due to the high molecularity of Equilibria (1) and (2) and the high concentration of oxygen-containing ions. It is more likely that W x O 3 x 1 2 + cations can be formed as a result of the dissociation of W x O 3 x + 1 2 anions. However, we must ensure that W x O 3 x 1 2 + cations are stable at the experimental temperatures (see Section 4) before modeling equilibria that involve these cations (see Section 5).

3.2. Density Measurement

The experimental dependences of the density on the WO3 mole fraction in (1 − ν)Na2WO4–νWO3 melts are shown in Figure 4. Figure 4a demonstrates that the decrease in density occurs both with an increase in temperature and with a decrease in ν. The measured dependences of density on temperature and concentration of tungsten trioxide can be described by the following equations:
ρ = −2.64·10−3·Tν −7.69·10−4·T + 3.781ν + 4.555, 0 < ν ≤ 0.3,
ρ = −1.20·10−3·Tν −1.01·10−3·T + 2.362ν + 4.816, 0.3 < ν ≤ 0.5.
Note that these dependences are much more accurate than the approximate equation proposed in [44] (see Figure 4a):
ρ = 4.808 +1.282 ν −0.962·10−3·T, 0 ≤ ν ≤ 0.2.
The dependences of the density of the tested melt on the WO3 mole fraction are presented in Figure 4b. These dependencies can be described by the equations
ρ(ν) = −5.4406ν3 + 4.0947ν2 + 0.4392ν + 3.8182, T = 983 K,
ρ(ν) = −3.8131ν3 + 3.3646ν2 + 0.3974ν + 3.7866, T = 1023 K,
ρ(ν) = −3.7622ν3 + 3.513ν2 + 0.2277ν + 3.7489, T = 1073 K.
These equations are suitable for calculating the dependence of the molar volume of the melt on the WO3 mole fraction.
V m ( ν ) = ( 1 ν ) M N a 2 WO 4 + ν M WO 3 / ρ ( ν ) ,
where M (g mol−1) is the molar mass.

4. DFT Calculation

To determine the stability of W x O 3 x 1 2 + and W x O 3 x + 1 2 compounds, the W–O bond energy, Eb (eV), was calculated according to the following formula:
Eb = − (EtotE1WNWE1ONO)/Nb,
where Etot (eV) is the total energy of the compound, E1W and E1O (eV) are the energies calculated for single W6+ and O2− ions, respectively, NW and NO are the amount of tungsten and oxygen ions in the system, respectively, and Nb is the amount W–O bonds in the compound.
The W–O binding energies in the NaW x O 3 x + 1 and   Na 2 W x O 3 x + 1 compounds were calculated, taking into account the interaction of Na+ with tungsten-containing anions.
Eb = (EtotEbNaE1WNWE1ONO)/Nb,
EbNa = EtotEWOENa,
where EbNa (eV) is the binding energy between Na+ and W x O 3 x + 1 2 , EWO (eV) is the energy calculated for the system without Na+, and ENa (eV) is the energy calculated for the sodium subsystem only.
The calculation results indicate that all the tested compounds are stable in the temperature range from 983 to 1073 K. The geometric structures of W x O 3 x 1 2 + and W x O 3 x + 1 2 at 983 K are shown in Figure 5a. An increase in temperature does not affect the geometric structure of the compounds in most cases. The only exception is the rearrangement of the [W4O13]2− compound at 1073 K (Figure 5b).
Figure 6a shows the dependences of the W–O bond energy on the number of tungsten atoms in W x O 3 x 1 2 + , W x O 3 x + 1 2 , NaW x O 3 x + 1 and Na 2 W x O 3 x + 1 compounds at 983 K. The Eb value in W x O 3 x 1 2 + is lower than in W x O 3 x + 1 2 by 54% at x = 1 and by 6.7% at x = 4. The W–O bond energy for the WO 2 2 + cation is minimal (Eb = −105.48 eV), which indicates its high stability. An increase in the number of tungsten atoms in W x O 3 x 1 2 + and W x O 3 x + 1 2 leads to an increase in the Eb value from −105.48 to −66.13 eV and from −68.41 to −62.01 eV, respectively. The addition of one or two sodium ions to the second coordination sphere contributes to a decrease in the W–O bond energy; namely, the Eb values are 2.5–8.3% (for NaW x O 3 x + 1 ) or 4.6–13.8% (for Na 2 W x O 3 x + 1 ) lower than those for W x O 3 x + 1 2 .
At 1073 K, the Eb value increases slightly in most cases (Figure 6b). The maximum increase is observed for Na2[W2O7] and does not exceed 0.7%. On the contrary, the W–O bond energy in the [W4O13]2− compound decreases from −62.01 eV at 983 K to −65.71 eV at 1073 K, which is associated with the breaking of one W–O bond (see Figure 5).

5. Model of Ionic Equilibria

The above results of DFT calculations and published data [29,30,31,32,33,34] indicate the presence of several species of ions in the Na2WO4–WO3 melt. The most probable set of ions includes Na+, WO 4 2 , W 2 O 7 2 , W 3 O 10 2 , W 4 O 13 2 , WO 2 2 + , and O2−. One can assume that the ratios between tungsten-containing ions depend on the WO3 mole fraction in the initial oxide–salt mixture before its melting. If the initial mixture consists of ν mol of WO3 and (1 − ν) mol of Na2WO4, then the melting process can be described as follows:
( 1 ν ) Na 2 WO 4 + ν WO 3 = 2 ( 1 ν ) Na + + a WO 2 2 + + b WO 4 2 + c W 2 O 7 2 + d W 3 O 10 2 + f W 4 O 13 2 + g O 2 ,
where 0 ≤ ν ≤ 0.55 and the coefficients a, b, c, d, f, and g are equal to the number of moles of the corresponding ion if the total amount of the initial mixture is one mol. The balances of tungsten atoms and electric charges for reaction (20) give the following equalities::
a + b + 2 c + 3 d + 4 f = 1 ,
b + c + d + f + g a = 1 ν .
The tungsten-containing ions interact and are in dynamic equilibrium with each other and with O2− ions.
WO 4 2     WO 2 2 + + 2 O 2 ,   k 0 = [ WO 2 2 + ]   [ O 2 ] 2 [ WO 4 2 ] = a g 2 b V m 2 ,
WO 4 2 + W 3 O 10 2     2 W 2 O 7 2 ,   k 1 = c 2 b d ,
W 2 O 7 2 + W 4 O 13 2     2 W 3 O 10 2 ,   k 2 = d 2 c f ,
2 WO 4 2     W 2 O 7 2 + O 2 ,   k 3 = c g b 2 ,
3 W 2 O 7 2     2 W 3 O 10 2 + O 2 ,   k 4 = d 2 g c 3 ,
4 W 3 O 10 2     3 W 4 O 13 2 + O 2 ,   k 5 = f 3 g d 4
where k0, k1, k2, k3, k4, and k5 are the equilibrium constants of the corresponding processes. Equilibria (24) and (25) were written under the assumption that all the tungsten trioxide is consumed for the formation of tungsten-containing anions as a result of the reactions WO 4 2 + WO 3 W 2 O 7 2 , W 2 O 7 2 + WO 3 W 3 O 10 2 , and W 3 O 10 2 + WO 3 W 4 O 13 2 .
It is easy to prove that the equilibrium constants k3, k4 and k5 are interdependent, as demonstrated by the following equations:
g = k 3 b 2 c = k 4 k 1 2 b 2 c = k 5 k 1 2 k 2 3 b 2 c ,
k 4 = k 3 k 1 2 ,   k 5 = k 3 k 1 2 k 2 3 .
Therefore, only the constants k0, k1, k2 and k3 are required to determine the coefficients in Equation (20) for any ν value. Accordingly, the problem is reduced to solving a system of non-linear equations.
a + b + 2 c + 3 d + 4 f = 1 b + c + d + f + g a = 1 ν log a + 2 log g log b = log k 0 + 2 log V m 2 log c log b log d = log k 1 2 log d log c log f = log k 2 log c + log g 2 log b = log k 3 .
The values of the equilibrium constants that provide the best agreement between the experimental and calculated ΔE(ν) dependences can now be found using a non-linear fit. As a result, a(ν), b(ν) c(ν), d(ν), f(ν), and g(ν) that correspond to the experimental conditions can be determined.

6. Results and Discussion

To calculate ΔE(ν), Equation (32) was used in the following form:
Δ E ( ν ) = k T 2 e ln g ( ν ) V m ( ν = 0.2 ) g ν = 0.2 V m ( ν ) .
Dependences Vm(ν) were found using Equation (16). To find g(ν), the numerical solution of System (31) was first performed with random values of k0, k1, k2 and k3, and then the equilibrium constants were iteratively changed according to the Levenberg–Marquardt algorithm, until the best agreement between the experimental and calculated ΔE(ν) dependences was achieved.
Figure 7 shows a comparison of the experimental and fitted ΔE(ν) dependences at 983, 1023 and 1073 K. Table 1 gives the values of the equilibrium constants for Reactions (23)–(26), i.e., k0, k1, k2 and k3, at which the optimal coincidence of the curves was achieved.
Figure 8 shows the dependences of the concentrations of WO 2 2 +   WO 4 2 , W 2 O 7 2 , W 3 O 10 2 , W 4 O 13 2 , and O2− ions on the WO3 mole fraction at 983, 1023 and 1073 K that correspond to the obtained values of the equilibrium constants. For some ν values, the concentrations of these ions and Na+ ions are presented in Table 2. The calculation results indicate the following. The Na+ and WO 4 2 ions predominate in molten sodium tungstate, which is consistent with the results of structural studies of this melt [29,31,32]. The fraction of other ions included in our model does not exceed 0.008%. The addition of even very small amounts of tungsten trioxide leads to a sharp decrease in the concentration of O2− ions, which is especially noticeable at lower temperatures (Figure 8f). The WO 4 2 concentration, [ WO 4 2 ] , also decreases as the WO3 mole fraction increases, and the decrease in [ WO 4 2 ] is directly proportional to ν and almost does not depend on temperature up to ν ≈ 0.3 (Figure 8b).
The W 2 O 7 2 concentration initially increases in proportion to ν (Figure 8c), and the slopes are almost equal at lower temperatures. At 1073 K, the slope is noticeably lower and the maximum value is reached at ν = 0.43. Interestingly, the melt contains not only Na+ cations and W 2 O 7 2 anions at ν = 0.5, contrary to the a priori assumption in [30,33,34,35]. Our calculations indicate that the fraction of W 2 O 7 2 in the 0.5Na2WO4–0.5WO3 melt (among other tungsten-containing ions) is only 72.7% at 1023 K and 51.2% at 1073 K.
The WO 2 2 + concentration (Figure 8a) increases dramatically as ν increases. For example, at T = 983 K, the [ WO 2 2 + ] = 2.57 × 106 cm−3 at ν = 0, while the [ WO 2 2 + ] = 2.72 × 1017 cm−3 at ν = 0.1. The WO 2 2 + concentration, as is the case for the W 2 O 7 2 concentration, increases more slowly at 1023 K. At a high mole fraction of WO3, the concentration of these cations is inferior to the concentrations of WO 4 2 and W 2 O 7 2 only by 1–2 orders of magnitude.
The concentrations of W 3 O 10 2 and W 4 O 13 2 ions increase as the WO3 mole fraction and temperature increase (Figure 8d,e). The largest values of d [ W 3 O 10 2 ] / d ν and d [ W 4 O 13 2 ] / d ν are observed at low ν values. When T = 1073 K and ν ≥ 0.4, the W 4 O 13 2 concentration exceeds those of W 3 O 10 2 and WO 2 2 + ions.
Thus, the proposed model allows us to calculate the equilibrium concentrations of ions, which are necessary for simulating electrodeposition from polytungstate melts and determining the composition of cathode products.

7. Conclusions

The problem of the ionic composition of Na2WO4–WO3 melts was studied within the framework of a model that considers the equilibria between tungsten-containing and O2− ions. The concentrations of O2−, WO 2 2 + , WO 4 2 , W 2 O 7 2 , W 3 O 10 2 , and W 4 O 13 2 ions were determined by fitting the EMF dependences calculated by Equation (32) to the experimental EMF values measured in the cell Pt ,   O 2 ( air ) | 0 . 8 Na 2 WO 4 0.2 WO 3 ¦ ( 1 ν ) Na 2 WO 4 ν WO 3 |   O 2 ( air ) , Pt with 0 ≤ ν ≤ 0.55 at 983, 1023 and 1073 K. Preliminarily, the adequacy of the approach was confirmed by measuring the EMF of the cell Pt ,   O 2 ( air + Ar ) | 0 . 8 Na 2 WO 4 0.2 WO 3 ¦ 0.8 Na 2 WO 4 0.2 WO 3 |   O 2 ( air ) , Pt , and the stability of tungsten-containing anions and cations at experimental temperatures was established using DFT calculations.
The calculated dependences of the ion concentrations on the mole fraction of tungsten trioxide, taking into account the change in the density of the melt, indicate the following. In molten sodium tungstate, Na+ and WO 4 2 ions predominate, and the fraction of other ions is less than 0.01%. The addition of tungsten trioxide causes a sharp decrease in the concentration of O2− ions and a sharp increase in the concentrations of WO 2 2 + , W 3 O 10 2 , and W 4 O 13 2 . In addition, an increase in the WO3 mole fraction leads to a decrease in the WO 4 2 concentration and an increase in the W 2 O 7 2 concentration. It was found that the W 2 O 7 2 ion does not dominate among other tungsten-containing ions, even in the 0.5Na2WO4–0.5WO3 melt, and its main competitors are WO 4 2 and WO 2 2 + at 1023 K or WO 4 2 , W 4 O 13 2 and W 3 O 10 2 at 1073 K.
In conclusion, we emphasize that the main advantage of the proposed approach is the ability to obtain quantitative information about the change in the ionic composition of the melt that is associated with a change in temperature and/or the ratio of components in the initial oxide–salt mixture, based on a simple procedure for analyzing experimental data.

Author Contributions

Conceptualization, A.V.K.; methodology, A.V.K., A.S.V. and A.O.K.; software, A.V.K. and A.S.V.; validation, S.V.V. and O.V.G.; formal analysis, A.V.K., O.V.G. and O.L.S.; investigation, A.V.K., O.L.S., A.S.V. and A.O.K.; resources, Y.P.Z.; writing—original draft preparation, O.V.G. and A.V.K.; writing—review and editing, O.V.G. and A.V.K.; supervision, Y.P.Z. and S.V.V.; project administration, Y.P.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

DFT calculations were performed on the hybrid cluster-type computer “URAN” at the N.N. Krasovskii Institute of Mathematics and Mechanics of the Ural Branch of the Russian Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bartha, L.; Kiss, A.B.; Szalay, T. Chemistry of tungsten oxide bronzes. Int. J. Refract. Met. Hard Mater. 1995, 13, 77–91. [Google Scholar] [CrossRef]
  2. El-Sayed, A.M.; Mousa, S.M.A. Some properties of sodium tungsten bronzes as a function of sodium concentration. Ind. J. Chem. Technol. 2005, 12, 304–308. [Google Scholar]
  3. Sano, K.; Nitta, Y.; Ōno, Y. Transition temperature of superconductivity in sodium tungsten bronze—Theoretical study based on first-principles calculations. J. Phys. Soc. Jpn. 2020, 89, 113704. [Google Scholar] [CrossRef]
  4. Lekshmi, I.C.; Hegde, M.S. Synthesis and electrical properties of cubic NaxWO3 thin films across the metal–insulator transition. Mater. Res. Bull. 2005, 40, 1443–1450. [Google Scholar] [CrossRef] [Green Version]
  5. Raj, S.; Matsui, H.; Souma, S.; Sato, T.; Takahashi, T.; Chakraborty, A.; Sarma, D.D.; Mahadevan, P.; Oishi, S.; McCarroll, W.H.; et al. Electronic structure of sodium tungsten bronzes NaxWO3 by high-resolution angle-resolved photoemission spectroscopy. Phys. Rev. B 2007, 75, 155116. [Google Scholar] [CrossRef] [Green Version]
  6. Bikberdina, N.Y.; Boronenko, M.P.; Gulyaev, P.Y.; Yunusov, R.D. Electrophysical properties of oxide bronze NaxWO3. Bull. Yugorsk State Univ. 2017, 3, 7–11. [Google Scholar] [CrossRef] [Green Version]
  7. Petrov, L.A.; Shishmakov, A.B.; Vakarin, S.V.; Semerikova, O.L.; Melyaeva, A.A.; Mikushina, Y.V.; Zaykov, Y.P.; Chupakhin, O.N. Behavior of nanosized tungsten oxide bronzes produced by high-temperature electrolysis in model processes of desulfurization of petroleum products. Russ. J. Inorgan. Chem. 2014, 59, 7–10. [Google Scholar] [CrossRef]
  8. Vakarin, S.V.; Melyaeva, A.A.; Semerikova, O.L.; Kondratuk, V.S.; Pankratov, A.A.; Plaksin, S.V.; Porotnikova, N.M.; Zaykov, Y.P.; Petrov, L.A.; Mikushina, Y.V.; et al. Catalase activity of coarse grained and nanosized oxide tungsten bronzes obtained by electrolysis of molten salts. Russ. Chem. Bull. 2011, 60, 1985–1988. [Google Scholar] [CrossRef]
  9. Semerikova, O.L.; Vakarin, S.V.; Kosov, A.V.; Plaksin, S.V.; Pankratov, A.A.; Grishenkova, O.V.; Zaykov, Y.P.; Shishmakov, A.B.; Mikushina, Y.V.; Petrov, L.A. Electrochemical synthesis of nanohybrid systems based on copper and the oxide tungsten bronzes. J. Electrochem. Soc. 2019, 166, D792–D797. [Google Scholar] [CrossRef]
  10. Wu, C.-M.; Naseem, S.; Chou, M.-H.; Wang, J.-H.; Jian, Y.-Q. Recent advances in tungsten-oxide-based materials and their applications. Front. Mater. 2019, 6, 49. [Google Scholar] [CrossRef] [Green Version]
  11. Yi, L.; Zhao, W.; Huang, Y.; Wu, X.; Wang, J.; Zhang, G. Tungsten bronze Cs0.33WO3 nanorods modified by molybdenum for improved photocatalytic CO2 reduction directly from air. Sci. China Mater. 2020, 63, 2206–2214. [Google Scholar] [CrossRef]
  12. Wang, L.; Zhan, J.; Fan, W.; Cui, G.; Sun, H.; Zhuo, L.; Zhao, X.; Tang, B. Microcrystalline sodium tungsten bronze nanowire bundles as efficient visible light-responsive photocatalysts. Chem. Commun. 2010, 46, 8833–8835. [Google Scholar] [CrossRef]
  13. Zimmer, A.; Gilliot, M.; Tresse, M.; Broch, L.; Tillous, K.E.; Boulanger, C.; Stein, N.; Horwat., D. Coloration mechanism of electrochromic NaxWO3 thin films. Opt. Lett. 2019, 44, 1104–1107. [Google Scholar] [CrossRef]
  14. Tegg, L.; Cuskelly, D.; Keast, V.J. The sodium tungsten bronzes as plasmonic materials: Fabrication, calculation and characterization. Mater. Res. Express 2017, 4, 065703. [Google Scholar] [CrossRef]
  15. Jie, S.; Guo, X.; Ouyang, Z. Tumor ablation using novel photothermal NaxWO3 nanoparticles against breast cancer osteolytic bone metastasis. Int. J. Nanomed. 2019, 14, 7353–7362. [Google Scholar] [CrossRef] [Green Version]
  16. Azimirad, R.; Khademi, A.; Akhavan, J.; Moshfegh, A.Z. Growth of Na0.3WO3 nanorods for the field emission application. J. Phys. D Appl. Phys. 2009, 42, 205405. [Google Scholar] [CrossRef]
  17. Wu, P.J.; Brahma, S.; Lu, H.H.; Huang, J.L. Synthesis of cesium tungsten bronze by a solution-based chemical route and the NIR shielding properties of cesium tungsten bronze thin films. Appl. Phys. A 2020, 126, 98. [Google Scholar] [CrossRef]
  18. Kaliev, K.A.; Baraboshkin, A.N. Electrocrystallization of transition metal oxide bronzes from molten salts. In Oxide Bronzes; Spitsyn, V.I., Ed.; Nauka: Moscow, Russia, 1982; pp. 137–175. [Google Scholar]
  19. Vakarin, S.V. Oriented Growth of Tungsten Bronzes at the Melts Electrolysis; UB RAS: Ekaterinburg, Russia, 2005. (In Russian) [Google Scholar]
  20. Kosov, A.V.; Semerikova, O.L.; Vakarin, S.V.; Grishenkova, O.V.; Trofimov, A.A.; Leonova, A.M.; Leonova, N.M.; Zaikov, Y.P. Photovoltaic response of silicon wafers treated in the K2WO4-Na2WO4-WO3 melt. J. Electrochem. Soc. 2021, 168, 126503. [Google Scholar] [CrossRef]
  21. Kosov, A.V.; Semerikova, O.L.; Vakarin, S.V.; Zaykov, Y.P. Electrochemical behavior of the nickel/oxide tungsten bronze system at cyclic potential sweep. Russ. Metall. 2019, 2019, 803–808. [Google Scholar] [CrossRef]
  22. Kosov, A.V.; Semerikova, O.L.; Vakarin, S.V.; Pankratov, A.A.; Plaksin, S.V.; Korzun, I.V.; Akashev, L.A.; Zaykov, Y.P. Electrochemical synthesis of tetragonal oxide tungsten bronze nanofilms on platinum. Russ. Metall. 2017, 2017, 152–157. [Google Scholar] [CrossRef]
  23. Kosov, A.V.; Semerikova, O.L.; Vakarin, S.V.; Pankratov, A.A.; Plaksin, S.V.; Zaykov, Y.P. Formation of nanocrystalline tetragonal oxide tungsten bronzes on platinum. Russ. Metall. 2017, 2017, 158–162. [Google Scholar] [CrossRef]
  24. Vakarin, S.V.; Semerikova, O.L.; Kosov, A.B.; Pankratov, A.A.; Plaksin, S.V.; Korzun, I.V.; Akashev, L.A.; Zaikov, Y.P. Electrochemical deposition of nanocrystalline tungsten bronze films on platinum. Int. J. Adv. Res. 2015, 3, 691–700. [Google Scholar]
  25. Shurdumov, B.K. On the mechanism of the formation of oxide tungsten bronzes of alkali metals during the electrolysis of high-viscosity melts of tungstate phosphate systems. Izvestiya Vuzov Khimiya Khimicheskaya Tekhnologiya 2001, 44, 152–155. (In Russian) [Google Scholar]
  26. Khubolov, B.M. Electrocrystallization of thin films of sodium—Tungsten bronzes. In Physical and Chemical Aspects of the Study of Clusters, Nanostructures and Nanomaterials; Tver State University: Tver, Russia, 2020; pp. 213–221. [Google Scholar]
  27. Khubolov, B.M. Electrochemical synthesis of oxide tungsten bronze powders. In Physical and Chemical Aspects of the Study of Clusters, Nanostructures and Nanomaterials; Tver State University: Tver, Russia, 2020; pp. 738–745. [Google Scholar]
  28. Drobasheva, T.I.; Snezhkov, V.I. Electrocrystallization and properties of alkali-metal tungsten and molybdenum bronzes. Inorg. Mater. 1998, 34, 1162–1165. [Google Scholar]
  29. Okada, K.; Miyake, M.; Iwai, S.; Ohno, H.; Furukawa, K. Structural analysis of molten Na2WO4. J. Chem. Soc. Faraday Trans. 2 1978, 74, 1141–1148. [Google Scholar] [CrossRef]
  30. Miyake, M.; Okada, K.; Iwai, S.; Ohno, H.; Furukawa, K. Structural analysis of molten Na2W2O7. J. Chem. Soc. Faraday Trans. 2 1978, 74, 1880–1884. [Google Scholar] [CrossRef]
  31. Voron’ko, Y.K.; Sobol’, A.A. Influence of cations on the vibrational spectra and structure of [WO4] complexes in molten tungstates. Inorg. Mater. 2005, 41, 420–428. [Google Scholar] [CrossRef]
  32. Voronko, Y.K.; Sobol, A.A.; Shukshin, V.E. Raman scattering study of molten alkali-metal molybdates and tungstates rich in basic oxides. Inorg. Mater. 2014, 50, 837–843. [Google Scholar] [CrossRef]
  33. Wang, J.; You, J.; Wang, M.; Lu, L.; Wan, S.; Sobol, A.A. In-situ studies on the micro-structure evolution of A2W2O7 (A = Li, Na, K) during melting by high temperature Raman spectroscopy and density functional theory. Spectrochim. Acta A 2017, 185, 188–196. [Google Scholar] [CrossRef]
  34. Wang, J.; You, J.L.; Sobol, A.A.; Lu, L.M.; Wang, M.; Wu, J.; Lv, X.M.; Wan, S.M. In-situ high temperature Raman spectroscopic study on the structural evolution of Na2W2O7 from the crystalline to molten states. J. Raman Spectrosc. 2017, 48, 298–304. [Google Scholar] [CrossRef]
  35. Wang, J.; You, J.L.; Wang, Y.Y.; Wang, M.; Jun, W.U. In-situ Raman spectroscopic study of the molten tungstates in Li2O-WO3 binary system. Chin. J. Light Scatt. 2016, 28, 149–152. [Google Scholar]
  36. Gunzi, K.; Kohsaka, S.; Yokokawa, T.E.m.f. measurements of molten oxide mixtures IV. Sodium oxide + molybdenum trioxide. J. Chem. Thermodyn. 1979, 11, 553–557. [Google Scholar] [CrossRef]
  37. Khvatov, A.Y.; Baraboshkin, A.N.; Tarasova, K.P. Study of the ionic composition of a tungstate melt by the EMF method. Elektrokhimiya 1985, 21, 1657–1660. (In Russian) [Google Scholar]
  38. Afonichkin, V.K.; Leontiev, V.N.; Komarov, V.E. Equilibrium electrode potentials of tungsten in melts of the Na2WO4–WO3 system. Elektrokhimiya 1993, 29, 341–347. (In Russian) [Google Scholar]
  39. Smirnov, M.; Stepanov, V. Density and surface tension of molten alkali halides and their binary mixtures. Electrochim. Acta 1982, 27, 1551–1563. [Google Scholar] [CrossRef]
  40. Khudorozhkova, A.O.; Isakov, A.V.; Kataev, A.A.; Red’kin, A.A.; Zaikov, Y.P. Density of KF–KCl–KI melts. Russ. Metall. 2020, 2020, 918–924. [Google Scholar] [CrossRef]
  41. Soler, J.M.; Artacho, E.; Gale, J.D.; Garcia, A.; Junquera, J.; Ordejon, P.; Sanchez-Portal, D. The SIESTA method for ab initio order-N materials simulation. J. Phys. Condens. Matter 2002, 14, 2745–2779. [Google Scholar] [CrossRef]
  42. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  43. Nose, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef]
  44. Bychin, V.P.; Antonov, A.A. Some physical and chemical properties of the Na2WO4–WO3. In Chemistry and Technology of Molybdenum and Tungsten. In Proceedings of the III All-Union Meeting, North Caucasus Mining and Metallurgical Institute, Ordzhonikidze, USSR, 7 October 1977; pp. 63–64. (In Russian). [Google Scholar]
Figure 1. Schemes of experimental setups (a) for determining the effect of oxygen partial pressure on the potential of the Pt electrode and (b) for measuring the dependence of ΔE on the mole fraction of WO3 in the melt. Designations: 1—shaft furnace, 2—quartz protective vessel, 3—platinum crucible with the 0.8Na2WO4–0.2WO3 melt (standard melt), 4—Pt electrodes, 5—quartz test tube, 6—alumina tubes, 7—Pt/Pt-Rh thermocouple, 8—stainless steel lid filled with kaolin wool, 9—thin-walled alumina crucible with the melt under study.
Figure 1. Schemes of experimental setups (a) for determining the effect of oxygen partial pressure on the potential of the Pt electrode and (b) for measuring the dependence of ΔE on the mole fraction of WO3 in the melt. Designations: 1—shaft furnace, 2—quartz protective vessel, 3—platinum crucible with the 0.8Na2WO4–0.2WO3 melt (standard melt), 4—Pt electrodes, 5—quartz test tube, 6—alumina tubes, 7—Pt/Pt-Rh thermocouple, 8—stainless steel lid filled with kaolin wool, 9—thin-walled alumina crucible with the melt under study.
Processes 10 02658 g001
Figure 2. (a) Experimental time dependence of the potential difference between two platinum–oxygen electrodes with a stepwise change in the volume fraction of air, Wair, above one of them. The Wair (vol%) values are shown in the figure. T = 983 K. Other details of the experiment are described in Section 2.2.1. (b) Dependence of the potential difference on the logarithm of the volume fraction of air.
Figure 2. (a) Experimental time dependence of the potential difference between two platinum–oxygen electrodes with a stepwise change in the volume fraction of air, Wair, above one of them. The Wair (vol%) values are shown in the figure. T = 983 K. Other details of the experiment are described in Section 2.2.1. (b) Dependence of the potential difference on the logarithm of the volume fraction of air.
Processes 10 02658 g002
Figure 3. Experimental dependences of the potential difference in Cell (4) on the WO3 mole fraction in the (1 − ν)Na2WO4–νWO3 melt at 983, 1023 and 1073 K.
Figure 3. Experimental dependences of the potential difference in Cell (4) on the WO3 mole fraction in the (1 − ν)Na2WO4–νWO3 melt at 983, 1023 and 1073 K.
Processes 10 02658 g003
Figure 4. (a) Experimental temperature dependences of the density of the (1 − ν)Na2WO4–νWO3 melts. The WO3 mole fraction values (ν) are shown in the figure. Our results are compared with the calculation by Equation (12) [44] at a ν equal to 0 (□), 0.1(△) and 0.2(○). (b) Dependences of density on the WO3 mole fraction at 983, 1023 and 1073 K.
Figure 4. (a) Experimental temperature dependences of the density of the (1 − ν)Na2WO4–νWO3 melts. The WO3 mole fraction values (ν) are shown in the figure. Our results are compared with the calculation by Equation (12) [44] at a ν equal to 0 (□), 0.1(△) and 0.2(○). (b) Dependences of density on the WO3 mole fraction at 983, 1023 and 1073 K.
Processes 10 02658 g004
Figure 5. (a) Geometric structures of W x O 3 x + 1 2 and W x O 3 x 1 2 + compounds after geometric optimization and ab initio molecular dynamics with a Nose–Hoover thermostat at 983 K. (b) Geometric structure of [W4O13]2− at 1073 K. Bond lengths are given in angstroms.
Figure 5. (a) Geometric structures of W x O 3 x + 1 2 and W x O 3 x 1 2 + compounds after geometric optimization and ab initio molecular dynamics with a Nose–Hoover thermostat at 983 K. (b) Geometric structure of [W4O13]2− at 1073 K. Bond lengths are given in angstroms.
Processes 10 02658 g005
Figure 6. Dependences of W–O bond energy on the number of W atoms in the compound at (a) 983 K and (b) 1073 K.
Figure 6. Dependences of W–O bond energy on the number of W atoms in the compound at (a) 983 K and (b) 1073 K.
Processes 10 02658 g006
Figure 7. Experimental (symbols) and fitted (lines) dependences of the redox potential on the mole fraction of WO3 in the Na2WO4–WO3 melt. The curves have been shifted relative to each other to show them in one figure (the actual relative position of the experimental points is shown in Figure 3). The points that correspond to the standard melt 0.8Na2WO4–0.2WO3E(ν = 0.2) = 0) are marked with blank symbols.
Figure 7. Experimental (symbols) and fitted (lines) dependences of the redox potential on the mole fraction of WO3 in the Na2WO4–WO3 melt. The curves have been shifted relative to each other to show them in one figure (the actual relative position of the experimental points is shown in Figure 3). The points that correspond to the standard melt 0.8Na2WO4–0.2WO3E(ν = 0.2) = 0) are marked with blank symbols.
Processes 10 02658 g007
Figure 8. Dependences of the concentrations of tungsten-containing ions and O2− ions on the mole fraction of tungsten trioxide in Na2WO4–WO3 melts at 983, 1023 and 1073 K.
Figure 8. Dependences of the concentrations of tungsten-containing ions and O2− ions on the mole fraction of tungsten trioxide in Na2WO4–WO3 melts at 983, 1023 and 1073 K.
Processes 10 02658 g008
Table 1. Equilibrium constants of Reactions (23)–(26).
Table 1. Equilibrium constants of Reactions (23)–(26).
T, Kk0, mol2cm−6k1k2k3
9839.82 × 10−338.40 × 1033.05 × 1011.77 × 10−13
10231.09 × 10−284.32 × 1021.653.45 × 10−11
10734.64 × 10−241.44 × 1015.82 × 10−21.44 × 10−8
Table 2. Calculated values of ionic concentrations in Na2WO4–WO3 melts.
Table 2. Calculated values of ionic concentrations in Na2WO4–WO3 melts.
WO3 Mole Fraction (ν)Concentrations of Ions (cm−3)
WO 2 2 + WO 4 2 W 2 O 7 2 W 3 O 10 2 W 4 O 13 2 O2−Na+
983 K
02.57 × 1067.82 × 10213.29 × 10151.65 × 1052.71 × 10−73.29 × 10151.56 × 1022
0.12.72 × 10176.52 × 10218.15 × 10201.21 × 10165.90 × 1099.25 × 1091.47 × 1022
0.22.42 × 10185.17 × 10211.72 × 10216.77 × 10168.77 × 10102.76 × 1091.38 × 1022
0.31.64 × 10193.70 × 10212.70 × 10212.35 × 10176.71 × 10118.96 × 1081.28 × 1022
0.41.22 × 10202.30 × 10213.62 × 10216.79 × 10174.17 × 10122.59 × 1081.16 × 1022
1023 K
01.47 × 1087.76 × 10214.55 × 10166.19 × 1085.114.55 × 10161.55 × 1022
0.18.02 × 10166.46 × 10218.06 × 10202.33 × 10174.10 × 10131.78 × 10121.45 × 1022
0.27.19 × 10175.10 × 10211.69 × 10211.30 × 10186.10 × 10145.29 × 10111.36 × 1022
0.35.04 × 10183.61 × 10212.67 × 10214.59 × 10184.78 × 10151.68 × 10111.26 × 1022
0.45.06 × 10192.06 × 10213.66 × 10211.51 × 10193.76 × 10164.01 × 10101.14 × 1022
0.53.73 × 10201.16 × 10214.17 × 10213.48 × 10191.77 × 10171.1 × 10109.96 × 1021
1073 K
01.52 × 10107.68 × 10219.21 × 10177.68 × 10121.10 × 1099.21 × 10171.54 × 1022
0.11.90 × 10166.37 × 10217.80 × 10206.62 × 10189.65 × 10177.51 × 10141.43 × 1022
0.21.53 × 10175.07 × 10211.57 × 10213.36 × 10191.24 × 10192.36 × 10141.34 × 1022
0.37.77 × 10173.76 × 10212.25 × 10219.38 × 10196.72 × 10199.02 × 10131.23 × 1022
0.43.46 × 10182.53 × 10212.63 × 10211.90 × 10202.35 × 10203.51 × 10131.12 × 1022
0.51.51 × 10191.50 × 10212.52 × 10212.92 × 10205.84 × 10201.30 × 10139.77 × 1021
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kosov, A.V.; Semerikova, O.L.; Vakarin, S.V.; Grishenkova, O.V.; Vorob’ev, A.S.; Khudorozhkova, A.O.; Zaikov, Y.P. Ionic Equilibria in Polytungstate Melts. Processes 2022, 10, 2658. https://doi.org/10.3390/pr10122658

AMA Style

Kosov AV, Semerikova OL, Vakarin SV, Grishenkova OV, Vorob’ev AS, Khudorozhkova AO, Zaikov YP. Ionic Equilibria in Polytungstate Melts. Processes. 2022; 10(12):2658. https://doi.org/10.3390/pr10122658

Chicago/Turabian Style

Kosov, Alexander V., Olga L. Semerikova, Sergey V. Vakarin, Olga V. Grishenkova, Alexey S. Vorob’ev, Anastasia O. Khudorozhkova, and Yury P. Zaikov. 2022. "Ionic Equilibria in Polytungstate Melts" Processes 10, no. 12: 2658. https://doi.org/10.3390/pr10122658

APA Style

Kosov, A. V., Semerikova, O. L., Vakarin, S. V., Grishenkova, O. V., Vorob’ev, A. S., Khudorozhkova, A. O., & Zaikov, Y. P. (2022). Ionic Equilibria in Polytungstate Melts. Processes, 10(12), 2658. https://doi.org/10.3390/pr10122658

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop