DOI QR코드

DOI QR Code

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

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

  • Published : 2008.06.20

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

Keywords

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.

References

  1. Jones, M. SIAM Rev. 1971, 13, 571
  2. Krishnamurthy, E. V. Int. J. Math. Educ. Sci. Technol. 1978, 9, 323 https://doi.org/10.1080/0020739780090310
  3. Das, S. C. Int. J. Math. Educ. Sci. Technol. 1986, 17, 191 https://doi.org/10.1080/0020739860170209
  4. Herndon, W. C. Int. J. Math. Educ. Sci. Technol. 1987, 18, 745 https://doi.org/10.1080/0020739870180516
  5. Das, S. C. Int. J. Math. Educ. Sci. Technol. 1987, 18, 746
  6. Sen, S. K; Agarwal, H; Sen, S. Math. & Comp. Modeling 2006, 44, 678 https://doi.org/10.1016/j.mcm.2006.02.004
  7. Risteski, I. B. SIAM Problems and Solutions, 2007, 1, http://www.siam.org/journals/problems/downloadfiles/71-025s.pdf
  8. Risteski, I. B. Int. J. Math. Manuscripts, 2007, 1, 180
  9. Brown, J. P; Brown, L. P; Redd, R. M. J. Chem. Educ. 1972, 49, 754 https://doi.org/10.1021/ed049p754
  10. Rosen, A. I. J. Chem. Educ. 1977, 54, 704 https://doi.org/10.1021/ed054p704
  11. Goldberg, S. Z. J. Chem. Educ. 1978, 55, 532 https://doi.org/10.1021/ed055p532
  12. Kennedy, J. H. J. Chem. Educ. 1982, 59, 523 https://doi.org/10.1021/ed059p523
  13. Blakley, G. R. J. Chem. Educ. 1982, 59, 728 https://doi.org/10.1021/ed059p728
  14. Alberty, R. A. J. Chem. Educ. 1983, 60, 102 https://doi.org/10.1021/ed060p102
  15. Loercher, W. J. Chem. Educ. 1986, 63, 74 https://doi.org/10.1021/ed063p74
  16. Cheluget, E. L; Missen, R. W. J. Phys. Chem. 1987, 91, 2428 https://doi.org/10.1021/j100293a044
  17. Missen, R. W; Smith, W. R. J. Chem. Educ. 1989, 66, 217 https://doi.org/10.1021/ed066p217
  18. Missen, R. W; Smith, W. R. J. Chem. Educ. 1990, 67, 876 https://doi.org/10.1021/ed067p876
  19. Alberty, R. A. J. Chem. Educ. 1981, 68, 984 https://doi.org/10.1021/ed068p984
  20. Filgueiras, C. A. L. J. Chem. Educ. 1992, 69, 276 https://doi.org/10.1021/ed069p276
  21. Riley, J; Richmond, Th. G. J. Chem. Educ. 1992, 69, 114
  22. Alberty, R. A. J. Chem. Educ. 1992, 69, 493 https://doi.org/10.1021/ed069p493
  23. Toby, S. J. Chem. Educ. 1994, 71, 270
  24. Weltin, E. J. Chem. Educ. 1994, 71, 295 https://doi.org/10.1021/ed071p295
  25. Wink, D. J. J. Chem. Educ. 1994, 71, 490 https://doi.org/10.1021/ed071p490
  26. Campanario, J. M. Computers Chem. 1995, 19, 85 https://doi.org/10.1016/0097-8485(95)00008-G
  27. Hernandes, M. C. Educ. Quim. 1997, 8, 56
  28. Smith, W. R; Missen, R. W. J. Chem. Educ. 1997, 74, 1369 https://doi.org/10.1021/ed074p1369
  29. Kumar, D. D. J. Sci. Educ. Technol. 2001, 10, 347 https://doi.org/10.1023/A:1012295119087
  30. Endslow, A. W. S. J. Chem. Educ. 1931, 8, 2453 https://doi.org/10.1021/ed008p2453
  31. Scheflan, L. J. Chem. Educ. 1932, 9, 358 https://doi.org/10.1021/ed009p358
  32. Kendall, J. J. Chem. Educ. 1932, 9, 360 https://doi.org/10.1021/ed009p360.1
  33. Clark, A. R. J. Chem. Educ. 1932, 9, 360 https://doi.org/10.1021/ed009p360.1
  34. Dietz, Jr. N. J. Chem. Educ. 1932, 9, 361 https://doi.org/10.1021/ed009p361.1
  35. Otto, E. J. Chem. Educ. 1932, 9, 361 https://doi.org/10.1021/ed009p361.1
  36. Menschutkin, B. J. Chem. Educ. 1932, 9, 560
  37. Slotterbeck, O. J. Chem. Educ. 1932, 9, 751 https://doi.org/10.1021/ed009p751
  38. Peet, B. W. J. Chem. Educ. 1932, 9, 752 https://doi.org/10.1021/ed009p752
  39. Lehrman, L. J. Chem. Educ. 1932, 9, 944 https://doi.org/10.1021/ed009p944
  40. Wernimont, G. J. Chem. Educ. 1932, 9, 1124 https://doi.org/10.1021/ed009p1124
  41. Lehrman, L. J. Chem. Educ. 1932, 9, 1125 https://doi.org/10.1021/ed009p1125
  42. Melville, J. H. J. Chem. Educ. 1932, 9, 1299 https://doi.org/10.1021/ed009p1299
  43. Smith, F. W. J. Chem. Educ. 1933, 10, 250
  44. McNeil, H. C. J. Chem. Educ. 1933, 10, 707 https://doi.org/10.1021/ed010p707
  45. Deming, H. G. J. Chem. Educ. 1934, 11, 125
  46. Crocker, C. J. Chem. Educ. 1968, 45, 731 https://doi.org/10.1021/ed045p731
  47. Copley, G. N. J. Chem. Educ. 1969, 46, 699
  48. Crocker, R. J. Chem. Educ. 1969, 46, 699
  49. Van Cleave, G. Chem 13 News 1974, 62, 12
  50. Robinson, D. A. Chem 13 News 1974, 64, 7
  51. Ludwig, O. G. J. Chem. Educ. 1996, 73, 507 https://doi.org/10.1021/ed073p507
  52. Hart, D. M. J. Chem. Educ. 1996, 73, A226
  53. Hart, D. M. J. Chem. Educ. 1997, 74, 1256 https://doi.org/10.1021/ed074p1256.1
  54. Zaugg, N. S. J. Chem. Educ. 1996, 73, A226
  55. Stout, R. J. Chem. Educ. 1996, 73, A227
  56. Balasubramanian, K. J. Math. Chem. 2001, 30, 219 https://doi.org/10.1023/A:1017979821326
  57. Ferguson, H. W. Chemistry 1967, 40(2), 18
  58. Fabishak, V. L. Chemistry 1967, 40(11), 18
  59. Copley, G. N. Chemistry 1967, 41(9), 22
  60. Swinehart, D. F. J. Chem. Educ. 1985, 62, 55 https://doi.org/10.1021/ed062p55
  61. Harjadi, W. J. Chem. Educ. 1986, 63, 978 https://doi.org/10.1021/ed063p978
  62. Olson, J. A. J. Chem. Educ. 1997, 74, 538 https://doi.org/10.1021/ed074p538
  63. Bottommley, J. Chem. News 1878, 37, 110
  64. Bennett, G. W. J. Chem. Educ. 1954, 31, 324 https://doi.org/10.1021/ed031p324
  65. Haas, R; Gayer, K. H. J. Chem. Educ. 1962, 39, 537 https://doi.org/10.1021/ed039p537
  66. Carrano, S. A. J. Chem. Educ. 1978, 55, 382 https://doi.org/10.1021/ed055p382
  67. Steinbach, O. F. J. Chem. Educ. 1944, 21, 66 https://doi.org/10.1021/ed021p66
  68. Hall, W. T. J. Chem. Educ. 1944, 21, 201 https://doi.org/10.1021/ed021p201
  69. Lehrman, A. J. Chem. Educ. 1944, 21, 202 https://doi.org/10.1021/ed021p202
  70. Porges, A. J. Chem. Educ. 1945, 22, 266 https://doi.org/10.1021/ed022p266
  71. McGavock, W. C. J. Chem. Educ. 1945, 22, 269 https://doi.org/10.1021/ed022p269
  72. Standen, A. J. Chem. Educ. 1945, 22, 461 https://doi.org/10.1021/ed022p461
  73. Swayze, D. R. J. Chem. Educ. 1963, 40, 269 https://doi.org/10.1021/ed040p269.1
  74. Swayze, D. R. Ind. Eng. Chem. Fundam. 1964, 3, 269 https://doi.org/10.1021/i160011a016
  75. Missen, R. W. J. Chem. Educ. 1970, 47, 785
  76. Stark, F. M. J. Chem. Educ. 1984, 61, 476
  77. Porter, S. K. J. Chem. Educ. 1985, 62, 507 https://doi.org/10.1021/ed062p507
  78. Willard, H. H. J. Am. Chem. Soc. 1912, 34, 1480 https://doi.org/10.1021/ja02212a006
  79. Jensen, W. B. J. Chem. Educ. 1987, 64, 646
  80. Aris, R; Math, R. H. S. Ind. Eng. Chem. Fundam. 1963, 2, 90 https://doi.org/10.1021/i160006a002
  81. Schneider, D. R; Reklaitis, G. V. Chem. Eng. Sci. 1975, 30, 243 https://doi.org/10.1016/0009-2509(75)80012-1
  82. Toth, Z. J. Chem. Educ. 1997, 74, 1363 https://doi.org/10.1021/ed074p1363
  83. Guo, Ch. J. Chem. Educ. 1997, 74, 1365 https://doi.org/10.1021/ed074p1365
  84. Griggs, M. A; Warner, M. E. J. Chem. Educ. 1926, 3, 425 https://doi.org/10.1021/ed003p425
  85. Simons, J. H. J. Chem. Educ. 1926, 3, 1305 https://doi.org/10.1021/ed003p1305
  86. Jette, E. R; La Mer V. K. J. Chem. Educ. 1927, 4, 1021 https://doi.org/10.1021/ed004p1021
  87. Jette, E. R; La Mer V. K. J. Chem. Educ. 1927, 4, 1158 https://doi.org/10.1021/ed004p1158
  88. Hall, W. T. J. Chem. Educ. 1929, 6, 479 https://doi.org/10.1021/ed006p479
  89. Reinmuth, O. J. Chem. Educ. 1929, 6, 527 https://doi.org/10.1021/ed006p527
  90. Brinkley, S. R. J. Chem. Educ. 1929, 6, 1894 https://doi.org/10.1021/ed006p1894
  91. Reinmuth, O. J. Chem. Educ. 1930, 7, 1181
  92. Meldrum, W. B. J. Chem. Educ. 1930, 7, 1688 https://doi.org/10.1021/ed007p1688
  93. O. Reinmuth, O. J. Chem. Educ. 1930, 7, 1689 https://doi.org/10.1021/ed007p1689
  94. Bennett, C. W. J. Chem. Educ. 1935, 12, 189 https://doi.org/10.1021/ed012p189
  95. Morris, K. B. J. Chem. Educ. 1938, 15, 538 https://doi.org/10.1021/ed015p538
  96. Wanderwerf, C. A; Davidson, A. W; Sisler, H. H. J. Chem. Educ. 1945, 22, 450 https://doi.org/10.1021/ed022p450
  97. Wanderwerf, C. A. J. Chem. Educ. 1945, 25, 547 https://doi.org/10.1021/ed025p547
  98. Bennett, G. W. J. Chem. Educ. 1954, 31, 157 https://doi.org/10.1021/ed031p157
  99. Loekwood, K. L. J. Chem. Educ. 1961, 38, 326 https://doi.org/10.1021/ed038p326
  100. Kolb, D. J. Chem. Educ. 1978, 55, 326 https://doi.org/10.1021/ed055p326
  101. Kolb, K. J. Chem. Educ. 1979, 56, 181 https://doi.org/10.1021/ed056p181
  102. Barbera, O. Ens. Cien. 1990, 8, 85
  103. Campanario, J. M; Ballesteros, R. Ens. Cien. 1991, 9, 205
  104. Garritz, A; Rincon, C. Educ. Quim. 1996, 7, 190
  105. Ten Hoor, M. J. J. Chem. Educ. 1997, 74, 1367 https://doi.org/10.1021/ed074p1367
  106. Glaister, P. J. Chem. Educ. 1997, 74, 1368 https://doi.org/10.1021/ed074p1368
  107. Waldbauer, L. J; Thurn, W. E. J. Chem. Educ. 1926, 3, 1430 https://doi.org/10.1021/ed003p1430
  108. Lochte, H. L. J. Chem. Educ. 1927, 4, 223 https://doi.org/10.1021/ed004p223
  109. Park, B. J. Chem. Educ. 1929, 6, 1136 https://doi.org/10.1021/ed006p1136
  110. Kolb, D. J. Chem. Educ. 1981, 58, 642-645 https://doi.org/10.1021/ed058p642
  111. Pernicone, E. J. Chem. Educ. 1981, 58, 966
  112. Yarroch, W. L. J. Res. Sci. Teach. 1985, 22, 449 https://doi.org/10.1002/tea.3660220507
  113. Freeman, W. A; Goh, A. K; Chia, L. S; Hart, D. M; Lucas, E. A; Peery, D. J; Subramanian, R; Ten Hoor, M. J. J. Chem. Educ. 1997, 74, 1271
  114. Karslake, W. J. Chem. News 1907, 66, 41
  115. García, A. J. Chem. Educ. 1987, 64, 247 https://doi.org/10.1021/ed064p247
  116. Johnson, O. C. Chem. News 1880, 42, 51
  117. Stout, R. J. Chem. Educ. 1995, 72, 1125 https://doi.org/10.1021/ed072p1125
  118. Nelson, R. J. Chem. Educ. 1997, 74, 1256 https://doi.org/10.1021/ed074p1256.1
  119. Sobkowiak, A. J. Chem. Educ. 1997, 74, 1256 https://doi.org/10.1021/ed074p1256.1
  120. Toth, Z. J. Chem. Educ. 1997, 74, 1270
  121. Woolf, A. A. J. Chem. Educ. 1997, 74, 1270
  122. Ludwig, O. J. Chem. Educ. 1997, 74, 1270
  123. Herndon, W. C. J. Chem. Educ. 1997, 74, 1359 https://doi.org/10.1021/ed074p1359
  124. Zimmermann, J. A. E. J. Chem. Educ. 1925, 2, 383 https://doi.org/10.1021/ed002p383
  125. Stone, Ch. H. J. Chem. Educ. 1944, 21, 550 https://doi.org/10.1021/ed021p550
  126. Brooks, D. W. J. Chem. Educ. 1987, 64, 53 https://doi.org/10.1021/ed064p53
  127. Kumar, D. D. The Agora 1993, 3, 8
  128. Brooks, D. W. J. Chem. Educ. 1995, 72, 233 https://doi.org/10.1021/ed072p233
  129. Cardinali, M. E. J. Chem. Educ. 1995, 72, 716 https://doi.org/10.1021/ed072p716
  130. Subramanian, R; Goh, N. K; Chia, L. S. J. Chem. Educ. 1995, 72, 894 https://doi.org/10.1021/ed072p894
  131. Ferguson, L. A. J. Chem. Educ. 1996, 73, 1129 https://doi.org/10.1021/ed073p1129
  132. Kumar, D. D; White, A. L; Helgeson, S. L. J. Sci. Educ. & Technol. 1994, 3, 187 https://doi.org/10.1007/BF01575179
  133. Kumar, D.D; Helgeson, S. L. Elect. J. Sci. Educ. 2000, 4(4), 1
  134. Niaz, M; Lawson, A. E. J. Res. Sci. Teach. 1985, 22, 41 https://doi.org/10.1002/tea.3660220104
  135. Moore, E. H. Bull. Amer. Math. Soc. 1920, 26, 394
  136. Penrose, R. Proc. Cambridge Phil. Soc. 1955, 51, 406 https://doi.org/10.1017/S0305004100030401
  137. Zeggeren, V. F; Storey, S. H. The Computation of Chemical Equlibria, Cambridge Univ. Press, London 1970
  138. Smith, W. R; Missen, R. W. Chemical Reaction Equilibrium Analysis: Theory and Algorithms, Wiley, New York 1982
  139. Risteski, I. B. Mat. & Technol. 2007, 41, 213
  140. Risteski, I. B. Math. J. Okayama Univ. 2008, 50, 1
  141. Risteski, I. B. 22nd Oct. Meeting of Min. & Metall. Collection of Papers, Oct. 1-2, 1990 Bor, 313

Cited by

  1. A New Topology of Solutions of Chemical Equations vol.57, pp.2, 2013, https://doi.org/10.5012/jkcs.2013.57.2.176
  2. Some Matrix Iterations for Computing Generalized Inverses and Balancing Chemical Equations vol.8, pp.4, 2015, https://doi.org/10.3390/a8040982
  3. A New Singular Matrix Method for Balancing Chemical Equations and Their Stability vol.56, pp.1, 2009, https://doi.org/10.1002/jccs.200900011