Ruthenium redox equilibria: 2. Thermodynamic analysis of disproportionation and comproportionation conditions doi:10.5599/jese.228 135 J. Electrochem. Sci. Eng. 6(1) (2016) 135-143; doi: 10.5599/jese.228 Open Access : : ISSN 1847-9286 www.jESE-online.org Original scientific paper Ruthenium redox equilibria 2. Thermodynamic analysis of disproportionation and comproportionation conditions Igor Povar, Oxana Spinu Institute of Chemistry of the Academy of Sciences of Moldova, 3 Academiei str., MD 2028, Chisinau, Moldova Corresponding Author: ipovar@yahoo.ca; Tel.: +373 22 73 97 36; Fax: +373 22 73 97 36 Received: September 30, 2015; Accepted: February 18, 2016 Abstract The key property of Frost diagram has been confirmed using thermodynamic and linear algebra methods. On the basis of the thermodynamic data, the areas of thermodynamic stability of ruthenium species of different valence states as a function of pH for each degree of oxidation have been determined. Subsequently, based on the diagrams calculated for several values of pH, a narrow ΔpH value is determined, in which the dismutation of appropriate form takes place. Based on thermodynamic analysis, the exact value of the beginning of disproportionation (or comproportionation) is found. Finally, the developed revised Frost diagrams of ruthenium heterogeneous chemical and redox equilibria, as a function of pH and the total concentration of metal ion in solution, have been built. Keywords Disproportionation and comproportionation equilibria; Revised Frost diagram; Reduced Gibbs energy of half reaction of oxidation; Ruthenium soluble and insoluble species. Introduction Predicting redox reactions is not necessarily a simple task. It involves a careful examination of the reaction experimental conditionsIn particular, the influence of a number of such factors, as the medium acidity, the complexation and precipitation processes of the redox species is necessary to examine. Predicting redox reactions with the aid of graphical means involves the knowledge of the predominance areas of the redox couples that participate in the equilibria. The superimposition of http://www.jese-online.org/ mailto:ipovar@yahoo.ca J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 RUTHENIUM REDOX EQUILIBRIA : 2. THERMODYNAMIC ANALYSIS 136 the predominance areas of the two different redox couples permits us to predict the evolution of the global redox reaction qualitatively. The Frost diagrams (FD) are used to represent the disproportionation conditions of ions and allow judging the possibility along with the extent of disproportionation (comproportionation) of the valence states of the element. In this diagram, there is represented the dependence of the reduced standard energy of formation of ions 0 0 f f( ) ( ) /G i G i F   (where F is the Faraday constant equal to 96485 C mol-1) on the degree of oxidation of an element (n) [1]. So, the Gibbs energy change is expressed in eV/mol. FD, being simple to accomplish, characterizes clearly the disproportionation (dismutation) processes of ions in solution. The disproportionation occurs if the value 0 f ( )G i of the analyzed ion is situated above the straight line joining the points 0 f ( )G i of neighboring valence forms on the diagram. In the case of oxygen-containing species (ions, molecules), on the ordinate axis the standard Gibbs free energy of formation of species ( 0 f ( )G i ) minus the Gibbs energy of a number x of water molecules, equal to the number of oxygen atoms in the examined species ( 0 0 0 f f f 2( ) - (H O)G i G x G    ), is placed. When the valence form contains more than one element, which is subject to redox transformations, the 0* f G values is calculated per one atom. The standard redox potential represents the negative slope of the straight line joining two points on the diagram, corresponding to any two valence forms, since the equality is valid: 0 0* f /E G n  Here the n quantity coincides with the degree of oxidation of the ion of interest. The angle of inclination of the curve for the given pair of the valence forms characterizes the ability of interaction with the formation of the products with a lower Gibbs energy. By means of this diagram it is easy to establish whether a particular ion is stable against disproportionation. Therefore, the Frost diagram is a convenient source of understanding complex redox chemistry of any element in its various oxidation states and sometimes under different reaction conditions like pH etc. Details and the mechanics of constructing Frost diagrams can be found in numerous basic chemistry books [1-17]. According to a rigorous line of argument, predictions based on Frost diagrams are not more accurate than those given by only considering the standard potentials of couples. In other words, the predictions given by this strategy are only accurate when the redox species involved in the reaction are in their standard states. But these predictions may sometimes be inverted in other experimental conditions. Because standard potentials involve pH 0, Frost diagrams are implicitly drawn for pH 0. For other pH values, the apparent standard potentials must be used in order to build the diagrams. Actually, they are often drawn for pH 0 and pH 14. Hitherto, the Frost diagrams can be used only under standard conditions, i.e. for semi-quantitative estimations. However, under non-standard, real conditions, a number of factors such as pH, concentrations of soluble forms, participation of each valence state of element in diverse auxiliary chemical reactions (precipitation, complex formation, hydrolysis etc.) can exert influence on the dispropor- tionation-comproportionation processes. The Frost diagrams, easily realizable, clearly represent the processes of ion dismutation. But FD drawn this way, suffer from a number of drawbacks, namely:  The diagrams can be only used under standard conditions. Under real conditions a decisive influence on the possibility of occurring dismutation processes can have such factors as acidity, concentration of soluble species and the formation of solid phases;  In function of the solution pH, metal ions are subject to hydrolysis with the formation of mononuclear and polynuclear hydroxocomplexes. In the second case the degree of disproportionation depends on the concentration of metal ions in the solution [18]; )()( ~ 0 nfiG f  I. Povar et al. J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 doi:10.5599/jese.228 137  Conditions of dismutation reactions may be influenced by complex formation reactions. The goal of this paper is to propose a method to build Frost diagrams under real, non-standard conditions (called here as modified Frost diagrams, MFD), based on rigorous thermodynamic analysis of chemical equilibria in the system Ru–H2O, taken into consideration all the factors mentioned above. Proof of the theorem for necessary conditions of disproportionation processes has been also completed. Theory and calculations Proof of the theorem for necessary conditions of disproportionation processes The disproportionation occurs if the value 0 f ( )G i of the analyzed ion is situated above the straight line joining the points 0 f ( )G i of neighboring valence forms on the 0 f ( ) ( )G i f n  diagram. In the case of oxygen-containing species (ions, molecules), on the ordinate axis the standard Gibbs free energy of formation of ion ( 0 f ( )G i ) minus the Gibbs energy of a number x of water molecules, equal to the number of oxygen atoms in the examined species ( 0* 0 0 f f f 2(H O)G G x G     ) is placed. When the valence form contains more than one element, which is subject to redox transformations, the 0* f G values is calculated per one atom. The standard redox potential represents the negative slope of the straight line joining two points on the diagram, corresponding to any two valence forms, since the equality is valid: 0 0 /fE G n    Here the n quantity coincides with the degree of oxidation of the ion of interest. The angle of inclination of the curve for the given pair of the valence forms characterizes the ability of interaction with the formation of the products with a lower Gibbs energy. By means of this diagram is easy to establish whether a particular ion is stable against disproportionation. The disproportionation occurs when 0* f G of ion lies above the straight line joining the Gibbs energies of two neighboring valence forms. This is explained by the fact that 0* f G of the products of disproportionation reaction corresponds to the point situated at the intersection of this line with the vertical passing through the point corresponding to this ion. In addition, the greater the gain in energy, the less stable the ion towards its disproportionation in solution. This property of diagram can be fundamented using thermodynamic and linear algebra methods. We present here a brief proof. We will examine the redox system formed from the ions of an element, which are in three valence states: a, b and c. Let the following reaction of disproportionation takes place in this system: (c-a)M b+ = (c-b)M a+ + (b-a)M c+ (1) and 0 0 0 0 r f f f (c-b) (a)+(b-a) (c) (c-a) (b)<0G G G G       (2) The reaction (1) is characterized by the following sequence of the standard redox potentials: 0 0 0 a/b a/c b/c> >E E E In the case of reaction (1) occurring, the points a ( 0fG (a), a), b ( 0 f G (b), b) and c ( 0fG (c), c) on the diagram 0 f ( ) ( )G i f n  correspond to the subsequent valence states (Fig. 1). J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 RUTHENIUM REDOX EQUILIBRIA : 2. THERMODYNAMIC ANALYSIS 138 Figure 1. The Frost diagram for the case analyzed in this paper. We will prove that if the point b ( 0 f G (b), b), corresponding to the intermediate valence, lies above the straight line ac joining the a and c points, which characterize respectively the forms with lower and larger valences, then the disproportionation reaction (1) occurs. If the straight line ac is drawn through the points a and c, then the point b is situated above the line. We will draw the line perpendicular to the axis n through the point b. We will denote the intersection point of the straight line perpendicular to straight line ac through x ( 0 f G (x), x). The equation of the straight line can be represented by the expression 0 0 0 0 f f f f ( ) ( ) ( ) ( )G x G a G c G a x b c b          (3) from where we get 0 0 0 f f f( ) ( ) ( ) ( ) ( ) ( )c x G a x a G c c a G x        (4) On the other hand, under the conditions needed for the reaction (1) to occur, from the inequality (2) it follows that 0 0 0 f f f( ) ( ) ( ) ( ) ( ) ( )c b G a b a G c c a G b        (5) Therefore, if the b point is located right above the ac straight line, the spontaneous disproportionation of M b+ to M a+ and M c+ takes place. So, the proof is done! Thermodynamic analysis of disproportionation - comproportionation conditions of Ru in different valence states Fig. 2 shows the disproportionation of Ru(II) in Ru(0) and Ru(III). The dotted line drawn through the points corresponding to Ru(0) and Ru(III) is well below the corresponding 0 f G point of the Ru(II) ion. Similarly, one deduces that species Ru(V), Ru(VI) and Ru(VII) are unstable with respect to disproportionation in Ru(IV) andRu (VII). For example, at disproportionation of Ru (VI), 0 f G comes down to the value corresponding to the cross point on the dotted line. Therefore, the disproportionation of Ru(VI) to Ru(IV) and Ru(VIII) is accompanied by the 0 f G decrease and, as a result, it occurs spontaneously. Higher is the energy gain, the more I. Povar et al. J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 doi:10.5599/jese.228 139 unstable is ion with respect to disproportionation in solution. But, as it was above-mentioned, FD drawn by this way, suffers from a number of weaknesses. The improved version of the Frost diagrams under conditions, different from standard ones, MFD, enables to use the Gibbs energy change of oxidation half reaction of metal to the corresponding valence state under real conditions instead of 0 f G . For example, in the case of Ru(IV) the change in standard Gibbs energy of half reaction: 2+ + 2 2Ru+ 2H O= Ru(OH) +2H + 4e (6) is equal to    0 0 2+ 0 0 2+r f 2 f 2 f 2Ru(OH) 2 (H O) Ru(OH)G G G G       (7) For the calculation of the change in standard Gibbs energy of half reaction (6) under real conditions the equation of isotherm reaction should be applied: 0 0 + 4 r r Rulog [H ]G G C    , where ln10RT F  , 0 Ru C is the ruthenium total concentration in mixture. Within the used approach of MFD, the slope of the line that connects two valence states, constitutes the formal (conditional) redox potential. In the case of formation of the solid phase, RuO2∙H2O(s) according to the equation of reaction + 2 2 2Ru + 4H O = RuO × 2H O (s) + 4H + 4e , the equation of isotherm is 0 + 4 0 r r rlog[H ] 4 pHG G G        . Results and discussion A preliminary step in the construction of MFD is to determine the thermodynamic stability of valence states, in function on the pH solution and 0 Ru C [19] (see Table 1). Table 1. Calculation of the MFD for Ru-H2O at different pH values Ru(II) Ru 2+ 0 < pH < 14.0 Ru(V) Ru2O5(am) 0 < pH < 14.0 Ru(III) Ru 3+ 0 < pH < 1.76 Ru(VI) 2- 4 RuO 0 < pH < 14.0 + 2 Ru(OH) 1.76 < pH < 2.12 Ru(VII) - 4 RuO 0 < pH < 14.0 Ru(OH)3(am) H2O 2.12 < pH < 14.0 Ru(VIII) 2 5H RuO 0 < pH < 11.53 Ru(IV) 2+ 2 Ru(OH) 0 < pH < 1.05 - 5 HRuO 11.5 < pH < 14.0 4+ 4 12 Ru (OH) 1.05 < pH < 3.93 RuO2∙H2O (s) 3.93 < pH < 14.0 Then, the change in Gibbs energy for the respective reactions is calculated. For example, at pH 13, we get 0 r1. Ru(0) : 0 eVG  J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 RUTHENIUM REDOX EQUILIBRIA : 2. THERMODYNAMIC ANALYSIS 140 2+ 0 r2. Ru(II) : Ru = Ru + 2e, 1.558 eVG  0 0 r r Rulog 1.321 eVG G C     + 0 2 3 2 r3. Ru(III): Ru + 4H O = Ru(OH) × H O (am) + 3H + 3e, = 1.894 eVG 0 r r 3 pH = -0.413 eVG G     + 0 2 2 2 r4. Ru(IV): Ru + 4H O = RuO × 2H O (s)+ 4H + 4e, = 2.671 eVG 0 r r 4 pH=-0.405 eVG G     + 0 2 2 5 r5. Ru(V): Ru + 5/2H O = 1/2Ru O (am) +5H + 5e, = 3.840 eVG 0 r r 5 pH = -0.005 eVG G     2- + 0 2 4 r6. Ru(VI): Ru + 4H O = RuO + 8H + 6e, Δ =6.655eVG 0 0 r r Rulog - 8 pH = 0.266 eG G C V     - + 0 2 4 r7. Ru(VII): Ru + 4H O = RuO + 8H + 7e, = 7.241 eVG 0 0 r r Ru+ log - 8 pH = 0.852 eVG G C    - + 0 2 5 r8. Ru(VIII): Ru + 5H O = HRuO + 9H + 8e, = 8.919 eVG 0 0 r r Rulog 9 pH = 1.761 eVG G C      MFD for 0 4 Ru 10C   and 10 -6 mol/L at different pH values are presented on Figs. 2 and 3. In all the cases MFD suffer essential changes with pH variation. On this base the following conclusions can be made: 1. In acidic medium (pH = 0 and 1) stable valence states are: for Ru(III) (Ru 3+ and + 2 Ru(OH) ), Ru(IV) ( 2+ 2 Ru(OH) or 4+ 4 12 Ru (OH) ) and Ru(VIII) (as 2 5 H RuO ). 2. Ru(II) disproportionates according to the equations: 2+ 3+ 3Ru = Ru + 2Ru 2+ + + 2 23Ru +4H O = Ru+ 2Ru(OH) + 4H . 3. Ru(V) at pH = 0 and pH = 1 disproportionates according to the scheme + 2+ 2 5 2 2 2 52Ru O (am) +H O+6H =3Ru(OH) +H RuO (aq) At the same pH values, Ru(VI) and Ru(VII) are instable: 2- + 2+ 4 2 5 2 22RuO + 6H = H RuO + Ru(OH) + H O (pH 1) 2- + 4+ 4 2 5 4 122RuO + 5H = H RuO + 1/4Ru (OH) (pH 3) - + 2+ 4 2 2 5 24RuO + H O+6H = 3H RuO + Ru(OH) (pH 1) - + 4+ 4 2 2 5 4 124RuO + 2H O+ 5H = 3H RuO + 1/4Ru (OH) (pH 3) The diagrams ΔGr(n) can also be used in the presence of complexing agent: 2+ - 2-n 2 2 nRu(IV): Ru(OH) + nCl = Ru(OH) Cl , n=0-5, 3+ - 3-n nRu(III): Ru + nCl = RuCl , n=0-6, 2+ - + Ru(II): Ru +Cl = RuCl . I. Povar et al. J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 doi:10.5599/jese.228 141 In this case ΔGr is a function of 3 variables: 0 0 Ru Cl ,C C and pH. For the construction of diagrams ΔGr = f(n) at different 0 Cl C , it is necessary to maintain constant other two parameters – pH and 0 Ru C . The distribution of different species of Ru(IV), Ru(III) and Ru(II) in function on log [Cl - ] (pH = 0, 0 -6 Ru =10 mol/LC ) is the following: Ru(III) Ru(IV) Ru 3+ log[Cl - ] <-2.17 2+ 2 Ru(OH) -log[Cl ]<-1.39 2+ RuCl --2.17-0.44 0 3RuCl - -0.540.48 3- 6 RuCl - log[Cl ]>0.40 The diagram ΔGr(n) for different 0 0 0 Cl Cl Ru ( ,C C C thus 0 Cl C = [Cl - ]) is shown in Fig.4. Figure 2. Modified Frost diagrams for the system Ru-H2O, 0 -4 C 10 / Ru mol L Figure 3. Change in Gibbs energy of the half reaction of oxidation of metallic Ru to respective valence state versus degree of oxidation (n) for the system Ru–H2O, 0 4 C 10 /Ru mol L   0 1 2 3 4 5 6 7 8 -1 0 1 2 3 4 5 6 7 8 pH = 15 pH = 10 pH = 7 pH = 3 pH = 1 pH = 0  G r, e V n 0 1 2 3 4 5 6 7 8 -1 0 1 2 3 4 5 6 7 8 pH = 13 pH = 10 pH = 7 pH = 8 pH = 1 pH = 0  G r, e V n J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 RUTHENIUM REDOX EQUILIBRIA : 2. THERMODYNAMIC ANALYSIS 142 Figure 4. The diagram ΔGr(n) for the system Ru–Cl - -H2O, 0 6 10 /RuC mol L   Conclusions 1. The development of Frost diagrams for the soluble and insoluble species of ruthenium has been performed. An essential point in the modified method is that instead of the standard Gibbs energy of formation of ions 0 f ( )G i , the use of the Gibbs energy change of the oxidation half reaction of element up to the respective valence state under non-standard conditions has been introduced. 2. A new type of diagrams for ruthenium species, based on the thermodynamic examination of the stability of solid phases and original conditions of mass balance in the heterogeneous systems, has been developed. The conditions of mass balance take into consideration all the chemical forms present in aqueous and solid phases, including inorganic and organic mononuclear and polynuclear forms. 3. Within the developed approach, the effect of major factors influencing on the distribution and concentration of soluble and insoluble species of key components in heterogeneous systems can be analyzed. The developed method of presentation of Frost diagrams removes the previous shortcomings and allows rapid and quantitative problem solving on the possibility of disproportionation of ions in real, non-standard conditions. References [1] D. Shriver, P. Atkins, C. Langford, Inorganic Chemistry, Oxford University Press, Oxford, UK, 1994. [2] J. Warren, T. Tronic, J. Mayer, Chemical reviews 110 (2010) 6961-7001. [3] J. M. Martinez de llarduya, F. Villafane, J. Chem. Ed. 71 (1994) 480. [4] A. Michałowska-Kaczmarczyk, A. Asuero, M. Toporek, T. Michałowski, Critical Reviews in Analytical Chemistry 45 (2015) 241-269. [5] W. Koppenol, D. Stanbury, P. Bounds, Free Radical Biology and Medicine 49 (2010) 317- 322. [6] F. Armstrong, J. Hirst, Proceedings of the National Academy of Sciences 108 (2011) 14049- 14054. [7] J. Burgot, Ionic Equilibria in Analytical Chemistry, Springer, New York, USA, 2012, p. 247- 263. [8] A. Michałowska-Kaczmarczyk, M. Toporek, T. Michałowski, Electrochimica Acta 155 (2015) 217-227. 0 1 2 3 4 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.9 0.2 0 -2.25 n lg  C l -/ m o l L -1) I. Povar et al. J. Electrochem. Sci. Eng. 6(1) (2016) 135-143 doi:10.5599/jese.228 143 [9] H. Eggenkamp, The Geochemistry of Stable Chlorine and Bromine Isotopes, Springer, Berlin Heidelberg, Germany, 2014, p. 21. [10] J. Schneider, H. Jia, J. Muckerman, E. Fujita, Chemical Society Reviews 41 (2012) 2036-2051. [11] J. Ibanez, M. Hernandez-Esparza, C. Doria-Serrano, A. Fregoso-Infante, M. Singh, Environmental Chemistry: Fundamentals, Springer Science & Business Media, New York, USA, 2010. [12] S. Chung, K. Chun, N. Kim, S. Jeong, H. Sim, J. Lee, H. Maeng, Match-Communications in Mathematical and in Computer Chemistry 72 (2014) 389-402. [13] J. Muckerman et al., Inorganic Chemistry 53 (2014) 6904-6913. [14] A. Lodge, M. Lacey, M. Fitt, N. Garcia-Araez, J. Owen, Electrochimica Acta 140 (2014) 168- 173. [15] L. Wang, Y. Hung, N. Shammas, Physicochemical Treatment Processes, Humana Press, New Jork, USA, 2005. [16] J. Kostera, J. McGarry, A. Pacheco, Biochemistry 49 (2010) 8546-8553. [17] D. Smith, Journal of Chemical Education 73 (1996) 1099. [18] I. Povar, Bulletin of Academy of Sciences of the MSSR, Series of Biological and chemical Sciences, 5 (1993) 64-67. [19] I. Povar, O. Spînu, Journal of Electrochemical Science and Engineering (2016) doi: 10.5599/jese.226 © 2016 by the authors; licensee IAPC, Zagreb, Croatia. This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/) http://creativecommons.org/licenses/by/4.0/