Open access peer-reviewed chapter

A Distinguishing Feature of the Balance 2∙f(O)−f(H) in Electrolytic Systems: The Reference to Titrimetric Methods of Analysis

Written By

Anna Maria Michałowska-Kaczmarczyk, Aneta Spórna-Kucab and Tadeusz Michałowski

Submitted: 20 March 2017 Reviewed: 14 April 2017 Published: 27 September 2017

DOI: 10.5772/intechopen.69249

From the Edited Volume

Advances in Titration Techniques

Edited by Vu Dang Hoang

Chapter metrics overview

1,899 Chapter Downloads

View Full Metrics

Abstract

The balance 2∙f(O)−f(H) provides a general criterion distinguishing between electrolytic redox and non-redox systems of any degree of complexity, in aqueous, non-aqueous and mixed-solvent media. When referred to redox systems, it is an equation linearly independent on charge (ChB) and elemental/core balances f(Yg) for elements/cores Yg ≠ H and O, whereas for non-redox systems, 2∙f(O)−f(H) is linearly dependent on these balances. The balance 2∙f(O)−f(H) formulated for redox systems is the primary form (pr-GEB) of the generalized electron balance (GEB) as the fundamental equation needed for resolution of these systems. Formulation of GEB for redox systems needs no prior knowledge of oxidation numbers for all elements of the system. Any prior knowledge of oxidation numbers for all elements in components forming a redox system and in the species of the system thus formed is not necessary within the Approach II to GEB. Oxidants and reductants are not indicated. Stoichiometry and equivalent mass are redundant concepts only. The GEB, together with charge balance and concentration balances for elements ≠ H and O, and the complete set of independent equations for equilibrium constants form an algorithm, resolvable with use of an iterative computer program. All attainable physicochemical knowledge can be included in the algorithm. Some variations involved with tests of possible reaction paths for metastable systems can also be made. The effects of incomplete physicochemical knowledge on the system can be also tested. One of the main purposes of this chapter is to provide the GEB formulation needed for resolution of redox systems and familiarize it to a wider community of chemists.

Keywords

  • electrolytes
  • redox systems
  • non-redox systems
  • generalized electron balance (GEB)
  • titration

Motto: ‘Nothing is too wonderful to be true if it be consistent with the laws of Nature’

(M. Faraday)

Advertisement

1. Introduction

Scientific theories describe particular units and rules governing the relationships between them. This description is considered as interpretation of reality. In particular, the thermodynamic description of any electrolytic system according to generalized approach to electrolytic systems (GATESs) [1] is based on fundamental, physical rules of conservation, expressed by charge balance (ChB) and elemental balances, for particular elements and/or cores in a closed system, separated from the environment by diathermal (freely permeable by heat) walls. The term ‘core’ is related to a cluster of elements with the same formula, structure and external charge. For example, HSO4−1∙n5H2O, SO4−2∙n6H2O and FeSO4∙n27H2O, in Eq. (3) (below), have the common core (SO4−2). In this context, the pairs of species: (i) C2O4−2 and CO3−2; (ii) C2O4−2 (from maleic acid) and C2O4−2 (from fumaric acid) and (iii) NO2−1 and NO2 have no common cores.

Chemical interactions in electrolytic systems, for example protonation, neutralization, hydration, hydrolysis or dilution phenomena, are usually accompanied by exothermic or endothermic effects. However, the mass change, Δm, resulting from these thermal phenomena, estimated according to the Einstein’s formula ΔE = Δm.c2 and put in context with their enthalpies ΔH (ΔE = ΔH), that is Δm = ΔH/c2, is negligibly small (not measurable). Therefore, the mass of a chemical system remains practically unchanged, regardless of whether the chemical reactions take place in it or not.

The heat exchange between the system and its environment through diathermal walls enables the temperature T of the system to be kept constant during the appropriate dynamic processes, such as titration, performed in a quasistatic manner. Stability of temperature T within the titrant (titrating solution, T), titrand (titrated solution, D) and D + T mixture, together with constancy of ionic strength (I) in D+T, is the preliminary condition ensuring stability of the corresponding equilibrium constants, Ki = Ki(T, I), related to the system in question. The diathermal walls separate condensed (liquid or liquid + solid) phases from its environment.

An open chemical system is an approximation of the closed system—provided that the matter (e.g., H2O, CO2 and O2) exchange between the system and its surroundings can be neglected, within a relatively short period of time needed to carry out the process considered, for example, titration.

On the initial stage of ChB and elemental/core balance formulation, it is advisable to start the quantitative considerations from the numbers of particular entities (components, species):

  • N0j for jth component constituting the system

  • Ni for ith species in the system thus formed

For example, H2O and gaseous HCl, as components, form an aqueous solution of HCl, with H2O and hydrates of H+1, OH−1 and Cl−1 as the species. Generally, when solid, liquid and/or gaseous solutes are introduced into water, a mono- or two-phase system is obtained. The resulting mixture is limited to the condensed (liquid or liquid + solid) phases. We refer mainly to aqueous media (W = H2O), where the physicochemical knowledge is relatively extensive, incomparably better than one for the system with non-aqueous, or mixed-solvent media [26], with amphiprotic co-solvents involved, also considered in this chapter. For such media, the elemental f(Eg)) or core f(coreg)) balances written in terms of numbers of individual entities containing the elements (Eg), or cores (coreg), are formulated.

In aqueous electrolytic systems, different entities Xizi exist as hydrated species, XiziniH2O; ni = niW = niH2O is the mean number of water (W = H2O) molecules attached to Xizi(ni0), zi is a charge of this species, expressed in elementary charge units and e = F/NA (F = Faraday constant, NA = Avogadro number). For these species present in static or dynamic systems, we apply the notation

Xizi(Ni,ni)E1

where Ni is the number of these entities (individual species). On this basis, the numbers of particular elements in these species are calculated; for example, in the solution II (see below), N04 molecules of FeSO4∙7H2O contain 14 N04 atoms of H, 11 N04 atoms of O and N11 atoms of Fe; N5 ions of HSO4−1∙n5H2O (N5,n5) in the set (2) of species specified below contain N5(1 + 2n5) atoms of H, N5(4 + n5) atoms of O and N5 atoms of S.

In further parts of this chapter, the terms linear combination and linear dependency/independency of equations are introduced. These terms, well known from the elementary algebra course, will be applied to elemental/core balances, as a system of algebraic equations. The elemental balances f(H) for hydrogen (H) and f(O) for oxygen (O) and the linear combination 2∙f(O)−f(H) are formulated and then combined with charge balance (ChB) and other elemental/core balances f(Yg) for other elements (Yg = Eg) or cores (Yg = coreg), Yg ≠ H and O. This way, the general properties of 2∙f(O)−f(H) in non-redox and redox systems are distinguished, see Refs. [718] and earlier references cited therein.

The 2f(O)−f(H), charge balance and elemental/core balances will be expressed first in terms of the numbers of particular entities. Next, the related balances will be presented in terms of molar concentrations, to be fully compatible with expressions for equilibrium constants that are also presented in terms of molar concentrations of the related species.

Static and dynamic systems are distinguished. A static system is obtained after disposable mixing with the respective components. For illustrative purposes, we consider first four solutions, as static non-redox systems, formed from the following components:

  • (I) N01 molecules of KMnO4, N02 molecules of CO2 and N03 molecules of H2O in V1 mL of the resulting solution

  • (II) N04 molecules of FeSO4∙7H2O, N05 molecules of H2SO4, N06 molecules of CO2 and N07 molecules of H2O in V2 mL of the resulting solution

  • (III) N08 molecules of H2C2O4∙2H2O, N09 molecules of H2SO4, N010 molecules of CO2 and N011 molecules of H2O in V3 mL of the resulting solution

  • (IV) N012 molecules of FeSO4∙7H2O, N013 molecules of H2C2O4∙2H2O, N014 molecules of H2SO4, N015 molecules of CO2 and N016 molecules of H2O in V4 mL of the resulting solution.

The CO2 in the respective solutions is primarily considered as one originated from ambient air, on the step of preparation of these solutions.

From these static systems, we prepare later different dynamic systems: (I) ⇒ (II), (I) ⇒ (III) and (I) ⇒ (IV), where (I) as titrant T is added into (II), (III) or (IV) as titrand D (T ⇒ D), and the D + T mixtures containing different species are formed.

To avoid a redundancy resulting from application of different subscripts within (Ni, ni) ascribed to the same species XiziniH2O in different solutions (I)–(IV), we apply the common basis of the species from which the components will be selected to the respective balances. The set of the species is as follows:

H2O(N1), H+1(N2, n2), OH1(N3, n3), K+1(N4, n4), HSO41(N5, n5), SO42(N6, n6),H2C2O4(N7, n7),and HC2O41(N8, n8),C2O42(N9, n9), H2CO3(N10, n10),HCO31(N11, n11), CO32(N12, n12), MnO41(N13, n13),MnO42(N14, n14), Mn+3(N15, n15),MnOH+2(N16, n16), MnC2O4+1(N17, n17), Mn(C2O4)21(N18, n18), Mn(C2O4)33(N19, n19),Mn+2(N20, n20), MnOH+1(N21, n21), MnSO4(N22, n22), MnC2O4(N23, n23), Mn(C2O4)22(N24, n24), Fe+2(N25, n25), FeOH+1(N26, n26)andFeSO4(N27, n27)Fe(C2O4)22(N28, n28),Fe(C2O4)34(N29, n29), Fe+3(N30, n30), FeOH+2(N31, n31), Fe(OH)2+1(N32, n32),Fe2(OH)2+4(N33, n33), FeSO4+1(N34, n34), Fe(SO4)21(N35, n35), FeC2O4+1(N36, n36),Fe(C2O4)21(N37, n37), Fe(C2O4)33(N38, n38),andFeC2O4(N39, n39),MnC2O4(N40, n40)E2
Advertisement

2. A short note

Referring to pure algebra, let us consider the set of G + 1 algebraic equations: fg(x) = φg(x)bg = 0, where g = 0,1,…,G, xT = (x1,…,xI), transposed (T) vector x, composed of independent (scalar) variables xi (i ϵ <1, I>); agi, bg ϵ R are independent (explicitly) on x. After multiplying the equations by the numbers ωg ϵ R, and addition of the resulting equations, we get the linear combination g=0Gωgfg(x)=0g=0Gωgφg(x)=g=0Gωgbg of the basic equations.

Formation of linear combinations is applicable to check the linear dependence or independence of the balances. A very useful/effective manner for checking/stating the linear dependence of the balances is the transformation of an appropriate system of equations to the identity 0 = 0 [2, 23]. For this purpose, we will try, in all instances, to obtain the simplest form of the linear combination. To facilitate these operations, carried out by cancellation of the terms on the left and right sides of equations after multiplying and changing sides of these equations, we apply the equivalent forms of the starting equations fg(x) = 0:

fg(x):φg(x)bg=0φg(x)=bgfg(x):φg(x)=bgbg=φg(x).E3

In this notation, fg(x) will be essentially treated not as the algebraic expression on the left side of the equation fg(x) = 0 but as an equation that can be expressed in alternative forms presented above.

Advertisement

3. Combination of elemental/core balances for non-redox systems

For the solution (I), we have the balances:

f1=f(H):2N1+ N2(1+2n2)+N3(1+2n3)+2N4n4+N10(2+2n10)+N11(1+2n11)+2N12n12+2N13n13=2N03E406
f2=f(O):N1+N2n2+N3(1+n3)+N4n4+N10(3+n10)+N11(3+n11)+N12(3+n12)+N13(4+n13)=4N01+2N02+N03E407
f12=2f(O)f(H):N2+N3+4N10+5N11+6N12+8N13=8N01+ 4N02E408
f0=ChB:    N2N3+ N4N112N12N13=0f3=f(K):N01=N47f4=7f(Mn):7N01=7N134f5=4f(CO3):4N02=4N10+4N11+4N12E409
f12+f0f37f44f5:0=0E4010

For the solution (II), we have the balances:

f1=f(H):2N1+N2(1+2n2)+N3(1+2n3)+N5(1+2n5)+2N6n6+N10(2+2n10)+N11(1+2n11)+2N12n12+2N25n25+N26(1+2n26)+2N27n27=14N04+2N05+2N07E4011
f2=f(O):N1+N2n2+N3(1+n3)+N5(4+n5)+N6(4+n6)+N10(3+n10)+N11(3+n11)+N12(3+n12)+N25n25+N26(1+n26)+N27(4+n27)=11N04+4N05+2N06+N07E4012
f12= 2f(O)f(H):N2+N3+7N5+8N6+4N10+5N11+6N12+N26+8N27=8N04+6N05+4N06E4013
f0=ChB:                     N2N3N52N6N112N12+2N25+N26=02f3=2f(Fe):2N04=2N25+2N26+2N276f4=6f(SO4):6N04+6N05=6N5+6N6+6N274f5=4f(CO3):4N06=4N10+4N11+4N12E4014
f12+f02f36f44f5:0=0E4015

For the solution (III), we have the balances:

f1=f(H):2N1+N2(1+2n2)+N3(1+2n3)+N5(1+2n5)+2N6n6+N7(2+2n7)+N8(1+2n8)+2N9n9+N10(2+2n10)+N11(1+2n11)+2N12n12=6N08+2N09+2N011E40016
f2=f(O):N1+N2n2+N3(1+n3)+N5(4+n5)+N6(4+n6)+N7(4+n7)+N8(4+n8)+N9(4+n9)+N10(3+n10)+N11(3+n11)+N12(3+n12)=6N08+4N09+2N010+N011E4017
f12= 2f(O)f(H):N2+N3+7N5+8N6+6N7+7N8+8N9+4N10+5N11+6N12=6N08+6N09+4N010E4018
f0=ChB:N2N3N52N6N82N9N112N12=06f3=6f(SO4):6N09=6N5+6N64f4=4f(CO3):4N010=4N10+4N11+4N126f5=6f(C2O4):6N08=6N7+6N8+6N9E4019
f12+f06f34f46f5:0=0E4020

For the solution (IV), we have the balances:

f1=f(H):2N1+N2(1+2n2)+N3(1+2n3)+N5(1+2n5)+2N6n6+N7(2+2n7)+N8(1+2n8)+2N9n9+N10(2+2n10)+ N11(1+2n11)+2N12n12+2N25n25+N26(1+2n26)+2N27n27+2N28n28+2N29n29+2N39n39=14N012+6N013+2N014+2N016E4021
f2=f(O):N1+N2n2+N3(1+n3)+N5(4+n5)+N6(4+n6)+N7(4+n7)+N8(4+n8)+N9(4+n9)+N10(3+n10)+N11(3+n11)+N12(3+n12)+N25n25+N26(1+n26)+N27(4+n27)+N28(8+n28) + N29(12+n29)+N39(4+n39)=11N012+ 6N013+ 4N014+ 2N015+ N016E4022
f12= 2f(O)f(H):N2+N3+7N5+8N6+6N7+7N8+8N9+4N10+5N11+6N12+N26+8N27+16N28+24N29+8N39=8N012+6N013+6N014+4N015E4023
f0=ChB:N2N3N52N6N82N9N112N12+2N26+N272N294N30=06f3=6f(SO4):6N012+6N014=6N5+6N6+6N284f4=4f(CO3):4N015=4N10+4N11+4N126f5=6f(C2O4):6N013=6N7+6N8+6N9+12N28+18N29+6N392f6=2f(Fe):2N012=2N25+2N26+2N27+2N28+2N29+2N39E4024
f12+f06f34f46f52f6:0=0E4025

Summarizing, for all the solutions (I)–(IV), we obtained the identities 0 = 0:

(I)f12+f0f37f44f5=0;(II)f12+f02f36f44f5=0;(III)f12+f06f34f46f5=0;(IV)f12+f06f34f46f52f6=0 E4

All the solutions are non-redox systems. Except protonation/hydrolytic effects in (I)–(IV), the complexation and precipitation occur in (II) and (IV); the precipitation of FeC2O4 does not occur there at sufficiently high concentrations of H2SO4.

The solutions can be mixed according to titrimetric mode. In particular, we refer to the D + T systems obtained in the titrations T ⇒ D indicated above, namely (I) ⇒ (II), (I) ⇒ (III) and (I) ⇒ (IV). According to the notation applied elsewhere, for example, in Refs. [1923], V0 mL of D is titrated with volume V mL of T, added up to a given point of the titration, and the D + T mixture with volume V0 + V mL is formed at this point if the assumption of the volumes additivity is valid.

We assume V1 = V, CV = 103∙N01/NA (NA: Avogadro’s number) and V2 = V0 and C0V0 = 103∙N04/NA for (I) ⇒ (II); V3 = V0 and C0V0 = 103∙N08/NA for (I) ⇒ (III); V4 = V0 and C01V0 = 103∙N013/NA and C02V0 = 103∙N012/NA for (I) ⇒ (IV). Concentrations of the species XiziniH2O in the related systems are defined by relation [Xizi](V0+V)=Ni/NA, where [Xizi] is the molar concentration of XiziniH2O for i ≥ 2. The progress of the titration in (I) ⇒ (II) and (I) ⇒ (III) can be defined by the fraction titrated [2429] value

Φ=CVC0V0E5

whereas V will be taken as a parameter varied on abscissa in the graphical presentation of the system (I) ⇒ (IV).

Advertisement

4. Formulation of dynamic redox systems

The D and T, formed out of particular components, are considered as subsystems of the D + T system thus obtained. The titration is considered as a quasistatic process carried out under isothermal conditions and perceived both from physicochemical and analytical viewpoints.

Let us consider four starting solutions composed of:

  • N01 molecules of KMnO4, N02 molecules of CO2 and N03 molecules of H2O in V1 mL of the resulting solution

  • N04 molecules of FeSO4∙7H2O, N05 molecules of H2SO4, N06 molecules of CO2 and N07 molecules of H2O in V2 mL of the resulting solution

  • N08 molecules of H2C2O4∙2H2O, N09 molecules of H2SO4, N010 molecules of CO2 and N011 molecules of H2O in V3 mL of the resulting solution

  • N012 molecules of FeSO4∙7H2O, N013 molecules of H2C2O4°2H2O, N014 molecules of H2SO4, N015 molecules of CO2 and N016 molecules of H2O in V4 mL of the resulting solution

We start our considerations from the most complex dynamic system (I) ⇒ (IV), where V mL KMnO4 (C) + CO2 (C1) is added into V0 mL FeSO4 (C01) + H2C2O4 (C02) + H2SO4 (C03) + CO2 (C04) at the defined point of the titration. The less complex dynamic systems (I) ⇒ (II) and (I) ⇒ (III) will be considered later as a particular case of the system (I) ⇒ (IV).

4.1. Formulation of GEB for the system (I) ⇒ (IV)

Referring to the set of species in Eq. (2), we apply a1 = 1 if pr1 = FeC2O4 is the equilibrium solid phase (precipitate) in the system, and a2 = 1 if pr2 = MnC2O4 is the equilibrium solid phase in the system; otherwise, we have a1 = 0 and/or a2 = 0. The elemental/core balances and ChB, formulated on the basis of the set of the species Eq. (2), are as follows:

f1=f(H):2N1+N2(1+2n2)+N3(1+2n3)+2N4n4+N5(1+2n5)+2N6n6+N7(2+2n7)+N8(1+2n8)+2N9n9+N10(2+2n10)+N11(1+2n11)+2N12n12+2N13n13+2N14n14+2N15n15+N16(1+2n16)+2N17n17+2N18n18+2N19n19+2N20n20+N21(1+2n21)+2N22n22+2N23n23+2N24n24+2N25n25+N26(1+2n26)+2N27n27+2N28n28+2N29n29+2N30n30+N31(1+2n31)+N32(2+2n32)+N33(2+2n33)+2N34n34+2N35n35+2N36n36+2N37n37+2N38n38+2a1N39n39+2a2N40n40=2N03+14N012+6N013+2N014+2N016E4028
f2=f(O):N1+N2n2+N3(1+n3)+N4n4+N5(4+n5)+N6(4+n6)+N7(4+n7)+N8(4+n8)+N9(4+n9)+N10(3+n10)+N11(3+n11)+N12(3+n12)+N13(4+n13)+N14(4+n14)+N15n15+N16(1+n16)+N17(4+n17)+N18(8+n18)+N19(12+n19)+N20n20+N21(1+n21)+N22(4+n22+ N23(4+n23)+N24(8+n24)+N25n25+N26(1+n26)+N27(4+n27)+N28(8+n28)+N29(12+n29)+N30n30+N31(1+n31)+N32(2+n32)+N33(2+n33)+N34(4+n34)+N35(8+n35)+N36(4+n36)+N37(8+n37)+N38(12+n38)+a1N39(4+n39)+a2N40(4+n40)=4N01+2N02+N03+11N012+6N013+4N014+2N015+N016E4029
f12=2f(O)f(H)N2+N3+7N5+8N6+6N7+7N8+8N9+4N10+ 5N11+6N12+8N13+8N14+N16+8N17+16N18+24N19N21+8N22+8N23+16N24+N26+8N27+16N28+24N29+N31+2N32+2N33+8N34+16N35+8N36+16N37+24N38+8a1N39+8a2N40=8N01+4N02+8N012+6N013+6N014+4N015E6
f0=ChBN2N3+N4N52N6N82N9N112N12N132N14+3N15+2N16+N17N183N19+2N20+N212N24+2N25+N262N284N29+3N30+2N31+N32+4N33+N34N35+N36N373N38=0E7
f3=f(K):N01=N4E8
6f4=6f(S)=6f(SO4):6N012+6N014=6N5+6N6+6N22+6N27+6N34+12N35E9
4f5=4f(C):4N02+8N013+4N015=8N7+8N8+8N9+4N10+4N11+4N12+8N17+16N18+24N19+8N23+16N24+16N28+24N29+8N36+16N37+24N38+8a1N39+8a2N40E10
3f6=3f(Fe) :3N012=3N25+3N26+3N27+3N28+3N29+3N30+3N31+3N32+6N33+3N34+3N35+3N36+3N37+3N38+3a1N39E11
2f7=2f(Mn):2N01=2N13+2N14+2N15+2N16+2N17+2N18+2N19+2N20+2N21+2N22+2N23+2N24+2a2N40E12
f12+f0f36f44f53f62f7:5N13+4N14+N15+N16+N012+2N013=5N01+2N7+2N8+2N9+N17+3N18+5N19+2N23+4N24+N25+N26+N27+5N28+7N29+2N36+4N37+6N38+3a1N39+2a2N40E13
5[MnO41]+4[MnO42]+[Mn+3]+[MnOH+2](2([H2C2O41]+[HC2O41]+[C2O42])+[MnC2O4+1]+3[Mn(C2O4)21]+5[Mn(C2O4)33]+2[MnC2O4]+4[Mn(C2O4)22]+[Fe+2]+[FeOH+1]+[FeSO4]+5[Fe(C2O4)22]+7[Fe(C2O4)34]+2[FeC2O4+1]+4[Fe(C2O4)21]+6[Fe(C2O4)33]+3a1[FeC2O4]+2a2[MnC2O4])=5CV/(V0+V)(C01+2C02)V0/(V0+V)2([H2C2O4]+[HC2O41]+[C2O42])+[MnC2O4+1]+3[Mn(C2O4)21]+5[Mn(C2O4)33]+2[MnC2O4]+4[Mn(C2O4)22]+[Fe+2]+[FeOH+1]+[FeSO4]+5[Fe(C2O4)22]+7[Fe(C2O4)34]+2[FeC2O4+1]+4[Fe(C2O4)21]+6[Fe(C2O4)33]+3a1[FeC2O4]+2a2[MnC2O4](5[MnO41]+4[MnO42]+[Mn+3]+[MnOH+2])=((C01+2C02)V05CV)/(V0+V)E14

Eq. (14) is the shortest/simplest form of GEB for the related system; it is, of course, different from the identity 0 = 0. On the basis of Eq. (14), one can also formulate the GEB for the system (I) ⇒ (II):

(KMnO4(C)+CO2(C1),V)(FeSO4(C01)+H2SO4(C03)+CO2(C04),V0)E15

and for the system (I) ⇒ (III):

(KMnO4(C)+CO2(C1),V)(H2C2O4(C02)+H2SO4(C03)+CO2(C04),V0)E16

Assuming C02 = 0, after omission of the related species involved with oxalates, from Eq. (14), we have the GEB valid for the system (I) ⇒ (II):

[Fe+2]+[FeOH+1]+[FeSO4](5[MnO41]+4[MnO42]+[Mn+3]+[MnOH+2])=(C01V05CV)(V0+V)E17

Assuming C01 = 0, after omission of the related Fe-species, from Eq. (14), we have the GEB valid for the system (I) ⇒ (III)

2([H202O4]+[HC2O41]+[C2O42])+[MnC2O4+1]+3[Mn(C2O4)21]+5[Mn(C2O4)33]+2[MnC2O4]+4[Mn(C2O4)22]+2a2[MnC2O4](5[MnO41]+4[MnO42]+[Mn+3]+[MnOH+2])=(2C02V05CV)(V0+V)E18

On the other hand, Eqs. (17) and Eqs. (18) are the simplest/shortest linear combinations for the related subsystems of the system (I) ⇒ (IV); both are also different from the identity, of course. For comparison, the linear combination f12 + f0f3 – 6f4 – 4f5 – 2f6 – 2f7 = 0, that is where −2f6 is put for −3f6, gives a more extended equation, where more components are involved. Anyway, we get here the equation, not the identity 0 = 0. It must be stressed that none of the linear combinations of these equations gives the identity. This is the general property of all redox systems, of any degree of complexity.

Advertisement

5. Confirmation of linear dependency of balances for non-redox systems

It can be stated that 2f(O)−f(H) is a linear combination of charge and elemental/core balances for non-redox systems of any degree of complexity; this means that 2f(O)−f(H) is not a new, independent balance in non-redox systems. From Eq. (4) we see, that for non-redox systems, f12 can be expressed as the linear combination of other balances of the system considered:

(I)f12=f3+7f4+4f5f0;(II)f12=2f3+6f4+4f5f0;(III)f12=6f3+4f4+6f5f0;(IV)f12=6f3+4f4+6f5+2f6f0E19

Eq. (4) can be also rewritten into equivalent forms:

(I). (+1)f1+(2)f2+(+1)f3+(+7)f4+(+4)f5f0=0;(II)(+1)f1+(2)f2+(+2)f3+(+6)f4+(+4)f5f0=0;(III)(+1)f1+(2)f2+2(+3)f3+(+4)f4+(+6)f5f0=0;(IV)(+1)f1+(2)f2+(+2)f6+2(+3)f3+(+4)f4+(+6)f5f0=0E20

As we see, the coefficient at the corresponding elemental/core balance in the related sum is equal to the oxidation number of the corresponding element. The linear dependence will be thus ascertained by multiplying the elemental/core balances by the appropriate oxidation numbers. After consecutive addition of the resulting balances to the sum of 2f(O)−f(H) and charge balance, followed by simplifications, the resulting sum is reduced to the identity 0 = 0. It is the simplest way of checking the linear dependency of the equations related to non-redox systems.

For redox systems, the appropriate linear combination of 2f(O)−f(H) with charge balance and elemental/core balances related to electron-non-active elements in the system in question leads to the simplest form of GEB named as generalized electron balance (GEB). This means that the GEB is a new balance, complementary/compatible with other (charge and elemental/core) balances related to the system in question.

Advertisement

6. Confirmation of linear independency of balances for redox systems

Applying a similar procedure, one can also state that 2f(O)−f(H) is not a linear combination of charge and elemental/core balances for redox systems of any degree of complexity; this means that 2f(O)−f(H) is a new/independent balance in redox systems.

The independency/dependency property of the balance 2f(O)−f(H) is the basis for the division of electrolytic systems into redox and non-redox systems [8, 9]. This rule is illustrated by the following examples, related to static and dynamic systems.

Advertisement

7. Confirmation of equivalency of approaches I and II to GEB for the system (I) ⇒ (IV)

We apply now the linear combination (algebraic sum) of Eqs. (69) for ChB and elemental/core balances, involving electron-non-active elements: H, O, K and S, perceived in terms of the Approach I to GEB as ‘fans’, we have:

f12+f0f36f4:6N7+6N8+6N9+4N10+4N11+4N12+7N13+6N14+3N15+3N16+9N17+15N18+21N19+2N20+2N21+2N22+8N23+14N24+2N25+2N26+2N27+14N28+20N29+3N30+3N31+3N32+6N33+3N34+3N35+9N36+15N37+21N38+8a1N39+8a2N40=7N01+6N013+4(N02+N015)6(N7+N8+N9)+4(N10+N11+N12)+7N13+6N14+3N15+3N16+9N17+15N18+21N19+2(N20+N21+N22)+8N23+14N24+2(N25+N26+N27)+14N28+20N29+3(N30+N31+N32+2N33+N34+N35)+9N36+15N37+21N38+8a1N39+8a2N40=7N01+6N013+4(N02+N015)E21

Denoting by ZC (= 6), ZFe (= 26) and ZMn (= 25), the atomic numbers for electron-active elements (‘players’) C, Fe and Mn, from Eqs. (10)Eqs. (12) and Eqs. (21), we have, by turns,

ZCf5+ZFef6+ZMnf7(f12+f0f36f4)(2ZC6)(N7+N8+N9)+(ZC4)(N10+N11+N12)+(ZMn7)N13+(ZMn6)N14+(ZMn3)(N15+N16)+N17(ZMn+2ZC9)+N18(ZMn+4ZC15)+N19(ZMn+6ZC21)+(ZMn2)(N20+N21+N22)+N23(ZMn+2ZC8)+N24(ZMn+4ZC14)+(ZFe2)(N25+N26+N27)(ZFe+4ZC12)N28+(ZFe+6ZC20)N29+(ZFe3)(N30+N31+N32)+2(ZFe3)N33+(ZFe3)(N34+N35)(ZFe+2ZC9)N36+(ZFe+4ZC15)N37+(ZFe+6ZC21)N38+a1(ZFe+2ZC8)N39+a2(ZMn+2ZC8)N40= (ZMn7)N01+(ZFe2)N012+2(ZC3)N013+(ZC4)(N02+N015)2(ZC3)(N7+N8+N9)+(ZC4)(N10+N11+N12)+(ZMn7)N13+(ZMn6)N14+(ZMn3)(N15+N16)+N17(ZMn+2ZC9)+N18(ZMn+4ZC15)+N19(ZMn+6ZC21)+(ZMn2)(N20+N21+N22)+N23(ZMn+2ZC8)+N24(ZMn+4ZC14)+(ZFe2)(N25+N26+N27)(ZFe+4ZC12)N28+(ZFe+6ZC20)N29+(ZFe3)(N30+N31+N32)+2(ZFe3)N33+(ZFe3)(N34+N35)(ZFe+2ZC9)N36+(ZFe+4ZC15)N37+(ZFe+6ZC21)N38+a1(ZFe+2ZC8)N39+a2(ZMn+2ZC8)N40=(ZMn7)N01+(ZFe2)N012+2(ZC3)N013+(ZC4)(N02+N015)2(ZC3)([H2C2O4]+[HC2O41]+[C2O42])+(ZC4)([H2CO3]+[HCO31]+[CO32])+(ZFe2)([Fe+2]+[FeOH+1]+[FeSO4])+(ZFe3)([Fe+3]+[FeOH+2]+[Fe(OH)2+1]+2[Fe2(OH)2+1]+[FeSO4+1]+[Fe(SO4)21])+(ZMn7)([MnO41]+(ZMn6)[MnO42]+(ZMn3)([Mn+3]+[MnOH+2])+(ZMn2)([Mn+2]+[MnOH+1]+[MnSO4])+(ZFe2+4(ZC3))[Fe(C2O4)22]+(ZFe2+6(ZC3))[Fe(C2O4)34]+a1(ZFe2+2(ZC3))[FeC2O4]+(ZFe3+2(ZC3))[FeC2O4+1]+(ZFe3+4(ZC3))[Fe(C2O4)21]+(ZFe3+6(ZC3))[Fe(C2O4)33]+(ZMn2+2(ZC3))[MnC2O4]+(ZMn2+4(ZC3))[Mn(C2O4)22]+(ZMn3+2(ZC3))[MnC2O4+1]+(ZMn3+4(ZC3))[Mn(C2O4)21]+(ZMn3+6(ZC3))[Mn(C2O4)33]+a2(ZMn2+2(ZC3))[MnC2O4]=2(ZC3)C01V0/(V0+V)+(ZFe2)C02V0/(V0+V)+(ZMn7)CV/(V0+V)E22

Eq. (22) is obtainable immediately according to the Approach I to GEB [1922]. Note, for example, that:

N17(ZMn+2ZC9)=N17(ZMn3+1°2(ZC3))N18(ZMn+4ZC15)=N18(ZMn3+2°2(ZC3))N19(ZMn+6ZC21)=N19(ZMn3+3°2(ZC3)E47

The equation for GEB thus obtained (according to the Approach II to GEB [1, 4, 718, 24, 25]) is then equivalent to GEB, obtained according to the Approach I to GEB, based on the principle of the common pool of electrons introduced by electron-active elements (‘players’) of the system in question. For redox systems, the GEB is the inherent part of the generalized approach to electrolytic systems (GATES) [1], denoted as GATES/GEB.

Advertisement

8. Some generalizing remarks on GEB

The linear combination 2f(O)−f(H) of elemental balances, f(H) for H and f(O) for O, is a keystone of the overall, thermodynamic knowledge on electrolytic systems. The 2f(O) −f(H) can be formulated both for non-redox and redox systems, with amphiprotic (co)solvent(s) involved. It is the basis for the Generalized Electron Balance (GEB) formulated according to the Approach II to GEB.

The principle of GEB formulation was presented for the first time in Refs. [30, 31] and then in Refs. [1922, 3136] as the Approach I to GEB. The GEB formulation according to the Approach I is based on the ‘card game’ principle, with a common pool of electrons as money, electron-active elements as players and electron-non-active elements as fans—not changing their oxidation degree, that is the fans’ accounts are intact in this convention [13, 23], see an illustration below. Electrons are considered as money, transferred between players; the knowledge of oxidation numbers of all elements in the system in question is needed there.

The Approach I to GEB, named also as the ‘short’ version of GEB, needs a knowledge of oxidation numbers for all elements in the species participating in the system that is considered. The equivalency of the Approaches I and II means that the equation obtained by a suitable linear combination of pr-GEB with charge balance and other elemental/core balances becomes identical with the one obtained directly from the Approach I to GEB.

Although derivation of GEB according to the Approach II is more extensive/laborious, it enables to formulate this balance without prior knowledge of oxidation numbers for the elements involved in the system. It is the paramount advantage of the Approach II to GEB, particularly when applied to more complex organic species, with radicals and ion-radicals involved. Moreover, within the Approaches I and II, the roles of oxidants and reducers are not ascribed a priori to particular components forming the redox system and the species formed in this system.

Ultimately, GEB, charge and elemental/core balances are expressed in terms of molar concentrations—to be fully compatible with expressions for equilibrium constants, interrelating molar concentrations of defined species on the basis of the mass action law applied to the correctly written reaction equation. The law of mass action is the one and only chemical law applied in GATES.

GEB is perceived as the law of matter conservation, as the general law of nature related to electrolytic (aqueous, non-aqueous or mixed-solvent media) redox systems and as a synthesis of physical and chemical laws [1, 14, 15, 23, 24, 27]. This law can also be extended on the systems with mixed (e.g., binary) solvents with amphiprotic (protophilic and protogenic) and aprotic properties. GEB is a rather unexpected consequence of the concentration balances for H and O and therefore the formulation of GEB, especially as the Approach II to GEB, is regarded as the scientific discovery and not as a confirmation of the obvious fact arising from other, fundamental laws of nature. This fact is emphasized in this chapter in the context of philosophical understanding of the scientific discoveries in the aspect of the laws of nature.

The GEB, together with charge and concentration/core balances and a set of independent equilibrium constants, provides a complete set of equations used for a thermodynamic description of a redox system taken for quantitative considerations within GATES/GEB ∈ GATES.

The roles of oxidants and reductants are not ascribed a priori to particular components forming the redox system and to the species formed in this system. In other words, full ‘democracy’ is established a priori within GATES/GEB.

The Approach II to GEB shows that the equivalent equations for GEB are derived from the common root of the elements conservation and then GEB is fully compatible with charge and concentration balances like ‘the lotus flower, lotus leaf and lotus seed come from the same root’ [13]. This compatibility is directly visible from the viewpoint of the Approach II to GEB. The GEB, based on a reliable law of the matter conservation, is equally robust as equations for charge and concentration balances. The complementarity of the GEB (Approaches I and II) to other balances is regarded as the expression of harmony of nature, and GATES/GEB is an example of excellent epistemological paradigm [27].

The number of electron-active elements (considered as players, in terms of Approach I to GEB) in a redox system, considered according to GATES/GEB principles, is practically unlimited; among others, the systems with three [24] or four [1] players were considered.

In the modeling of real systems, it is assumed that an effect of the matter (such as H2O, CO2 and O2) exchange with the environment is negligibly small within the period designed for certain chemical operations made on the system.

Advertisement

9. Completion of balances

The set of balances for the system (I) ⇒ (IV) is composed of GEB (e.g., 14 or 22) and equations obtained from the balances Eqs. (7)Eqs. (12) are expressed in terms of molar concentrations, namely:

[H+1][OH1]+[K+1][HSO41]2[SO42][HC2O41]2[C2O42][HCO31] 2[CO32][MnO41]2[MnO42]+3[Mn+3]+2[MnOH+2]+[MnC2O4+1][Mn(C2O4)21]3[Mn(C2O4)33]+2[Mn+2]+[MnOH+1]2[Mn(C2O4)22]+2[Fe+2]+[FeOH+1]2[Fe(C2O4)22]4[Fe(C2O4)34]+3[Fe+3]+2[FeOH+2]+[Fe(OH)2+1]+4[Fe2(OH)2+4]+[FeSO4+1][Fe(SO4)21]+[FeC2O4+1][Fe(C2O4)21]3[Fe(C2O4)33]=0E23
[K+1]=CV/(V0+V)E24
CB(S):[HSO41]+[SO42]+[MnSO4]+[FeSO4]+[FeSO4+1]+2[Fe(SO4)21]=(C01+C03)V0/(V0+V)E25
CB(C):2([H2C2O4]+[HC2O41]+[C2O42])+[H2CO3]+[HCO31]+[CO32]+2[MnC2O4]+4[Mn(C2O4)22]+4[Fe(C2O4)22]+6[Fe(C2O4)34]+2[FeC2O4+1]+4[Fe(C2O4)21]+6[Fe(C2O4)33]+2a1[FeC2O4]+2a2[MnC2O4]=((2C02+C04)V0+C1V)/(V0+V)E26
CB(Fe):[Fe+2]+[FeOH+1]+[FeSO4]+[Fe(C2O4)22]+[Fe(C2O4)34]+[Fe+3]+[FeOH+2]+[Fe(OH)2+1]+2[Fe2(OH)2+4]+[FeSO4+1]+[Fe(SO4)21]+[FeC2O4+1]+[Fe(C2O4)21]+[Fe(C2O4)33]+a1[FeC2O4]=C01V0/(V0+V)E27
CB(Mn):[MnO41]+[MnO42]+[Mn+3]+[MnOH+2]+[MnC2O4+1]+[Mn(C2O4)21]+[Mn(C2O4)33] +[Mn+2]+[MnOH+1]+[MnSO4]+[MnC2O4]+[Mn(C2O4)22]+a2[MnC2O4]=CV/(V0+V)E28

The balances can be specified as equations or equalities. The equality is represented here by relation Eqs. (24), where only one species is involved. In Eqs. (23) and Eqs. (25)Eqs. (28), we have concentrations of more species, interrelated in expressions for equilibrium constants, formulated on the basis of the proper stoichiometric reaction notations. As with results from the above, we have seven balances: six equations and one equality for the system (I) ⇒ (IV); the equality Eqs. (24) can enter immediately the charge balance (23). The current volume V of titrant T that is added is a parameter (steering variable) in these balances.

10. The arrangement of relations for equilibrium constants

The balances written in terms of molar concentrations are congruent with a complete set of independent (non-contradictory [9, 16]) relations for the equilibrium constants, interrelating concentrations of some species in the balances.

The complete set of equilibrium constants, namely ionic product of water, dissociation constants, standard potentials, stability constants of complexes and interrelated concentrations of the species (except K+1) is presented in Eq. (2).

{1}[H+1][OH1]=1014;{2}[HSO41]=101.8[H+1][SO42];{3}[H2C2O4]=105.2[H+1][C2O42] ;{4}[HC2O41]=103.8[H+1][C2O42] ;E54
{5}[H2CO3]=1016.4[H+1]2[CO32];{6}[HCO31]=1010.1[H+1][CO32];{7}[H2CO3]=[H2C2O40.5°10A(E+0.396)+pH;E55
{8}[MnO41]=[Mn+2]105A(E1.507)+8pH;{9}[MnO42]=[Mn+2]104A(E1.743)+8pH;{10}[Mn+3]=[Mn+2]10A(1.509);E56
[11][Fe+3]=[Fe+2]10A(E0.771);{12}[FeOH+1]=104.5[Fe+2][OH1] ;{13}[FeSO4]=102.3[Fe+2][SO42] ;E57
{14}[Fe(C2O4)22]=104.52[Fe+2][C2O42]2;{15}[Fe(C2O4)34]=105.22[Fe+2][C2O42]3;{16}[FeOH+2]=1011.0[Fe+3][OH1] ;E58
{ 17 }[Fe(OH)2+1]=1021.7[Fe+3][OH1]2;{ 18 }[Fe2(OH)2+4]=1025.1[Fe+3]2[OH1]2;{ 19 }[FeSO4+1]=104.18[Fe+3][SO42] ;E59
{ 20 }[Fe(SO4)21]=107.4[Fe+3][SO42]2;{ 21 }[FeC2O4+1]=107.53[Fe+3][C2O42];{ 22 }[Fe(C2O4)21]=1013.64[Fe+3][C2O42]2;E60
{23}[Fe(C2O4)33]=1018.46[Fe+3][C2O42]3;{24}[MnOH+1]=103.4[Mn+2][OH1];{25}[MnSO4]=102.28[Mn+2][SO42] ;E61
{26}[MnC2O4]=103.82[Mn+2][C2O42];{27}[Mn(C2O4)22]=105.25[Mn+2][C2O42]2;{28}[MnOH+2]=1014.2[Mn+3][OH1] ;E62
{29}[MnC2O4+1]=109.98[Mn+3][C2O42] ;{30}[Mn(C2O4)21]=1016.57[Mn+3][C2O42]2;{31}[Mn(C2O4)33] =1019.42[Mn+3][C2O42]3;E63
{32}[Fe+2][C2O42]=106.7for pr1=FeC2O4and{33}[Mn+2][C2O42]=105.3for pr2=MnC2O4E29

11. Relation between the numbers of species, balances and equilibrium constants

For any electrolytic system, one can define the relationship between the numbers of (i) kinds of species (P) (with free H2O molecules included), (ii) independent balances (Q) and (iii) independent equilibrium constant (R) values.

We refer first to non-redox systems discussed in Section 3, and we have:

  • for the system I: P = 8, Q = 4 (from f0, f3, f4, f5) and R = 3 ({1,5,6})

  • for the system II: P = 11, Q = 4 (from f0, f3, f4, f5) and R = 6 ({1,2,5,6,12,13})

  • for the system III: P = 11, Q = 4 (from f0, f3, f4, f5) and R = 6 ({1–6})

  • for the system IV: P = 17, Q = 5 (from f0, f3, f4, f5, f6) and R = 11 ({1–-6,12–15,32})

Referring now to the redox systems, we have:

  • for the system (I) ⇒ (II): P = 25, Q = 7 (GEB,ChB,CB(K),CB(CO3),CB(S),CB(Fe),CB(Mn)) and R = 18 ({1,2,3,4,5,6,12,13,14,15,32})

  • for the system (I) ⇒ (III): P = 25, Q = 6 (GEB,ChB,CB(K),CB(C),CB(S),CB(Mn)) and R = 19 ({1–10,24–31,33})

  • for the system (I) ⇒ (IV): P = 40 (collected in 2), Q = 7 (GEB,ChB,CB(K),CB(C),CB(S),CB(Fe),CB(Mn)) and R = 33 (collected in 29)

On this basis, one can state the relationships:

  • P = Q + R + 1—for non-redox systems;

  • P = Q + R—for redox systems.

To standardize this relationship, it is (informally) assumed that the electron is one of the species in the redox systems. In this way, the number of species is increased by 1, and we can suggest the relationship

P=Q+R+1E30

as common for both redox and non-redox electrolytic systems, regardless of their degree of complexity. The relation (30), applicable in resolution of electrolytic systems, was first presented in Ref. [37]; it can be perceived as a useful counterpart of the Gibbs’ phase rule (of a similar ‘degree of complexity’) in this area of the knowledge.

It should be noted that the total number, P = 40, of kinds of species involved in the system (I) ⇒ (IV) is relatively high.

12. The steps realized within GATES/GEB

Modeling the electrolytic systems according to GATES/GEB consists of several interacting steps [1]: (1) the collection of preliminary data; (2) preparation of computer programs; (3) calculations and data handling and (4) the knowledge gaining; all the steps are indicated in Figure 1.

Figure 1.

The steps of modeling any electrolytic system [1].

12.1. Collection of preliminary data

The necessary physicochemical knowledge is mainly attainable in tables of physicochemical data, exemplified by monographs [3840]. It should be noted that the period of interest in this field of research is currently the past. On the other hand, the physicochemical constants originate from works released over several decades, which are clearly seen in Ref. [40], where the relevant information is included. The point is that these physicochemical constants were determined using models that had been adapted to current computing capabilities, especially in the pre-computer era; these calculations were based (exclusively, in principle) on the reaction stoichiometry. For example, the solubility products were determined on the basis of molar solubility, see for example, the remark in Refs. [41, 42], without checking whether the precipitate is the equilibrium solid phase in the system [4348]. Acid dissociation constants were mainly determined on the basis of the Ostwald’s formula (see e.g., [49, 50]). Conditional (‘formal’, not normal) potentials were determined for many redox pairs [51]. In the computer era, some new models resolved with use of iterative methods and were elaborated. The assumptions and implementation of these models in relevant experimental studies aroused a number of concerns expressed, among others, in Refs. [52, 53]. Despite these circumstances, GATES, and GATES/GEB in particular, provides a new, reliable tool, applicable for physicochemical knowledge gaining. Thanks to this tool, it will be possible, in the immediate future, for a renaissance of interest in this—so important, after all!—field of fundamental research, which cannot be creatively developed on the basis of the previous ‘paradigm’ [27] based on the stoichiometry concept.

12.2. Preparation of computer program

Modeling of electrolytic systems can be realized with the use of iterative computer programs, for example MATLAB, perceived as a universal and high-level programming language [54, 55]. From the viewpoint of the form of mathematical models, MATLAB is focused on the matrix algebra procedure. MATLAB allows to make a quick and accurate computation and visualization of numerical data.

The iterative computer programs, written in MATLAB language, are exemplified in Refs. [11, 12, 16].

12.3. Calculations and data handling

The calculations can be made at different levels of the preliminary, physicochemical knowledge about the system in question. What is more, some ‘variations on the subject’ can also be done for this purpose; it particularly refers to metastable and non-equilibrium systems. A special emphasis will be put on complex redox systems, where all types of elementary chemical reactions proceed simultaneously and/or sequentially. In all instances, one can follow measurable quantities (potential E, pH) in dynamic and static processes and gain the information about many details not measurable in real experiments; it particularly refers to dynamic speciation.

We refer here to dynamic (titration) redox systems, represented by the system of 2 + k nonlinear equations composed of GEB, charge balance and k (≥ 1) concentration balances. The results of calculations, made with the use of an iterative computer program, are presented graphically. Thus, the plots E = E(Φ) and pH = pH(Φ) for potential E and pH of the solution and log[Xizi] versus Φ (Eq. (5)) relationships (speciation curves) will be drawn.

The Φ concept is used for simpler systems, providing a kind of normalization (independence on V0) in the systems considered. Φ plays a key role in the formulation of the generalized equivalence mass (GEM) concept, introduced also by Michałowski [25]. In more complex systems, the volume V is put on the abscissa.

The numerical data can be visualized in the form of two- or three-dimensional graphs (2D, 3D), see for example Ref. [12].

12.4. Computer program for the system (I) ⇒ (III)

The set of independent equilibrium constants is involved in the algorithm needed for calculation purposes [14], realized in the system (I) ⇒ (IV) as specified above. An algorithm is a well-defined procedure, expressed by a sequence of unambiguous instructions, which allows a computer to solve a problem according to a computer program implemented for this purpose. The term ‘unambiguous’ indicates that there is no room for subjective interpretation.

The system (I) ⇒ (IV) and its subsystems (I) ⇒ (II) and (I) ⇒ (III) were simulated using an iterative computer program (MATLAB software, included in the optimization toolbox™). In particular, the computer program for the system (I) ⇒ (III) is as follows.

function F = Function_KMnO4_Na2C2O4(x)

global V Vmin Vstep Vmax V0 C C1 C0 C01 C02 fi H OH pH E Kw pKw A aa

global H2C2O4 HC2O4 C2O4 H2CO3 HCO3 CO3 K

global logH2C2O4 logHC2O4 logC2O4 logH2CO3 logHCO3 logCO3 logK

global Mn7O4 Mn6O4 HSO4 SO4 Na

global logMn7O4 logMn6O4 logHSO4 logSO4 logNa

global Mn3 Mn3OH Mn3C2O4 Mn3C2O42 Mn3C2O43

global logMn3 logMn3OH logMn3C2O4 logMn3C2O42 logMn3C2O43

global Mn2 Mn2OH Mn2SO4 Mn2C2O4 Mn2C2O42

global logMn2 logMn2OH logMn2SO4 logMn2C2O4 logMn2C2O42

global pr logpr q logq

pH=x(1);

E=x(2);

if aa==0

 Mn2=10.^-x(3);

 pr=0;

else

 pr=10.^-x(3);

end;

H2C2O4=10.^-x(4);

SO4=10.^-x(5);

H=10.^-pH;

pKw=14;

Kw=10.^-14;

OH=Kw./H;

A=16.92;

ZMn=25;

ZC=6;

Ksp=10.^-5.3;

HC2O4=10.^(pH-1.25).*H2C2O4;

C2O4=10.^(pH-4.27).*HC2O4;

H2CO3=10.^(A.*(E+0.386)).*H2C2O4.^0.5;

HCO3=10.^(pH-6.3).*H2CO3;

CO3=10.^(pH-10.1).*HCO3;

if aa==1

 Mn2=Ksp./C2O4;

end;

HSO4=10.^(1.8-pH).*SO4;

Mn7O4=Mn2.*10.^(5.*A.*(E-1.507)+8.*pH);

Mn6O4=10.^(A.*(0.56-E)).*Mn7O4;

Mn2OH=10.^3.4.*Mn2.*OH;

Mn2SO4=10.^2.28.*Mn2.*SO4;

Mn2C2O4=10.^3.82.*Mn2.*C2O4;

Mn2C2O42=10.^5.25.*Mn2.*C2O4.^2;

Mn3=Mn2.*10.^(A.*(E-1.509));

Mn3OH=10.^(pH-0.2).*Mn3;

Mn3C2O4=10.^9.98.*Mn3.*C2O4;

Mn3C2O42=10.^16.57.*Mn3.*C2O4.^2;

Mn3C2O43=10.^19.42.*Mn3.*C2O4.^3;

K=C.*V./(V0+V);

Na=C0.*V0./(V0+V);

 %Charge balance

F=[(H-OH+K+Na-HSO4-2.*SO4-HC2O4-2.*C2O4-HCO3-2.*CO3-Mn7O4-2.*Mn6O4…

 +3.*Mn3+2.*Mn3OH+Mn3C2O4-Mn3C2O42-3.*Mn3C2O43+2.*Mn2+Mn2OH…

 -2.*Mn2C2O42);

 %Concentration balance of Mn

 (Mn7O4+Mn6O4+Mn3+Mn3OH+Mn3C2O4+Mn3C2O42+Mn3C2O43…

 +Mn2+Mn2OH+Mn2SO4+Mn2C2O4+Mn2C2O42+aa.*pr-C.*V./(V0+V));

 %Concentration balance of C

 (2.*H2C2O4+2.*HC2O4+2.*C2O4+H2CO3+HCO3+CO3+2.*Mn2C2O4…

 +4.*Mn2C2O42+2.*Mn3C2O4+4.*Mn3C2O42+6.*Mn3C2O43…

 +2.*aa.*pr-(2.*C0.*V0+C02.*V0+C1.*V)./(V0+V));

 %Concentration balance of S

 (HSO4+SO4+Mn2SO4-C01.*V0./(V0+V));

 %Electron balance

 ((ZMn-7).*Mn7O4+(ZMn-6).*Mn6O4+(ZMn-3).*(Mn3+Mn3OH)…

 +(ZMn-2).*(Mn2+Mn2OH+Mn2SO4)+(ZC-4).*(H2CO3+HCO3+CO3)…

 +2.*(ZC-3).*(H2C2O4+HC2O4+C2O4)+(ZMn-3+2.*ZC-6).*Mn3C2O4…

 +(ZMn-3+4.*ZC-12).*Mn3C2O42+(ZMn-3+6.*ZC-18).*Mn3C2O43…

 +(ZMn-2+2.*ZC-6).*Mn2C2O4+(ZMn-2+4.*ZC-12).*Mn2C2O42+…

 +(ZMn-2+2.*ZC-6).*aa.*pr…

 -((2.*ZC-6).*C0.*V0+(ZC-4).*C02.*V0+(ZC-4).*C1.*V…

 +(ZMn-7).*C.*V)./(V0+V))];

q=Mn2.*C2O4./Ksp;

logMn2=log10(Mn2);

logMn2OH=log10(Mn2OH);

logMn2SO4=log10(Mn2SO4);

logq=log10(q);

logpr=log10(pr);

logMn2C2O4=log10(Mn2C2O4);

logMn2C2O42=log10(Mn2C2O42);

logMn3=log10(Mn3);

logMn3OH=log10(Mn3OH);

logMn3C2O4=log10(Mn3C2O4);

logMn3C2O42=log10(Mn3C2O42);

logMn3C2O43=log10(Mn3C2O43);

logMn6O4=log10(Mn6O4);

logMn7O4=log10(Mn7O4);

logH2CO3=log10(H2CO3);

logHCO3=log10(HCO3);

logCO3=log10(CO3);

logH2C2O4=log10(H2C2O4);

logHC2O4=log10(HC2O4);

logC2O4=log10(C2O4);

logHSO4=log10(HSO4);

logSO4=log10(SO4);

logNa=log10(Na);

logK=log10(K);

%The end of program

13. Graphical presentation of the data

The results of calculations in the system (I) ⇒ (IV) are presented in Figures 24. More detailed, numerical data are specified in Ref. [14].

Figure 2.

The (a) E versus V and (b) pH versus V relationships plotted at V0 = 100, C = 0.02, C03 = 0.5, C1 = C04 = 0.001 and indicated pairs of C01 and C02 values.

Figure 3.

The speciation curves plotted for (a) Mn and (b) Fe species at V0 = 100, C = 0.02, C03 = 0.5, C1 = C04 = 0.001 and C01 = C02 = 0.002.

Figure 4.

The log(qi) versus V relationships (see Eq. (31)) plotted for (a) Fe (i = 1) and (b) Mn (i = 2) oxalates.

There are valid relationships: CVeq1 = 0.2∙C01V0 for iron and CVeq2 = 0.4∙C02V0 for oxalate. For V0 = 100, C = 0.02 we have, in particular, Veq1 = 10 and Veq2 = 20 at C01 = 0.01, C02 = 0.01 and Veq1 = 20 and Veq2 = 40 at C01 = 0.02 and C02 = 0.02. This agrees exactly with the position of the points where jumps of potential E occur and should be compared with the plots of titration curves for individual analytes: FeSO4 (Figure 5a) and H2C2O4 (Figure 6a), where abscissas are expressed in terms of the fraction titrated Φ (Eq. (5)). The related pH versus Φ relationships are presented in Figures 5b and 6b. To explain/formulate the reactions occurred in the systems together with their relative efficiencies, the speciation diagrams depicted in Figures 3a, b, and 5c, d are used. From Figure 4a and b, we see that

Figure 5.

The relationships: (a) E = E(Φ), (b) pH = pH(Φ) and the speciation curves log[Xizi] versus Φ for (c) manganese and (d) iron species plotted for titration of V0 = 100 mL of C0 = 0.01 mol/L FeSO4 + C01 mol/L H2SO4 as D with V mL of C = 0.02 mol/L KMnO4 as T.

Figure 6.

The relationships: (a) E = E(Φ), (b) pH = pH(Φ) and the speciation curves log[Xizi] for (c) Mn-, (d) C-species for titration of V0 = 100 mL of C0 = 0.01 mol/L H2C2O4 + C01 mol/L H2SO4 as D with V mL of C = 0.02 mol/L KMnO4 as T.

q1=[Fe+2][C2O42]/Ksp1<1andq2=[Mn+2][C2O42]/Ksp2<1E31

that is the precipitates FeC2O4 (Ksp1) and MnC2O4 (Ksp2) do not exist in this system as the equilibrium solid phases at the pre-assumed, sufficiently high concentration C03 of H2SO4; MnO2 is not formed there as well, that is, a1 = a2 = 0 in Eqs. (26)Eqs. (28).

14. Deficiency and veracity of equilibrium data

In some ‘variations on the subject’, we try to know what would happen if some constraints put on the metastable system are removed and the reaction is conducted in a thermodynamic manner, in accordance with the conditions imposed by the equilibrium constants [1, 36]. One can also analyze the data resulting from (intentional) omission or (factual, presupposed) incomplete physicochemical knowledge on the system studied.

Some computer simulations can be used to check some effects involved with complexation phenomena. For example, we intend to check the effect involved with formation of sulfate complexes FeSO4, MnSO4 and (particularly) FeSO4+1, Fe(SO4)2−1 in the system (I) ⇒ (II). The shapes of the titration curves E = E(Φ) are compared in Figure 7.

Figure 7.

The E versus Φ relationships plotted for the system (I) ⇒ (II): (1) at pre-assumed physicochemical knowledge and (2) after intentional omission of all sulfate complexes; C0 = 0.01, C01 = 1.0 and C = 0.02.

Some equilibrium constants used in calculations may be unknown/doubtful on the stage of collection of equilibrium data. In such instances, the pre-assumed/virtual data can be introduced for comparative purposes, and the effects involved with omission/inclusion of some types of complexes can be checked.

The possible a priori complexes of Mn(SO4)i+3−2i are unknown in literature. To check the effect of formation of these complexes on the shape of the titration curve E = E(Φ) in the system (I) ⇒ (II), the pre-assumed stability constants K3i of the complexes, [Mn(SO4)i+3−2i] = K3i[Mn3+][SO42−]i, specified in legend for Figure 8, were applied in the related algorithm, where concentrations of MnSO4+1 and Mn(SO4)2−1 with the corresponding multipliers were inserted in electron (GEB) and charge balances and in concentration balances for Mn and sulfate. As we see, at higher K3i values (comparable to ones related to Fe(SO4)i+3−2i (i=1,2) complexes [39]), the new inflection point appears at Φ = 0.25 and disappears at lower K3i values assumed in the simulating procedure. Comparing the simulated curves with the one obtained experimentally [25], one can conclude that the complexes Mn(SO4)i+3−2i do not exist at all, or the K3i values are small, when compared with those for Fe(SO4)i+3−2i [33].

Figure 8.

Fragments of hypothetical titration curves for V0 = 100 mL of FeSO4 (C0 = 0.01 mol/L) + H2SO4 (Ca0 = 0.1 mol/L) titrated with C = 0.02 mol/L KMnO4, plotted at different pairs of stability constants (K31, K32) of the sulfate complexes Mn(SO4)i+3–2i: (1) (104, 107), (2) (103, 106), (3) (102.5, 105), (4) (102, 104), (5) (104, 0), (6) (103, 0), (7) (102, 0) and (8) (0, 0).

Other interesting examples involved with ‘variations on the subject’ are presented in Ref. [1], and other references cited therein.

15. Advantages of GATES and usefulness of chemical processes simulation

Mathematical formalism of electrolytic systems tested and resolved according to GATES principles formulated by Michałowski (1992) arises from synthesis of three laws: (1o) law of charge conservation, (2o) law of elements conservation and (3o) law of mass action. All other chemical laws result from conjunction of those laws; it particularly refers to the stoichiometry and equivalent mass concepts.

GATES, based on physical, physicochemical and chemical laws, is considered as the best thermodynamic approach to equilibrium, non-equilibrium and metastable, mono- and poly-phase, static and dynamic and redox and non-redox systems, of any degree of complexity, with liquid-liquid extraction systems included [31].

GATES related to redox systems is denoted as GATES/GEB. All these systems are resolvable with use of iterative computer programs, for example, MATLAB. The complexity of chemical systems is here of a secondary importance from the point of view of the computational capabilities inherent in iterative computer programs.

GATES is a confirmation of the thesis that ‘everything brilliant is simple’. GATES/GEB is the unique tool to obtain information about the thermodynamics of redox systems on the basis of balances and equilibrium constants values.

GATES enables to simulate all possible (from a thermodynamic point of view) processes obtained after pre-assumed crossing of one or more reaction paths in metastable systems.

GATES enables to simulate the processes impossible to track experimentally; for example, dissolving a solid phase in the electrolytic system of a pre-established composition.

GATES relies on the assumption that the chemistry involved with such systems is predictable on the basis of knowledge of physicochemical properties of the species involved in the system in question. A complete set of non-contradictory relations for the equilibrium constants must be used in calculations; this ‘iron rule’ of mathematics is then obligatory also in calculations related to electrolytic systems.

GATES is the intrinsically consistent theory, joining fundamental laws of physics and chemistry [10, 28]. The knowledge gaining from redox systems is the most comprehensive way for studying such systems. Note that equations-based simulations are most commonly used in physics and related sciences.

GATES joins, on the thermodynamic basis, four kinds of chemical interactions, named as acid-base, redox, precipitation and complexation reactions, extended on a liquid-liquid extraction in mono- and poly-phase systems. To a certain degree—one can perceive GATES as a spitting image of theory of everything (ToE), as the main, still unresolved issue in physics, aiming to elaborate the consistent theory, that links together four: strong, weak, electromagnetic and gravitational interactions.

GATES referred to electrolytic non-redox and redox systems and is considered as the best thermodynamic approach to such systems. GATES, based on physical (charge conservation), physicochemical (conservation of elements) and chemical (mass action) laws, is the best tool applicable for computer simulation of equilibrium, non-equilibrium and metastable and mono- and poly-phase electrolytic systems. GATES is the basis for the Generalized Equivalence Mass (GEM) concept, with no relevance to the chemical reaction notation.

One can also express a conviction that the discovery of the Approach II to GEB in context with GATES will lead to gradual elimination of the stoichiometry concept from the consciousness of chemists.

16. Final comments

This chapter provides comprehensive, compatible and consistent knowledge on modeling electrolytic redox and non-redox systems and further steps applied to gain the thermodynamic knowledge on the systems, referred mainly to aqueous media.

The Generalized Electron Balance (GEB) concept, related to electrolytic redox systems, is put in context with the principle of conservation of all elements in electrolytic redox systems, in aqueous, non-aqueous or mixed-solvent media. The GEB is fully compatible with charge and concentration balances and completes the set of 2 + k equations needed for quantitative description of a redox system, with 2 + k independent/scalar variables xT=(E,pH,pX1,,pXk). Two equivalent approaches (I and II) to GEB were proposed (1992, 2006) by Michałowski. The Approach I to GEB is based on a card-game principle, with electron-active elements as gamblers, electron-non-active elements as fans and common pool of electrons introduced by electron-active elements as money. The Approach II to GEB is based on the linear combination 2f(O)−f(H) of elemental balances: f(H) for H and f(O) for O. The linear independency/dependency of 2f(O)−f(H) from charge and other elemental/core balances referred to the system in question provides the general criterion distinguishing between redox and non-redox systems. For non-redox systems, 2f(O)−f(H) is the linear combination of those balances, that is, it is not a new, independent equation in such systems. In redox systems, 2f(O)−f(H) is the independent equation, considered as the primary form of GEB and denoted as pr-GEB. The balances for elements/cores ≠ H, O are the basis for k concentration balances, forming—with GEB and charge balance—the set of 2+k independent balances, expressed in terms of concentrations. The Approach I to GEB, considered as the ‘short’ version of GEB, can be applied if the oxidation numbers for all elements in components forming a system and in the species of the system are known beforehand. The Approach II to GEB needs none prior information on oxidation numbers of all elements in components and the species in the system. Within the Approaches I and II to GEB, the roles of oxidants and reducers are not ascribed to the components and particular species. The GEB is put in context with the Generalized Approach to Electrolytic Systems (GATES) as GATES/GEB, where all quantitative thermodynamic knowledge on the redox system is involved in the complete set of independent equilibrium constants, where standard potentials E0i are involved. The GATES/GEB provides the best thermodynamic formulation of electrolytic redox systems of any degree of complexity, namely: equilibrium, non-equilibrium and metastable, mono- and poly-phase and static and dynamic electrolytic systems, resolvable with the use of iterative computer programs, applied to the set of nonlinear equations, with no simplifying assumptions needed. The GATES/GEB can also be referred to as redox systems in mixed-solvent media, provided that the related thermodynamic knowledge is attainable. This chapter is referred to dynamic systems, realized according to the titrimetric mode. The results obtained from the calculations can be presented graphically on 2D or 3D diagrams. The speciation diagrams obtained according to GATES/GEB have indisputable advantage over Pourbaix predominance diagrams. The GEB concept, unknown before 1992, is perceived as an emanation of the matter/elements conservation, as the general law of nature. The redox systems are formulated on simple principles, unknown in earlier literature. Earlier approaches to electrolytic redox systems, based on stoichiometric principles, are thus invalidated.

References

  1. 1. Michałowski T. Application of GATES and MATLAB for resolution of equilibrium, metastable and non-equilibrium electrolytic systems. In: Michałowski T, editor. Applications of MATLAB in Science and Engineering. Rijeka: InTech—Open Access publisher in the fields of Science, Technology and Medicine; 2011. pp. 1–34
  2. 2. Pilarski B, Dobkowska A, Foks H, Michałowski T. Modelling of acid–base equilibria in binary-solvent systems: A comparative study. Talanta. 2010;80(3):1073–1080
  3. 3. Michałowski T, Pilarski B, Asuero AG, Dobkowska A, Wybraniec S. Determination of dissociation parameters of weak acids in different media according to the isohydric method. Talanta. 2011;86:447–451
  4. 4. Michałowska-Kaczmarczyk AM, Michałowski T. Generalized electron balance for dynamic redox systems in mixed-solvent media. Journal of Analytical Sciences, Methods and Instrumentation. 2014;4(4):102–109
  5. 5. Michałowska-Kaczmarczyk AM, Michałowski T. Linear dependence of balances for non-redox electrolytic systems. American Journal of Analytical Chemistry. 2014;5:1285–1289
  6. 6. Michałowski T, Pilarski B, Asuero AG, Michałowska-Kaczmarczyk AM. Modeling of acid-base properties in binary-solvent systems. In: Wypych G, editor. Handbook of Solvents. Vol. 1, Properties. Toronto: ChemTec Publishing; 2014. pp. 623–648
  7. 7. Michałowski T, Toporek M, Michałowska-Kaczmarczyk AM, Asuero AG. New trends in studies on electrolytic redox systems. Electrochimica Acta. 2013;109:519–531
  8. 8. Michałowski T, Michałowska-Kaczmarczyk AM, Toporek M. Formulation of general criterion distinguishing between non-redox and redox systems. Electrochimica Acta. 2013;112:199–211
  9. 9. Michałowska-Kaczmarczyk AM, Michałowski T. Comparative balancing of non-redox and redox electrolytic systems and its consequences. American Journal of Analytical Chemistry. 2013;4(10):46–53
  10. 10. Michałowska-Kaczmarczyk AM, Michałowski T. Compact formulation of redox systems according to GATES/GEB principles. Journal of Analytical Sciences, Methods and Instrumentation. 2014;4(2):39–45
  11. 11. Michałowska-Kaczmarczyk AM, Toporek M, Michałowski T. Speciation diagrams in dynamic iodide + dichromate system. Electrochimica Acta. 2015;155:217–227
  12. 12. Toporek M, Michałowska-Kaczmarczyk AM, Michałowski T. Symproportionation versus disproportionation in bromine redox systems. Electrochimica Acta. 2015;171:176–187
  13. 13. Michałowska-Kaczmarczyk AM, Asuero AG, Toporek M, Michałowski T. “Why not stoichiometry” versus “Stoichiometry—why not?” part II. GATES in context with redox systems. Critical Reviews in Analytical Chemistry. 2015;45(3):240–268
  14. 14. Michałowska-Kaczmarczyk AM, Michałowski T, Toporek M, Asuero AG. “Why not stoichiometry” versus “Stoichiometry—why not?” part III, extension of GATES/GEB on complex dynamic redox systems. Critical Reviews in Analytical Chemistry. 2015;45(4):348–366
  15. 15. Michałowska-Kaczmarczyk AM, Michałowski T, Toporek M, Pietrzyk A. Solubility and dissolution in terms of generalized approach to electrolytic systems principles. Journal of Analytical Sciences, Methods and Instrumentation. 2015;5:47–58
  16. 16. Michałowska-Kaczmarczyk AM, Michałowski T, Toporek M. Formulation of dynamic redox systems according to GATES/GEB principles. International Journal of Electrochemical Science. 2016;11:2560–2578
  17. 17. Meija J, Michałowska-Kaczmarczyk AM, Michałowski T. Redox titration challenge. Analytical and Bioanalytical Chemistry. 2017;409(1):11–13
  18. 18. Michałowski T, Michałowska-Kaczmarczyk AM, Meija J. Solution of redox titration challenge. Analytical and Bioanalytical Chemistry. 2017;409(17):4113–4115
  19. 19. Michałowski T, Kupiec K, Rymanowski M. Numerical analysis of the Gran methods. A comparative study. Analytica Chimica Acta. 2008;606(2):172–183
  20. 20. Ponikvar M, Michałowski T, Kupiec K, Wybraniec S, Rymanowski M. Experimental verification of the modified Gran methods applicable to redox systems. Analytica Chimica Acta. 2008;628(2):181–189
  21. 21. Michałowski T, Toporek M, Rymanowski M. Overview on the Gran and other linearization methods applied in titrimetric analyses. Talanta. 2005;65(5):1241–1253
  22. 22. Michałowski T, Baterowicz A, Madej A, Kochana J. An extended Gran method and its applicability for simultaneous determination of Fe(II) and Fe(III). Analytica Chimica Acta. 2001;442(2):287–293
  23. 23. Michałowski T, Stępak R. Evaluation of the equivalence point in potentiometric titrations with application to traces of chloride. Analytica Chimica Acta. 1985;172:207–214
  24. 24. Michałowski T. The generalized approach to electrolytic systems: I. Physicochemical and analytical implications. Critical Reviews in Analytical Chemistry. 2010;40(1):2–16
  25. 25. Michałowski T, Pietrzyk A, Ponikvar-Svet M, Rymanowski M. The generalized approach to electrolytic systems: II. The generalized equivalent mass (GEM) concept. Critical Reviews in Analytical Chemistry. 2010;40(1):17–29
  26. 26. Asuero AG, Michałowski T. Comprehensive formulation of titration curves referred to complex acid-base systems and its analytical implications. Critical Reviews in Analytical Chemistry. 2011;41(2):151–187
  27. 27. Michałowski T, Asuero AG. New approaches in modeling the carbonate alkalinity and total alkalinity. Critical Reviews in Analytical Chemistry. 2012;42(3):220–244
  28. 28. Michałowski T, Asuero AG, Ponikvar-Svet M, Toporek M, Pietrzyk A, Rymanowski M. Principles of computer programming applied to simulated pH–static titration of cyanide according to a modified Liebig-Denigès method. Journal of Solution Chemistry. 2012;41(7):1224–1239
  29. 29. Michałowski T. Some remarks on acid-base titration curves. Chemia Analityczna. 1981;26:799–813
  30. 30. Michałowski T, Lesiak A. Euroanalysis VIII Book of Abstracts, PJ 35, PM 10, Edinburgh, Scotland; Young RA. Analyst. 1994;119:75N-76N
  31. 31. Michalowski T. Calculation of pH and potential E for bromine aqueous solutions. Journal of Chemical Education. 1994;71(7):560–562
  32. 32. Michałowski T, Lesiak A. Acid-base titration curves in disproportionating redox systems. Journal of Chemical Education. 1994;71(8):632–636
  33. 33. Michałowski T, Lesiak A. Formulation of generalized equations for redox titration curves. Chemia Analityczna (Warsaw). 1994;39:623–637
  34. 34. Michałowski T, Wajda N, Janecki D. A unified quantitative approach to electrolytic systems. Chemia Analityczna (Warsaw). 1996;41:667–685
  35. 35. Michałowski T. Calculations in analytical chemistry with elements of computer programming (in Polish). Cracow: PK; 2001.
  36. 36. Michałowski T, Rymanowski M, Pietrzyk A. Non-typical Brönsted’s acids and bases. Journal of Chemical Education. 2005;82(3):470–472
  37. 37. Michałowska-Kaczmarczyk AM, Michałowski T. General relation valid for electrolytic systems. Journal of Analytical Sciences, Methods and Instrumentation. 2015;5:74–85
  38. 38. Sillén LG, Martell AE. Stability Constants of Metal-ion Complexes. London: The Chemical Society; 1964. Suppl. No. 1, 1971 Other options of the notation met in literature: Sillén, L.G., Martell, A.E. Stability Constants of Metal-Ion Complexes, The Chemical Society, London, 1964, Suppl. No.1, 1971.
  39. 39. Lurie J. Handbook of Analytical Chemistry. Moscow: Mir Publishers; 1975
  40. 40. Inczédy J. Analytical Applications of Complex Equilibria. New York: Wiley; 1976
  41. 41. Michałowska-Kaczmarczyk AM, Michałowski T. Solubility product challenge. Analytical and Bioanalytical Chemistry. 2015;407(6):1789–1791
  42. 42. Michałowska-Kaczmarczyk AM, Michałowski T. Solution to the solubility product challenge. Analytical and Bioanalytical Chemistry. 2015;407(17):4877–4878
  43. 43. Michałowski T, Pietrzyk A. A thermodynamic study of struvite+water system. Talanta. 2006;68:594–601
  44. 44. Michałowski T, Pietrzyk A. Quasistatic processes in non-equilibrium two-phase systems with ternary salts: I. struvite + aqueous solution (CO2 + KOH). Chemia Analityczna (Warsaw). 2008;53:33–46
  45. 45. Michałowski T, Borzęcka M, Toporek M, Wybraniec S, Maciukiewicz P, Pietrzyk A. Quasistatic processes in non-equilibrium two-phase systems with ternary salts: II. Dolomite + aqueous media. Chemia Analityczna (Warsaw). 2009;54:1203–1217
  46. 46. Michałowski T, Asuero AG. Thermodynamic modelling of dolomite behaviour in aqueous media. Journal of Thermodynamics. 2012;2012:12. Article ID 723052
  47. 47. Michałowska-Kaczmarczyk AM, Michałowski T. Calculation of solubility of oxyquinolinates. Journal of Analytical Sciences, Methods and Instrumentation. 2014;4(3):71–79
  48. 48. Michałowska-Kaczmarczyk AM, Michałowski T. Evaluation of transition points between different solid phases in aqueous media. Journal of Analytical Sciences, Methods and Instrumentation. 2014;4(3):87–94
  49. 49. Michałowska-Kaczmarczyk AM, Michałowski T. Ostwald’s dilution law challenge. Analytical and Bioanalytical Chemistry. 2014;406(12):2741–2742
  50. 50. Michałowska-Kaczmarczyk AM, Michałowski T. Solution to Ostwald’s dilution law challenge. Analytical and Bioanalytical Chemistry. 2014;406(22):5251–5251
  51. 51. http://goldbook.iupac.org/C01240.html [accessed on August 5, 2017]
  52. 52. Michałowska-Kaczmarczyk AM, Asuero AG, Michałowski T. “Why not stoichiometry” versus “Stoichiometry—why not?” Part I, general context. Critical Reviews in Analytical Chemistry. 2015;45(2):166–188
  53. 53. Beck MT. Chemistry of Complex Equilibria. New York: Van Nostrand Reinhold; 1970
  54. 54. Michałowski T, editor. Applications of MATLAB in Science and Engineering. Published by InTech Janeza Trdine 9, 51000 Rijeka, Croatia, 2011
  55. 55. Katsikis VN, editor. MATLAB—A Fundamental Tool for Scientific Computing and Engineering Applications. Vol. 2. Published by InTech Janeza Trdine 9, 51000 Rijeka, Croatia, 2012

Written By

Anna Maria Michałowska-Kaczmarczyk, Aneta Spórna-Kucab and Tadeusz Michałowski

Submitted: 20 March 2017 Reviewed: 14 April 2017 Published: 27 September 2017