Open access peer-reviewed chapter

Generalized Electron Balance (GEB) as the Law of Nature in Electrolytic Redox Systems

Written By

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

Submitted: 13 December 2016 Reviewed: 20 April 2017 Published: 06 September 2017

DOI: 10.5772/intechopen.69341

Chapter metrics overview

1,769 Chapter Downloads

View Full Metrics


This chapter refers to fundamental/general/obligatory regularities of electrolytic systems. The linear combination 2·f(O) − f(H) of elemental balances, f(H) for H and f(O) for O, provides a rigorous criterion distinguishing between redox and non-redox systems is presented as the general relation distinguishing between electrolytic redox and non-redox systems in aqueous media. As the linearly independent equation for a redox system, 2·f(O) − f(H) is considered as the primary form of the generalized electron balance (GEB), perceived as a law of nature, as the hidden connection of physicochemical laws and the breakthrough in thermodynamic theory of electrolytic redox systems. GEB completes the set of 2+K equations necessary for thermodynamic resolution of redox systems according to generalized approach to electrolytic systems (GATES) applying all relevant, physicochemical knowledge available. GATES/GEB, perceived as an example of excellent paradigm, provides the best thermodynamic approach to electrolytic redox systems of any degree of complexity, in aqueous, non-aqueous, and mixed-solvent media. The formulation of GEB does not need prior knowledge of oxidation numbers for all elements in components forming any electrolytic system, within GATES/GEB, the stoichiometry, oxidation number, oxidant, reductant and equivalent mass are as derivative concepts.


  • electrolytic redox systems
  • approach I to GEB
  • approach II to GEB

1. Introduction

According to the principles assumed in the generalized approach to electrolytic systems (GATES) introduced/formulated by Michałowski in 1992 [1], a balancing of any electrolytic system is based on the rules of conservation of particular elements/cores Yg (g = 1,…, G), and on a charge balance (ChB) expressing the rule of electroneutrality of this system. The closed systems, separated from its environment by diathermal walls, are considered for modelling purposes. The elemental f(Eg) or/and core f(coreg) balances are denoted briefly as f(Yg) (f(Eg) or f(coreg). The balance for the gth element (Eg) or core (coreg) is expressed by the equation interrelating the numbers of gth atom or core in components as units composing the system with the numbers of atoms or cores of gth kind in the species of the system thus formed. For ordering purposes, we assume E1 = H (hydrogen) and E2 = O (oxygen); then, we have f(H) for Y1 = E1 = H, f(O) for Y2 = E2 = O, etc. Free water particles and water bound in hydrates are included in balances for f(H) and f(O).

The ChB interrelates the charged species (ions) in this system. A core is a cluster of elements with defined composition, expressed by its chemical formula, structure and external charge, which remains unchanged in a system considered, e.g. SO4−2 in Eq. (54).

The rules of conservation, formulated according to GATES principles [136], have the form of algebraic equations related to closed systems, composed of condensed (e.g. liquid, liquid + solid, etc.) phases separated from its environment by diathermal (freely permeable by heat) walls; it enables the heat exchange between the system and its environment. Any chemical process, such as titration, is carried out under isothermal conditions, in a quasistatic manner; constant temperature (T = constant) is one of the conditions securing constancy of equilibrium constant values. Any exchange of the matter (H2O, CO2, O2,…) between the system and its environment is thus forbidden for modelling purposes.

We refer first to aqueous media, where the species Xizi exist as hydrates XiziniW; zi = 0, ±1, ±2,… is a charge, expressed in terms of elementary charge unit, e = F/NA (F = 96485.333 C/mol, Faraday’s constant; NA = 6.022141·1023 mol−1, Avogadro’s number; ni (≡ niWniH2O) ≥ 0 is a mean number of water (W = H2O) molecules attached to Xizi; the case niW = 0 is then also admitted. For ordering purposes, we assume X2z2n2W= H+1n2W, X3z3n3W=OH1n3W, …, i.e. z2 = 1, z3 = −1,…. Molar concentration of the species XiziniW is denoted as [Xizi]. The ni = niW values are virtually unknown, even for X2z2=H+1 [37] in aqueous media, and depend on ionic strength (I) of the solution. The Xizis with different numbers of H2O molecules involved, e.g. H+1, H3O+1, H9O4+1; H4IO6−1, IO4−1; H2BO3−1, B(OH)4−1; AlO2−1 and Al(OH)4−1 are considered equivalently [27], i.e. as the same species in this medium. The ChB interrelates charged species (ions, zj ≠ 0) in the system.

From f(H) and f(O), the linear combination 2∙f(O) – f(H) is formulated and termed as the primary form of generalized electron balance (GEB), pr-GEB = 2∙f(O) – f(H), obtained according to approach II to GEB; this leitmotiv will be extended in further parts of this chapter. The GEB is the immanent part of GATES; the computer software applied to redox systems is denoted as GATES/GEB [1]. When related to redox systems, GATES is based on the generalized electron balance (GEB) [136] concept, perceived as a law of nature [1, 2, 14, 15, 22], as the hidden connection of physicochemical laws, as a breakthrough in the theory of electrolytic redox systems. The GATES refers to mono- and poly-phase, redox and non-redox, equilibrium and metastable [3842] static and dynamic systems, in aqueous, non-aqueous and mixed-solvent media [26, 29] and in liquid-liquid extraction systems [16, 43].

The generalized electron balance (GEB) concept, discovered by Michałowski as the approach I in 1992 and approach II in 2006 to GEB, plays the key role in redox systems; both approaches are equivalent:

approachII to GEB aapproach I to GEB

The GEB is fully compatible with charge balance (ChB) and concentration balances f(Yg) (g = 3,…, G), formulated for different elements and cores Yg. The elemental f(Eg) or/and core f(coreg) balances are transformed into concentration balances: CB(Yg) for g = 3,…, G.

To avoid redundant terms/naming, the acronyms ChB and GEB are applied both to equations, expressed in terms of particular units (Ni, N0j), or in terms of molar concentrations. On the basis of Eqs. (1a) and (1b) exemplified in Eqs. (2a) and (2b), this should not cause any misinterpretations. Then, i ∈ < 1, I > enumerate species, j ∈ < 1, J > enumerate components, g ∈ < 0, G > enumerate equations for ChB (g = 0) and elements/cores: g = 1 for H, g = 2 for O, g ∈ < 3, G > for other elements/cores, i.e. Y1 = H, Y2 = O,…, for ordering purposes.

The terms components and species are distinguished. In the notation applied here, N0j (j = 1, 2,…, J) is the number of molecules of components of jth kind composing the static or dynamic D + T system, whereby the D and T are composed separately, from defined components, including water. The mono- or two-phase electrolytic system thus obtained involve N1 molecules of H2O and Ni species of ith kind, XiziniW (i = 2, 3,…, I), specified briefly as Xizi (Ni, ni), where niniWniH2O; then, we have H+1 (N2, n2), OH−1 (N3, n3),…. Thus, the components form a (sub)system, and the species XiziniW enter the system thus formed. A solid XiniW (precipitate, zi = 0), as a species in a two-phase system, is marked by bold letters, e.g. I2(s) and AgCl.

In Example 1 (Section 4.1), CuSO4⋅5H2O is one of components, and Cu(OH)3−1n9H2O is one of the species in the system. N01 molecules of CuSO4⋅5H2O involve 10N01 atoms of H, 9N01 atoms of O and N01 atoms of Cu, and N01 atoms of S. The notation Cu(OH)3−1 (N9, n9) refers to N9 ions of Cu(OH)3−1n9H2O involving N9(3 + 2n9) atoms of H, N9(3 + n9) atoms of O and N9 atoms of Cu.

Molar concentration [Xizi] of XiziniW is as follow:

(a)[Xizi] =103(Ni/NA)/V0for a static system,or (b)[Xizi] =103(Ni/NA)/(V0+V)for a dynamic D+T systemE1

In a static or dynamic system, the balances are ultimately expressed in terms of molar concentrations of compounds and species, like the expressions for equilibrium constants. In particular, the charge balance (ChB) formulated as follows:


interrelates charged (zi ≠ 0) species of this system. In terms of molar concentrations [mol/L] (Eq. (1a) or (1b)), the charge balance has the form


where z1 = 0 for X1z1=H2O, z2 = +1 for X2z2=H+1, z3 = –1 for X3z3=OH1,…

In non-aqueous and mixed-solvent media, with amphiprotic (co)solvent(s) involved, we assume/allow the formation of mixed solvates Xizini1ni2niS, where nis=niAs(≥0) are the mean numbers of As (s = 1,…, S) molecules attached to Xizi. We apply the notation Xizi(Ni,niA1,niA2,,niAS), where Ni is a number of entities of these species in the system [25, 27, 28, 4446].


2. Preliminary information

The balances f(Yg) and ChB will be related to some dynamic redox and non-redox systems for comparative purposes. The balances for a given system are combined according to linear combination principles, and some general properties of the resulting balances are indicated.

The components and species in redox systems are involved in GEB, charge (ChB) and elemental balances. All the balances are founded on the well-established physical, physicochemical and chemical rules involved with charge and all elements conservation, and on the so-named mass action law (MAL), with its ‘old-fashioned’ principles. However, to be fully compatible with the GATES idea, introduced to avoid the stoichiometric reasoning, the equilibrium law (EL), suggested by Michałowski in 2016, was put instead of MAL; the principle of EL formulation is based on the idea of Lagrange multipliers for searching the local extrema of a function that subjects some constraints, expressed by GEB, ChB and concentration balances. The Gibbs free energy G is applied here as a function of the measurable, intensive properties: p and T and the numbers of constituents are convenient in the study of chemical reaction equilibrium. The problem in question is then consistent with the GATES formulation.

The generalized equivalent mass (GEM) concept, suggested by Michałowski in 1979 [11, 19, 47] is also based on the GATES principles, contrary to IUPAC recommendations [48] based on the stoichiometry of chemical reactions. Within GATES, the stoichiometric reaction notations are used only to formulate the expressions for the set of equilibrium constants related to the system in question, not to one ‘representative’ reaction, selected arbitrarily.

Within GATES, the mass conservation law is not limited only to components/species of one, specified/ ‘responsible’ chemical reaction equation but relates to all components/species of an electrolytic system. The mass change involved with an (exo- or endothermal) effect is negligible when compared with the total mass of the system. For example, the mass change, Δm, involved with enthalpy ΔH° of the reaction H2(g) + 0.5O2(g) = H2O(l)H° = – 286 kJ/mol H2O), equal Δm = ΔH°/c2 = –3.18∙10−9 g, is negligible (not measurable) when compared with 18 g of H2O; c = 299,792,458 m/s is the speed of light in vacuum. The neutralization or dilution gives much smaller heat effects. The resulting law of mass preservation is then fulfilled, irrespectively on whether stoichiometric or non-stoichiometric chemical reactions occur (or do not occur) in the system.

All considerations made in this chapter refer, in principle, to the systems where the elements are formed by stable/non-radioactive isotopes [2], i.e. where no nuclear (α, β+, β, or electron capture) transformations occur [7], with emission of γ and/or X-ray radiation. However, it is possible to extend the description of redox systems on the systems with radioactive elements, in which the concentrations of the respective components are expressed by dependencies, identified as Bateman’s system of linear differential equations [4952], binding quantitatively the radioactive elements with products of their decay. Radioactive elements and their decay products are also included in the balances for other elements and dependencies for the equilibrium constants.


3. Linear dependence of algebraic equations and transformation into identity

Linear combination of algebraic equations plays a fundamental/decisive role in thermodynamics of electrolytic systems, considered according to GATES [1]. An elementary information related to linear combination, perceived (mainly) from mathematical/algebraic viewpoint, can be found in Ref. [23]. It should be noted that the results of simple addition, i.e. f1 + f2 and simple subtraction, i.e. f1f2 or f2f1, are also linear combinations of any equations f1 and f2.

We refer here to the problem of linear dependency of balances, analogous to the problem of dependency of linear equations, considered in elementary algebra, see the picture below. In this context, the general property of linear independency, inherently involved with redox systems, will be emphasized.

(1) For the beginning, let us take the set of linear equations:


completed by the linear combination of these equations, i.e.


Applying the matrix algebra, we see that the determinant


has zero value


irrespectively on the c1 and c2 values; at D = 0, calculation of x1, x2 and x3 is then impossible.

(2) Let us consider the set of G + 1 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 dependency or independency 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 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)=bg  fg(x):φg(x)=bg  bg=φg(x)E8

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.

We refer now to the set of linear, algebraic equations


where agi are the coefficients and bg is the free terms. Multiplying Eq. (9) by ωg, after subsequent summation and rearrangement, we have




from Eqs. (10) and (11) we have


Referring to the problem in question, and placing xi = Ni, x0j = N0j in Eq. (11), we write


Then for ChB (Eq. (2a)), where the right side equals zero, we have


The charge balance (ChB) is expressed by Eqs. (2a) and (2b), where Xizi is defined by Eq. (1a) (for a static system) or Eq. (1b) (for a dynamic system). The elemental/core balances: f(H), f(O) and f(Yg) (Yg ≠ H, O, g = 3,…, G) are written as follows:

f1=f(H)=i=1I(a1i+2ni)Nij=1Jb1jN0j=0 forY1=H,f2=f(O)=i=1I(a2i+ni)Nij=1Jb2jN0j=0 forY2= O,,fg=i=1IagiNij=1JbgjN0j=0,,fG=i=1IaGiNij=1JbGjN0j=0E15

where agi and bgj are the numbers of atoms/cores of gth kind in the ith species as a constituent of the system and in the component of jth kind, respecively. Denoting, for a moment, ω1 = –1, ω2 = 2, we transform the balance


In Eq. (17), f1* and f2* have the shape similar to the general expression for fg (g = 3,…, G) in Eq. (15).

In the balances related to aqueous media, the terms involved with water, i.e. N0j (for j related to H2O, as the component), N1, and all ni = niW are not involved (in f0, f3,…, fG) or are cancelled within f12 (Eq. (17)). Other species, such as CH3COOH transformable (mentally, purposefully) into C2H4O2 ≡ C2(H2O)2, and the species in which H and O are not involved in Xizi are also cancelled within f12.

On the basis of relations (13) and (17), the linear combination of G + 1 balances f0, f12, f3,…, fG expressed by Eq. (14) can be presented in equivalent forms:


All multipliers at Ni and N0j in Eq. (18a) are cancelled simultaneously, if we have


for all i = 1,…, I and j = 1,…, J; then Eq. (18) is transformed into identity

i=1INi0=j=1JN0j0  0 = 0E20

Transformation of a set of the equations into the identity, 0 = 0, proves the linear dependence between the equations considered. Then from Eq. (18a), we have


i.e. f12 is the dependent balance.

Briefly, from G + 1 starting balances: f0, f1, f2, f3, …, fG we obtain G balances: f0, f12, f3,…, fG. If f12 is the dependent balance, we have G − 1 independent balances: f0, f3,…, fG; it is the case related to non-redox systems. If f12 is the independent balance, we have G independent balances: f0, f12, f3,…, fG, that will be rearranged, optionally, as the set (f12, f0, f3,…, fG) related to GEB, ChB and f(Yg) (g = 3,…, G), respectively. The number of elemental/core balances f(Yg) (Yg ≠ H, O, g = 3,…, G) and then the number of concentration balances CB(Yg) in redox systems equals K = G − 2.

As stated above, the linear combination of 2∙f(O) − f(H) with ChB and elemental/core balances f(Yg) (g = 3,…, G) provides the criterion distinguishing between non-redox and redox systems of different complexity [3, 23]. To obtain the simplest form of the linear combination, some useful/general rules allowing to select multipliers ωg of the corresponding balances are suggested. Namely, the ωg values in Eq. (10) are equal to the oxidation numbers of the electron-non-active elements in the related system. The proper linear combination of the balances is reducible (a) to the identity, 0 = 0, for non-redox systems or (b) to the simplest equation different from (not reducible to) the identity, but involving some terms related to the species and components of the system. This way we state that 2∙f(O) − f(H) is linearly dependent on the other balances in non-redox systems, i.e. it is not a new/independent balance in such systems. In redox systems, the balance 2∙f(O) − f(H) and its linear combinations with ChB and f(Yg) (g = 3,…, G) are the new equation, completing, as the generalized electron balance (GEB), the set of equations needed for resolution of electrolytic redox systems.

Static and dynamic non-redox and redox systems, with water as the main component, are considered. A static system of volume V0 mL is obtained by disposable mixing different components. The dynamic system is realized according to the titrimetric mode. At defined point of the titration B (C, V) ⇒ A (C0, V0), denoted briefly as T ⇒ D, V mL of titrant T containing the component B (C mol/L) is added into V0 mL of titrand D containing the component A (C0 mol/L), and V0 + V mL of D + T system/mixture is thus obtained, if the assumption of volume additivity is valid/tolerable; then D and T are sub-systems of the D + T system. The progress of the titration is represented by the fraction titrated Φ [11, 19, 47, 53, 54] value


For comparative purposes, we first refer to non-redox systems, where the derivation of relevant formulas will be carried out in detail also for training purposes. The regularities resulting from these examples, namely linear dependency of (linear) algebraic equations stated for non-redox electrolytic systems, will provide the reference point for further parts of this chapter, where redox systems are considered. Further examples in this chapter, related only to redox systems, will be presented in a synthetic manner, indicating the similarities and differences in the respective equations. The conclusions drawn here will provide the basis for further, important generalizations.


4. Linear combination of balances for non-redox systems

4.1. Examples

Example 1 (static): V0 mL of CuSO4 solution is prepared from N01 molecules of CuSO4⋅5H2O and N02 molecules of H2O. The resulting solution consists of the following species: H2O (N1), H+1 (N2, n2), OH−1 (N3, n3), HSO4−1 (N4, n4), SO4−2 (N5, n5), Cu+2 (N6, n6), CuOH+1 (N7, n7), Cu(OH)2 (N8, n8), Cu(OH)3−1 (N9, n9), Cu(OH)4−2 (N10, n10) and CuSO4 (N11, n11). The components and species are involved in the balances for particular elements, Yg: H, O, Cu, S, i.e.

f2= f(O): N1+N2n2+N3(1+n3)+N4(4+n4)+N5(4+n5)+N6n6+N7(1+n7) +N8(2+n8)+N9(3+n9)+N10(4+n10)+N11(4+n11)=9N01+N02E23
f12=2f(O)f(H): N2+N3+7N4+8N5+N7+2N8+3N9+4N10+8N11= 8N01E24
f0=ChB:N2N3N42N5+ 2N6+N7N92N10= 0E25
2f3=2f(Cu):2N01= 2N6+2N7+2N8+2N9+2N10+ 2N11E26
6f4=6f(S)=6f(SO4): 6N01= 6N4+6N5+6N11E27

Simple addition of the elemental/core balances (Eqs. (24)–(27)) gives the identity


As we see, the multipliers ωg in the transformed identity are equal to oxidation numbers of the indicated elements.

Example 2 (dynamic): V0 mL of CuSO4 solution, as a titrand D composed of N01 molecules of CuSO4⋅5H2O and N02 molecules of H2O, is titrated with V mL of titrant T, composed of N03 molecules of NaOH and N04 molecules of H2O, added up to a defined point of the titration. In the resulting D + T mixture (two-phase system of volume ca. V0 + V mL), we have the species as in Example 1 and, additionally: Cu(OH)2 (N12, n12), Na+1 (N13, n13). From the balances:

f1=f(H):   ()+N12(2+2n12)+2N13n13=()+N03+2N04E29
f2=f(O):   ()+N12(2+n12)+N13n13=()+N03+N04E30
f12=2f(O)f(H):   N2+N3+7N4+8N5+N7+2N8+3N9+4N10+8N11+2N12=8N01+N03E31
f0=ChB:   N2N3N42N5+2N6+N7N9 2N10+N13=0E32
2f3= 2f(Cu):    2N01= 2N6+2N7+2N8+2N9+2N10+2N11+2N12 E33
6f4=6f(S)=6f(SO4):       6N01=6N4+6N5+6N11E34
f5=f(Na):         N03=N13E35

That means we obtain the identity, 0 = 0, for

2f(O)f(H)+ChBf(Na)2f(Cu)6f(S)=0 (+1)f(H)+(2)f(O)+(+2)f(Cu)+(+6)f(S)ChB=0E36

The terms within (…)’s in Eqs. (29) and (30) are the same as in Eqs. (22) and (23).

4.2. A comment

The equations for f(H) and f(O) in involve the numbers ni = niW of hydrating water molecules, attached to particular species. In the linear combinations 2∙f(O) – f(H), all the niW are cancelled, together with N1, i.e. the niW and N1 do not enter the equation for 2∙f(O) – f(H). The niW is not involved in ChB and in f(Yj) for Yj ≠ H, O and then they are not introduced in the linear combinations of 2∙f(O) – f(H) with these equations. The multipliers ωg in Eq. (10) are equal to oxidation numbers of particular elements involved in the related elemental/core balances. Then, the linear combination of the balances related to a non-redox system, when put in context with Eqs. (19) and (20) gives the identity 0 = 0. This regularity is obligatory for all non-redox systems considered above (Examples 1 and 2), and elsewhere, e.g. [22, 37, 55]. Note that 2∙f(O) – f(H) involves –N2 + N3, whereas ChB involves N2N3 = – (–N2 + N3), and the related terms are cancelled within the combination 2∙f(O) – f(H) + ChB. Moreover, all components related to the species XiziniW, not involving H and/or O in Xizi, are cancelled too within 2∙f(O) – f(H).

In further sections, we formulate the balances for redox systems containing one, two or three electron-active elements in redox systems (aqueous media), where the complete set of expressions for independent equilibrium constants interrelating concentrations of different species is involved. We start our considerations from the systems with one electron-active element involved with disproportionation and symproportionation [5, 13, 15, 16, 19]. These properties are appropriate for the elements that form compounds and species at three or more oxidation degrees. In particular, bromine (Br) forms the species on five oxidation degrees (–1, –1/3, 0, 1, 5), see Figure 1, e.g. 0 (for Br2) and 1 (for HBrO or NaBrO) ϵ (–1, 5). There are possible transitions between different bromine species, associated with changes of the oxidation states of this element, see Figure 1.

Figure 1.

Schemes of disproportionation and symproportionation within bromine species [5].


5. Disproportionating systems

5.1. Formulation of balances for NaOH (C, V) ⇒ Br2 (C0, V0) system

5.1.1. Approach II to GEB

Note that V0 mL of the titrand D composed of Br2 (N01) + H2O (N02) is titrated with V mL of titrant T composed of NaOH (N03) + H2O (N04). For the V0 + V mL of the D + T mixture, with the set of the species completed by Na+1 (N11, n11).

The balances can be formulated as follows:

f0=ChB:N2N3N5N7N9N10+N11= 0E39
f4=f(Br): N4+N5+N6+N7+ 2N8+ 3N9+N102N01= 0E41

Then we formulate the linear combinations:

f0123=f12+f0f3:5(N4+N5)+N6+N7N9N10= 0E44

For comparative purposes, we formulate also the linear combination ZBrf3f0123

(ZBr5)(N4+N5)+(ZBr5)(N6+N7)+2ZBrN8+ (3ZBr+1)N9+ (ZBr+1)N10= 2ZBrN01E46

where ZBr = 35 is the atomic number for bromine (Br).

All the linear combinations (42)–(46) are not reducible to identity, 0 = 0; in the simplest case, we have six constituents (species, components) involved in f0123 (Eq. (44)). Eqs. (42)–(46) are equivalent forms of GEB, expressed in terms of particular constituents (components, species). To express them in terms of molar concentrations, we apply the relations: (1a) and C0V0 = 103N01/NA, CV = 103∙N03/NA. In this way, we consider the titration of V0 mL of Br2 (C0 mol/L) solution with V mL of C mol/L NaOH solution added up to a defined point of the titration. In particular, from Eq. (46), we have the balance for GEB

(ZBr5)([HBrO3]+[BrO31])+(ZBr1)([HBrO]+[BrO1])+2ZBr[ Br2 ]+(3ZBr+1)[Br31]+(ZBr+1)[Br1] = 2ZBrC0V0/(V0+V)E46a

completed by the balances for ChB, CB(Br) and CB(Na):

[ HBrO3 ]+[BrO31]+[ HBrO ]+[BrO1]+2[ Br2 ]+3[Br31]+[Br1]2C0V0/(V0+V)=0E46c

5.1.2. Approach I to GEB

The (optional) formulas for GEB, expressed by Eqs. (42)–(46), were obtained according to approach II to GEB, based on pr-GEB = 2∙f(O) – f(H). In this section, we apply the approach I to GEB, known also as the ‘short’ version of GEB, where prior knowledge of oxidation numbers of all elements in the system is assumed/required.

In the NaOH (C, V) ⇒ Br2 (C0, V0) system considered in Section 4.2, bromine is the only electron-active element, considered as the carrier of its own, bromine electrons. One atom of Br has ZBr bromine electrons, and then one molecule of Br2 has 2ZBr bromine electrons, i.e. N01 molecules of Br2 involve 2ZBrN01 bromine electrons. The oxidation degree x of an atom in simple species, such as ones formed here by bromine, is calculated on the basis of known oxidation degrees: +1 for H and –2 for O and external charge of this species. Then for HBrO3, we have, by turns: 1∙1 + 1∙x + 3∙(–2) = 0 ⇒ x = 5; for BrO3−1: 1∙x + 3∙(–2) = –1 ⇒ x = 5; for HBrO: 1∙1 + 1∙x + 1∙(–2) = 0 ⇒ x = 1, etc.

The oxidation number is the net charge resulting from the presence of charge carriers, inherently involved in an atom: protons in nuclei and orbital electrons, expressed in elementary charge units—as +1 for protons and –1 for electrons. The number y of bromine electrons in one molecule of HBrO3 is calculated from the formula: ZBr∙(+1) + y∙(–1) = 5, i.e., bromine in HBrO3 involves y = ZBr – 5 bromine electrons, etc. On this basis, we state that [15]

N4 species HBrO3n4H2O involve(ZBr – 5)∙N4bromine electrons;
N5 species BrO3−1n5H2O involve(ZBr– 5)∙N5bromine electrons;
N6 species HBrO∙n6H2O involve(ZBr – 1)∙N6bromine electrons;
N7 species BrO−1n7H2O involve(ZBr – 1)∙N7bromine electrons;
N8 species Br2n8H2O involve2ZBrN8bromine electrons;
N9 species Br3−1n9H2O involve(3ZBr + 1)∙N9bromine electrons;
N10 species Br−1∙n10H2O involve(ZBr + 1)∙N10bromine electrons;
N01 molecules of Br2 involved2ZBrN01bromine electrons.

Balancing the bromine electrons, we get Eq. (46), and then Eq. (46a), considered as the GEB obtained immediately from the approach I to GEB; this proves the equivalency of approaches I and II to GEB.

5.1.3. Preparation of an algorithm and computer program

Let us refer again to balances (46a)–(46d) interrelating concentrations [Xizi] of the species Xizi with the total concentrations C0, C of indicated components of the system. In this place, one should distinguish between the terms equation and equality. The term equation is related here to the balance where at least two species are involved; these species are interrelated by expressions for the corresponding equilibrium constants values (Table 1). Such a requirement is fulfilled by Eqs. (46a)–(46c), whereas in Eq. (46d) we have concentration of one species; at any V value, [Na+1] is a number (not variable) and, as such, it enters immediately ChB (Eq. (46b)).

No.ReactionEquilibrium equationEquilibrium data
1BrO3–1+6H+1+6e–1 = Br–1+3H2O[Br–1] = Ke1⋅[BrO3–1][H+1]6[e–1]6E04 = 1.45 V
2BrO–1+H2O+2e–1 = Br–1+2OH–1[Br–1] = Ke2⋅[BrO–1][H+1]2[e–1]2/KW2E03 = 0.76 V
3Br2 + 2e–1 = 2Br–1[Br–1]2 = Ke3⋅[Br2][e–1]2E03 = 1.087 V
4Br3–1 + 2e–1 = 3Br–1[[Br–1]3 = Ke4⋅[Br3-1][e–1]2E04 = 1.05 V
5HBrO3 = H+1 + BrO3–1[H+1][BrO3–1] = K51⋅[HBrO3]pK51 = 0.7
6HBrO = H+1 + BrO–1[H+1][BrO–1] = K11⋅[HBrO]pK11 = 8.6
7H2O = H+1 + OH–1[H+1][OH–1] = KWpKW = 14.0

Table 1.

Equilibrium data related to different bromine species.

where logKe1 = 6AE01, logKei = 2AE0i (i = 2, 3, 4); A = 16.92; pK51 = −logK51, pK11 = −logK11, pKW = −logKW.

From the interrelations obtained on the basis of equilibrium data [56] collected in Table 1 we have

[H+1]=10pH;[OH1]=10pH14;[BrO31]=106A(E1.45)pBr+6pH;[BrO1]=102A(E0.76)pBr+2pH28;[Br2]=102A(E1.087)2pBr;[Br31]=102A(E1.05)2pBr;[ HBrO3 ]=100.7pH[BrO31];[ HBrO ]=108.6pH[BrO1]E47

where the uniformly defined (scalar) variables: E, pH and pBr, forming a vector x = (E, pH, pBr)T, are involved; AE = – log[e–1], pH = –log[H+1], pBr = –log[Br–1]. All the variables are in the exponents of the power for 10 in [e–1] = 10-AE, [H+1] = 10-pH, [Br–1] = 10-pBr. The number of the (independent) variables equals to the number of equations, K = 2 + 1 = 3; this ensures a unique solution of the equations, at a pre-set C0, C and V0 values, and the V-value at which the calculations are realized, at a defined step of the calculation procedure, according to iterative computer program presented below.

Computer program for the NaOH_Br2 system

function F = NaOH_Br2(x)

global V C0 V0 C yy

  E = x(1);

  pH = x(2);

  pBr = x(3);

  H = 10^(-pH);

  Kw = 10^-14;

  pKw = 14;

  OH = Kw/H;

  A = 16.92;

  Br = 10^-pBr;

  ZBr = 35;








  F = [%Charge balance


   %Concentration balance for Br


   %Electron balance












5.1.4. Preparation of the computer program for NaOH (C, V) ⇒ HBrO (C0, V0) system

Modification of the computer program makes it possible to carry out calculations for other systems associated with bromine, e.g. for the system NaOH (C, V) ⇒ HBrO (C0, V0). In this case, the changes are made in the following lines:

Computer program for the NaOH_HBrO system

Function F = NaOH_HBrO(x)






5.1.5. Graphical presentation of the data and discussion

The systems considered here are characterized by distinctly marked jumps in E and pH values (Figures 2A and B) at the vicinity of the equivalence points, occurring at Φeq1 = 1 for NaOH (C, V) → HBrO (C0, V0) (system I), and at Φeq2 = 2 for NaOH (C, V) → Br2 (C0, V0) (system II), see Table 2 [15, 16].

System IISystem I
NaOH (C, V) → HBrO (C0, V0)NaOH (C, V) → Br2 (C0, V0)

Table 2.

The sets of rounded (Φ, pH, E) values taken from the vicinity of the equivalence points; V0 = 100, C0 = 0.01, C = 0.1.

Figure 2.

(A) pH = pH(Φ) and (B) E = E(Φ) relationships plotted for the systems I and II, and the related speciation diagrams for the systems: I (C) and II (D).

Reactions occurred in the systems I and II can be formulated from the analysis of the respective speciation diagrams. As results from Figure 2C, the disproportionation reaction in the system I occurs according to the scheme


i.e. Φeq = 2 ⇒ CVeq1 = 2C0V0. After crossing the equivalent point, the main products of this disproportionation reaction are Br−1 and BrO3−1, not Br−1 and BrO−1 corresponding to the reaction


with the same stoichiometry, 3 : 6 = 1 : 2. From detailed calculations at Φ = 2.5, i.e. at an excess of NaOH added, we find [BrO3−1]/[BrO−1] = 10−2.574/10−6.7584 = 1.53∙104, i.e. efficiency of the reaction (48) is more than 4 orders of magnitude greater than the efficiency of reaction (49).

Scheme of Br2 disproportionation affected by NaOH (C) can also be calculated [16]. At Φ = 2.5, we have [BrO3−1]/[Br−1] = 10−2.574/10−1.875 = 10−0.699 = 0.2 = 1:5, i.e. at an excess of NaOH added (see Figure 2c), the disproportionation occurs mainly according to the scheme indicated by reaction (48) from the products side.

In the system II, disproportionation of HBrO affected by NaOH (C) added according to the titrimetric mode is presented in Figure 2D [16]. The [Br−1]/[BrO3−1] ratio equals: 10−2.2553/10−2.5563 at Φ = 2.0; 10−2.2730/10−2.5740 at Φ = 2.5, i.e. 100.3010 = 2 = 2:1, corresponding to the stoichiometric ratio of products of this reaction. As results from Figure 2C, the disproportionation, at an excess of NaOH added, occurs mainly according to reaction 3HBrO + 3OH−1 = 2Br−1 + BrO3−1 + 3H2O (stoichiometry 3:3 = 1:1), resulting from half reactions: HBrO + 2e−1 + H+1 = Br−1 + H2O, HBrO – 4e−1 + 2H2O = BrO3−1 + 5H+1, and 3H+1 + 3OH−1 = 3H3O. The (Φ, pH, E) values from the close vicinity of the corresponding equivalence points on the curves in Figures 2A and B are collected in Table 2. The Br2 solution is acidic, as results, e.g. from the ChB (Eq. (29a)): at [Na+1] = 0 (Φ = 0) we have [H+1] – [OH−1] = [BrO3−1]+[BrO−1]+[Br3−1]+[Br−1] > 0, i.e. [H+1] > [OH−1]; this inequality is also obtained from Eq. (42), at N03 = 0: [H+1] – [OH−1] = 5[HBrO3]+6[BrO3−1] + [HBrO] + 2[BrO−1] > 0. Br2 is an acid with a strength comparable to that of acetic acid; at C0 = 0.01, pH equals 3.40 for Br2 and 3.325 for CH3COOH (pK1 = 4.65). Disproportionation of Br2 occurs initially to a small extent (several %), according to the scheme Br2 + OH−1 = HBrO + Br−1, compare with [57].

In C0 = 0.01 mol/L HBrO solution, more than 90% HBrO disproportionates occur according to the reaction 5HBrO = BrO3−1 + 2Br2 + 2H2O + H+1; at V = 0, we have [Br2] = 10−2.4406, [BrO3−1] = 10−2.7442, i.e. [Br2]/[BrO3−1] = 100.3036 ≈ 2, which confirms this stoichiometry of the reaction. The H+1 ions formed in this reaction acidify the solution significantly: at C0 = 0.01 and V = 0, we have pH = 2.74, although HBrO itself is a relatively weak acid.

The numerical values of the concentrations given here are taken from the corresponding files with results of iterative calculations.

5.2. Titration in NaOH (C, V) ⇒ HIO (C0, V0) system

The curves plotted in Figure 3a–c are related to titration of V0 = 10 mL of HIO (C0 = 0.1 mol/L) with V mL of C = 0.1 mol/L NaOH. At the initial part of the titration, we have the reactions:

Figure 3.

Plots of (a) E = E(Φ), (b) pH = pH(Φ), and (c) speciation curves for NaOH (C, V) ⇒ HIO (C0, V0) system at V0 = 10, C0 = 0.01, C = 0.1.


In the following, at Φ ca. 0.20–0.22, a pronounced increase in [I−1] occurs, as a result of reaction


The increase in [I−1] is accompanied by an increase in [I3−1]

(I2(s),I2) +I1=I31E52

This leads to the gradual disappearance of I2(s) (which is ultimately ended at Φ = 0.5347) and lowering of [I2] and [I3−1]; all them disproportionate

3(I2(s),I2,I31)+6OH1=IO31+(5, 5, 8)I1+3H2OE53

At [I2(s)] > 0, we have [I2] = s = const; s = 1.33⋅10−3 mol/L is the solubility of I2(s) in water, at 20°C. Finally, the disproportionation of HIO, affected by NaOH, can be expressed by the equation


Note that the stoichiometry of the reaction (54) is 3 : 3 = 1 : 1, which corresponds to the jump on the curves presented in Figures 3a,b, occurring at Φ = 1. For Φ > 1, we have [I−1]/[IO3−1] = 2, i.e. the stoichiometry of the products of reaction (54) equals to 1 : 2. The jumps on the curves E = E(Φ) and pH = pH(Φ) (Figures 3a,b) occur at Φ ca. 0.2 (which corresponds to the stoichiometry 1:5 of the reaction (50)) and at Φ ca. 1 (which corresponds to the stoichiometry 3:3 = 1:1 of the reaction (54)); the maxima on the corresponding derivative curves in Figures 4a,b fit the stoichiometric ratios.

Figure 4.

The (a) ΔpH/ΔΦ = (pHj+1 − pHj)/(Φj+1 − Φj), (b) ΔE/ΔΦ = (Ej+1Ej)/(Φj+1 − Φj) vs. Φ = (Φj + Φj+1)/2 relationships for the NaOH (C, V) ⇒ HIO (C0, V0) system.

5.3. Titration in HCl (C, V) ⇒ NaIO (C0, V0) system

V0 = 10 mL of C0 = 0.01 mol/L NaIO is titrated with C = 0.1 mol/L HCl. The related curves are presented in Figures 5a‐d. Initially, the reaction

Figure 5.

The speciation curves for indicated (a) iodine and (b) chlorine species Xizi and (c) E = E(Φ), (d) pH = pH(Φ) functions for HCl (C, V) ⇒ NaIO (C0, V0) system; V0 = 10, C0 = 0.01, C = 0.1.


and then, the reaction


occurs. Then I−1 from Eq. (55) and I2 from Eq. (56) form I3−1 in the reaction I2 + I−1 = I3−1 and [I3−1] increases. At Φ = 0.4654, I2(s) appears as the solid phase

5IO1+ 4H+1=2I2(s)+ IO31+2H2OE57

At [I2(s)] > 0, we have [I2] = const. The increase in [Cl−1], resulted from addition of HCl, causes an increase in [I2Cl−1], and—to a lesser extent—the increase in [ICl2−1] and [ICl]. The addition of HCl lowers pH of the solution, and then [HIO] becomes larger than [IO−1]; [HIO3] also increases. In effect, the summary concentration [HIO] + [IO−1] after addition of an excess of HCl is higher than in the starting NaIO solution.

In the algorithm, we have allowed the participation of Cl−1 ions from HCl solution in the redox reaction. However, the concentration of Cl2 and HClO as the main products of Cl−1 oxidation (Figure 5b) is quite negligible. This way, one can state that the Cl−1 ions practically do not participate in the redox reaction as a reducing agent. From linear combination of reactions (55) and


(multiplication by 5 and 2 respectively), cancellations and division by 3, we get the reaction


with stoichiometry 4:5 = 0.8, which corresponds to Φ = 4:5 = 0.8, where the inflection point on the curves in Figures 5c and d are observed. The I−1 and I3−1 ions are consumed in reactions (58) and


See Figure 5a.

Computer program for HCl (C, V) ⇒ NaIO (C0, V0) system

function F = System_NaIO_HCl(x);


% Titration of V0 mL of NaIO (C0) with V mL HCl (C).

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

global I I3 I2 I2s HIO IO HI5O3 I5O3 H5I7O6 H4I7O6 H3I7O6 Na

global logI logI3 logI2 logI2s logHIO logIO logHI5O3 logI5O3 logH5I7O6

global logH4I7O6 logH3I7O6 logNa

global Cl Cl2 HClO ClO HCl3O2 Cl3O2 Cl4O2 Cl5O3 Cl7O4 I2Cl ICl ICl2

global logCl logCl2 logHClO logClO logHCl3O2 logCl3O2 logCl4O2 logCl5O3

global logCl7O4 logI2Cl logICl logICl2 pI pCl



































if I2>1.33e-3








 %Charge balance



 %Concentration balance for I



 %Concentration balance for Cl



 %Electron balance































% The end of program

A remark: Some notations applied in this program are as follows: [HIO3] → HI5O3, [H3IO6−2] → H3I7O6, [ClO2−1] → Cl3O2, [ClO2] → Cl4O2, etc.

5.4. Disproportionation in static redox systems

In disproportionating redox systems, it is advisable to check an effect of dilution of the corresponding solutes on values of the corresponding variables. It is advisable to plot the desired relationships in the figures with the values pC = –logC on the abscissa, where C is the concentration [mol/L] of the solute considered. Formally, the related plots correspond to presentation, in extended logarithmic scale, the results of titration of the initial solution of this solute with pure water as a titrant T.

The static systems with C mol/L solutions of (1) HIO and (2) NaIO are shown graphically in Figures 6a‐c and 7a‐c, where pH, E and log [Xizi] values related to different concentrations C of the corresponding solutes with the values pC = −logC on the abscissa.

Figure 6.

The plots of (a) E vs. pC, (b) pH vs. pC relationships and (c) speciation curves for indicated iodine species Xizi in C mol/L HIO.

Figure 7.

The functions: (a) E = E(pC), (b) pH = pH(pC) and (c) speciation curves for indicated iodine species Xizi plotted for C mol/L NaIO.

5.4.1. C mol/L HIO

As results from the speciation diagram in Figure 6c, in more concentrated HIO solutions, i.e. at lower pC values, the predominating reactions are as follows: 5HIO = 2(I2(s), I2) + IO3−1 + 2H2O + H+1 and 5HIO = 2(I2(s), I2) + HIO3 + 2H2O; solubility of I2(s) s = 1.33∙10−3 mol/L (25°C). At further dilution of HIO, the reaction 3HIO = 2I−1 + IO3−1 + 3H+1 occurs in an increasing degree. This change in disproportionation scheme, more significant at pC 4–5, resulted in a change of the shapes of the plots: E = E(pC) (Figure 6a) and pH = pH(pC) (Figure 6b).

5.4.2. C mol/L NaIO

The disproportionation of IO−1 introduced by NaIO proceeds mainly according to the scheme 3HIO = IO3−1 + 2I−1 + 3H+1 (see Figure 7c), where [I−1]/[IO3−1] ≅ 2. Concentration of I2 is lower than 10−6 mol/L and then solid I2(s) is not formed.


6. Symproportionating systems

6.1. Titration in KBrO3 (C, V) ⇒ NaBr (C0, V0) system

In this case, symproportionation practically does not occur (Figure 8a); concentration of HBrO, as the major product formed in the symproportionation reaction


Figure 8.

Plots for the KBrO3 (C, V) ⇒ NaBr (C0, V0) system: (a) speciation diagram, (b) E = E(Φ), and (c) pH = pH(Φ) curves; V0 = 10, C0 = 0.01, C = 0.1.

is ca. 10−6 mol/L. The potential E increases monotonically (Figure 8b), whereas pH first increases, passes through maximum and then decreases (Figure 8c). The relevant pH and E changes are small. Binding the H+1 ions in reaction (61) causes a weakly alkaline reaction (Figure 8b).

6.2. Titration in KBrO3 (C, V) ⇒ NaBr (C0) + H2SO4 (C03) V0 system

The stoichiometry 1:5, i.e. Φeq = 0.2, stated for C03 values indicated at the curves plotted in Figure 9 (column a), results from reaction

Figure 9.

Plots for the KBrO3 (C, V) ⇒ NaBr (C0) + H2SO4 (C03) V0 system: speciation diagrams (column a); E = E(Φ) (column b) and pH = pH(Φ) (column c), at V0 = 100, C0 = 0.01, C = 0.1 and indicated C03 [mol/L] values for H2SO4.


For Φ > 0.2, an increase of efficiency of the competitive reaction


is noted. A growth of C03 value causes a small extension of the potential range in the jump region, on the side of higher E values (Figure 9, column b). With an increase in the C03 value, the graphs of pH vs. Φ resemble two almost straight-line segments intersecting at Φeq = 0.2 (Figure 9, column c). However, the pH ranges covered by the titration curves are gradually narrowed (Figure 9, column c).


7. Redox systems with two electron-active elements

In the redox systems considered above, one electron-active element was involved in disproportionation or symproportionation reactions affected by NaOH (in dynamic systems) or water (in static systems). In the HCl (C, V) ⇒ NaIO (C0, V0) system, where HCl was used as disproportionating reagent, the possibility of oxidation of Cl−1 ions was allowed a priori, in which HCl was used as a disproportionating reagent, i.e. chlorine in HCl was treated as an electron-active element. It turned out, however, that the oxidation of Cl−1 occurred only in an extremely small extent, ca. 10−16/10−2 ≈ 10−14 part of Cl−1 ions was oxidized. So, we can consider that it is virtually the disproportionation of IO−1 ions originating from the NaIO.

In this section, we consider the systems where two electron-active elements are factually present.

7.1. Titration in K2Cr2O7 (C, V) ⇒ KI (C0) + H2SO4 (C01) V0 system

Note that V0 mL of titrand D composed of KI (N03) + H2SO4 (N04) + H2O (N05) is titrated with V mL of K2Cr2O7 (N01) + H2O (N02) as titrant T added up to a defined point of the titration. In V0 + V mL of D + T mixture, we have the following species:

H2O(N1),H+1(N2,n2)OH1(N3,n3), K+1(N4,n4),HSO41(N5,n5),SO42(N6,n6),I1(N7,n7),I31(N8,n8),I2(N9,n9),I2(s)(N10,n10),HIO(N11,n11), IO1(N12,n12),HIO3(N13,n13),IO31(N14,n14), H5IO6(N15,n15),H4IO61(N16,n16),H3IO62(N17,n17),HCr2O71(N18,n18),Cr2O72(N19,n19), H2CrO4(N20,n20),HCrO41(N21,n21),CrO42(N22,n22),Cr+3(N23,n23), CrOH+2(N24,n24),Cr(OH)2+1(N25,n25),Cr(OH)41(N26,n26),CrSO4+1(N27,n27).

On this basis, we formulate f1 = f(H), f2 = f(O) and then:

f12=2f(O)f(H)N2+N3+7N5+8N6+N11+2N12+5N13+6N14+7N15+8N16+9N17+13N18+14N19+6N20+7N21+8N22+N24+2N25+4N26+8N27=14N01+ 6N04E64
f0=ChB:N2N3+N4N52N6N7N8N12N14N162N17N182N19N212N22+3N23+2N24+N25N26+N27= 0E65
f3=f(K):  N4=2N01+N03)E66
f4=f(S)=f(SO4):  N5+N6+N27=N04E67
(1)N7+(1/3)3N8+0(N9+N10)+(+1)(N11+N12)+(+5)(N13+N14)+(+7)(N15+N16+N17)+(+6)(2N18+2N19+N20+N21+N22)+(+3)·(N23+N24+N25+N26+N27)=(+6)2N01+(1)N03  E68b
(1)[I1]+(1/3)3[I31]+0(2[I2]+2[I2(s)])+(+1)([HIO][IO1])+(+5)([HIO3]+[IO31])+(+7)([H5IO6]+[H4IO61]+[H3IO62])+(+6)(2[HCr2O71])+2[Cr2O72]+[H2CrO4]+[HCrO41]+[CrO42]) +(+3)([Cr+3]+[CrOH+2]+[Cr(OH)2+1]+[Cr(OH)41]+[CrSO4+1])=((+6)2CV+(1)C0V0)/(V0+V)E69

where C0V0 = 103N03/NA, CV = 103N01/NA. Note that Eq. (69) was obtained only from linear combination of electron-non-active elements (fans) in this system. The balances for electron-active elements are as follows:

f5=f(I): N7+3N8+2N9+2N10+N11+N12+N13+N14+N15+N16+N17=N03E70
f6=f(Cr): 2N18+2N19+N20+N21+N22+N23+N24+N25+N26+N27=2N01E72
2[Cr2O72]+[H2CrO4]+[HCrO41]+[CrO42])+[Cr+3]+[CrOH+2]+[Cr(OH)2+1] + [Cr(OH)41]+[CrSO4+1]=CV/(V0+V)E73

Subtraction of 3∙f(Cr) (Eq. (72)) from Eq. (68a) and further operations give


The simplest/shortest form of GEB, obtained from Eqs. (68), (70) and (72) is the relation

2f(O)f(H)+ ChBf(K)6f(S)+f(I)6f(Cr)=0E75
N8+N9+N10+N11+N12+ 3(N13+N14)+4(N15+N16+N17) =1.5(N23+N24+N25+N26)E75a
[I31]+[I2]+[I2(s)]+[HIO]+[IO1]+3([HIO3]+[IO31])+4([H5IO6]+[H4IO61]+[H3IO62]) =1.5([Cr+3]+[CrOH+2]+[Cr(OH)2+1][Cr(OH)41]+[CrSO4+1])E75b

Note that the numbers (N01, N03) of components forming the system are not involved in Eqs. (75a) and (75b).

Applying the atomic numbers: ZI = 53 for I and ZCr = 24 for Cr, we obtain the linear combination

ZIf(I)+ZCrf(Cr)(2f(O)f(H)+ChBf(K) 6f(S))=0  (ZI+1)N7+(3ZI+1)N8+2ZI(N9+N10)+(ZI1)(N11+N12)+(ZI5)(N13+N14)+(ZI7)(N15+N16+N17)+(ZCr6)(2N18+2N19+N20+N21+N22)+(ZCr3)(N23+N24+N25+N26)=(ZI+1)N03+2(ZCr6)N01E76
(ZI+1)[I1]+(3ZI+1)[I31]+2ZI([I2]+[I2])+(ZI1)([HIO]+[IO1])+(ZI5)([HIO3]+[IO31]) +(ZI7)([H5IO6]+ [H4IO61]+[H3IO62])+(ZCr6)(2[HCr2O71]+2[Cr2O72]+[H2CrO4]+[HCrO41]+[CrO42])+(ZCr3)([Cr+3]+[CrOH+2]+[Cr(OH)2+1]+[Cr(OH41)]+[CrSO4+1])=((ZI+1)C0V0+2(ZCr6)CV)/(V0+V)E77

Eqs. (69), (74), (75b) and (77) (and other linear combinations, as well) are equivalent forms of GEB for this system. Note that Eq. (77) is identical with the one obtained immediately on the basis of the approach I to GEB [4]. The E = E(Φ) and pH = pH(Φ) and some speciation curves for iodine and chromium species are plotted in Figures 10 and 11, where Φ is the fraction titrated.

Figure 10.

The plots K2Cr2O7 (C, V) ⇒ KI (C0) + H2SO4 (C01) V0 system for (a) E = E(Φ) and (b) pH = pH(Φ) functions, at V0 = 100, C0 = 0.01 and C01 values indicated at the corresponding curves.

Figure 11.

Speciation diagrams plotted for iodine (a1, a2, a3) and chromium (b1, b2, b3) species, at indicated C01 values.

For C01 = 0.01, I2, I2(s) and I3−1 are formed in reactions:

Cr2O72+(6,6,9)I1+12H+1+2HSO41=2Cr+3+3(I2(s),I2,I31) +2CrSO4+1+7H2OE79
Cr2O72+(6,6,9)I1+14H+1+ 2SO42= 2Cr+3+3(I2(s),I2,I31)+2CrSO4+1+7H2OE80

where predominating products are involved. Binding the H+1 ions corresponds to the pH increase, which is largest for low C01 value (relatively low buffer capacity of the solution, compare with [53, 54, 5860]). In reactions (78)–(80), protons are consumed and then dpH/dΦ > 0; moreover, dE/dΦ > 0 for Φ > 0.

At C01 = 0.02, E = E(Φ) and [I2(s)] pass through maximum at Φ ca. 0.2. The plot of pH = pH(Φ) shows a slight distortion of the course at Φ = 1/6, and DpH/DΦ > 0 for Φ > 0. The [IO3−1] is comparable with [I2(s)] and [I2]; [I3−1] is small here because [I] < 10−6. The reactions

3(I2,I2(s)) + 5Cr2O72+ 34H+1= 6IO31+10Cr+3+17H2O E81
3(I2,I2(s)) + 5Cr2O72+ 24H+1+10HSO41=6IO31+10CrSO4+1+17H2OE82

are clearly indicated.

At C01 = 0.05, the function E = E(Φ) has a complex course: it first increases, reaches a maximum at Φ = 0.295, then decreases, reaches a minimum at Φ = 1, increases again, passes through a flat maximum at Φ ca. 1.6 and decreases. The curve pH = pH(Φ) breaks at Φ = 1/6 and Φ = 1. The I2(s) exists as the equilibrium solid phase at 0.0934 < Φ < 0.762. A decrease in [I2] value is more expressed at Φ > 1. The ratio [HIO3]/[IO3−1] grows with growth of C01 value. The stoichiometry at Φ = Φeq2 = 1 is described by the reaction


and by reactions where HSO4−1, SO4−2, CrSO4+1 and HIO3 are involved.

At C01 ≥ 0.1, a jump on the E = E(Φ) curve at Φ = 1/6 is clearly marked. The growth of jump at Φ = 1 results from a more significant decrease in the [I2] value at Φ > 1. The E-range covered by the jump at Φ = 1 extends with an increase in the C01 value (Figure 10a). For more details, see Ref. [4].

7.2. A comment

Redox systems with two electron-active elements were widely discussed from the GATES/GEB viewpoint. From the earlier literature, one can recall the simulated systems: KMnO4 (C, V) ⇒ VSO4 (C0) + H2SO4 (C01) V0 [18]; K2Cr2O7 (C, V) ⇒ VSO4 (C0) + H2SO4 (C01) V0 [18]; KMnO4 (C, V) ⇒ FeSO4 (C0) + H2SO4 (C01) [18]; Ce(SO4)2 (C) + H2SO4 (C1) V ⇒ FeSO4 (C0) + H2SO4 (C01) V0 [18]; Cl2 (C, V) ⇒ KI (C0, V0) [19]; KMnO4 (C, V) ⇒ KBr (C0), H2SO4 (C01) V0 [19]. Except the NaOH ⇒ Br2 and NaOH ⇒ HBrO systems, in [16] were simulated also the systems: NaOH (C) + CO2 (C1) V ⇒ I2 (C0) + KI (C01) + CO2 (C02); NaOH (C) + CO2 (C1) V ⇒ I2 (C0) + KI (C01) + CO2 (C02) V0 + CCl4 (V*) (dynamic liquid-liquid extraction system). In [16], an interesting qualitative reaction where MgSO4 solution was added into I2 + KI solution previously alkalinized with an excess of NaOH. This addition causes precipitation of Mg(OH)2 and a decrease of pH value, a growth of solubility of this precipitate and a shift of equilibrium (symproportionation) involved with formation of iodide (I2) that adsorbes on the Mg(OH)2 precipitate turning it red-brown (similarity with starch action in visual iodometric titrations). The dynamic solubility curve was plotted, i.e. the MgSO4 addition is formally treated as titration. The speciation diagrams in all the systems considered were plotted. Among others, it is stated that the disproportionation of iodine (I2(s), I2) occurs according to the scheme (A) 3(I2(s), I2) + 6OH−1 = IO3−1 + 5I−1 + 3H2O, not the scheme (B) (I2(s), I2) + 2OH−1 = IO−1 + I−1 + H2O, as stated in almost all contemporary textbooks. One can calculate [16] that the yield of the reaction (A) is 2.5∙109 times higher than the yield of reaction (B) of the same stoichiometry, 3:6 = 1:2. It is interesting to add that 80 years ≈ 2.5∙109 s. Thus, the yield of the reaction B relative to A is in such a ratio as 1 s in relation to the average length of human’s life.

It should also be added that GATES/GEB enables to verify experimental results. Such a case was mentioned in Ref. [18] to the system, where KIO3 (C0) + HCl (C01) + H2SeO3 (C02) + HgCl2 (C03) was titrated with ascorbic acid (C6H8O6). Authors of the work [61] cited in there committed, among other things, a simple mistake resulting from improper recalculation of numerical values of potentials; for more details, see Refs. [1, 17, 18].

In all instances, full attainable (quantitative+qualitative) physicochemical knowledge was involved in the related algorithms. The quantitative knowledge was related to the complete set of equilibrium constants values, whereas qualitative knowledge was helpful in aspect of metastable and kinetic systems.

7.3. An illustrative presentation of the approach I to GEB

A redox reaction can also participate two or more electron-active elements. In this convention, the approach I to GEB can be perceived as the card game, with electron-active elements as ‘players’, electron-non-active elements—as ‘fans’, and electrons—as ‘money’ transferred between ‘players’, see the picture from Ref. [62]. The ‘players’ provide a common pool of their own electrons into the system in question. The ‘money’ is transferred between ‘players’. As a result of the game, the players’ accounts are changed, while the fans’ accounts remain intact/unchanged, in this convention.

In C0 mol/L aqueous solution of Br2, bromine as one electron-acitve element in Br2 can be perceived as a ‘distributor’ of its own electrons, and H, O as ‘fans’. For example, Eq. (42) involves H and O as ‘fans’, Eq. (43) involves Na as ‘fan’, whereas Eq. (44) does not involve ‘fans’. Generally, ‘fans’ are eliminated after the proper (i.e. indicated above) combination of 2∙f(O) – f(H) with the balances for elements considered as ‘fans’. Further simplification can also result from further combination with the balance for ‘player(s)’. But, in any case, a linear combination of balances related to a redox system does not lead to the identity 0 = 0.

As usually happens in the ‘card game’ practice, the players devote to the game only a part of their cash resources. Similarly, in redox reactions, electrons may participate from the valence shells of atoms of electron-active elements; the electrons from the valence shell of the reductant atoms are transferred onto the valence shell of the oxidant atoms. However, this restriction to the valence electrons is not required here. For example, replacing ZBr in Eq. (46a) by ζBrBr < ZBr), we have


Replacing ZI by ζI and ZCr by ζCrI < ZI, ζCr < ZCr) in Eq. (77), we have


In particular, we can put ζCr = ζCr = 0 in Eq. (86). Obviously, we get the relation


identical to the one obtained from Eq. (68a). If we put ζBr = 0 in Eq. (85), we get the relation

5([HBrO3]+[BrO31])([ HBrO ]+[BrO1] )+ [Br31] + [Br1]=05([HBrO3]+[BrO31])+([ HBrO ]+[BrO1])[Br31][Br1]=0E88

identical to the ones obtained from Eq. (44). This way, we recall the card game without ‘live cash’ but with ‘debt of honor’—in not accidental reference to the title of the thriller novel by T. Clancy; btw, nota bene, the “Debt of honor” was published in 1994, like the papers [15, 16, 18].


8. Conclusions

The GEB formulated according to approach I was named first as electron prebalance and presented, in totally mature form, in three papers issued in 1994 [15, 16, 18] and then followed by further articles and other communications. Currently, in context with the approach II to GEB, it is named as the approach I to GEB, fully equivalent to the approach II to GEB, explicitly related to the law of conservation of H and O.

The linear combination 2⋅f(O) – f(H) is a keystone for the overall thermodynamic knowledge on electrolytic systems. It can be formulated both for non-redox and redox systems, in aqueous, non-aqueous and mixed-solvent systems, with amphiprotic (co)solvent(s) involved. The 2⋅f(O) – f(H) is linearly independent on ChB and other balances, for elements/cores f(Ym) ≠ H, O, in any redox system. For any non-redox system, 2⋅f(O) – f(H) is linearly dependent on those balances. Then, the linear independency/dependency of 2⋅f(O) – f(H) on the other balances is the general criterion distinguishing between redox and non-redox systems. The equation for 2⋅f(O) – f(H), considered as the primary form of GEB, prGEB=2f(O)f(H), is the basis of GEB formulation for redox systems according to approach II to GEB.

Ultimately, GEB, ChB and elemental/core balances are expressed in terms of molar concentrations—to be fully compatible/congruent with expressions for equilibrium constants, interrelating molar concentrations of defined species, on the basis of the mass action law applied to correctly written reaction equation. The mass action law is the one and only chemical law applied in GATES. A complete set of independent (non-contradictory [22]) relations for the equilibrium constants is needed for this purpose.

When compared with the approach I, the approach II to GEB offers several advantages. Although derivation of GEB according to the approach II is more laborious, it enables to formulate this balance without prior knowledge of oxidation numbers for the elements involved frequently in complex components and species of the system. Only the composition (expressed by chemical formula) of components forming the system and composition of the species formed in the system, together with their external charges, are required, i.e. it provides an information sufficient to formulate the GEB; it is the paramount advantage of the approach II to GEB. Anyway, the oxidation number, representing the degree of oxidation of an element in a compound or a species, is commonly perceived as a contractual concept. In this regard, formulation of GEB according to approach II is more useful than the approach I when applied to complex organic species in redox systems of biological origin [6367]. The approach II to GEB is advantageous/desired, inter alia, for redox systems where radical and ion-radical species are formed. What is more, the ‘players’ and ‘fans’, as ones perceived from the approach I to GEB viewpoint, are not indicated a priori within the approach II to GEB. The approach I, considered as a ‘short’ version of GEB, is more convenient when oxidation numbers for all elements of the system are known beforehand. Within the approaches I and II to GEB, 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, where the oxidation number, oxidant, reductant, equivalent mass and stoichiometric reaction notation are the redundant concepts only.

The 1 + K balances composed of ChB and K concentration balances related to equations (not equalities, see e.g. Eq. (40)) for elemental balances form a complete set of equations related to a non-redox system. The set of 1 + K variables in x = (pH, pX1,…, pXK)T is involved in the algorithm applied for calculation purposes. For example, 1 + K = 1 + 2 = 3 balances are related to non-redox systems presented in Examples 1 and 2, with independent variables represented by components of the vector x = (pH, pCu, pSO4)T. In this context, Eq. (35) is considered as equality, not equation.

The 2 + K balances composed of GEB, ChB and K concentration balances related to equations (not equalities) for elemental balances form a complete set of equations related to the redox system. The set of 2 + K variables in x = (E, pH, pX1,…, pXK)T is involved in the algorithm applied for calculation purposes. In this chapter, some results of the simulations of electrolytic redox systems, made according to the GATES/ GEB principles, with use of iterative computer programs, are graphically presented and discussed. The computer simulation realized within GATES with use of iterative computer programs, e.g. MATLAB, provides quite a new quality in knowledge gaining. It enables to follow the details of the process, registered with use of measurable quantities, such as pH and/or potential E. All these calculations are made under assumption that the relevant reactions take place in quasi-static manner under isothermal conditions. The reactions proceeding in the respective systems were formulated under assumption that all equilibrium constants found in the relevant tables and then used in the calculations are correct.

The number of electron-active elements (considered as ‘players’, in terms of the approach I to GEB) in a redox system, is practically unlimited and adapted according to current needs; among others, the systems with one, two, three [13] or four [1, 18] ‘players’ were considered.

The GATES and GATES/GEB (in particular) can be applied for thermodynamic resolution of systems of any degree of complexity. An example is the four-step process involved with iodometric method of copper(+2) determination, considered in detail in Refs. [14, 68], where 47 species are involved in the system with three electron-active elements, eight equations and two equalities, interrelated in 35 independent expressions for equilibrium constants. The systems with similar complexity were resolved. The complexity is limited, however, by factors of the qualitative and quantitative nature, that is, the knowledge of species the knowledge of equilibrium constants interrelating concentration of complex species with concentrations of their constituting parts. The physicochemical data are incomplete or not reliable, in many cases. This problem was raised lately in Ref. [12].

Each of the components of x = (pH, pX1,…, pXK)T is informally ascribed to the corresponding balance: ChB, and the corresponding concentration balances for elements/cores ≠ H, O. Each of the components of x = (E, pH, pX1,…,pXk)T is informally ascribed to the corresponding balance: GEB, ChB, and the corresponding concentration balances for elements/cores ≠ H, O. A unequivocal solution of equations is obtained when the number of independent equations equals to the number of independent variables; it is the ‘iron rule’ obligatory in mathematics; pe = – log[e−1] = A∙E, where 1/A = RT/F∙ln10, pH = –log[H+1], pXi = –log[Xizi].

GATES avoids the necessity of quantitative inferences based on fragile/rachitic chemical reaction notation, involving only some of the species existing in the system; it is only a faint imitation of a true, algebraic notation. From the GATES viewpoint, the ‘stoichiometry’ can be perceived as a mnemonic term only. In calculations, the metastable state is realized by omission of potential products in the related balances, whereas ‘opening’ a reaction pathway in metastable state is based on insertion of possible (from equilibrium viewpoint) products in the related balances [17, 18]. One can also test the interfering effects of different kinds.

All the inferences made within GATES/GEB are based on firm, mathematical (algebraic) foundations. The proposed approach allows us to understand far better all physicochemical phenomena occurring in the system in question and improve some methods of analysis. All the facts testify very well about the potency of simulated calculations made, according to GATES, on the basis of all attainable physicochemical knowledge. Testing the complex redox and non-redox systems with use of iterative computer programs deserves wider popularization among physicochemists and chemists-analysts.

The GATES/GEB is put in context with constructivistic and deterministic principles, and GEB is perceived as the general law of nature, referred to as electrolytic (aqueous media) redox systems. It is proved that stoichiometry of reactions is not a primary concept in chemistry, and its application provides false results, for obvious reasons. From the GATES viewpoint, the stoichiometric reactions are only the basis to formulate the related equilibrium constants. GATES/GEB referred to modeling of redox titration curves in context with earlier approaches to this problem. The GATES/GEB is also presented in three other chapters issued in 2017 within InTech [6870].

The dependency/independency criteria ascribed to 2∙f(O) – f(H) distinguishing between the relevant (non-redox and redox) systems are the properties of the equations obtained from the linear combination of the balances for H and O. Namely, the resulting equation is not independent of non-redox systems, since it is a linear combination of the remaining (charge and concentration) balances, whereas in the case of redox systems, this equation is linearly independent of those balances. This is a general property of nature, independent of the complexity of the system under consideration, which is the electrolytic system. GATES and GATES/GEB, in particular, are clear confirmation of the fact that the nature is mathematically designed and the true laws of nature are mathematical. In other words, the quantitative, mathematical method became the essence of science. To paraphrase a Chinese proverb, one can state that ‘the lotus flower, lotus leaf and lotus seed come from the same root’ [2]. Similarly, the three kinds of balances: GEB, charge and concentration balances come from the same family of fundamental laws of preservation. This compatibility is directly visible from the viewpoint of the approach II to GEB. The equivalent equations for GEB, based on a reliable law of the matter conservation, are 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 as an example of excellent epistemological paradigm.

All earlier (dated from the 1960s) efforts made towards formulation of electrolytic redox systems were only clumsy attempts of resolution of the problem in question, as stated in review papers [2, 1014]. These approaches were slavishly related to the stoichiometric reaction notations, involving only two pairs of indicated species participating in redox reaction; there were usually minor species of the system considered. The species different from those involved in the reaction notation were thus omitted in considerations. Moreover, the charge balance and concentration balances for accompanying substances were also omitted. Theoretical considerations were related to virtual cases, not to real, electrolytic redox systems.


9. Afterword: reductionism

The principles of science gaining are based on a conviction that complex phenomena occurred in nature, e.g. in electrolytic systems, that can be explained in terms of some general laws of the matter conservation; it is the basic assumption of reductionism [2]. These laws are expressed in terms of mathematical equations, valid for the systems of any degree of complexity. Reductionism is very similar to and has its roots in Occam’s razor principle that gives precedence to simplicity, i.e. the explanation which requires the fewest assumptions. The power of reductionism lies in prediction and formulation; it is perceived as a good approximation of the macroscopic world. The knowledge thus obtained is verifiable and based on logical premises. This way the quantitative knowledge gained from the study of relatively simple systems can be synthesized in the knowledge obtainable from more complex systems. From this viewpoint, the knowledge obtained from physicochemical analysis involved, e.g. with determination of the stability constants of complex species formed in a particular system, can be perceived as a ‘stone’ used in construction of more complex systems. In this place, one can recall H. Poincaré who stated that ‘Science is facts; just as houses are made of stones, so science is made of facts’. Any complex species must be equipped with its equilibrium constant value; a qualitative knowledge only (e.g. chemical formula) is insufficient in this respect. However, to construct the knowledge on more complex systems, these stones should be arranged according to a defined scheme (design), based on a set of compatible balances. Closely associated with reductionism is determinism—the philosophy that everything has a cause, and that a particular cause leads to a unique effect.

In any complex system, many particular reactions occur; the resultant reaction is the combination (superposition) of these elementary reactions which occur with different efficiencies that are known only after thorough physicochemical examination of the system in question. Similarly, the vibration in a polyatomic molecule is a superposition of normal vibrations. The physicochemical/thermodynamic knowledge on electrolytic systems is based on equilibrium constant values, referred to as the equilibrium system and all information about possible paths of particular reactions occurred in the system in question. Not all paths of chemical reactions are accessible, under defined conditions of analysis, involved with temperature and/or the presence of catalytic agents. This problem has been raised in Ref. [1], in context with GATES/GEB.

The approach II to GEB, preceded by the approach I to GEB, indicated new insight into redox systems and electrolytic systems in general. This way, the thermodynamic knowledge about redox systems was built in 1992 practically from scratch. The combination 2∙f(O) – f(H) of elemental balances for H and O (approach II to GEB) is the quintessence of the generalized electron balance (GEB) that is the link needed for mathematical/algebraic description of electrolytic redox systems of any degree of complexity within the generalized approach to electrolytic systems as GATES/GEB.

Within GATES, a resolution of electrolytic systems is realized with use of iterative computer programs. Formally, the manner of resolution of this task is even easier than the one based on the formulation of some functional dependencies that require some simplifications, as a rule. The simplifications are not necessary in iterative methods.

Summarizing, the generalized electron balance (GEB), perceived within GATES as GATES/GEB, fulfils all the requirements imposed on reductionism. Its formulation is possible for any electrolytic redox systems—static and dynamic (titration), mono- and multiphase, equilibrium, metastable and non-equilibrium systems, of any degree of complexity, provided that all necessary physicochemical knowledge is attainable. Simply, it is the best possible thermodynamic approach to electrolytic redox systems.

GEB, within the context of GATES/GEB, confirms the validity of some statements expressed in the past:

  • A hidden connection is stronger than an obvious one (Heraclitus)

  • All truths are easy to understand once they are discovered; the point is to discover them (Galileo Galilei)

  • Everything should be made as simple as possible, but not simpler (A. Einstein)

  • Minds are like parachutes: they function only when they are open (J. Dewar)

  • A good theory is the best practical tool (J.C. Maxwell)

  • Nothing is too wonderful to be true if it be consistent with the laws of nature (M. Faraday).

  • A new scientific truth does not triumph by convincing its opponents (…), but rather because its opponents eventually die, and a new generation grows up that is familiar with it (M. Planck).

  • Art of reasoning is nothing more than a well-ordered language

  • Telling the truth in a time of universal deceit is a revolutionary act (G. Orwell)

  • Man learns throughout his life, with the exception of school years

  • Do not believe in any messages just because they are in force for a long time, in many countries. Do not believe in something just because a lot of people say/write it for a long time. Do not accept anything just because someone else said/wrote it. Do not believe in something just because it sounds plausible. Have confidence in what you have accepted as true after a long experience/checking, what brings success to you and others (own apprenticeship).



ChBcharge balance
ELequilibrium law
G (≠ G)Gibbs’ function
Gnumber of balances
GEBgeneralized electron balance
Knumber of concentration balances
NAAvogadro’s number
Rgas constant
T (≠ T)temperature (in K)
Ttransposition sign for a vector
Vvolume [mL] of T
V0volume [mL] of D


  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. InTech-Open Access publisher in the fields of Science, Technology and Medicine; Rijeka, Croatia, 2011. pp. 1–34. Chapter 1.
  2. 2. 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
  3. 3. 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
  4. 4. Michałowska-Kaczmarczyk AM, Toporek M, Michałowski T. Speciation diagrams in dynamic iodide + dichromate system. Electrochimica Acta. 2015;155:217–227
  5. 5. Toporek M, Michałowska-Kaczmarczyk AM, Michałowski T. Symproportionation versus disproportionation in bromine redox systems. Electrochimica Acta. 2015;171:176–187
  6. 6. Michałowski T, Kupiec K, Rymanowski M. Numerical analysis of the Gran methods. A comparative study. Analytica Chimica Acta. 2008;606(2):172–183.
  7. 7. 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
  8. 8. 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.
  9. 9. 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.
  10. 10. Michałowski T. The generalized approach to electrolytic systems: I. Physicochemical and analytical implications. Critical Reviews in Analytical Chemistry. 2010;40(1):2–16
  11. 11. 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
  12. 12. 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.
  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):241–269
  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. Michalowski T. Calculation of pH and potential E for bromine aqueous solution. Journal of Chemical Education. 1994;71:560–562. DOI: 10.1021/ed071p560
  16. 16. Michałowski T, Lesiak A. Acid-Base titration curves in disproportionating redox ystems. Journal of Chemical Education. 1994;71:632–636. DOI: 10.1021/ed071p632
  17. 17. Michałowski T, Rymanowski M, Pietrzyk A. Non-typical Brönsted’s acids and bases. Journal of Chemical Education. 2005;82(3):470–472.
  18. 18. Michałowski T, Lesiak A. Formulation of generalized equations for redox titration curves. Chemia Analityczna (Warsaw). 1994;39:623–637.
  19. 19. Michałowski T, Wajda N, Janecki D. A unified quantitative approach to electrolytic systems. Chemia Analityczna (Warsaw). 1996;41:667–685.
  20. 20. Michałowski T, Ponikvar-Svet M, Asuero AG, Kupiec K. Thermodynamic and kinetic effects involved with pH titration of As(III) with iodine in a buffered malonate system. Journal of Solution Chemistry. 2012;41(3):436–446. DOI: 10.1007/s10953-012-9815-6
  21. 21. 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
  22. 22. 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.
  23. 23. Toporek M, Michałowska-Kaczmarczyk AM, Michałowski T. Disproportionation reactions of HIO and NaIO in static and dynamic systems. American Journal of Analytical Chemistry. 2014;5:1046–1056.
  24. 24. Michałowska-Kaczmarczyk AM, Rymanowski M, Asuero AG, Toporek M, Michałowski T. Formulation of titration curves for some redox systems. American Journal of Analytical Chemistry. 2014;5:861–878.
  25. 25. 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.
  26. 26. Michałowska-Kaczmarczyk AM, Michałowski T, General relation valid for electrolytic systems. Journal of Analytical Sciences, Methods and Instrumentation. 2015;5:74–85.;
  27. 27. 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.
  28. 28. Michałowski T, Pilarski B, Asuero AG, Michałowska-Kaczmarczyk AM. Modeling of acid-base properties in binary-solvent systems. In: Wypych G, editors. Handbook of Solvents, Vol. 1. Properties. Toronto: ChemTec Publishing; 2014. pp. 623–648. Chap. 9.4
  29. 29. Michałowska-Kaczmarczyk AM, Michałowski T. GATES as the unique tool for simulation of electrolytic redox and Non-Redox systems. Journal of Analytical & Bioanalytical Techniques.;;i3915
  30. 30. Meija J, Michałowska-Kaczmarczyk AM, Michałowski T. Redox titration challenge. Analytical and Bioanalytical Chemistry. 2017;409(1):11–13
  31. 31. Meija J, Michałowska-Kaczmarczyk AM, Michałowski T. Solution to redox titration challenge. Analytical and Bioanalytical Chemistry. 2017;409(17)
  32. 32. Michałowski T, Pietrzyk A. Complementarity of physical and chemical laws of preservation in aspect of electrolytic systems (in Polish). Wiadomości Chemiczne. 2007;61:625–640
  33. 33. Michałowski T. Calculations in Analytical Chemistry with Elements of Computer Programming [Internet]. Cracow, 2001. Available from:;
  34. 34. Michałowski T, Pietrzyk A. The generalized electron balance concept. Derivation based on elementary rules of the matter conservation. In: Zuba D, Parczewski A, editors. Chemometrics: Methods and Applications. Publishing House of the Institute of Forensic Research, Cracow, 2006. pp. 415–422.
  35. 35. Michałowski T. Electron balance as the basis of calculations in redox systems (in Polish). In: Use of Information Technology in Academic Teaching of Chemistry. Department of Chemistry, Jagiellonian University, Drukarnia Cyfrowa Eikon Plus,, 2007. pp. 162–169.
  36. 36. Michałowski T. Generalized electron balance (GEB) as a law of preservation for electrolytic redox systems. In: 65th Annual Meeting of the International Society of Electrochemistry. Lausanne, Switzerland. 2014.
  37. 37. Marx D, Tuckerman ME, Hutter J, Parrinello M. The nature of the hydrated excess proton in water. Nature. 1999;397:601–604. DOI: 10.1038/17579
  38. 38. Michałowski T, Pietrzyk A. A thermodynamic study of struvite+water system, Talanta. 2006;68:594–601.
  39. 39. Pietrzyk A, Michałowski T. Computer simulation of the process of struvite MgNH4PO4 6H2O dissolution in aqueous solutions (in Polish). In: Use of Information Technology in Academic Teaching of Chemistry. Department of Chemistry, Jagiellonian University, Drukarnia Cyfrowa Eikon Plus,, 2007. pp.137–146.
  40. 40. 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.
  41. 41. Pietrzyk A, Michałowski T. Conditions of struvite formation in aqueous systems. In: Zuba D, Parczewski A, editors. Chemometrics: Methods and Applications. Department of Chemistry, Jagiellonian University, Drukarnia Cyfrowa Eikon Plus,, 2006. Krak?w; pp. 423–427.
  42. 42. Michałowski T, Asuero AG. Thermodynamic modelling of dolomite behavior in aqueous media. Journal of Thermodynamics. 2012;2012:12. Article ID 723052. DOI: 10.1155/2012/723052.
  43. 43. Michałowski T. Effect of mutual solubility of solvents in multiple extraction. Journal of Chemical Education. 2002;79:1267–1268. DOI: 10.1021/ed079p1267.
  44. 44. Michałowski T, Pilarski B, Dobkowska A, Młodzianowski J. Mathematical modeling and physicochemical studies on acid-base equilibria in binary-solvent systems. Wiadomości Chemiczne. 2010;54:124–154.;
  45. 45. 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.
  46. 46. Asuero AG, Pilarski B, Dobkowska A, Michałowski T. On the isohydricity concept—Some comments. Talanta. 2013;112:49–54.
  47. 47. Michałowski T. Some remarks on acid-base titration curves. Chemia Analityczna. 1981;26:799–813
  48. 48. Inczédy J, Lengyel T, Ure AM. Compendium of Analytical Nomenclature, Definitive Rules, 1997. International Union of Pure and Applied Chemistry. Blackwell Scientific Publications, Oxford 1998 (the Orange Book).
  49. 49. Segré E. Nuclei and Particles. 2nd ed. Reading, MA: Benjamin/Cummings; 1977
  50. 50. Cetnar J. General solution of Bateman equations for nuclear transmutations. Annals of Nuclear Energy. 2006;33:640–645
  51. 51. Massey F, Prentis J. Power law approximations for radioactive decay chains. Applied Mathematics and Computation. 2014;245:135–144
  52. 52. Thibes R, de Oliveira SL. General solution to Bateman’s differential equations with direct index notation. International Journal of Pure and Applied Mathematics. 2014;93(6):879–883
  53. 53. 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:151–187.
  54. 54. Michałowski T, Asuero AG. New approaches in modeling the carbonate alkalinity and total alkalinity. Critical Reviews in Analytical Chemistry. 2012;42:220–244.
  55. 55. 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.
  56. 56. Yu L. Handbook of Analytical Chemistry. Moscow: Mir Publishers; 1975
  57. 57. Wang H, Qiu Y, Schaefer 3rd HF. Pathways for the OH + Br2 → HOBr + Br and HOBr + Br → HBr + BrO reactions. Journal of Physical Chemistry A. 2016;120(5):805–816. DOI: 10.1021/acs.jpca.5b11524
  58. 58. Michałowska-Kaczmarczyk AM, Michałowski T. Dynamic buffer capacity in acid–base systems. Journal of Solution Chemistry. 2015;44(6):1256–1266. DOI: 10.1007/s10953-015-0342-0.
  59. 59. Michałowska-Kaczmarczyk AM, Michałowski T, Asuero AG. Formulation of dynamic buffer capacity for phytic acid. American Journal of Chemistry and Applications. 2015;2(1):5–9
  60. 60. Michałowska-Kaczmarczyk AM, Michałowski T. Application of simms constants in modeling the titrimetric analyses of fulvic acids and their complexes with metal ions. Journal of Solution Chemistry. 2016;45:200–220. DOI: 10_1007-S10953-015-0430-1?locale=pl.
  61. 61. Erdey L, Bodor A, Buzás H. Ascorbinometrische Bromatbestimmung. Fresenius' Zeitschrift für Analytische Chemie. 1952;134(6):412-419. DOI:10.1007/BF00433426
  62. 62. Alajalov C. Commuter Card Game [Internet]. Available from: [Accessed: 15 March 1947]
  63. 63. Nemzer B, Pietrzkowski Z, Spórna A, Stalica P, Thresher W, Michałowski T, Wybraniec S. Betalainic and nutritional profiles of Pigment-Enriched root (Beta vulgaris L.) dried extracts. Food Chemistry. 2011;127(1):42–53.
  64. 64. Wybraniec S, Michałowski T. New pathways of betanidin and betanin enzymatic oxidation. Journal of Agricultural and Food Chemistry. 2011;59(17):9612–9622.
  65. 65. Wybraniec S, Stalica P, Spórna A, Nemzer B, Pietrzkowski Z, Michałowski T. Antioxidant activity of betanidin: Electrochemical study in aqueous media. Journal of Agricultural and Food Chemistry. 2011;59(22):12163–12170.
  66. 66. Wybraniec S, Starzak K, Skopińska A, Szaleniec M, Słupski J, Mitka K, Kowalski P, Michałowski T. Effects of metal cations on betanin stability in aqueous-organic solutions. Food Science and Biotechnology. 2013;22(2):353–363.
  67. 67. Wybraniec S, Starzak K, Skopińska A, Nemzer B, Pietrzkowski Z, Michalowski T. Studies on non-enzymatic oxidation mechanism in neobetanin, betanin and decarboxylated betanins. Journal of Agricultural and Food Chemistry. 2013;61(26):6465–6476.
  68. 68. Michałowska-Kaczmarczyk AM, Spórna-Kucab A, Michałowski T. Solubility Products and Solubility Concepts. in: Descriptive Inorganic Chemistry Researches of Metal Compounds. InTech, Rijeka, Croatia; 2017
  69. 69. Michałowska-Kaczmarczyk AM, Spórna-Kucab A, Michałowski T. A Distinguishing Feature of the Balance 2∙f(O) – f(H) in Electrolytic Systems. The Reference to Titrimetric Methods of Analysis. In: Advances in Titration Techniques. InTech, Rijeka, Croatia, 2017
  70. 70. Michałowska-Kaczmarczyk AM, Spórna-Kucab A, Michałowski T. Principles of Titrimetric Analyses According to Generalized Approach to Electrolytic Systems (GATES). In: Advances in Titration Techniques. InTech, Rijeka, Croatia, 2017

Written By

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

Submitted: 13 December 2016 Reviewed: 20 April 2017 Published: 06 September 2017