A New Pseudoinverse Matrix Method For Balancing Chemical Equations And Their Stability

화학반응식 균형과 안정성을 위한 새로운 유사 역행렬법

Abstract

this work is given a new pseudoniverse matrix method for balancing chemical equations. Here offered method is founded on virtue of the solution of a Diophantine matrix equation by using of a Moore-Penrose pseudoinverse matrix. The method has been tested on several typical chemical equations and found to be very successful for the all equations in our extensive balancing research. This method, which works successfully without any limitations, also has the capability to determine the feasibility of a new chemical reaction, and if it is feasible, then it will balance the equation. Chemical equations treated here possess atoms with fractional oxidation numbers. Also, in the present work are introduced necessary and sufficient criteria for stability of chemical equations over stability of their extended matrices.

이 논문에 화학반응식 균형을 맞추기 위한 새로운 유사 역행렬법이 기술되었다. 여기에 제공된 방법은 Moore-Penrose 유사 역행렬을 사용한 Diophantin 행렬식의 해에 기초를 둔다. 방법은 전형적인 여러 화학반응식에 시험적용되었고 폭넓은 균형연구에서 모든 반응식에 매우 성공적이었다. 이 방법은 아무 제한없이 성공적으로 적용되고, 또한 새로운 화학반응식의 타당성에 대한 검증력도 있고, 만일 새식이 타당하다면 화학식 균형을 이룰 것이다. 여기서 다루어진 화학반응식들은 소수산화수를 지닌 원자를 포함하고 있다. 또한, 화학반응식의 확장된 행렬의 안정성에 대한 화학반응식의 안정성의 필요충분조건을 이 연구에 소개하였다.

1. INTRODUCTION

Chemical equations are the most exploited part of theoretical as well as industrial chemistry. Every chemistry instructor doubtless has his favorite technique of balancing equations. Chemical equations are commonly written for two purposes: to show the probable reaction products, and to show the quantities of materials that enter into reaction. It is fortunate indeed that an extensive and comprehensive knowledge of sub-atomic phenomena is not essential in the balancing of equations. Every student, which has general chemistry as an obligatory subject, is bound to come across balancing chemical equations. It is highly important for chemistry student to be able to write and balance chemical equations readily and accurately. Actually, balancing chemical equations provides an excellent demonstrative and pedagogical example of interconnection between stoichiometrical principles and linear algebra.

A chemical equation is a symbolic representation of a chemical reaction. The substances initially involved in a chemical reaction are called reactants, but the newly formed substances are called the products. The products are new substances with properties that are different from those of reactants. Classically, chemical reactions encompass changes that strictly involve the motion of electrons in the forming and breaking of chemical bonds, although the general concept of a chemical reaction, in particular the notion of a chemical equation, is applicable to transformations of elementary particles.

A chemical equation is not only the shorthand writing of the chemist, but it should be a mental picture of an actual reaction. To the observant investigator, the equation should immediately remind him as to the physical nature and properties of the reactants, viz., color, state, etc., as well as the chemical result and its physical nature. Thus, a great deal of significance should be attached to the writing of chemical equations. The part of chemical mathematics called Stoichiometry deals with the weight relations determined by chemical equations and formulas. According to it, the balancing of chemical equations is important in this area. Since a chemical reaction, when it is feasible, is a natural process, the consequent equation is always consistent. Therefore, we must have a nontrivial solution and we should be able to obtain it assuming its existence. Such an assumption is absolutely valid and does not introduce any error. If the reaction is infeasible, then the only solution is a trivial one, i. e., the all coefficients are equal to zero.

 

2. HISTORICAL SURVEY

During the past years, numerous articles have appeared in scientific literature relative to the balancing of chemical equations. In this section is made a short survey of selected articles for balancing chemical equations. The selection criteria for references were intentionally wide, in order to include a large variety of topics. Balancing chemical equations in the scientific literature is considered from four points of view: mathematical, computational, chemical and pedagogical.

Now, shortly we will describe these glances.

· Jones for the first time in mathematics proposed the general problem for balancing chemical equations.1 Krishnamurthy in his article2 gave an algebraic approach for balancing chemical equations founded on virtue of a generalized matrix inverse. Little bit late Das3 offered a simple scholarly technique, which was discussed in.4,5 One other technique for balancing chemical equation over an integer programming approach is given in,6 while in7 by using of a reflexive generalized inverse matrix is solved the general problem of balancing chemical equations given in.1 The newest mathematical results for balancing chemical equations and their stability over a nonsingular matrix method are obtained in.8 Actually, to date in mathematics and chemistry there are only three strictly formalized consistent mathematical methods for balancing chemical equations, particularly they are the methods given in7,8 and right now presented method in this work, while other so called methods in chemical sense have a limited usage, and they are useful only for particular cases, especially for balancing chemical equations which possess atoms with integer oxidation numbers.

· In chemistry there are many published articles,9-29 which consider the use of computers for balancing chemical equations. All of these computational methods use some commercial softer packet.

· University textbooks of general chemistry generally support the ion-electron technique as basic procedure for balancing chemical equation, because it makes the best use of fundamental chemical principles. Also, some authors advocated other techniques, which involve less algebraic manipulation that may deserve attention - particularly in classes of chemistry and chemical engineering majors.30-56

Several simple chemical equations are solved by elementary algebraic techniques in.45,57-62 Bottomley published the earliest article that makes use of the linear algebra method.63 A set of various modifications, which implement this approach, is documented in.31,44,45,64-66 The case when the chemical equation has no unique solution received considerable attention in the education articles.18,67-79 The equation represents two or more independently occurring reactions can be combined in varying stoichiometric ratios.80,81 Fixed ratios of reagents, observed experimentally in particular cases, are equivalent to a restriction on the coefficients that make a unique solution.18

It is necessary to stress out that balancing chemical equations by inspection is equivalent to using the algebraic method or a computerized matrix algebra approach.82,83 The valence change method30,84-106 and the ion-electron method85-88,92,100,107-113 are also simple algebraic inspection techniques, subjected to exactly the same controls and limitations as the algebraic and matrix methods. Here it is good to emphasis that first Karslake in114 considered balancing of ionic chemical equations. Actually, the technique suggested by Garcia115 can reduces the number of algebraic steps for ion-electron method. Above both mentioned methods - the valence change method and ion-electron method begin by establishing the relative proportions of reagents taking part in separate oxidation and reduction components of a redox reaction. Then, each technique uses a lowest common multiplier to enforce a principle of conservancy - for instance, conservation of oxidation number change in the case of the oxidation number method. Johnson in his article116 defined the equivalent term oxidation stage change on this subject.

Stout in117 presented three redox reactions as puzzles. Each one can be shown as simple redox system, which may easily be balanced using here offered method. After this article was published, the followed other debatable articles with critical accent.118-123

· Balancing chemical equations through the pedagogical point of view is given in the articles.112,124-133 This approach is very interesting for the education of chemical research. A check of the hypothesis that formal reasoning and a sufficiently large mental capacity are required to balance more complex many-step equations is made over a test to determine level of intellectual development, mental capacity, and degree of field dependence/field independence of the students.134

 

3. PRELIMINARIES

Now we will introduce some well known results from the matrix algebra. Throughout, the set of m×n matrices over a field will be denoted by ℝ m×n.

Definition 3. 1. The Moore-Penrose pseudoinverse A+ of a matrix A ∈ℝm×n is the unique matrix satisfying the following criteria: AA+A = A, A+AA+ = A+, (AA+)T = AA+, (A+A)T = A+A.

Suppose the matrix A has full rank, then

A = [aij]n×n and rank A = n⇒ A+ = A-1,

A = [aij]m×n (m ＞ n) and rankA = n⇒A+ = (ATA)-1AT,

A = [aij]m×n (m＜ n) and rankA = m⇒A+ = AT (AAT )-1.

If the matrix A does not have full rank, i.e., A is an m×n matrix and rank A = r＜min(m, n), then pseudoinverse A+ can be constructed from the singular value decomposition A = UDVT by A+ = VD+UT.

Also the Moore-Penrose pseudoinverse A+ of a matrix A, we can derive on this way. There do always exist two matrices C = [cij]m×r and D = [dij]r×n of rank r, such that A = CD. Using these matrices it holds that A+ = DT(DDT)-1(CTC)-1CT.

Pseudoinverse matrix A+ of a matrix A was independently defined by Moore135 and Penrose136. If the matrix A+ satisfies first two equalities of definition 3. 1, then it is called a reflexive generalized inverse of A, denoted by A–. This matrix is not unique.

Definition 3. 2. The characteristic equation of an n×n matrix A is the equation in one variable λ

where det(·) denotes a determinant and I is an n×n identity matrix.

Definition 3. 3. The polynomial

which results from evaluating the determinant (3.1) is the characteristic polynomial of the matrix A.

The above polynomial of degree n≥1 with real coefficients ai (0 ≤ i ≤ n), where an≠0, by the fundamental theorem of algebra has n (not necessarily distinct) roots λ1, λ2, … , λn.

Definition 3. 4. The roots of the characteristic polynomial (3. 2) are precisely the eigenvalues of the matrix A.

Definition 3. 5. The reaction matrix has the following format A¦B = [A¦B]m×(r + s), where A = [aij]m×r(r ＜ m) and B = [bij]m×s (s ＜ m).

Definition 3. 6. An extended matrix

of (r+s )×(r+s) dimension, is singular if det(ExtA¦B) = 0, where A = [aij]m×r (r ＜ m), B = [bij]m×s (s ＜ m), I is the identity matrix of (r + s - m)×(r + s - m) dimension and O is a zero matrix of (r + s - m)×m dimension.

In8 is treated the non-singular case.

Let σ(ExtA¦B) = {λi, 1≤i ≤ k} be the spectrum of ExtA¦B. Let | · | denotes a vector norm in ℝk.

Definition 3. 7. The LozinskiĬ measure μ on ℝn with respect to | · | is defined by

Definition 3. 8. The LozinskiĬ measures of ExtA¦B = [aij]p×p with respect to the three common norms

are

where

stab(ExtA¦B) = max{λ, λ∈σ(ExtA¦B)}

is the stability modulus of ExtA¦B and T denoting transpose.

Definition 3. 9. The extended matrix ExtA¦B is stable if stab(ExtA¦B) ＜ 0.

 

4. MAIN RESULTS

In this section we will give a completely new method for balancing and stability of chemical equations. Given analysis is done for arbitrary chemical equation presented in its general form.

Proposition 4. 1. Any chemical equation may be presented in this form

where xj (1≤ j ≤r) and yj (1≤ j≤s) are unknown rational coefficients, Ψi and Ωi (1≤i≤m) are chemical elements in reactants and products, respectively, aij (1≤i≤m; 1≤ j≤r) and bij (1≤i≤m; 1≤ j≤s) are numbers of atoms of elements Ψi and Ωi, respectively, in j-th molecule.

Proof. Let there exists an arbitrary chemical equation from s distinct elements and m + n molecules

where Ξj=Ψ1a1jΨ2a2j···Ψmamj (1≤ j≤r) and Θj = Ω1b1jΩ2b2j···Ωmbmj (1≤ j≤ s). Then previous expression becomes

Now, if we write the above equation in a compact form, then immediately follows (4.1).

The coefficients xj, yj satisfy three basic principles (corresponding to a closed input-output static model137,138)

· the low of conversation of atoms,

· the low of conversation of mass, and

· the time-independence of the reaction.

Theorem 4. 2. The chemical equation (4.1) can be reduced to the following Diophantine matrix equation

where A=[aij]m×r is a matrix of the reactants, B=[bij]m×s is a matrix of the products, xT=(x1, x2, …, xr) and yT=(y1, y2,…,ys) are column vectors of unknown coefficients and T denotes transpose.

Proof. If we develop the molecules of the reaction (4.1) in an explicit form, then we obtain a matrix of the reactants A

and a matrix of the products B

From the above development we obtain that

and

If we substitute (4.5) and (4.6) into (4.2), follows

or

i.e., the above expression holds if only if

Last equation if we present in a matrix form, actually we obtain (4. 4).

Now we will prove the following result.

Theorem 4. 3. Let A∈ℝm×r and B∈ℝm×s. The matrix Diophantine equation (4. 4) is consistent if and only if for some A+

where

If this is the case, a representation of the general solution is

and

with arbitrary vectors u∈ ℝs×1 and v∈ ℝr×1.

Proof. Matrix equation (4. 4) is consistent if and only if there exist vectors x∈ℝr×1 and y∈ℝs×1 such that

From (4. 13) we have Ax = AA+By, since A(I - A+A) = 0. Thus, from (4. 14) follows

On the other hand, from (4. 12) it follows that

Immediately, from (4.16) and (4.15) follows (4.4).

Remark 4. 4. This theorem generalizes the theorem proved in7. By the above theorem, a century old chemistry problem of balancing chemical equations in a general form is completely solved.

Definition 4. 5. Chemical equation (4.1) is stable if stab(ExtA¦B) ＜ 0.

Lemma 4. 6. For any nonsingular matrix U and any vector norm │ · │, with the induced Lozinski measure μ, │Ux│ defines another vector norm and its induced matrix measure μU is given by

Proof. The proof of this lemma follows directly from the definition 3.7.

Theorem 4. 7. For any matrix ExtA¦B∈ℝ it holds

Proof. The relation (4.18) obviously holds for diagonizable matrices in view of (4.17) and the first two relations in (3.6). Furthermore, the infimum in (4.18) can be achieved if the matrix ExtA¦B is diagonizable. The general case can be shown based on this observation, the fact that ExtA¦B can be approximated by diagonizable matrices in and the continuity of μ(·), which is implied by the property

|μ(ℑ) - μ(ℵ)| ≤ |ℑ - ℵ|.

Remark 4. 8. From the above proof it follows that

stab(ExtA¦B)=inf{μ∞[U(ExtA¦B)U-1], detU≠0}.

The same relation holds if μ is replaced by μ1.

Corollary 4. 9. Let ExtA¦B∈ℝ. Then stab(ExtA¦B)＜0⇔μ(ExtA¦B) ＜ 0 for some LozinskiĬ measure μ on ℝn.

There are more results for stability criteria obtained in the works.139,140

For the matrices A and B given in the theorem 4.2, let rankA¦B=ρ and kerA¦B=κ, where ker denotes the nullity or kernel of matrix A¦B. According to141, the deterministic approach is important, since it enables us to classify the chemical reaction as:

1˚ infeasible when the nullity of the reaction matrix A¦B is zero;

2˚ unique (within relative proportions) when the nullity of the reaction matrix A¦B is one; or

3˚ non-unique when the nullity of the reaction matrix A¦B is bigger than one.

Possible cases of balancing chemical equations are the following

1. If m≥ r + s = ρ, then κ = r + s - ρ = 0, i. e., trivial solution x, y = 0, the reaction is infeasible.

2. If m≥r + s, ρ = r + s - 1, then κ = r + s - ρ = 1, unique solution x, y ＞ 0, i. e., the reaction is feasible and is unique.

In practical terms this means that the general procedure for obtaining these coefficients is to solve the system of linear equations derived from the principles of conservation of matter and charge, applied to the reaction element-by-element.

3. If m≥r + s, ρ ＜ r + s - 1, then κ = r + s - ρ ＞ 1, κ (＞1) linearly independent solutions x, y ＞ 0, i. e., the reaction is feasible and is non-unique.

4. If m ＜ r + s, ρ ≤ m, then κ = r + s - ρ≥1, κ (≥1) linearly independent solutions x, y ＞ 0, i. e., the reaction is feasible and is non-unique.

Last two kinds of the reactions are puzzling in that they exhibit infinite linearly independent solution all of which satisfy the chemical balance, and yet they are not all chemically feasible solutions for a given set of experimental conditions. A unique solution is obtained by imposing a chemical constraint, namely, that reactants have to react only in certain proportions.

 

5. AN APPLICATION OF THE MAIN RESULTS

In this section will be applied above method on many chemical equations for their balancing. All chemical equations balanced here appear first time in professional literature and they are chosen with an intention to be avoided to date all well know chemical equations which were repeated many times in the chemical journals for explanation of certain particular techniques for balancing of some chemical equations using only atoms with integer oxidation numbers.

1˚ First we will consider an infeasible reaction, i. e., the case when the nullity of the reaction matrix is zero.

Example 5. 1. Consider this equation

By the schemes

are determined the matrix of reactants A and the matrix of the products B, an according to it the reaction matrix A¦B obtains the form

The rank of the above matrix is ρ = 11. Since the nullity of the reaction matrix is κ = r + s - ρ = 5 + 6 - 11 = 0, then we have only a trivial solution x, y = 0, that means that the reaction is infeasible.

2˚ Next, we will consider the case when the chemical reaction is feasible and is unique, i. e., the nullity of its reaction matrix is one. Here we will balance many special chemical equations with a goal to show the power of the offered mathematical method.

Example 5. 2. As a next example to illustrate above method, we will balance this chemical equation

From the schemes given below

are determined the matrix of reactants A and the matrix of the products B, i. e.,

The reaction matrix A¦B, according to the definition 3. 5 obtains the form

The rank of the above matrix is ρ = 9. Since the nullity of the reaction matrix is κ = r + s - ρ = 3 + 7 - 9 = 1, then we have only a unique solution, that means that the reaction is feasible. Next, we will determine the solution x, y.

The Moore-Penrose pseudoinverse A+ = (ATA)-1AT of the matrix A is

and the matrix G has this form

688 G = (688 I - AA+) B =

Required vectors y and x, according to (4.12) and (4.13) are

Now immediately follows balanced equation in its conventional form

0.1356805781 Pt(NH3)BrNO3+0.5130378112 Cu(NH3)KNO3+1.5338181820 BeCO3= 0.0348595041 C44H50O14.98+0.5130378112 Cu(NO3)2+0.1356805781 PtO3+0.0681811950 Br1.99NO2+0.2604252849 K1.97O+1.5338181820 BeO+0.2031799611 HNO3.

Since the reaction matrix A¦B has a rectangular format, according to definition 3. 6, we will extend it to a square matrix by adding a new row. To do that, we will add the following vector-row (0, 0, 0, 0, 0, 0, 0, 0, 0, 1) as a 10th row, such that ExtA¦B obtains this form

The eigenvalues of the matrix [ExtA¦B + (ExtA¦B)T]/2 are λ1=35.111, λ2=-33.595, λ3=7.769, λ4=1.630, λ5= 1.140, λ6=0.919, λ7=0.205, λ8=-1.285, λ9=-1.087, λ10= -0.807.

The LozinskiĬ measures of ExtA¦B given by (3.6) with respect to the three common norms (3.5) are

μ∞(ExtA¦B) = max (2, 8, 57, 2.99, 36.98, 2, 2.97, 2, 45, 1) = 57,

μ1(ExtA¦B)=max (10, 10, 5, 108.98, 9, 4, 4.99, 2.97, 2, 6)=108.98,

μ2{[ExtA¦B + (ExtA¦B)T]/2} = 35.111.

Since μ2{[ExtA¦B+(ExtA¦B)T]/2}>0, according to the definition 4.5 the chemical equation (5.2) is unstable.

Example 5. 3. The chemical equation

like a characteristic case will be balanced too.

From the schemes given below

are determined the matrix of reactants A and the matrix of the products B.

The reaction matrix A¦B, according to the definition 3. 5 obtains this form

The rank of the above matrix is ρ = 10. Since the nullity of the reaction matrix is κ = r + s - ρ = 5 + 6 - 10 = 1, then we have only a unique solution, that means that the reaction is feasible. Next, we will determine the solution x, y.

The Moore-Penrose pseudoinverse A+=(ATA)-1AT of the matrix A is

3821829615 A+ =

and the matrix G has this form

3821829615 G = (3821829615 I - AA+)B =

Required vectors y and x, according to (4.12) and (4.13) are

Now immediately follows balanced equation in its conventional form

1.0653080872 Pt(NH3)BrNO3+0.0208499800 [4Fe(CN)3·3Fe(CN)2] + 0.0320088416 C44H50O15 + 1.2365580746 Cu(NO3)2 + 0.9079895101 K2BeO2 =1.2365580746 Cu(NH3)KNO3+0.9079895101 BeCO3 + 0.1459498604 K3.97Fe(CN)6 + 1.0653080872 PtO2 + 0.5435245343 Br1.96NO2 + 1.0866921189 HNO3.

According to the definition 3. 6, extended matrix obtains the form

The eigenvalues of the matrix [ExtA¦B + (ExtA¦B)T]/2 are λ1=59.824, λ2=21.786, λ3=-10.767, λ4=5.648, λ5=-4.030, λ6=-2.045, λ7=-0.940, λ8=-0.044, λ9=1.678, λ10=0.831, λ11=1.059.

The LozinskiĬ measures of ExtA¦B given by (3.6) with respect to the three common norms (3.5) are

μ∞(ExtA¦B) = max (2, 32, 57, 2.96, 39, 8, 69, 2, 6.97, 2, 1) = 69,

μ1(ExtA¦B) = max (10, 43, 109, 9, 5, 10, 5, 16.97, 3, 4.96, 6) = 109,

μ2{[ExtA¦B + (ExtA¦B)T]/2} = 59.824.

Since μ2{[ExtA¦B + (ExtA¦B)T]/2} > 0, according to the definition 4.5 the chemical equation (5.3) is unstable.

Example 5. 4. We will now balance this equation

The matrix of reactants A is

and the matrix of the products B is

The reaction matrix A¦B, according to the definition 3. 5 obtains this form

The rank of the above matrix is ρ = 15. Since the nullity of the reaction matrix is κ = r + s - ρ = 8 + 8 - 15 = 1, then we have only a unique solution, that means that the reaction is feasible. Next, we will determine the solution x, y.

The Moore-Penrose pseudoinverse A+=(ATA)-1AT of the matrix A is

94083728 A+ =

and the matrix G has this form

94083728 G = (94083728 I - AA+)B =

Required vectors y and x, according to (4.12) and (4.13) are

Now immediately follows balanced equation in its conventional form

3.0772746861 KAu(CN)2+0.0710721984 AgRuAuTe8 + 0.3865922972 Fe2(SO4)3 + 0.0339665957 N2Se4 + 0.1358663829 WO3 + 0.1623687648 Na2CO3 + 1.1170850360 H2CO3 + 0.8214763117 HCl = 0.0710721984 [Ru(C10H8N2)3]Cl2·6H2O + 0.1656824131 C4H3AuNa1.96OS7 + 0.0679331915 [WCl4(NSeCl)]2 + 0.7731845945 K3.98Fe(CN)6 + 1.4913322357 Au2O3 + 0.5685775876 TeO3 + 0.0710721984 AgO + 1.0210754230 NO2.

According to the definition 3. 6, extended matrix obtains the form

The eigenvalues of the matrix [ExtA¦B + (ExtA¦B)T]/2 are λ1=28.656, λ2=-24.826, λ3=9.278, λ4=-6.650, λ5=-6.847, λ6=5.301, λ7=4.806, λ8=3.495, λ9=-3.389, λ10=-2.805, λ11=1.297, λ12=-0.937, λ13=0.759, λ14=-0.312, λ15=0.040, λ16=0.134.

The LozinskiĬ measures of ExtA¦B given by (3.6) with respect to the three common norms (3.5) are

μ∞(ExtA¦B) = max (4.98, 5, 44, 19, 2, 2, 9, 3, 10, 37, 6, 3, 3.96, 42, 13, 1) = 44,

μ1(ExtA¦B) = max (6, 11, 17, 6, 4, 6, 6, 2, 81, 17.96, 16, 16.98, 5, 4, 2, 4) = 81,

μ2{[ExtA¦B + (ExtA¦B)T]/2} = 28.656.

Since μ2{[ExtA¦B + (ExtA¦B)T]/2} > 0, according to the definition 4. 5 the chemical equation (5.4) is unstable.

Example 5. 5. Also, we balanced this very interesting chemical equation

According to the definition 3. 6, extended matrix obtains the form

The eigenvalues of the matrix [ExtA¦B + (ExtA¦B)T]/2 are λ1=143.364, λ2=-41.450, λ3=14.630, λ4=-11.150, λ5=11.154, λ6=-6.730, λ7=4.060, λ8=-3.424, λ9=-3.357, λ10=3.432, λ11=-2.112, λ12=2.481, λ13=1.348, λ14=1.453, λ15=0.544, λ16=-0.614, λ17=-0.575, λ18=-0.137, λ19=0.090, λ20=-0.006.

The LozinskiĬ measures of ExtA¦B given by (3.6) with respect to the three common norms (3.5) are

μ∞(ExtA¦B) = max (377, 296, 4, 35, 85, 9, 6, 3, 2, 3, 9.986, 15, 2, 4, 2, 2, 13.97, 4, 4, 1) = 377,

μ1(ExtA¦B) = max (168, 86, 95, 88, 90, 29, 2, 5, 3, 3, 85, 85, 98, 6, 11.986, 3, 3, 7.97, 3, 6) = 168,

μ2{[ExtA¦B + (ExtA¦B)T]/2} = 143.364.

Since μ2{[ExtA¦B+(ExtA¦B)T]/2}>0, according to the definition 4. 5 the chemical equation (5.5) is unstable.

Example 5. 6. An other chemical equation which we balanced is the equation

The reaction matrix A¦B, according to the definition 3.5 has the form

The eigenvalues of the matrix [A¦B+(A¦B)T]/2 are λ1=131.840, λ2=-24.264, λ3=15.402, λ4=-13.779, λ5=-8.097, λ6=4.189, λ7=3.564, λ8=-3.485, λ9=-2.396, λ10=1.459, λ11=-1.347, λ12=-1.023, λ13=0.923, λ14=0.748, λ15=0.449, λ16=-0.722, λ17=-0.513, λ18=0.116, λ19=0.011, λ20=-0.075.

The LozinskiĬ measures of A¦B with respect to the three common norms are

μ∞(A¦B) = max (239, 195, 6, 8, 86, 12, 8, 3, 2, 2, 3, 2.99, 10, 4, 8, 3, 7.96, 3, 3, 3) = 239,

μ1(A¦B) = max (158, 88, 13, 20, 7, 2, 5, 3, 4, 7, 3, 152, 83, 41.96, 5, 2, 3, 4.99, 4, 3) =158,

μ2{[A¦B + (A¦B)T]/2} = 131.840.

Because μ2{[A¦B+(A¦B)T]/2}>0, that means that the chemical equation (5. 6) is unstable.

Example 5. 7. It will be benefit for us if we balance some little bit bigger chemical equation, for instance like this

According to the definition 3.6, extended matrix obtains the form

The eigenvalues of the matrix [ExtA¦B+(ExtA¦B)T]/2 are λ1=112.534, λ2=-45.894, λ3=22.466, λ4=-22.572, λ5=-11.987, λ6=11.234, λ7=-5.757, λ8=-5.367, λ9=5.011, λ10=3.199, λ11=-2.849, λ12=-2.010, λ13=1.948, λ14=-1.549, λ15=1.571, λ16=1.269, λ17=1.042, λ18=-1.070, λ19=0.814, λ20=0.614, λ21=-0.809, λ22=-0.516, λ23=-0.398, λ24=0.154, λ25=-0.091, λ26=0.012.

The LozinskiĬ measures of ExtA¦B given by (3.6) with respect to the three common norms (3.5) are

μ∞(ExtA¦B) = max (222, 210, 3, 26, 84, 3, 39, 5.97, 3, 11, 2, 13, 2, 6, 5, 3, 7, 2, 2, 4, 2, 2, 2, 4, 8, 1) = 222,

μ1(ExtA¦B) = max (137, 16, 5, 6, 3, 68, 5, 4, 8, 17, 4, 5, 4, 13, 138, 1, 151, 12, 6, 4, 14, 8, 8, 23, 7.97, 4) = 151,

μ2{[ExtA¦B + (ExtA¦B)T]/2} = 112.534.

Since μ2{[ExtA¦B+(ExtA¦B)T]/2}>0, according to the definition 4. 5 the chemical equation (5.7) is unstable.

Example 5. 8. Also we balanced this chemical equation

Also, the above equation is unstable.

Example 5. 9. This balanced chemical equation

is unstable too.

3˚ Next, we will consider the case when the chemical reaction is non-unique, i. e., when the nullity of its reaction matrix is bigger than one.

Example 5. 10. Last considered equation in this work will be this equation

From (5. 10) follows this system of linear equations

The general solution of the above system is

where x1, x2, x3, y5 and y6 are arbitrary real numbers.

Now, the balanced equation has a form

where x1, x2, x3, y5 and y6 are arbitrary real numbers.

A particular case of equation (5. 10) for y4 = 0 is solved in8, but there the author balanced this equation on completely different way by using of nonsingular matrix method. So, it is the mathematical consideration of this equation, but for chemical purposes stoichiometric coefficients must be unique.

Now, if we employ here offered method we obtain that only possible chemical equation is this equation

which actually is a particular case of the chemical equation (5.10) for x1=1, x2=-1/3, x3=-3/2, y5=9/10 and y6=-2/5.

Also, in the article78, is balanced a particular case for x1=34, x2=36, x3=8, and y2=y6=0.

 

6. CONCLUSION

The practical power of the matrix procedure as the most general tool for balancing chemical equations is demonstrable. By this method are balanced completely new classes of chemical equations with atoms, which have fractional oxidation numbers. Research shown that employed pseudoinverse matrix method founded on virtue of the solution of the Diophantine matrix equation works perfectly for the all chemical equations.

Actually, this method is unique method both in mathematics and chemistry, which generalize to date the all known methods and particular techniques for balancing chemical equations cited in the bibliography given below. In other words speaking, the mathematical method given here is applicable for the all possible cases for balancing chemical equations, does not matter what kind of atoms they possess – fractional or integer oxidation numbers. Also, this method determine and r + s - ρ linear independent solutions when κ > 1, i.e., when the reaction is feasible and is non-unique. This is the main advantage of the method in relation of other known particular methods and techniques.

For all considered chemical equations which have a unique solution is made stability analysis, and as shown results all of them are unstable. This stability analysis is founded on virtue of the Lozinski measures of extended matrix.