Next Article in Journal
Dark Energy, QCD Axion, and Trans-Planckian-Inflaton Decay Constant
Previous Article in Journal
Extra-Dimensional “Metamaterials”: A Model of Inflation Due to a Metric Signature Transition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Unique Mathematical Derivation of the Fundamental Laws of Nature Based on a New Algebraic-Axiomatic (Matrix) Approach †

Logic and Philosophy of Science Research Group, Hokkaido University, Hokkaido 060-0808, Japan
Invited and presented article at The 21st International Conference on General Relativity and Gravitation, Columbia University, New York, NY, USA, 2016. ‡ This work has been done and published during my research fellowship, 2007–2016.
Universe 2017, 3(4), 67; https://doi.org/10.3390/universe3040067
Submission received: 10 July 2017 / Revised: 25 August 2017 / Accepted: 28 August 2017 / Published: 22 September 2017

Abstract

:
In this article, as a new mathematical approach to origin of the laws of nature, using a new basic algebraic axiomatic (matrix) formalism based on the ring theory and Clifford algebras (presented in Section 2), “it is shown that certain mathematical forms of fundamental laws of nature, including laws governing the fundamental forces of nature (represented by a set of two definite classes of general covariant massive field equations, with new matrix formalisms), are derived uniquely from only a very few axioms.” In agreement with the rational Lorentz group, it is also basically assumed that the components of relativistic energy-momentum can only take rational values. In essence, the main scheme of this new mathematical axiomatic approach to the fundamental laws of nature is as follows: First, based on the assumption of the rationality of D-momentum and by linearization (along with a parameterization procedure) of the Lorentz invariant energy-momentum quadratic relation, a unique set of Lorentz invariant systems of homogeneous linear equations (with matrix formalisms compatible with certain Clifford and symmetric algebras) is derived. Then by an initial quantization (followed by a basic procedure of minimal coupling to space-time geometry) of these determined systems of linear equations, a set of two classes of general covariant massive (tensor) field equations (with matrix formalisms compatible with certain Clifford, and Weyl algebras) is derived uniquely as well.
Each class of the derived general covariant field equations also includes a definite form of torsion field appearing as the generator of the corresponding field’ invariant mass. In addition, it is shown that the (1 + 3)-dimensional cases of two classes of derived field equations represent a new general covariant massive formalism of bispinor fields of spin-2, and spin-1 particles, respectively. In fact, these uniquely determined bispinor fields represent a unique set of new generalized massive forms of the laws governing the fundamental forces of nature, including the Einstein (gravitational), Maxwell (electromagnetic) and Yang-Mills (nuclear) field equations. Moreover, it is also shown that the (1 + 2)-dimensional cases of two classes of these field equations represent (asymptotically) a new general covariant massive formalism of bispinor fields of spin-3/2 and spin-1/2 particles, corresponding to the Dirac and Rarita–Schwinger equations.
As a particular consequence, it is shown that a certain massive formalism of general relativity—with a definite form of torsion field appeared originally as the generator of gravitational field’s invariant mass—is obtained only by first quantization (followed by a basic procedure of minimal coupling to space-time geometry) of a certain set of special relativistic algebraic matrix equations. It has been also proved that Lagrangian densities specified for the originally derived new massive forms of the Maxwell, Yang-Mills and Dirac field equations, are also gauge invariant, where the invariant mass of each field is generated solely by the corresponding torsion field. In addition, in agreement with recent astronomical data, a new particular form of massive boson is identified (corresponding to the U(1) gauge symmetry group) with invariant mass: mγ ≈ 4.90571 × 10−50 kg, generated by a coupled torsion field of the background space-time geometry.
Moreover, based on the definite mathematical formalism of this axiomatic approach, along with the C, P and T symmetries (represented basically by the corresponding quantum operators) of the fundamentally derived field equations, it is concluded that the universe could be realized solely by the (1 + 2) and (1 + 3)-dimensional space-times (where this conclusion, in particular, is based on the T-symmetry). It is proved that ‘CPT’ is the only (unique) combination of C, P, and T symmetries that could be defined as a symmetry for interacting fields. In addition, on the basis of these discrete symmetries of derived field equations, it has also been shown that only left-handed particle fields (along with their complementary right-handed fields) could be coupled with the corresponding (any) source currents. Furthermore, it has been shown that the metric of background space-time is diagonalized for the uniquely derived fermion field equations (defined and expressed solely in (1 + 2)-dimensional space-time), where this property generates a certain set of additional symmetries corresponding uniquely to the SU(2)L⊗U(2)R symmetry group for spin-1/2 fermion fields (representing “1 + 3” generations of four fermions, including a group of eight leptons and a group of eight quarks), and also the SU(2)L⊗U(2)R and SU(3) gauge symmetry groups for spin-1 boson fields coupled to the spin-1/2 fermionic source currents. Hence, along with the known elementary particles, eight new elementary particles including four new charge-less right-handed spin-1/2 fermions (two leptons and two quarks), a spin-3/2 fermion, and also three new spin-1 (massive) bosons, are predicted uniquely by this mathematical axiomatic approach. As a particular result, based on the definite formulation of derived Maxwell (and Yang-Mills) field equations, it has also been concluded that magnetic monopoles cannot exist in nature.1

1. Introduction and Summary

Why do the fundamental forces of nature (i.e., the forces that appear to cause all the movements and interactions in the universe) manifest in the way, shape, and form that they do? This is one of the greatest ontological questions that science can investigate. In this article, we’ll consider this basic and crucial question (and a number of relevant issues) via a new axiomatic mathematical formalism. By definition, a basic law of physics (or a scientific law in general) is: “A theoretical principle deduced from particular facts, applicable to a defined group or class of phenomena, and expressible by the statement that a particular phenomenon always occurs if certain conditions be present” [1]. Eugene Wigner’s foundational paper—“On the Unreasonable Effectiveness of Mathematics in the Natural Sciences”—famously observed that purely mathematical structures and formalisms often lead to deep physical insights, in turn serving as the basis for highly successful physical theories [2]. However, all the known fundamental laws of physics (and corresponding mathematical formalisms which are used for their representations), are generally the conclusions of a number of repeated experiments and observations over years and have become accepted universally within the scientific community [3,4]. It should be noted that all the current fundamental field equations that comprise the general formulatory framework of physics are formulated and obtained from the results of experiments (such as the Einstein, Maxwell, Dirac, Schrödinger equations and so on, using in all modern classical and quantum field equations). The extent of their validity is merely the extent to which they correctly predict and agree with experimental results. From a historical point of view, these field equations can be also equivalently written and formulated on the basis of Lagrangian, Hamiltonian, and similar field theoretic formulations in connection with the principle of stationary action [5,6,7,8,9,10].
This article is based on my earlier publications (Refs. [11,12,13,14], Springer, 1996–1998). In this article, as a new mathematical approach to origin of the laws of nature, using a new basic algebraic axiomatic (matrix) formalism based on the ring theory and Clifford algebras (presented in Section 2), “it is shown that certain mathematical forms of fundamental laws of nature, including laws governing the fundamental forces of nature (represented by a set of two definite classes of general covariant massive field equations, with new matrix formalisms), are derived uniquely from only a very few axioms”and, in agreement with the rational Lorentz group, it is also basically assumed that the components of relativistic energy-momentum can only take rational values.
Concerning the basic assumption of the rationality of relativistic energy-momentum, it is necessary to note that the rational Lorentz symmetry group is not only dense in the general form of Lorentz group, but also is compatible with the necessary conditions required basically for the formalism of a consistent relativistic quantum theory [15]. In essence, the main scheme of this new mathematical axiomatic approach to fundamental laws of nature is as follows. First based on the assumption of rationality of D-momentum, by linearization (along with a parameterization procedure) of the Lorentz invariant energy-momentum quadratic relation, a unique set of Lorentz invariant systems of homogeneous linear equations (with matrix formalisms compatible with certain Clifford, and symmetric algebras) is derived. Then by initial quantization (followed by a basic procedure of minimal coupling to space-time geometry) of these determined systems of linear equations, a set of two classes of general covariant massive (tensor) field equations (with matrix formalisms compatible with certain Clifford, and Weyl algebras) is derived uniquely as well. Each class of the derived general covariant field equations also includes a definite form of torsion field appearing as the generator of the corresponding field’s invariant mass. In addition, it is shown that the (1 + 3)-dimensional cases of two classes of derived field equations represent a new general covariant massive formalism of bispinor fields of spin-2, and spin-1 particles, respectively. In fact, these uniquely determined bispinor fields represent a unique set of new generalized massive forms of the laws governing the fundamental forces of nature, including the Einstein (gravitational), Maxwell (electromagnetic) and Yang-Mills (nuclear) field equations. Moreover, it is also shown that the (1 + 2)-dimensional cases of two classes of these field equations represent (asymptotically) a new general covariant massive formalism of bispinor fields of spin-3/2 and spin-1/2 particles, respectively, corresponding to the Dirac and Rarita–Schwinger equations.
Of particular consequence, it is shown that a certain massive formalism of general relativity—with a definite form of torsion field appeared originally as the generator of gravitational field’s invariant mass—is obtained only by initial quantization (followed by a basic procedure of minimal coupling to space-time geometry) of a certain set of special relativistic algebraic matrix equations. It has been also proved that Lagrangian densities specified for the originally derived new massive forms of the Maxwell, Yang-Mills and Dirac field equations, are also gauge invariant, where the invariant mass of each field is generated solely by the corresponding torsion field. In addition, in agreement with recent astronomical data, a new particular form of massive boson is identified (corresponding to U(1) gauge group) with invariant mass: mγ ≈ 4.90571 × 10−50 kg, generated by a coupled torsion field of the background space-time geometry.
Moreover, based on the definite mathematical formalism of this axiomatic approach, along with the C, P and T symmetries (represented basically by the corresponding quantum operators) of the fundamentally derived field equations, it has been concluded that the universe could be realized solely with the (1 + 2) and (1 + 3)-dimensional space-times (where this conclusion, in particular, is based on the T-symmetry). It is proved that ‘CPT’ is the only (unique) combination of C, P, and T symmetries that could be defined as a symmetry for interacting fields. In addition, on the basis of these discrete symmetries of derived field equations, it has been also shown that only left-handed particle fields (along with their complementary right-handed fields) could be coupled to the corresponding (any) source currents. Furthermore, it has been shown that the metric of background space-time is diagonalized for the uniquely derived fermion field equations (defined and expressed solely in (1 + 2)-dimensional space-time), where this property generates a certain set of additional symmetries corresponding uniquely to the SU(2)L⊗U(2)R symmetry group for spin-1/2 fermion fields (representing “1 + 3” generations of four fermions, including a group of eight leptons and a group of eight quarks), and also the SU(2)L⊗U(2)R and SU(3) gauge symmetry groups for spin-1 boson fields coupled to the spin-1/2 fermionic source currents. Hence, along with the known elementary particles, eight new elementary particles, including: four new charge-less right-handed spin-1/2 fermions (two leptons and two quarks, represented by “ze, zn and zu, zd”), a spin-3/2 fermion, and also three new spin-1 massive bosons (represented by W ˜ + , W ˜ , Z , where in particular, the new boson Z is complementary right-handed particle of ordinary Z boson), have been predicted uniquely and expressly by this new mathematical axiomatic approach.
As a particular result, presented in Section 3.6, based on the definite and unique formulation of the derived Maxwell’s equations (and also determined Yang-Mills equations, represented uniquely with two specific forms of gauge symmetries, in Section 3.15), it has also been concluded generally that magnetic monopoles cannot exist in nature.

1.1. The Main Results Obtained in This Article are Based on the Following Three Basic Assumptions (as Postulates)

(1)- “A new definite axiomatic generalization of the axiom of “no zero divisors” of integral domains (including the ring of integers ℤ);”
This algebraic postulate (as a new mathematical concept) is formulated as follows:
“Let A = [ a i j ] be a n × n matrix with entries expressed by the following linear homogeneous polynomials in s variables over the integral domain ℤ: a i j = a i j ( b 1 , b 2 , b 3 , , b s ) = k = 1 s H i j k b k ; suppose also “   r ℕ: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n ”, where F ( b 1 , b 2 , b 3 , , b s ) is a homogeneous polynomial of degree r ≥ 2, and I n is n × n identity matrix; Then the following axiom is assumed (as a new axiomatic generalization of the ordinary axiom of “no zero divisors” of integral domain ℤ):
( A r = 0 ) ( A × M = 0 , M 0 )
where M is a non-zero arbitrary n × 1 column matrix”.
The axiomatic relation (1) is a logical biconditional, where ( A r = 0 ) and ( A × M = 0 , M 0 ) are respectively the antecedent and consequent of this biconditional. In addition, based on the initial assumption   r ℕ: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n , the axiomatic biconditional (1) could be also represented as follows:
[ F ( b 1 , b 2 , b 3 , , b s ) = 0 ] ( A × M = 0 , M 0 )
where the homogeneous equation F ( b 1 , b 2 , b 3 , , b s ) = 0 , and system of linear equations ( A × M = 0 , M 0 ) are respectively the antecedent and consequent of biconditional (1-1). The axiomatic biconditional (1-1), defines a system of linear equations of the type A × M = 0 ( M 0 ) , as the algebraic equivalent representation of r t h degree homogeneous equation F ( b 1 , b 2 , b 3 , , b s ) = 0 (over the integral domain ℤ). In addition, according to the Ref. [16], since F ( b 1 , b 2 , b 3 , , b s ) = 0 is a homogeneous equation over ℤ, it is also concluded that homogeneous equations defined over the field of rational numbers ℚ, obey the axiomatic relations (1) and (1-1) as well. As particular outcome of this new mathematical axiomatic formalism (based on the axiomatic relations (1) and (1-1), including their basic algebraic properties), in Section 3.4, it is shown that using a unique set of general covariant massive (tensor) field equations (with new matrix formalism compatible with Clifford, and Weyl algebras), corresponding with the fundamental field equations of physics, are derived—where, in agreement with the rational Lorentz symmetry group, it has been basically assumed that the components of relativistic energy-momentum can only take the rational values. In Section 3.2, Section 3.3, Section 3.4, Section 3.5, Section 3.6, Section 3.7, Section 3.8, Section 3.9, Section 3.10, Section 3.11, Section 3.12, Section 3.13, Section 3.14 and Section 3.15, we present in detail the main applications of this basic algebraic assumption (along with the following basic assumptions (2) and (3)) in fundamental physics.
(2)- “In agreement with the rational Lorentz symmetry group, we assume basically that the components of relativistic energy-momentum (D-momentum) can only take the rational values.”
Concerning this assumption, it is necessary to note that the rational Lorentz symmetry group is not only dense in the general form of the Lorentz group, but is also compatible with the necessary basic conditions required for the formalism of a consistent relativistic quantum theory [15]. Moreover, this assumption is clearly also compatible with any quantum circumstance in which the energy-momentum of a relativistic particle is transferred as integer multiples of the quantum of action “h” (Planck constant).
Before defining the next basic assumption, it should be noted that from basic assumptions (1) and (2), it follows directly that the Lorentz invariant energy-momentum quadratic relation (represented by Formula (52), in Section 3.1) is a particular form of homogeneous quadratic equation (represented by Formula (18-2) in Section 2.2). Hence, using the set of systems of linear equations that are determined uniquely as equivalent algebraic representations of the corresponding set of quadratic homogeneous equations (given by Equation (18-2) in various numbers of unknown variables, respectively), a unique set of the Lorentz invariant systems of homogeneous linear equations (with matrix formalisms compatible with certain Clifford, and symmetric algebras) are also determined, representing equivalent algebraic forms of the energy-momentum quadratic relation in various space-time dimensions, respectively. Subsequently, we have shown that by an initial quantization (followed by a basic procedure of minimal coupling to space-time geometry) of these determined systems of linear equations, a unique set of two definite classes of general covariant massive (tensor) field equations (with matrix formalisms compatible with certain Clifford, and Weyl algebras) is also derived, corresponding with various space-time dimensions, respectively. In addition, it is also shown that this derived set of two classes of general covariant field equations represent new tensor massive (matrix) formalism of the fundamental field equations of physics, corresponding with fundamental laws of nature (including the laws governing the fundamental forces of nature). Following these essential results, in addition to the basic assumptions (1) and (2), it would also be assumed that:
(3)- “We assume that the mathematical formalism of the fundamental laws of nature, are defined solely by the axiomatic matrix constitution formulated uniquely on the basis of postulates (1) and (2)”.
In addition to this basic assumption, in Section 3.11, the C, P and T symmetries of the uniquely derived general covariant field equations (that are field Equations (3) and (4) in Section 1.2), would be represented by their corresponding quantum matrix operators.

1.2. A Summary of the Main Consequences of Basic Assumptions (1)–(3)

In the following, we present a summary description of the main consequences of basic assumptions (1)–(3) (mentioned in Section 1.1) in fundamental physics. In this article, the metric signature (+ − ... −), the geometrized units [17] and also the following sign conventions have been used in the representations of the Riemann curvature tensor R σ μ ν ρ , Ricci tensor R μ ν and Einstein tensor G μ ν :
R σ μ ν ρ = ( ν Γ σ μ ρ + Γ λ ν ρ Γ σ μ λ ) ( μ Γ σ ν ρ + Γ λ μ ρ Γ σ ν λ ) , σ R μ ν ρ σ = ν R μ ρ μ R ν ρ , G μ ν = 8 π Τ μ ν + .
On the basis of assumptions (1)–(3), two sets of the general covariant field equations (compatible with the Clifford algebras) are derived solely as follows:
( i α μ μ m 0 ( R ) α ˜ μ k μ ) Ψ R = 0
( i α μ D μ m 0 ( F ) α ˜ μ k μ ) Ψ F = 0
where
α μ = β μ + β μ , α ˜ μ = β μ β μ
i μ and i D μ are the general relativistic forms of energy-momentum quantum operator (where μ is the general covariant derivative and D μ is gauge covariant derivative, defining in Section 3.4, Section 3.5 and Section 3.6), m 0 ( R ) and m 0 ( F ) are the fields’ invariant masses, k μ = ( c / g 00 , 0 , , 0 ) is the general covariant velocity in stationary reference frame (that is a time-like covariant vector), β μ and β μ are two contravariant square matrices (given by Formula (6) and (7)), Ψ R is a column matrix given by Formula (6) and (7), which contains the components of field strength tensor R μ ν ρ σ (equivalent to the Riemann curvature tensor), and also the components of a covariant quantity which defines the corresponding source current (by relations (6) and (7)), Ψ F is also a column matrix given by Formula (6) and (7), which contains the components of tensor field F μ ν (defined as the gauge field strength tensor), and also the components of a covariant quantity that defines the corresponding source current (by relations (6) and (7)). In Section 3.11, based on a basic class of discrete symmetries of general covariant field Equations (3) and (4), it would be concluded that these equations could be defined solely in (1 + 2) and (1 + 3) space-time dimensions, where the (1 + 2) and (1 + 3)-dimensional cases these field equations are given uniquely as follows (in terms of the above-mentioned quantities), respectively:
For (1 + 2)-dimensional space-time we have:
β 0 = [ 0 0 0 ( σ 0 + σ 1 ) ] , β 0 = [ σ 0 + σ 1 0 0 0 ] , β 1 = [ 0 σ 2 σ 2 0 ] , β 1 = [ 0 σ 3 σ 3 0 ] ,
β 2 = [ 0 σ 1 σ 0 0 ] , β 2 = [ 0 σ 0 σ 1 0 ] , σ 0 = [ 1 0 0 0 ] , σ 1 = [ 0 0 0 1 ] , σ 2 = [ 0 1 0 0 ] , σ 3 = [ 0 0 1 0 ] ,
Ψ R = [ R ρ σ 10 0 R ρ σ 21 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 0 F 21 ϕ ( F ) ] , J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ;
For (1 + 3)-dimensional space-time of we get:
β 0 = [ 0 0 0 ( γ 0 + γ 1 ) ] , β 0 = [ ( γ 0 + γ 1 ) 0 0 0 ] , β 1 = [ 0 γ 2 γ 3 0 ] , β 1 = [ 0 γ 3 γ 2 0 ] ,
β 2 = [ 0 γ 4 γ 5 0 ] , β 2 = [ 0 γ 5 γ 4 0 ] , β 3 = [ 0 γ 6 γ 7 0 ] , β 3 = [ 0 γ 7 γ 6 0 ] ,
γ 0 = [ 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 ] , γ 1 = [ 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 ] , γ 2 = [ 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 ] , γ 3 = [ 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 ] ,
γ 4 = [ 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 ] , γ 5 = [ 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 ] , γ 6 = [ 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 ] , γ 7 = [ 0 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 ] ,
Ψ R = [ R ρ σ 10 R ρ σ 20 R ρ σ 30 0 R ρ σ 23 R ρ σ 31 R ρ σ 12 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 F 20 F 30 0 F 23 F 31 F 12 ϕ ( F ) ] , J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ;
In Formula (6) and (7), J ν ρ σ ( R ) and J ν ( F ) are the covariant source currents expressed necessarily in terms of the covariant quantities ϕ ρ σ ( R ) and ϕ ( F ) (as initially given quantities). Moreover, in Section 3.4, Section 3.5, Section 3.6, Section 3.7, Section 3.8, Section 3.9, Section 3.10, Section 3.11, Section 3.12, Section 3.13, Section 3.14 and Section 3.15, it has been shown that the field equations in (1 + 2) dimensions, are compatible with the matrix representation of Clifford algebra Cℓ1,2, and represent (asymptotically) a new general covariant massive formalism of bispinor fields of spin-3/2 and spin-1/2 particles, respectively. It has also been shown that these field equations in (1 + 3) dimensions are compatible with the matrix representation of Clifford algebra Cℓ1,3, and represent solely a new general covariant massive formalism of bispinor fields of spin-2 and spin-1 particles, respectively.
In addition, from the field Equations (3) and (4), the following field equations (with ordinary tensor formulations) could be also obtained, respectively:
λ R ρ σ μ ν + μ R ρ σ ν λ + ν R ρ σ λ μ = T λ μ τ R ρ σ τ ν + T μ ν τ R ρ σ τ λ + T ν λ τ R ρ σ τ μ ,
μ R ρ σ μ ν ( i m 0 ( R ) / ) k μ R ρ σ μ ν = J ρ σ ( R ) ν ;
R σ μ ν ρ = ( ν Γ σ μ ρ + Γ λ ν ρ Γ σ μ λ ) ( μ Γ σ ν ρ + Γ λ μ ρ Γ σ ν λ ) ,
J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , T τ μ ν = i m 0 ( R ) 2 ( g τ μ k ν g τ ν k μ ) .
and
D λ F μ ν + D μ F ν λ + D ν F λ μ = 0 ,
D μ F μ ν = J ν ( F ) ;
F μ ν = D ν A μ D μ A ν , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ,   Z τ μ ν = i m 0 ( F ) 2 ( g τ μ k ν g τ ν k μ ) .
where in Equations (3-1)–(3-2), Γ σ μ ρ is the affine connection given by: Γ σ μ ρ = Γ ¯ σ μ ρ K σ μ ρ , Γ ¯ σ μ ρ is the Christoffel symbol (or the torsion-free connection), K σ μ ρ is a contorsion tensor defined by: K ρ σ μ = ( i m 0 ( R ) / 2 ) g ρ μ k σ (that is anti-symmetric in the first and last indices), T ρ σ μ is its corresponding torsion tensor given by: T ρ σ μ = K ρ μ σ K ρ σ μ (as the generator of the gravitational field’s invariant mass), μ is general covariant derivative defined with torsion T ρ σ μ . In Equations (4-1)–(4-3), D μ is the general relativistic form of gauge covariant derivative defined with torsion field Z τ μ ν (which generates the gauge field’s invariant mass), and A μ denotes the corresponding gauge (potential) field.
In Section 3.11, on the basis of definite mathematical formalism of this axiomatic approach, along with the C, P and T symmetries (represented basically by the corresponding quantum operators, in Section 3.11) of the fundamentally derived field equations, it has been concluded that the universe could be realized solely with the (1 + 2) and (1 + 3)-dimensional space-times (where this conclusion, in particular, is based on the T-symmetry). It is proved that ‘CPT’ is the only (unique) combination of C, P, and T symmetries that could be defined as a symmetry for interacting fields. In addition, on the basis of these discrete symmetries of derived field equations, it has been also shown that only left-handed particle fields (along with their complementary right-handed fields) could be coupled to the corresponding (any) source currents. Furthermore, it has been shown that the metric of background space-time is diagonalized for the uniquely derived fermion field equations (defined and expressed solely in (1 + 2)-dimensional space-time), where this property generates a certain set of additional symmetries corresponding uniquely to the SU(2)L⊗U(2)R symmetry group for spin-1/2 fermion fields (representing “1 + 3” generations of four fermions, including a group of eight leptons and a group of eight quarks), and also the SU(2)L⊗U(2)R and SU(3) gauge symmetry groups for spin-1 boson fields coupled to the spin-1/2 fermionic source currents. Hence, along with the known elementary particles, eight new elementary particles, including four new charge-less right-handed spin-1/2 fermions (two leptons and two quarks, represented by “ze, zn and zu, zd”), a spin-3/2 fermion, and also three new spin-1 massive bosons (represented by W ˜ + , W ˜ , Z , where in particular the new boson Z is complementary right-handed particle of ordinary Z boson), have been predicted uniquely by this new mathematical axiomatic approach (as shown in Section 3.15).
As a consequence, in Section 3.6 it is shown that a certain massive formalism of the general theory of relativity—with a definite torsion field which generates the gravitational field’s mass—is obtained only by first quantization (followed by a basic procedure of minimal coupling to space-time geometry) of a set of special relativistic algebraic matrix relations. In Section 3.9, it is also proved that Lagrangian densities specified for the derived unique massive forms of Maxwell, Yang-Mills and Dirac equations, are gauge-invariant as well, where the invariant mass of each field is generated by the corresponding torsion field. In addition, in Section 3.10, in agreement with recent astronomical data, a new massive boson is identified (corresponding to U(1) gauge group) with invariant mass: mγ ≈ 4.90571 × 10−50 kg, generated by a coupling torsion field of the background space-time geometry. Furthermore, in Section 3.6, based on the definite and unique formulation of the derived Maxwell’s equations (and also determined Yang-Mills equations, represented unqiely with two specific forms of gauge symmetries), it is also concluded that magnetic monopoles cannot exist in nature.
As shown in Section 3.7, if the Ricci curvature tensor R μ ν is defined by the following relation in terms of Riemann curvature tensor (which is determined by field Equations (3-1)–(3-3)):
( σ + i m 0 ( R ) k σ ) R μ ν ρ σ = ( ν + i m 0 ( R ) k ν ) R μ ρ ( μ + i m 0 ( R ) k μ ) R ν ρ ,
then from this expression for the current in terms of the stress-energy tensor T μ ν :
J ρ σ ν ( R ) = 8 π [ ( σ + i m 0 ( R ) k σ ) T ρ ν ( ρ + i m 0 ( R ) k ρ ) T σ ν ] + 8 π B [ ( σ + i m 0 ( R ) k σ ) T g ρ ν ( ρ + i m 0 ( R ) k ρ ) T g σ ν ]
where T = T μ μ , the gravitational field equations (including a cosmological constant Λ emerged naturally in the course of derivation process) could be equivalently derived from the massless case of tensor field Equations (3-1)–(3-3) in (1 + 3) space-time dimensions, as follows:
R μ ν = 8 π T μ ν + 4 π T g μ ν Λ g μ ν
Let us emphasize again that the results obtained in this article are the direct outcome of a new algebraic-axiomatic approach,2 presented in Section 2. This algebraic approach—in the form of a basic linearization theory—has been constructed on the basis of a new single axiom (that is the axiom (17) in Section 2.1) proposed to replace with the ordinary axiom of “no zero divisors” of integral domains (that is the axiom (16) in Section 2). In fact, as noted in Section 1.1 and Section 2.1, the new proposed axiom is a definite generalized form of ordinary axiom (16), which has been formulated in terms of square matrices (using basically as primary objects for representing the elements of underlying algebra, i.e., integral domains including the ring of integers). In Section 3, based on this new algebraic axiomatic formalism, as a new mathematical approach to origin of the laws of nature, “it is shown that certain mathematical forms of fundamental laws of nature, including laws governing the fundamental forces of nature (represented by a set of two definite classes of general covariant massive field equations, with new matrix formalisms), are derived uniquely from only a very few axioms” in which, in agreement with the rational Lorentz group, it is also assumed that the components of relativistic energy-momentum can only take rational values.

2. Theory of Linearization: A New Algebraic-Axiomatic (Matrix) Formalism Based on the Ring Theory and Clifford Algebras

In this Section, a new algebraic theory of linearization (including the simultaneous parameterization) of the homogeneous equations is presented that is formulated on the basis of ring theory and matrix representation of the generalized Clifford algebras (associated with homogeneous forms of degree r ≥ 2 defined over the integral domain ℤ).
Mathematical models of physical processes include certain classes of mathematical objects and the relations between these objects. The models of this type, which are most commonly used, are groups, rings, vector spaces, and linear algebras. A group is a set G with a single operation (multiplication) a × b = c ; a , b , c G which obeys the known conditions [18,19]. A ring is a set of elements R, where two binary operations, namely, addition and multiplication, are defined. With respect to addition this set is a group, and multiplication is connected with addition by the distributivity laws: a × ( b + c ) = ( a × b ) + ( a × c ) , ( b + c ) × a = ( b × a ) + ( c × a ) ; a , b , c R . The rings reflect the structural properties of the set R. As distinct from the group models, those connected with rings are not frequently applied, although in physics various algebras of matrices, algebras of hyper-complex numbers, Grassman and Clifford algebras are widely used. This is due to the intricacy of finding a connection between the binary relations of addition and multiplication and the element of the rings [12,22,23]. This Section is devoted to the development of a rather simple approach of establishing such a connection and an analysis of concrete problems on this basis.
I have found that if the algebraic axiom of “no zero divisors” of integral domains is generalized expressing in terms of the square matrices (as it has been formulated by the axiomatic relation (17)), fruitful new results hold. In this Section, first on the basis of the matrix representation of the generalized Clifford algebras (associated with homogeneous polynomials of degree r ≥ 2 over the integral domain ℤ), we have presented a new generalized formulation of the algebraic axiom of “no zero divisors” of integral domains. Subsequently, a linearization theory has been constructed axiomatically that implies (necessarily and sufficiently) any homogeneous equation of degree r ≥ 2 over the integral domain ℤ, should be linearized (and parameterized simultaneously), and its solution investigated systematically via its equivalent linearized-parameterized formulation (representing as a certain type of system of linear homogeneous equations). In Section 2.2 and Section 2.4, by this axiomatic approach a class of homogeneous quadratic equations (in various numbers of variables) over ℤ has been considered explicitly.

2.1. The Basic Properties of the Integral Domain ℤ

The ordinary basic properties of the integral domain ℤ with binary operations ( + , × ) are represented as follows, respectively [22,23] ( a i , a j , a k ,   ℤ):
Closure:
a k + a l ,   a k × a l
Associativity:
a k + ( a l + a p ) = ( a k + a l ) + a p ,   a k × ( a l × a p ) = ( a k × a l ) × a p
Commutativity:
a k + a l = a l + a k ,   a k × a l = a l × a k
Existence of identity elements:
a k + 0 = a k ,   a k × 1 = a k
Existence of inverse element (for operator of addition):
a k + ( a k ) = 0
Distributivity:
a k × ( a l + a p ) = ( a k × a l ) + ( a k × a p ) ,   ( a k + a l ) × a p = ( a k × a p ) + ( a l × a p )
No zero divisors (as a logical bi-conditional for operator of multiplication):
a k = 0 ( a k × a l = 0 , a l 0 )
Axiom (16), equivalently, could be also expressed as follows,
( a k = 0 a l = 0 ) a k × a l = 0
In this article, as a new basic algebraic property of the domain of integers, we present the following new axiomatic generalization of the ordinary axiom of “no zero divisors” (16), which particularly has been formulated on the basis of matrix formalism of Clifford algebras (associated with homogeneous polynomials of degree r ≥ 2, over the integral domain ℤ):
“Let A = [ a i j ] be a n × n matrix with entries expressed by the following linear homogeneous polynomials in s variables over the integral domain ℤ: a i j = a i j ( b 1 , b 2 , b 3 , , b s ) = k = 1 s H i j k b k ; suppose also “   r ℕ: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n ”, where F ( b 1 , b 2 , b 3 , , b s ) is a homogeneous polynomial of degree r ≥ 2, and I n is n × n identity matrix; Then the following axiom is assumed (as a new axiomatic generalization of the ordinary axiom of “no zero divisors” of integral domain ℤ):
( A r = 0 ) ( A × M = 0 , M 0 )
where M is a non-zero arbitrary n × 1 column matrix”.
The axiomatic relation (17) is a logical biconditional, where ( A r = 0 ) and ( A × M = 0 , M 0 ) are respectively the antecedent and consequent of this biconditional. In addition, based on the initial assumption   r ℕ: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n , the axiomatic biconditional (17) could be also represented as follows:
[ F ( b 1 , b 2 , b 3 , , b s ) = 0 ] ( A × M = 0 , M 0 )
where the homogeneous equation F ( b 1 , b 2 , b 3 , , b s ) = 0 , and system of linear equations ( A × M = 0 , M 0 ) are respectively the antecedent and consequent of biconditional (17-1). The axiomatic biconditional (17-1), defines a system of linear equations of the type A × M = 0 ( M 0 ) , as the algebraic equivalent representation of r t h degree homogeneous equation F ( b 1 , b 2 , b 3 , , b s ) = 0 (over the integral domain ℤ). The axiom (17) (or (17-1)) for n = 1 , is equivalent to the ordinary axiom of “no zero divisors” (16). In fact, the axiom (16), as a particular case, can be obtained from the axiom (17) (or (17-1)), but not vice versa.
Moreover, according to the Ref. [16], since F ( b 1 , b 2 , b 3 , , b s ) = 0 is a homogeneous equation over ℤ, it is also concluded that homogeneous equations defined over the field of rational numbers ℚ, obey the axiomatic relations (17) and (17-1).
As a crucial additional issue concerning the axiom (17), it should be noted that the condition “   r ℕ: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n ” which is assumed initially in the axiom (17), is also compatible with matrix representation of the generalized Clifford algebras [24,25,26,27,28,29,30,31,32,33] associated with the rth degree homogeneous polynomials F ( b 1 , b 2 , b 3 , , b s ) . In fact, we may represent uniquely the square matrix A (with assumed properties in the axiom (17)) by this homogeneous linear form: A = k = 1 s b k E k , then the relation: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n implies that the square matrices E k (which their entries are independent from the variables b k ) would be generators of the corresponding generalized Clifford algebra associated with the rth degree homogeneous polynomial F ( b 1 , b 2 , b 3 , , b s ) . However, in some particular cases and applications, we may also assume some additional conditions for the generators E k , such as the Hermiticity or anti-Hermiticity (see Section 2.2, Section 2.4 and Section 3.1). In Section 3, we use these algebraic properties of the square matrix A (corresponding with the homogeneous quadratic equations), where we present explicitly the main applications of the axiomatic relations (17) and (17-1) in foundations of physics (where we also use the basic assumptions (2) and (3) mentioned in Section 1.1).
It is noteworthy that since the axiom (17) has been formulated solely in terms of square matrices, in Ref. [34]3 we have shown that all the ordinary algebraic axioms (10)–(15) of the integral domain ℤ (except the axiom of “no zero divisors” (16)) in addition to the new axiom (17) could be also reformulated uniformly in terms of the set of square matrices. Hence, we may conclude that the square matrices, logically, are the most elemental algebraic objects for representing the basic properties of set of integers (as the most fundamental set of mathematics).
In the following, based on the axiomatic relation (17-1), we have constructed a corresponding basic algebraic linearization (including a parameterization procedure) approach applicable to the all classes of homogeneous equation. Hence, it could be also shown that for any given homogeneous equation of degree r ≥ 2 over the ring ℤ (or field ℚ), a square matrix A exists that obeys the relation (17-1). In this regard, for various classes of homogeneous equations, their equivalent systems of linear equations would be derived systematically. As a particular crucial case, in Section 2.2 and Section 2.4, by derivation of the systems of linear equations equivalent to a class of quadratic homogeneous equations (in various number of unknown variables) over the integral domain ℤ (or field ℚ), these equations have been analyzed (and solved) thoroughly by this axiomatic approach. In the following, the basic schemes of this axiomatic linearization-parameterization approach are described.
First, it should be noted that since the entries a i j of square matrix A are linear homogeneous forms expressed in terms of the integral variables b p , i.e., a i j = k = 1 s H i j k b k , we may also represent the square matrix A by this linear matrix form: A = k = 1 s b k E k , then (as noted above) the relation: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n implies that the square matrices E k (the entries of which are independent from the variables b k ) would be generators of the corresponding generalized Clifford algebra associated with the rth degree homogeneous polynomial F ( b 1 , b 2 , b 3 , , b s ) [28,29,30,31,32,33]. However, for some particular cases of the rth degree homogeneous forms F ( b 1 , b 2 , b 3 , , b s ) (for r ≥ 2, such as the standard quadratic forms defined in the quadratic Equation (18) in Section 2.2), without any restriction in the existence and procedure of derivation of their corresponding square matrices A = k = 1 s b k E k (with the algebraic properties assuming in axiom (17)) obeying the Clifford algebraic relation: A r =   F ( b 1 , b 2 , b 3 , , b s ) I n , we may also assume certain additional conditions for the matrix generators E k , such as the Hermiticity (or anti-Hermiticity), and so on (see Section 2.2, Section 2.4 and Section 3.1). In fact, these conditions could be required, for example if a homogeneous invariant relation (of physics) is represented by a homogeneous algebraic equation of the type: F ( b 1 , b 2 , b 3 , , b s ) = 0 , with the algebraic properties as assumed in the axiom (17), where the variables b k denote the components of corresponding physical quantity (such as the relativistic energy-momentum, as it has been assumed in Section 3.1 of this article based on the basic assumption (2) noted in Section 1.1 and Section 3.1).
In Section 2.2, as one of the main applications of the axiomatic relations (17) and (17-1), we derive a unique set the square matrices A n × n (by assuming a minimum value for n , i.e., the size of the corresponding matrix A n × n ) corresponding to the quadratic homogeneous equations of the type: i = 0 s e i f i = 0 , for s = 0, 1, 2, 3, 4, ..., respectively. Subsequently, in Section 2.4, by solving the corresponding systems of linear equations A × M = 0 , we obtain the general parametric solutions of the quadratic homogeneous equations i = 0 s e i f i = 0 , for s = 0, 1, 2, 3, 4,.., respectively. In addition, in Section 2.3 using this systematic axiomatic approach, for some particular forms of homogeneous equations of degrees 3, 4 and 5, their equivalent systems of linear equations have been derived as well. It is noteworthy that using this general axiomatic approach (on the basis of the logical biconditional relations (17) and (17-1)), for any given rth degree homogeneous equation in s unknown variables over the integral domain over ℤ, its equivalent system(s) of linear equations A × M = 0 is derivable (with a unique size, if in the course of the derivation, we also assume a minimum value for n , i.e., the size of corresponding square matrix A n × n ). Furthermore, for a given homogeneous equation of degree r in s unknown variables, the minimum value for n , i.e., the size of the corresponding square matrix A n × n in its equivalent matrix equation: A × M = 0 , is: n min = r s 1 × r s 1 for r = 2 , and n min = r s × r s for r > 2 . For additional detail concerning the general methodology of the derivation of square matrix A n × n and the matrix equation A × M = 0 equivalent to a given homogeneous equation of degree r in s unknown variables—on the basis of the axiomatic relations (17) and (17-1)—see also the preprint versions of this article in Ref. [34].

2.2. The Applications of Axiom (17-1) to Homogeneous Quadratic Equations

In this section, on the basis of axiomatic relation (17-1) and the general methodological notes (mentioned above), for the following general form of homogeneous quadratic equations their equivalent systems of linear equations are derived (uniquely):
Q ( e 0 , f 0 , e 1 , f 1 , , e s , f s ) = i = 0 s e i f i = 0
Equation (18) for s = 0, 1, 2, 3, 4, … is represented by, respectively:
i = 0 0 e i f i = e 0 f 0 = 0 ,
i = 0 1 e i f i = e 0 f 0 + e 1 f 1 = 0 ,
i = 0 2 e i f i = e 0 f 0 + e 1 f 1 + e 2 f 2 = 0 ,
i = 0 3 e i f i = e 0 f 0 + e 1 f 1 + e 2 f 2 + e 3 f 3 = 0 ,
i = 0 4 e i f i = e 0 f 0 + e 1 f 1 + e 2 f 2 + e 3 f 3 + e 4 f 4 = 0 .
It is necessary to note that quadratic Equation (18) is isomorphic to the following ordinary representations of homogeneous quadratic equations:
i , j = 0 s G i j c i c j = 0 ,
i , j = 0 s G i j c i c j = i , j = 0 s G i j d i d j ,
using the linear transformations:
[ e 0 e 1 e 3 . . . e s ] = [ G 00 G 01 G 02 . . . G 0 s G 10 G 11 G 12 . . . G 1 s G 20 G 21 G 22 . . . G 2 s . . . G s 0 G s 1 G s 2 . . . G s s ] [ c 0 + d 0 c 1 + d 1 c 2 + d 2 . . . c s + d s ] , [ f 0 f 1 f 3 . . . f s ] = [ c 0 d 0 c 1 d 1 c 2 d 2 . . . c s d s ]
where [ G i j ] is a symmetric and invertible square matrix, i.e, G i j = G j i and det [ G i j ] 0 , and the quadratic form i , j = 0 s G i j c i c j in Equation (18-1) could be obtained via transformations (18-3), only by taking d i = 0 .
We note here, and show below, that the reason for choosing Equation (18) as the standard general form for representing the homogeneous quadratic equations (that could also be transformed to the ordinary representations of homogeneous quadratic Equations (18-1) and (18-2), by linear transformations (18-3)) is not only its very simple algebraic structure, but also the simple linear homogeneous forms of the entries of square matrices A (expressed in terms of variables e i , f i ) in the corresponding systems of linear equations A × M = 0 obtained as the equivalent form of quadratic Equation (18) in various number of unknown variables.
Moreover, as shown in the following, we may also assume certain Hermiticity and anti-Hermiticity conditions for the deriving square matrices A (in the corresponding systems of linear equations A × M = 0 equivalent to the quadratic Equation (18)), without any restriction in the existence and procedure of derivation of these matrices. By adding these particular conditions, for a specific number of variables in Equation (18), its equivalent matrix equation A × M = 0 could be determined uniquely. In Section 3, where we use the algebraic results obtained in Section 2.2 and Section 2.4 on the basis of axiomatic relations (17) and (17-1) in fact the assumption of these Hermiticity and anti-Hermiticity properties is a necessary issue. These Hermiticity and anti-Hermiticity additional conditions are defined as follows:
“First, by supposing: e 0 = f 0 and e i = f i (for i = 1 , 2 , , s ), the quadratic Equation (18) would be represented as: e 0 2 i = 1 s e i 2 = 0 , and consequently the corresponding square matrix A in the deriving system of linear equations A × M = 0 (which equivalently represent the quadratic Equation (18), based on the axiomatic relation (17-1)) could be also expressed by the homogeneous linear matrix form: A = i = 0 s e i E i , where the real matrices E i are generators of the corresponding Clifford algebra associated with the standard quadratic form e 0 2 i = 1 s e i 2 .
Now for defining the relevant Hermiticity and anti-Hermiticity conditions, we assume that any square matrix A in the deriving matrix equation: A × M = 0 (as the equivalent representation of quadratic Equation (18)), should also has this additional property that by supposing: e 0 = f 0 and e i = f i by which the square matrix A could be represented as: A = i = 0 s e i E i , the matrix generator E 0 be Hermitian: E 0 = E 0 , and matrix generators E i (for i = 1 , 2 , , s ) be anti-Hermitian: E i = E i
As noted previously and shown below, by assuming the above additional Hermiticity and anti-Hermiticity conditions, the system of linear equations A × M = 0 corresponding to quadratic Equation (18), is determined uniquely for any specific number of variables e i , f i . Hence, starting from the simplest (or most trivial) case of quadratic Equation (18)—i.e., Equation (19)—its equivalent system of linear equations is given uniquely as follows:
A × M = [ 0 e 0 f 0 0 ] [ m 1 m 2 ] = 0
where it is assumed M 0 , and in agreement with (17-1) we also have:
A 2 = [ 0 e 0 f 0 0 ] × [ 0 e 0 f 0 0 ] = ( e 0 f 0 ) I 2
For Equation (20), the corresponding equivalent system of linear equations is determined as:
A × M = [ 0 A A 0 ] [ M M ] = [ 0 0 e 0 f 1 0 0 e 1 f 0 f 0 f 1 0 0 e 1 e 0 0 0 ] [ m 1 m 2 m 3 m ] = 0
where we have:
M = [ M M ] 0 , A 2 = [ 0 0 e 0 f 1 0 0 e 1 f 0 f 0 f 1 0 0 e 1 e 0 0 0 ] × [ 0 0 e 0 f 1 0 0 e 1 f 0 f 0 f 1 0 0 e 1 e 0 0 0 ] = ( e 0 f 0 + e 1 f 1 ) I 4
Notice that matrix Equation (25) could be represented by two matrix equations, as follows:
A × M = [ e 0 f 1 e 1 f 0 ] [ m 3 m ] = 0 ,
A × M = [ f 0 f 1 e 1 e 0 ] [ m 1 m 2 ] = 0
The matrix Equations (25-2) and (25-3) are equivalent (due to the assumption of arbitrariness of parameters m 1 , m 2 , m 3 , m ), so we may choose the matrix Equation (25-2) as the system of linear equations equivalent to the quadratic Equation (20)—where for simplicity in the indices of parameters m i , we may simply replace arbitrary parameter m 3 with arbitrary parameter m 1 , as follows (for [ m 1 m ] 0 ):
[ e 0 f 1 e 1 f 0 ] [ m 1 m ] = 0 .
The system of linear equations corresponding to the quadratic Equation (21) is obtained as:
A × M = [ 0 A A 0 ] [ M M ] = [ 0 0 0 0 e 0 0 e 2 f 1 0 0 0 0 0 e 0 e 1 f 2 0 0 0 0 f 2 f 1 f 0 0 0 0 0 0 e 1 e 2 0 f 0 f 0 0 e 2 f 1 0 0 0 0 0 f 0 e 1 f 2 0 0 0 0 f 2 f 1 e 0 0 0 0 0 0 e 1 e 2 0 e 0 0 0 0 0 ] [ m 1 m 2 m 3 m 4 m 5 m 6 m 7 m ] = 0 ,
where in agreement with (17) we have:
A 2 = [ 0 0 0 0 e 0 0 e 2 f 1 0 0 0 0 0 e 0 e 1 f 2 0 0 0 0 f 2 f 1 f 0 0 0 0 0 0 e 1 e 2 0 f 0 f 0 0 e 2 f 1 0 0 0 0 0 f 0 e 1 f 2 0 0 0 0 f 2 f 1 e 0 0 0 0 0 0 e 1 e 2 0 e 0 0 0 0 0 ] × [ 0 0 0 0 e 0 0 e 2 f 1 0 0 0 0 0 e 0 e 1 f 2 0 0 0 0 f 2 f 1 f 0 0 0 0 0 0 e 1 e 2 0 f 0 f 0 0 e 2 f 1 0 0 0 0 0 f 0 e 1 f 2 0 0 0 0 f 2 f 1 e 0 0 0 0 0 0 e 1 e 2 0 e 0 0 0 0 0 ] = ( e 0 f 0 + e 1 f 1 + e 2 f 2 ) I 8
In addition, similar to Equation (25), the obtained matrix Equation (27) is equivalent to a system of two matrix equations, as follows:
A × M = [ e 0 0 e 2 f 1 0 e 0 e 1 f 2 f 2 f 1 f 0 0 e 1 e 2 0 f 0 ] [ m 5 m 6 m 7 m ] = 0 ,
A × M = [ f 0 0 e 2 f 1 0 f 0 e 1 f 2 f 2 f 1 e 0 0 e 1 e 2 0 e 0 ] [ m 1 m 2 m 3 m 4 ] = 0
The matrix Equations (27-2) and (27-3) are equivalent (due to the assumption of arbitrariness of parameters m 1 , m 2 , , m 7 , m ), so we may choose the Equation (27-2) as the system of linear equations corresponding to the quadratic Equation (21)—where for simplicity in the indices of parameters m i , we may simply replace the arbitrary parameters m 5 , m 6 , m 7 with parameters m 1 , m 2 , m 3 , as follows:
[ e 0 0 e 2 f 1 0 e 0 e 1 f 2 f 2 f 1 f 0 0 e 1 e 2 0 f 0 ] [ m 1 m 2 m 3 m ] = 0 ,   [ m 1 m 2 m 3 m ] 0 .
Similarly, for the quadratic Equations (22) the corresponding system of linear equations is obtained uniquely as follows:
[ e 0 0 0 0 0 e 3 e 2 f 1 0 e 0 0 0 e 3 0 e 1 f 2 0 0 e 0 0 e 2 e 1 0 f 3 0 0 0 e 0 f 1 f 2 f 3 0 0 f 3 f 2 e 1 f 0 0 0 0 f 3 0 f 1 e 2 0 f 0 0 0 f 2 f 1 0 e 3 0 0 f 0 0 e 1 e 2 e 3 0 0 0 0 f 0 ] [ m 1 m 2 m 3 m 4 m 5 m 6 m 7 m ] = 0
where the column parametric matrix M in (29) is non-zero M 0 .
In a similar manner, the uniquely obtained system of linear equations corresponding to the quadratic Equation (23), is given by:
[ e 0 0 0 0 0 0 0 0 0 0 0 e 4 0 e 3 e 2 f 1 0 e 0 0 0 0 0 0 0 0 0 e 4 0 e 3 0 e 1 f 2 0 0 e 0 0 0 0 0 0 0 e 4 0 0 e 2 e 1 0 f 3 0 0 0 e 0 0 0 0 0 e 4 0 0 0 f 1 f 2 f 3 0 0 0 0 0 e 0 0 0 0 0 e 3 e 2 e 1 0 0 0 f 4 0 0 0 0 0 e 0 0 0 e 3 0 f 1 f 2 0 0 f 4 0 0 0 0 0 0 0 e 0 0 e 2 f 1 0 f 3 0 f 4 0 0 0 0 0 0 0 0 0 e 0 e 1 f 2 f 3 0 f 4 0 0 0 0 0 0 f 4 0 f 3 f 2 f 1 f 0 0 0 0 0 0 0 0 0 0 f 4 0 f 3 0 e 1 e 2 0 f 0 0 0 0 0 0 0 0 f 4 0 0 f 2 e 1 0 e 3 0 0 f 0 0 0 0 0 0 f 4 0 0 0 f 1 e 2 e 3 0 0 0 0 f 0 0 0 0 0 0 f 3 f 2 e 1 0 0 0 e 4 0 0 0 0 f 0 0 0 0 f 3 0 f 1 e 2 0 0 e 4 0 0 0 0 0 0 f 0 0 0 f 2 f 1 0 e 3 0 e 4 0 0 0 0 0 0 0 0 f 0 0 e 1 e 2 e 3 0 e 4 0 0 0 0 0 0 0 0 0 0 f 0 ] [ m 1 m 2 m 3 m 4 m 5 m 6 m 7 m 8 m 9 m 10 m 11 m 12 m 13 m 14 m 15 m ] = 0
where we have assumed the parametric column matrix M in (30) is non-zero, M 0 .
In a similar manner, the systems of linear equations (written in matrix forms similar to the matrix Equations (24), (26), (28)–(30)) with larger sizes are obtained for the quadratic Equation (18) in more variables (i.e., for s = 6 , 7 , 8 , ), where the size of the square matrices of the corresponding matrix equations is 2 s × 2 s (which could be reduced to 2 s 1 × 2 s 1 for s 2 ). In general (as it has been also mentioned in Section 2.1), the size of the n × n square matrices A (with the minimum value for n ) in the matrix equations A × M = 0 corresponding to the homogeneous polynomials F ( b 1 , b 2 , b 3 , , b s ) of degree r defined in axiom (17) is r s × r s (which for r = 2 this size, in particular, could be reduced to 2 s 1 × 2 s 1 ). Moreover, based on the axiom (17), by solving the obtained system of linear equations corresponding to a homogeneous equation of degree r we may systematically show (and decide) whether this equation has the integral solution.

2.3. The Applications of Axiom (17-1) to Higher Degree Homogeneous Equations

Similar to the uniquely obtained systems of linear equations corresponding to the homogeneous quadratic equations (in Section 2.2), in this section in agreement with the axiom (17), we present the obtained systems of linear equations, i.e., A × M = 0 (by assuming the minimum value for n , i.e., the size of square matrix A n × n ), corresponding to some homogeneous equations of degrees 3, 4 and 5, respectively. For the homogeneous equation of degree three of the type:
F ( e 0 , f 0 , e 1 , f 1 , e 2 , f 2 ) = e 0 2 f 0 e 0 f 0 2 + e 2 2 f 2 e 2 f 2 2 + e 1 f 1 g 1 = 0 ,
the corresponding system of linear equations is given as follows:
A × M = [ 0 0 A 1 A 2 0 0 0 A 3 0 ] [ m 1 m 2 : : : m 27 ] = 0 ,
where A is a 27 × 27 square matrix written in terms of the square 9 × 9 matrices A 1 , A 2 and A 3 , given by:
A 1 = [ e 2 + f 2 0 0 0 0 0 e 0 + f 0 e 1 0 0 e 2 + f 2 0 0 0 0 0 e 0 g 1 0 0 e 2 + f 2 0 0 0 f 1 0 f 0 f 0 e 1 0 e 2 0 0 0 0 0 0 e 0 + f 0 g 1 0 e 2 0 0 0 0 f 1 0 e 0 0 0 e 2 0 0 0 0 0 0 e 0 e 1 0 f 2 0 0 0 0 0 0 f 0 g 1 0 f 2 0 0 0 0 f 1 0 e 0 + f 0 0 0 f 2 ] , A 2 = [ f 2 0 0 0 0 0 e 0 + f 0 e 1 0 0 f 2 0 0 0 0 0 e 0 g 1 0 0 f 2 0 0 0 f 1 0 f 0 f 0 e 1 0 e 2 + f 2 0 0 0 0 0 0 e 0 + f 0 g 1 0 e 2 + f 2 0 0 0 0 f 1 0 e 0 0 0 e 2 + f 2 0 0 0 0 0 0 e 0 e 1 0 e 2 0 0 0 0 0 0 f 0 g 1 0 e 2 0 0 0 0 f 1 0 e 0 + f 0 0 0 e 2 ] , A 3 = [ e 2 0 0 0 0 0 e 0 + f 0 e 1 0 0 e 2 0 0 0 0 0 e 0 g 1 0 0 e 2 0 0 0 f 1 0 f 0 f 0 e 1 0 f 2 0 0 0 0 0 0 e 0 + f 0 g 1 0 f 2 0 0 0 0 f 1 0 e 0 0 0 f 2 0 0 0 0 0 0 e 0 e 1 0 e 2 + f 2 0 0 0 0 0 0 f 0 g 1 0 e 2 + f 2 0 0 0 0 f 1 0 e 0 + f 0 0 0 e 2 + f 2 ]
The uniquely obtained system of linear equations (i.e., A × M = 0 , by assuming the minimum size for the square matrix A n × n ) corresponds to the well-known homogeneous equation of degree three:
F ( a , b , c ) = 2 ( a 3 c 3 + B b 3 ) = 0
and has the following form (in compatible with the new axiom (17)):
A × M = [ 0 0 A 1 A 2 0 0 0 A 3 0 ] [ m 1 m 2 : : : m 27 ] = 0 ,
where A is a 27 × 27 square matrix written in terms of the 9 × 9 matrices A 1 , A 2 and A 3 given by:
A 1 = [ a 0 0 0 0 0 2 c b 0 0 a 0 0 0 0 0 c 2 b 0 0 a 0 0 0 b 0 c c b 0 a 0 0 0 0 0 0 2 c 2 b 0 a 0 0 0 0 b 0 c 0 0 a 0 0 0 0 0 0 c b 0 2 a 0 0 0 0 0 0 c 2 b 0 2 a 0 0 0 0 b 0 2 c 0 0 2 a ] , A 2 = [ 2 a 0 0 0 0 0 2 c b 0 0 2 a 0 0 0 0 0 c 2 b 0 0 2 a 0 0 0 b 0 c c b 0 a 0 0 0 0 0 0 2 c 2 b 0 a 0 0 0 0 b 0 c 0 0 a 0 0 0 0 0 0 c b 0 a 0 0 0 0 0 0 c 2 b 0 a 0 0 0 0 b 0 2 c 0 0 a ] , A 3 = [ a 0 0 0 0 0 2 c b 0 0 a 0 0 0 0 0 c 2 b 0 0 a 0 0 0 b 0 c c b 0 2 a 0 0 0 0 0 0 2 c 2 b 0 2 a 0 0 0 0 b 0 c 0 0 2 a 0 0 0 0 0 0 c b 0 a 0 0 0 0 0 0 c 2 b 0 a 0 0 0 0 b 0 2 c 0 0 a ] .
For the 4th degree homogeneous equation of the type:
F ( e 1 , e 2 , f 1 , f 2 , f 3 , f 4 ) = e 1 e 2 3 + e 1 3 e 2 + f 1 f 2 f 3 f 4 = 0
the corresponding system of linear equations is given as:
A × M = [ 0 0 0 A 1 A 2 0 0 0 0 A 3 0 0 0 0 A 4 0 ] [ m 1 m 2 : : : m 16 ] 0 ,
where A is a 16 × 16 square matrix represented in terms of the 4 × 4 matrices A 1 , A 2 , A 3 , A 4 :
A 1 = [ e 1 + e 2 0 0 f 1 f 2 e 1 0 0 0 f 3 e 1 + e 2 0 0 0 f 4 e 2 ] , A 2 = [ e 2 0 0 f 1 f 2 e 1 + e 2 0 0 0 f 3 e 1 0 0 0 f 4 e 1 + e 2 ] , A 3 = [ e 1 + e 2 0 0 f 1 f 2 e 2 0 0 0 f 3 e 1 + e 2 0 0 0 f 4 e 1 ] , A 4 = [ e 1 0 0 f 1 f 2 e 1 + e 2 0 0 0 f 3 e 2 0 0 0 f 4 e 1 + e 2 ] .
In addition, the system of linear equations corresponding to 5th degree homogeneous equation of the type,
F ( e 1 , e 2 , f 1 , f 2 , f 3 , f 4 , f 5 ) = e 1 4 e 2 e 1 3 e 2 2 e 1 2 e 2 3 + e 1 e 2 4 + f 1 f 2 f 3 f 4 f 5 = 0
is determined as:
A × M = [ 0 0 0 0 A 1 A 2 0 0 0 0 0 A 3 0 0 0 0 0 A 4 0 0 0 0 0 A 5 0 ] [ m 1 m 2 m 3 : : : m 25 ] = 0 ,
where A is a 25 × 25 square matrix expressed in terms of the following 5 × 5 matrices A 1 , A 2 , A 3 , A 4 , A 5 :
A 1 = [ e 1 e 2 0 0 0 f 1 f 2 e 1 0 0 0 0 f 3 e 2 0 0 0 0 f 4 e 1 + e 2 0 0 0 0 f 5 e 1 e 2 ] , A 2 = [ e 1 e 2 0 0 0 f 1 f 2 e 1 e 2 0 0 0 0 f 3 e 1 0 0 0 0 f 4 e 2 0 0 0 0 f 5 e 1 + e 2 ] , A 3 = [ e 1 + e 2 0 0 0 f 1 f 2 e 1 e 2 0 0 0 0 f 3 e 1 e 2 0 0 0 0 f 4 e 1 0 0 0 0 f 5 e 2 ] , A 4 = [ e 2 0 0 0 f 1 f 2 e 1 + e 2 0 0 0 0 f 3 e 1 e 2 0 0 0 0 f 4 e 1 e 2 0 0 0 0 f 5 e 1 ] , A 5 = [ e 1 0 0 0 f 1 f 2 e 2 0 0 0 0 f 3 e 1 + e 2 0 0 0 0 f 4 e 1 e 2 0 0 0 0 f 5 e 1 e 2 ] .

2.4. Obtaining the General Parametric Solutions of Homogeneous Quadratic Equations by Solving their Corresponding Systems of the Linear Equations (Derived in Section 2.2)

In this section by solving the derived systems of the linear Equations (26), (28)–(30) that correspond with the quadratic homogeneous Equations (20)–(23) in Section 2.2, the general parametric solutions of these equations are obtained for unknowns e i and f i . There are the standard methods for obtaining the general solutions of the systems of homogeneous linear equations in integers [35,36]. Using these methods, for the system of linear Equation (26) (and consequently, its corresponding quadratic Equation (20)) we get directly the following general parametric solutions for unknowns e 0 , e 1 and f 0 , f 1 :
e 0 = l k 0 m , f 0 = l k 1 m 1 , e 1 = l k 1 m , f 1 = l k 0 m 1
where k 0 , k 1 , m 1 , m , l are arbitrary parameters. In the matrix representation, the general parametric solution (43) has the following form:
[ e 0 e 1 ] = l M e K = l m [ 1 0 0 1 ] [ k 0 k 1 ] , [ f 0 f 1 ] = l M f K = l [ 0 m 1 m 1 0 ] [ k 0 k 1 ]
where M e = m I 2 , K is a column parametric matrix and M f is also a parametric anti-symmetric matrix.
For the system of linear Equation (28) (and, consequently, for its corresponding quadratic homogeneous Equation (21)), the following general parametric solution is obtained directly:
e 0 = l k 0 m , f 0 = l ( k 1 m 1 k 2 m 2 ) , e 1 = l k 1 m , f 1 = l ( k 2 m 3 k 0 m 1 ) , e 2 = l k 2 m , f 2 = l ( k 0 m 2 k 1 m 3 )
where k 0 , k 1 , k 2 , m 1 , m 2 , m 3 , m , l are arbitrary parameters. In matrix representation the general parametric solution (44) could be also written as follows:
[ e 0 e 1 e 2 ] = l M e K = l m [ 1 0 0 0 1 0 0 0 1 ] [ k 0 k 1 k 2 ] , [ f 0 f 1 f 2 ] = l M f K = l [ 0 m 1 m 2 m 1 0 m 3 m 2 m 3 0 ] [ k 0 k 1 k 2 ]
where M e = m I 3 , K is a column parametric matrix and M f is also a parametric anti-symmetric matrix.
In addition, it could be simply shown that by adding two particular solutions of the types { e i , f i } and { e i , f i } of homogeneous quadratic Equation (18), the new solution { e i + e i , f i } is also obtained, as follows:
( i = 0 s e i f i = 0 , i = 0 s e i f i = 0 ) i = 0 s ( e i f i + e i f i ) = 0 i = 0 s ( e i + e i ) f i = 0
Using the general basic property (44-2) in addition to the general parametric solution (44) of quadratic Equation (21) (which has been obtained directly from the system of linear Equation (28) corresponding to quadratic Equation (21)), exceptionally, the following equivalent general parametric solution is also obtained for quadratic Equation (21):
e 0 = l ( k 0 m k m 3 ) , f 0 = l ( k 1 m 1 k 2 m 2 ) , e 1 = l ( k 1 m k m 2 ) , f 1 = l ( k 2 m 3 k 0 m 1 ) , e 2 = l ( k 2 m k m 1 ) , f 2 = l ( k 0 m 2 k 1 m 3 )
where k 0 , k 1 , k 2 , k , m 1 , m 2 , m 3 , m , l are arbitrary parameters.
Moreover, the parametric solution (45) by the direct bijective replacements of six unknown variables ( e i , f i ) (where i = 0, 1, 2) with the six new variables of the type h μ ν , is given by: e 0 h 23 , e 1 h 20 , e 2 h 21 , f 0 h 10 , f 1 h 31 , f 2 h 03 , in addition to the replacements of nine arbitrary parameters u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , w , with new nine parameters of the types u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , w , given as: k 0 u 3 , k 1 u 0 , k 2 u 1 , k u 2 , k u 2 , m 1 v 1 , m 2 v 0 , m 3 v 3 , m v 2 , l w , exceptionally, could be also represented as follows:
h 23 = w ( u 3 v 2 u 2 v 3 ) , h 10 = w ( u 0 v 1 u 1 v 0 ) , h 20 = w ( u 0 v 2 u 2 v 0 ) , h 31 = w ( u 1 v 3 u 3 v 1 ) , h 21 = w ( u 1 v 2 u 2 v 1 ) , h 03 = w ( u 3 v 0 u 0 v 3 ) ;
where it could be expressed by a single uniform formula as well (for μ, ν = 0,1,2,3):
h μ ν = w ( u ν v μ u μ v ν )
A crucial and important issue concerning the algebraic representation (45-2) (as the differences of products of two parametric variables u μ and v ν ) for the general parametric solution (45), is that it generates a symmetric algebra Sym(V) on the vector space V , where ( u μ , v ν ) V [37]. This essential property of the form (45-2) could be used for various purposes in the following and also in Section 3, in which we show the applications of this axiomatic linearization-parameterization approach and the results obtained in this Section and Section 2.4, to the foundations of physics.
In addition, as it is also shown in the following, it should be mentioned again that the algebraic form (45-2) (representing the symmetric algebra Sym(V)), exceptionally, is determined solely by the parametric solution (44-1) (obtained from the system of Equation (28)) by using the identity (44-2). In fact, from the parametric solutions obtained directly from the subsequent systems of linear equations. i.e., Equations (29), (30) and so on (corresponding to the quadratic Equations (22) and (23), …, and subsequent equations, i.e., i = 0 s e i f i = 0 for s 3 ), the expanded parametric solutions of the type (45) (equivalent to the algebraic form (45-2)) are not derived.
In the following (also see Ref. [34]), we present the parametric solutions that are obtained directly from the systems of linear Equations (29) and (30) and so on, which also would be the parametric solutions of their corresponding quadratic Equation (18) in various number of unknowns (on the basis of axiom (17)). Meanwhile, the following obtained parametric solutions for the systems of linear Equations (29) and (30) and so on, similar to the parametric solutions (43) and (44), include one parametric term for each of unknowns e i , and the sum of s parametric terms for each of unknowns f i (where i = 0 , 1 , 2 , 3 , , s ).
Hence, the following parametric solution is derived directly from the system of linear Equation (29) (that would be also the solution of its corresponding quadratic Equation (22)):
e 0 = l k 0 m , f 0 = l ( k 1 m 1 + k 2 m 2 + k 3 m 3 ) , e 1 = l k 1 m , f 1 = l ( k 0 m 1 + k 3 m 6 k 2 m 7 ) , e 2 = l k 2 m , f 2 = l ( k 0 m 2 k 3 m 5 + k 1 m 7 ) , e 3 = l k 3 m , f 3 = l ( k 0 m 3 + k 2 m 5 k 1 m 6 ) .
where k 0 , k 1 , k 2 , k 3 , l are arbitrary parameters. In the matrix representation, the parametric solution (46) is represented as follows:
[ e 0 e 1 e 2 e 3 ] = l M e K = l m [ 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ] [ k 0 k 1 k 2 k 3 ] , [ f 0 f 1 f 2 f 3 ] = l M f K = l [ o m 1 m 2 m 3 m 1 0 m 7 m 6 m 2 m 7 0 m 5 m 3 m 6 m 5 0 ] [ k 0 k 1 k 2 k 3 ]
where M e = m I 4 , K is a column parametric matrix and M f is also a parametric anti-symmetric matrix.
However, in solutions (46) or (46-1) the parameters m 1 , m 2 , m 3 , m 4 , m 5 , m 6 , m 7 , m are not arbitrary and in fact, in the course of obtaining the solution (46) from the system of linear Equation (29), a condition appears for these parameters as follows:
m 4 m + m 1 m 5 + m 2 m 6 + m 3 m 7 = 0
The condition (47) is also a homogeneous quadratic equation that should be solved first, in order to obtain a general parametric solution for the system of linear Equation (29). Since the parameter m 4 has not appeared in the solution (46), it could be assumed that m 4 = 0 , and the condition (47) is reduced to the following homogeneous quadratic equation, which is equivalent to the quadratic Equation (20) (corresponding to the system of linear Equation (28)):
m 4 = 0 , m 1 m 5 + m 2 m 6 + m 3 m 7 = 0
where the parameter m would be arbitrary. The condition (47-1) is equivalent to the quadratic Equation (21). Hence by using the general parametric solution (45-1) (as the most symmetric solution obtained for quadratic Equation (21) by solving its corresponding system of linear Equation (28)), the following general parametric solution for the condition (47-1) is obtained:
m 1 = w ( u 0 v 1 u 1 v 0 ) , m 2 = w ( u 0 v 2 u 2 v 0 ) , m 3 = w ( u 0 v 3 u 3 v 0 ) , m 5 = w ( u 3 v 2 u 2 v 3 ) , m 6 = w ( u 1 v 3 u 3 v 1 ) , m 7 = w ( u 2 v 1 u 1 v 2 )
where u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , w , m are arbitrary parameters. Now by replacing the solutions (48) (obtained for m 1 , m 2 , m 3 , m 5 , m 6 , m 7 in terms of the new parameters u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , w ) in the relations (46), the general parametric solution of the system of linear Equation (29) (and its corresponding quadratic Equation (22)) is obtained in terms of the arbitrary parameters k 0 , k 1 , k 2 , k 3 , u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , w , m , l .
For the system of linear Equation (30) and its corresponding quadratic Equation (23), the following parametric solution is obtained:
e 0 = l k 0 m , f 0 = l ( k 1 m 1 k 2 m 2 + k 3 m 3 k 4 m 5 ) , e 1 = l k 1 m , f 1 = l ( k 0 m 1 + k 4 m 12 + k 3 m 14 + k 2 m 15 ) , e 2 = l k 2 m , f 2 = l ( k 0 m 2 + k 4 m 11 + k 3 m 13 k 1 m 15 ) , e 3 = l k 3 m , f 3 = l ( k 0 m 3 + k 4 m 10 k 2 m 13 k 1 m 14 ) , e 4 = l k 4 m , f 4 = l ( k 0 m 5 k 3 m 10 k 2 m 11 k 1 m 12 )
where k 0 , k 1 , k 2 , k 3 , k 4 , l are arbitrary parameters. In the matrix representation, the solution (49) can be also written as follows:
[ e 0 e 1 e 2 e 3 e 4 ] = l M e K = l m [ 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 ] [ k 0 k 1 k 2 k 3 k 4 ] , [ f 0 f 1 f 2 f 3 f 4 ] = l M f K = l [ 0 m 1 m 2 m 3 m 5 m 1 0 m 15 m 14 m 12 m 2 m 15 0 m 13 m 11 m 3 m 14 m 13 0 m 10 m 5 m 12 m 11 m 10 0 ] [ k 0 k 1 k 2 k 3 k 4 ]
where M e = m I 5 , K is a column parametric matrix and M f is also a parametric anti-symmetric matrix.
However, similar to the system of Equation (29), in the course of obtaining the solutions (49) or (49-1) from the system of linear Equation (30), the following conditions appear for parameters m 1 , m 2 , m 3 , m 5 , m 10 , m 11 , m 12 , m 13 , m 14 , m 15 :
m 4 m = m 1 m 13 m 2 m 14 m 3 m 15 , m 6 m = m 1 m 11 + m 2 m 12 m 5 m 15 , m 7 m = m 1 m 10 m 3 m 12 m 5 m 14 , m 8 m = m 2 m 10 + m 3 m 11 + m 5 m 13 , m 9 m = m 10 m 15 m 11 m 14 + m 12 m 13 .
that are similar to condition (47). Here also by the same approach, since the parameters m 4 , m 6 , m 7 , m 8 , m 9 have not appeared in the solution (49), it can be assumed that m 4 = m 6 = m 7 = m 8 = m 9 = 0 , and the set of conditions (50) are reduced to the following system of homogeneous quadratic equations which are similar to the quadratic Equation (20) (corresponding to the system of linear Equation (28)):
m 4 = m 6 = m 7 = m 8 = m 9 = 0 , m 1 m 10 m 3 m 12 m 5 m 14 = 0 , m 1 m 11 + m 2 m 12 m 5 m 15 = 0 , m 1 m 13 + m 2 m 14 + m 3 m 15 = 0 , m 2 m 10 + m 5 m 13 + m 3 m 11 = 0 , m 10 m 15 + m 12 m 13 m 11 m 14 = 0 ;
The conditions (50-1) are also similar to the quadratic Equation (21). Hence using again the general parametric solution (45-1), the following general parametric solutions for the system of homogeneous quadratic Equation (50-1) are obtained directly:
m 1 = w ( u 0 v 1 u 1 v 0 ) , m 2 = w ( u 2 v 0 u 0 v 2 ) , m 3 = w ( u 0 v 3 u 3 v 0 ) , m 4 = 0 , m 5 = w ( u 4 v 0 u 0 v 4 ) , m 6 = 0 , m 7 = 0 , m 8 = 0 , m 9 = 0 , m 10 = w ( u 3 v 4 u 4 v 3 ) , m 11 = w ( u 2 v 4 u 4 v 2 ) , m 12 = w ( u 1 v 4 u 4 v 1 ) , m 13 = w ( u 2 v 3 u 3 v 2 ) , m 14 = w ( u 1 v 3 u 3 v 1 ) , m 15 = w ( u 1 v 2 u 2 v 1 ) .
where u 0 , u 1 , u 2 , u 3 , u 4 , v 0 , v 1 , v 2 , v 3 , v 4 and w are arbitrary parameters. Now by replacing the solution (51) (that has been obtained for m 1 , m 2 , m 3 , m 5 , m 10 , m 11 , m 12 , m 13 , m 14 , m 15 in terms of the new parameters u 0 , u 1 , u 2 , u 3 , u 4 , v 0 , v 1 , v 2 , v 3 , v 4 , w ) in the relations (49), the general parametric solution of the system of linear Equation (30) (and its corresponding quadratic Equation (23)) is obtained in terms of the arbitrary parameters k 0 , k 1 , k 2 , k 3 , k 4 , u 0 , u 1 , u 2 , u 3 , u 4 , v 0 , v 1 , v 2 , v 3 , v 4 , w , m , l .
Meanwhile, similar to relations (48) and (51), it should be noted that arbitrary parameter m 1 in the general parametric solution (43) and the arbitrary parameters m 1 , m 2 , m 3 in the general parametric solution (44) (which have been obtained as the solutions of quadratic Equations (20) and (21), respectively, by solving their equivalent systems of linear Equations (26) and (28)), by keeping their arbitrariness property, could be expressed in terms of new arbitrary parameters u 0 , u 1 , v 0 , v 1 and u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , as follows, respectively:
m 1 = w ( u 0 v 1 u 1 v 0 ) ;
m 1 = w ( u 0 v 1 u 1 v 0 ) , m 2 = w ( u 0 v 2 u 2 v 0 ) , m 3 = w ( u 0 v 3 u 3 v 0 ) .
In fact, as a particular common algebraic property of both parametric relations (43-2) and (44-2), it could be shown directly that by choosing appropriate integer values for parameters u 0 , u 1 , v 0 , v 1 , w in the relation (43-2), the parameter m 1 (defined in terms of arbitrary parameters u 0 , u 1 , v 0 , v 1 , w ) could take any given integer value, and similarly, by choosing appropriate integer values for parameters u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w in the relation (43-2), the parameters m 1 , m 2 , m 3 (defined in terms of arbitrary parameters u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w ) could also take any given integer values.
Therefore, using this common algebraic property of the parametric relations (43-2) and (44-2), the arbitrary parameter m 1 in general parametric solutions (43), and arbitrary parameters m 1 , m 2 , m 3 in general parametric solutions (44), could be equivalently replaced by new arbitrary parameters u 0 , u 1 , v 0 , v 1 , w and u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w , respectively. In addition, for the general quadratic homogeneous Equation (18) with a larger number of unknowns, the general parametric solutions could be obtained by the same approaches used above for quadratic Equations (20)–(23), i.e., by solving their corresponding systems of linear equations (defined on the basis of axiom (17)). Moreover, using the isomorphic transformations (18-3) and the above general parametric solutions obtained for quadratic Equations (20)–(23), …, (via solving their corresponding systems of linear Equations (26), (28)–(30), …), the general parametric solutions of quadratic equations of the regular type (18-2) (in various number of unknown) are also obtained straightforwardly. All the parametric solutions that are obtained by this new systematic matrix approach for the homogeneous quadratic equations and also higher degree homogeneous equations of the type F ( x 1 , x 2 , x 3 , , x s ) = 0 (defined in the axiom (17)), are fully compatible with the solutions and conclusions that have been obtained previously for various homogeneous equations by different and miscellaneous methods and approaches [16,35,36]. In Section 3, we have used the uniquely specified systems of homogeneous linear equations (and also their general parametric solutions) corresponding with the homogeneous quadratic equations. It has been assumed that the components of the relativistic energy-momentum vector (as one of the most basic physical quantities) in the Lorentz invariant energy-momentum (homogeneous) quadratic relation, can only take the rational values.

3. A Unique Mathematical Derivation of the Laws Governing the Fundamental Forces of Nature: Based on a New Algebraic-Axiomatic (Matrix) Approach

In this Section, as a new mathematical approach to the origin of the laws of nature, using the new basic algebraic axiomatic (matrix) formalism as presented in Section 2 “it is shown that certain mathematical forms of fundamental laws of nature, including laws governing the fundamental forces of nature (represented by a set of two definite classes of general covariant massive field equations, with new matrix formalisms), are derived uniquely from only a very few axioms” where in agreement with the rational Lorentz group, it is also basically assumed that the components of relativistic energy-momentum can only take rational values. Concerning the basic assumption of the rationality of relativistic energy-momentum, it is necessary to note that the rational Lorentz symmetry group is not only dense in the general form of Lorentz group, but also is compatible with the necessary conditions required basically for the formalism of a consistent relativistic quantum theory [15]. In essence, the main scheme of this new mathematical axiomatic approach to fundamental laws of nature is as follows. First in Section 3.1, based on the assumption of rationality of D-momentum, by linearization (along with a parameterization procedure) of the Lorentz invariant energy-momentum quadratic relation, a unique set of Lorentz invariant systems of homogeneous linear equations (with matrix formalisms compatible with certain Clifford, and symmetric algebras) is derived.
Then in Section 3.4, by initial quantization (followed by a basic procedure of minimal coupling to space-time geometry) of these determined systems of linear equations, a set of two classes of general covariant massive (tensor) field equations (with matrix formalisms compatible with certain Clifford, and Weyl algebras) is derived. Each class of the derived general covariant field equations also includes a definite form of torsion field appeared as generator of the corresponding field’ invariant mass. In addition, in Section 3.4 and Section 3.11, it is shown that the (1 + 3)-dimensional cases of two classes of derived field equations represent a new general covariant massive formalism of bispinor fields of spin-2, and spin-1 particles, respectively. In fact, these uniquely determined bispinor fields represent a unique set of new generalized massive forms of the laws governing the fundamental forces of nature, including the Einstein (gravitational), Maxwell (electromagnetic) and Yang-Mills (nuclear) field equations. Moreover, it is also shown that the (1 + 2)-dimensional cases of two classes of these field equations represent (asymptotically) a new general covariant massive formalism of bispinor fields of spin-3/2 and spin-1/2 particles, respectively, corresponding with the Dirac and Rarita–Schwinger equations.
As a particular consequence, in Section 3.6, it is shown that a certain massive formalism of general relativity—with a definite form of torsion field appearing originally as the generator of the gravitational field’s invariant mass—is obtained only by the initial quantization (followed by a basic procedure of minimal coupling to space-time geometry) of a certain set of special relativistic algebraic matrix equations. In Section 3.9, it has been proven that Lagrangian densities specified for the originally derived new massive forms of the Maxwell, Yang-Mills and Dirac field equations, are also gauge-invariant, where the invariant mass of each field is generated solely by the corresponding torsion field. In addition, in Section 3.10, in agreement with recent astronomical data, a particular new form of massive boson is identified (corresponding to U(1) gauge group) with invariant mass: mγ ≈ 4.90571 × 10−50 kg, generated by a coupled torsion field of the background space-time geometry.
Moreover, in Section 3.12, based on the definite mathematical formalism of this axiomatic approach, along with the C, P and T symmetries (represented basically by the corresponding quantum operators) of the fundamentally derived field equations, it has been concluded that the universe could be realized solely with the (1 + 2) and (1 + 3)-dimensional space-times (where this conclusion, in particular, is based on the time-reversal symmetry). In Section 3.13 and Section 3.14, it is proved that ‘CPT’ is the only (unique) combination of C, P, and T symmetries that could be defined as a symmetry for interacting fields. In addition, in Section 3.14, on the basis of these discrete symmetries of derived field equations, it has been also shown that only left-handed particle fields (along with their complementary right-handed fields) could be coupled to the corresponding (any) source currents. Furthermore, in Section 3.15, it has been shown that metric of the background space-time is diagonalized for the uniquely derived fermion field equations (defined and expressed solely in (1 + 2)-dimensional space-time), where this property generates a certain set of additional symmetries corresponding uniquely to the SU(2)L⊗U(2)R symmetry group for spin-1/2 fermion fields (representing “1 + 3” generations of four fermions, including a group of eight leptons and a group of eight quarks), and also the SU(2)L⊗U(2)R and SU(3) gauge symmetry groups for spin-1 boson fields coupled to the spin-1/2 fermionic source currents. Hence, along with the known elementary particles, eight new elementary particles, including: four new charge-less right-handed spin-1/2 fermions (two leptons and two quarks, represented by “ze, zn and zu, zd”), a spin-3/2 fermion, and also three new spin-1 massive bosons (represented by W ˜ + , W ˜ , Z , where in particular, the new boson Z is complementary right-handed particle of ordinary Z boson), are uniquely predicted by this new mathematical axiomatic approach.
Furthermore, as a particular result, it is generally concluded in Section 3.6—based on the definite and unique formulation of the derived Maxwell’s equations (and also determined Yang-Mills equations, represented uniquely with two specific forms of gauge symmetries, in Section 3.15, Formula (114-4)–(114-9)— that magnetic monopoles cannot exist in nature.

3.1. The Main Results Obtained in This Article are Based on the Following Three Postulates

As noted in Section 1.1, the main results obtained in this article are based on the following three basic assumptions (as postulates):
(1)- “A new definite axiomatic generalization of the axiom of “no zero divisors” of integral domains (including the integr ring ℤ) is assumed (represented by Formula (17), in Section 2.1);”
This basic assumption (as a postulate) is a new mathematical concept. In Section 2.1, based on this new axiom, a general algebraic axiomatic (matrix) approach (in the form of a basic linearization-parameterization theory) to homogeneous equations of degree r ≥ 2 (over the integer domain, extendable to field of rational numbers), has been formulated. A summary of the main results obtained from this axiomatic approach have been presented in Section 1.1. As particular outcome of this new mathematical axiomatic formalism (based on the axiomatic relations (17) and (17-1), including their basic algebraic properties presented in detail, in Section 2.1, Section 2.2, Section 2.3 and Section 2.4), in Section 3.4, it is shown that using, a unique set of general covariant massive (tensor) field equations (with new matrix formalism compatible with Clifford, and Weyl algebras), corresponding with the fundamental field equations of physics, are derived—where, in agreement with the rational Lorentz symmetry group, it has been assumed that the components of relativistic energy-momentum can only take the rational values. In Section 3.2, Section 3.3, Section 3.4, Section 3.5, Section 3.6, Section 3.7, Section 3.8, Section 3.9, Section 3.10, Section 3.11, Section 3.12, Section 3.13, Section 3.14 and Section 3.15, we present in detail the main applications of this basic algebraic assumption (along with the following basic assumptions (2) and (3)) to fundamental physics.
(2)- “In agreement with the rational Lorentz symmetry group, we assume basically that the components of relativistic energy-momentum (D-momentum) can only take the rational values;”
Concerning this assumption, it is necessary to note that the rational Lorentz symmetry group is not only dense in the general form of Lorentz group, but also is compatible with the necessary conditions required basically for the formalism of a consistent relativistic quantum theory [15]. Moreover, this assumption is clearly also compatible with any quantum circumstance in which the energy-momentum of a relativistic particle is transferred as integer multiples of the quantum of action “h” (Planck constant).
Before defining the next basic assumption, it should be noted that from basic assumptions (1) and (2), it follows directly that the Lorentz invariant energy-momentum quadratic relation (represented by Formula (52)) is a particular form of homogeneous quadratic equation (18-2). Hence, using the set of systems of linear equations that have been determined uniquely as equivalent algebraic representations of the corresponding set of quadratic homogeneous equations (given by equation (18-2) in various number of unknown variables, respectively), a unique set of the Lorentz invariant systems of homogeneous linear equations (with matrix formalisms compatible with certain Clifford, and symmetric algebras) are also determined, representing equivalent algebraic forms of the energy-momentum quadratic relation in various space-time dimensions, respectively. Subsequently, we have shown that by first quantization (followed by a basic procedure of minimal coupling to space-time geometry) of these determined systems of linear equations, a unique set of two definite classes of general covariant massive (tensor) field equations (with matrix formalisms compatible with certain Clifford, and Weyl algebras) is also derived, corresponding to various space-time dimensions, respectively. In addition, it is also shown that this derived set of two classes of general covariant field equations represent new tensor massive (matrix) formalism of the fundamental field equations of physics, corresponding to fundamental laws of nature (including the laws governing the fundamental forces of nature). Following these essential results, in addition to the basic assumptions (1) and (2), it would be also basically assumed that:
(3)- “We assume that the mathematical formalism of the fundamental laws of nature, are defined solely by the axiomatic matrix constitution formulated uniquely on the basis of postulates (1) and (2)”.
In addition to this basic assumption, in Section 3.11, the C, P and T symmetries of uniquely derived general covariant field equations (that are Equations (71) and (72), in Section 3.4), are also represented basically by their corresponding quantum matrix operators.
As the next step, in the following, based on the basic assumption (2), i.e., the assumption of rationality of the relativistic energy-momentum, the following Lorentz invariant quadratic relations (expressed in terms of the components of D-momentums p μ , p μ of a relativistic massive particle (given bytwo reference frames), and also components of quantity p μ s t = m 0 k μ , where m0 is the invariant mass of particle and k μ is its covariant velocity in the stationary reference frame):
g μ ν p μ p ν = g μ ν p μ p ν ,
g μ ν p μ p ν = g μ ν p μ s t p ν s t = g μ ν ( m 0 k μ ) ( m 0 k ν ) = g 00 ( m 0 k 0 ) 2 = ( m 0 c ) 2 .
would be particular cases of homogeneous quadratic Equation (18-2) in Section 2.2, and hence, they would be necessarily subject to the process of linearization (along with a parameterization procedure) using the systematic axiomatic approach presented Section 2, Section 2.2 and Section 2.4 (formulated based on the basic assumption (1)).
The Lorentz invariant relations (51) and (52) (as the norm of the relativistic energy-momentum) have been defined in the D-dimensional space-time, where m0 is the invariant mass of the particle, p μ and p μ are its relativistic energy-momentums (i.e., D-momentums) given respectively in two reference frames, k μ is a time-like covariant vector given by: k μ = ( k 0 , 0 , , 0 ) = ( c / g 00 , 0 , , 0 ) , “c” is the speed of light, and the components of metric have constant values. As noted in Section 1.2, in this article the sign conventions (2) (including the metric signature (+ – – … – )) and geometric units would be used (where in particular “c = 1”). However, for the clarity, in some of relativistic Formula (such as the relativistic matrix relations), the speed of light “c” is indicated formally.
As a crucial issue here, it should be noted that in the invariant quadratic relations (51) and (52), the components of metric which have the constant values (as assumed), have necessarily been written by their general representations g μ ν (and not by the Minkowski metric η μ ν , and so on). This follows from the fact that by axiomatic approach of linearization-parameterization (presented in Section 2.1, Section 2.2, Section 2.3 and Section 2.4) of quadratic relations (51) and (52) (as particular forms of homogeneous quadratic Equation (18-2) which could be expressed equivalently by quadratic equations of the types (18) via the linear transformations (18-3)), their corresponding algebraic equivalent systems of linear equations could be determined uniquely. In fact, based on the formulations of systems of linear equations obtained uniquely for the quadratic Equation (18) in Section 2.2, Section 2.3 and Section 2.4, it is concluded directly that the algebraic equivalent systems of linear equations corresponding to the relations (51) and (52), are determined uniquely if and only if these quadratic relations be expressed in terms of the components g μ ν represented by their general forms (and not in terms of any special presentation of the metric’s components, such as the Minkowski metric, and so on). However, after the derivation of corresponding systems of linear equations (representing uniquely the equivalent algebraic matrix forms of the quadratic relations (51) and (52) in various space-time dimensions), the Minkowski metric could be used in these equations (and the subsequent relativistic equations and relations as well).
Hence, using the systems of linear Equations (24), (26), (28)–(30), …, obtained uniquely on the basis of the axiom (17) by linearization (along with a parameterization procedure) of the homogeneous quadratic Equations (19)–(23), … (which could be transformed directly to the general quadratic Equation (18-2), by the isomorphic linear transformations (18-3)) and also using the parametric relations (43-2), (44-3), (48) and (52) (expressed in terms of the arbitrary parameters u μ and v μ ), as the result of linearization (along with a parameterization procedure) of the invariant quadratic relations (51) and (52), the following systems of linear equations are also derived uniquely corresponding with various space-time dimensions, respectively:
For (1 + 0)-dimensional case of the invariant relation (51), we obtain:
[ g 0 ν ( p ν + p ν ) ] [ s ] = 0
where ν = 0 and parameter s is arbitrary;
For (1 + 1)-dimensional space-time we have:
[ g 0 ν ( p ν + p ν ) p 1 p 1 g 1 ν ( p ν + p ν ) ( p 0 p 0 ) ] [ ( u 0 v 1 u 1 v 0 ) w s ] = 0
where ν = 0,1 and u 0 , u 1 , v 0 , v 1 , w , s are arbitrary parameters;
For (1 + 2)-dimensional space-time we have:
[ g 0 ν ( p ν + p ν ) 0 g 2 ν ( p ν + p ν ) p 1 p 1 0 g 0 ν ( p ν + p ν ) g 1 ν ( p ν + p ν ) ( p 2 p 2 ) ( p 2 p 2 ) ( p 1 p 1 ) ( p 0 p 0 ) 0 g 1 ν ( p ν + p ν ) g 2 ν ( p ν + p ν ) 0 ( p 0 p 0 ) ] [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 1 v 2 u 2 v 1 ) w s ] = 0
where ν = 0, 1, 2 and u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w , s are arbitrary parameters;
For (1 + 3) -dimensional space-time we obtain:
[ e 0 0 0 0 0 e 3 e 2 f 1 0 e 0 0 0 e 3 0 e 1 f 2 0 0 e 0 0 e 2 e 1 0 f 3 0 0 0 e 0 f 1 f 2 f 3 0 0 f 3 f 2 e 1 f 0 0 0 0 f 3 0 f 1 e 2 0 f 0 0 0 f 2 f 1 0 e 3 0 0 f 0 0 e 1 e 2 e 3 0 0 0 0 f 0 ] [ ( u 0 v 1 u 1 v 0 ) w ( u 0 v 2 u 2 v 0 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 3 v 2 u 2 v 3 ) w ( u 1 v 3 u 3 v 1 ) w ( u 2 v 1 u 1 v 2 ) w s ] = 0
where ν = 0, 1, 2, 3 and w , u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , s are arbitrary parameters, and we also have:
e 0 = g 0 ν ( p ν + p ν ) , f 0 = ( p 0 p 0 ) 0 , e 1 = g 1 ν ( p ν + p ν ) , f 1 = ( p 1 p 1 ) , e 2 = g 2 ν ( p ν + p ν ) , f 2 = ( p 2 p 2 ) , e 3 = g 3 ν ( p ν + p ν ) , f 3 = ( p 3 p 3 ) ;
For (1 + 4)-dimensional case, the system of linear equations corresponding to the invariant quadratic relation (51) is specified as follows:
[ e 0 0 0 0 0 0 0 0 0 0 0 e 4 0 e 3 e 2 f 1 0 e 0 0 0 0 0 0 0 0 0 e 4 0 e 3 0 e 1 f 2 0 0 e 0 0 0 0 0 0 0 e 4 0 0 e 2 e 1 0 f 3 0 0 0 e 0 0 0 0 0 e 4 0 0 0 f 1 f 2 f 3 0 0 0 0 0 e 0 0 0 0 0 e 3 e 2 e 1 0 0 0 f 4 0 0 0 0 0 e 0 0 0 e 3 0 f 1 f 2 0 0 f 4 0 0 0 0 0 0 0 e 0 0 e 2 f 1 0 f 3 0 f 4 0 0 0 0 0 0 0 0 0 e 0 e 1 f 2 f 3 0 f 4 0 0 0 0 0 0 f 4 0 f 3 f 2 f 1 f 0 0 0 0 0 0 0 0 0 0 f 4 0 f 3 0 e 1 e 2 0 f 0 0 0 0 0 0 0 0 f 4 0 0 f 2 e 1 0 e 3 0 0 f 0 0 0 0 0 0 f 4 0 0 0 f 1 e 2 e 3 0 0 0 0 f 0 0 0 0 0 0 f 3 f 2 e 1 0 0 0 e 4 0 0 0 0 f 0 0 0 0 f 3 0 f 1 e 2 0 0 e 4 0 0 0 0 0 0 f 0 0 0 f 2 f 1 0 e 3 0 e 4 0 0 0 0 0 0 0 0 f 0 0 e 1 e 2 e 3 0 e 4 0 0 0 0 0 0 0 0 0 0 f 0 ] [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 4 v 0 u 0 v 4 ) w 0 0 0 0 ( u 3 v 4 u 4 v 3 ) w ( u 2 v 4 u 4 v 2 ) w ( u 1 v 4 u 4 v 1 ) w ( u 2 v 3 u 3 v 2 ) w ( u 1 v 3 u 3 v 1 ) w ( u 1 v 2 u 2 v 1 ) w s ] = 0
where ν = 0, 1, 2, 3, 4, u 0 , u 1 , u 2 , u 3 , u 4 , v 0 , v 1 , v 2 , v 3 , v 4 , w , s are arbitrary parameters, and we have:
e 0 = g 0 ν ( p ν + p ν ) , f 0 = p 0 p 0 , e 1 = g 1 ν ( p ν + p ν ) , f 1 = p 1 p 1 , e 2 = g 2 ν ( p ν + p ν ) , f 2 = p 2 p 2 , e 3 = g 3 ν ( p ν + p ν ) , f 3 = p 3 p 3 , e 4 = g 4 ν ( p ν + p ν ) , f 4 = p 4 p 4 ;
The systems of linear equations that are obtained for (1 + 5) and higher dimensional cases of the invariant quadratic relation (51), have also the formulations similar to the obtained systems of linear Equations (53)–(57), and would be expressed by the matrix product of a 2 N × 2 N square matrix and a 2 N × 1 column matrix in (1 + N)-dimensional space-time. For (1 + 5)-dimensional case of the invariant relation (51), the column matrix of the corresponding system of linear equations (expressed by the matrix product of a 32 × 32 square matrix and a 32 × 1 column matrix) are given by (where u 0 , u 1 , u 2 , u 3 , u 4 , u 5 , v 0 , v 1 , v 2 , v 3 , v 4 , v 5 , w , s are arbitrary parameters):
S = [ S S ] ,   S = [ ( u 0 v 1 u 1 v 0 ) w ( u 0 v 2 u 2 v 0 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 0 v 4 u 4 v 0 ) w 0 0 0 ( u 0 v 5 u 5 v 0 ) w 0 0 0 0 0 0 0 ] ,   S = [ 0 0 0 ( u 5 v 4 u 4 v 5 ) w 0 ( u 3 v 5 u 5 v 3 ) w ( u 5 v 2 u 2 v 5 ) w ( u 1 v 5 u 5 v 1 ) w 0 ( u 4 v 3 u 3 v 4 ) w ( u 2 v 4 u 4 v 2 ) w ( u 4 v 1 u 1 v 4 ) w ( u 2 v 3 u 3 v 2 ) w ( u 3 v 1 u 1 v 3 ) w ( u 1 v 2 u 2 v 1 ) w s ] .
In a similar manner, using the axiomatic approach presented in Section 2, the systems of linear equations corresponding to the energy-momentum invariant relation (52) in various space-time dimensions are obtained uniquely as follows, respectively (note that by using the geometric units, we would take c = 1 ):
For (1 + 0)-dimensional space-time we obtain:
[ g 0 ν p ν g 00 m 0 c g 00 ] [ s ] = 0
where ν = 0 and parameter s is arbitrary;
For (1 + 1)-dimensional space-time we have:
[ g 0 ν p ν g 00 ( m 0 c g 00 ) p 1 g 1 ν p ν + g 10 ( m 0 c g 00 ) ( p 0 + ( m 0 c g 00 ) ) ] [ ( u 0 v 1 u 1 v 0 ) w s ] = 0
where ν = 0, 1 and u 0 , u 1 , v 0 , v 1 , w , s are arbitrary parameters;
For (1 + 2) dimensions we have (where ν = 0, 1, 2 and u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w , s are arbitrary parameters):
[ g 0 ν p ν g 00 ( m 0 c g 00 ) 0 g 2 ν p ν + g 20 ( m 0 c g 00 ) p 1 0 g 0 ν p ν g 00 ( m 0 c g 00 ) g 1 ν p ν + g 10 ( m 0 c g 00 ) p 2 p 2 p 1 ( p 0 + ( m 0 c g 00 ) ) 0 g 1 ν p ν + g 10 ( m 0 c g 00 ) g 2 ν p ν g 20 ( m 0 c g 00 ) 0 ( p 0 + ( m 0 c g 00 ) ) ] [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 1 v 2 u 2 v 1 ) w s ] = 0
For (1 + 3)-dimensional space-time we obtain:
[ e 0 0 0 0 0 e 3 e 2 f 1 0 e 0 0 0 e 3 0 e 1 f 2 0 0 e 0 0 e 2 e 1 0 f 3 0 0 0 e 0 f 1 f 2 f 3 0 0 f 3 f 2 e 1 f 0 0 0 0 f 3 0 f 1 e 2 0 f 0 0 0 f 2 f 1 0 e 3 0 0 f 0 0 e 1 e 2 e 3 0 0 0 0 f 0 ] [ ( u 0 v 1 u 1 v 0 ) w ( u 0 v 2 u 2 v 0 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 3 v 2 u 2 v 3 ) w ( u 1 v 3 u 3 v 1 ) w ( u 2 v 1 u 1 v 2 ) w s ] = 0
where ν = 0, 1, 2, 3 and u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , w , s are arbitrary parameters, and we also have:
e 0 = g 0 ν p ν g 00 ( m 0 c / g 00 ) , f 0 = p 0 + ( m 0 c / g 00 ) , e 1 = g 1 ν p ν g 10 ( m 0 c / g 00 ) , f 1 = p 1 , e 2 = g 2 ν p ν g 20 ( m 0 c / g 00 ) , f 2 = p 2 , e 3 = g 3 ν p ν g 30 ( m 0 c / g 00 ) , f 3 = p 3 ;
For (1 + 4)-dimensional space-time, the system of linear equations corresponding to the invariant quadratic relation (52) is derived as follows:
[ e 0 0 0 0 0 0 0 0 0 0 0 e 4 0 e 3 e 2 f 1 0 e 0 0 0 0 0 0 0 0 0 e 4 0 e 3 0 e 1 f 2 0 0 e 0 0 0 0 0 0 0 e 4 0 0 e 2 e 1 0 f 3 0 0 0 e 0 0 0 0 0 e 4 0 0 0 f 1 f 2 f 3 0 0 0 0 0 e 0 0 0 0 0 e 3 e 2 e 1 0 0 0 f 4 0 0 0 0 0 e 0 0 0 e 3 0 f 1 f 2 0 0 f 4 0 0 0 0 0 0 0 e 0 0 e 2 f 1 0 f 3 0 f 4 0 0 0 0 0 0 0 0 0 e 0 e 1 f 2 f 3 0 f 4 0 0 0 0 0 0 f 4 0 f 3 f 2 f 1 f 0 0 0 0 0 0 0 0 0 0 f 4 0 f 3 0 e 1 e 2 0 f 0 0 0 0 0 0 0 0 f 4 0 0 f 2 e 1 0 e 3 0 0 f 0 0 0 0 0 0 f 4 0 0 0 f 1 e 2 e 3 0 0 0 0 f 0 0 0 0 0 0 f 3 f 2 e 1 0 0 0 e 4 0 0 0 0 f 0 0 0 0 f 3 0 f 1 e 2 0 0 e 4 0 0 0 0 0 0 f 0 0 0 f 2 f 1 0 e 3 0 e 4 0 0 0 0 0 0 0 0 f 0 0 e 1 e 2 e 3 0 e 4 0 0 0 0 0 0 0 0 0 0 f 0 ] [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 4 v 0 u 0 v 4 ) w 0 0 0 0 ( u 3 v 4 u 4 v 3 ) w ( u 2 v 4 u 4 v 2 ) w ( u 1 v 4 u 4 v 1 ) w ( u 2 v 3 u 3 v 2 ) w ( u 1 v 3 u 3 v 1 ) w ( u 1 v 2 u 2 v 1 ) w s ] = 0
where ν = 0, 1, 2, 3, 4 and u 0 , u 1 , u 2 , u 3 , u 4 , v 0 , v 1 , v 2 , v 3 , v 4 , w , s are arbitrary parameters, and we have:
e 0 = g 0 ν p ν g 00 ( m 0 c / g 00 ) , f 0 = p 0 + ( m 0 c / g 00 ) , e 1 = g 1 ν p ν g 10 ( m 0 c / g 00 ) , f 1 = p 1 , e 2 = g 2 ν p ν g 20 ( m 0 c / g 00 ) , f 2 = p 2 , e 3 = g 3 ν p ν g 30 ( m 0 c / g 00 ) , f 3 = p 3 , e 4 = g 4 ν p ν g 40 ( m 0 c / g 00 ) , f 4 = p 4 .
The systems of linear equations that are obtained for (1 + 5) and higher dimensional cases of the energy-momentum quadratic relation (52), have also the formulations similar to the obtained systems of linear Equations (58)–(62), and would be expressed by the matrix product of a 2 N × 2 N square matrix and a 2 N × 1 column matrix in (1 + N)-dimensional space-time. For the (1 + 5)-dimensional case of energy-momentum relation (52), the column matrix of the corresponding system of linear equations (expressed by the matrix product of a 32 × 32 square matrix and a 32 × 1 column matrix, similar to (57-2)) is given by:
S = [ S S ] ,   S = [ ( u 0 v 1 u 1 v 0 ) w ( u 0 v 2 u 2 v 0 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 0 v 4 u 4 v 0 ) w 0 0 0 ( u 0 v 5 u 5 v 0 ) w 0 0 0 0 0 0 0 ] ,   S = [ 0 0 0 ( u 5 v 4 u 4 v 5 ) w 0 ( u 3 v 5 u 5 v 3 ) w ( u 5 v 2 u 2 v 5 ) w ( u 1 v 5 u 5 v 1 ) w 0 ( u 4 v 3 u 3 v 4 ) w ( u 2 v 4 u 4 v 2 ) w ( u 4 v 1 u 1 v 4 ) w ( u 2 v 3 u 3 v 2 ) w ( u 3 v 1 u 1 v 3 ) w ( u 1 v 2 u 2 v 1 ) w s ] .
where u 0 , u 1 , u 2 , u 3 , u 4 , u 5 , v 0 , v 1 , v 2 , v 3 , v 4 , v 5 , w , s are arbitrary parameters.

3.2. Deivation of the Rational Lorentz Transformations

From the derived systems of linear Equations (54)–(57) corresponding to the (1 + 1)–(1 + 4)-dimensional cases of the invariant relation (51), and also using the general parametric solutions (43)–(51) (obtained for systems of linear Equations (26)–(30)), the rational Lorentz transformations (which are completely dense in the standard group of Lorentz transformations [15], as noted in Section 3.1) are derived for momentums p μ and p μ . For instance, assuming the Minkowski metric from the system of linear Equation (55), a parametric form of rational Lorentz transformations for three-momentums p μ and p μ in (1 + 2) dimensional space-time is derived as follows:
[ 1 + z 0 2 + z 1 2 + z 2 2 1 z 0 2 z 1 2 + z 2 2 2 ( z 0 + z 1 z 2 ) 1 z 0 2 z 1 2 + z 2 2 2 ( z 1 z 0 z 2 ) 1 z 0 2 z 1 2 + z 2 2 2 ( z 0 z 1 z 2 ) 1 z 0 2 z 1 2 + z 2 2 1 + z 0 2 z 1 2 z 2 2 1 z 0 2 z 1 2 + z 2 2 2 ( z 2 z 0 z 1 ) 1 z 0 2 z 1 2 + z 2 2 2 ( z 1 + z 0 z 2 ) 1 z 0 2 z 1 2 + z 2 2 2 ( z 2 + z 0 z 1 ) 1 z 0 2 z 1 2 + z 2 2 1 z 0 2 + z 1 2 z 2 2 1 z 0 2 z 1 2 + z 2 2 ] [ p 0 p 1 p 2 ] = [ p 0 p 1 p 2 ]
where the parameters s μ in (63) are given by the formulas: z 0 = ( u 0 v 1 u 1 v 0 ) w , z 2 = ( u 2 v 0 u 0 v 2 ) w , z 3 = ( u 1 v 2 u 2 v 1 ) w , that are expressed in terms of the arbitrary parameters u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w . These parameters would be also determined and expressed in terms of the initially given physical variables (such as the relative velocity between the reference frames). However, as it has been also noted in Section 2.4 concerning a particular common algebraic property of parametric relations (43-2) and (44-3) which are equivalent to the above expressions, by choosing appropriate integer values for parameters u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w , the parameters z 0 , z 1 , z 2 could take any given integer values. Thus, we may directly determine the relevant expressions for parameters s μ in terms of the initially given physical values and variables. Hence, as a particular case, from the isomorphic transformations (63), in addition to these determined expressions for the parameters s μ (in terms of the relative velocity between the reference frames in x-direction and the speed of light): z 0 = β / ( 1 + γ ) , z 1 = z 2 = 0 , γ = 1 / 1 β 2 , β = v / c , we obtain the equivalent form of Lorentz transformations in the standard configuration [38]:
[ 1 + z 0 2 1 z 0 2 2 z 0 1 z 0 2 2 z 0 1 z 0 2 1 + z 0 2 1 z 0 2 ] [ p 0 p 1 ] = [ γ β γ β γ γ ] [ p 0 p 1 ] = [ p 0 p 1 ]
Similar to the derived transformations (63-1), the Lorentz transformations (in standard configuration) are derived by the same approach for higher-dimensional space-times.

3.3. Matrix Formalism of the Lorentz Invariant Systems of Linear Equations (59)–(62),…, as Equivalent Forms of the Relativistic Energy-Momentum Quadratic Relation in Various Space-Time Dimensions

The Lorentz invariant systems of linear Equations (59)–(62), …, (obtained on the basis of the axiom (17) and relevant general results presented in Section 2.2 and Section 2.4 for homogeneous quadratic equations) as equivalent forms of the Lorentz invariant energy-momentum quadratic relation (52), could be expressed generally by the following matrix formulation in (1 + N)-dimensional space-time:
( α μ p μ m 0 α ˜ μ k μ ) S = 0 ,
where
α μ = β μ + β μ , α ˜ μ = β μ β μ ,
m 0 is the invariant mass of a relativistic particle and k μ = ( c / g 00 , 0 , , 0 ) is its covariant velocity (that is a time-like covariant vector) in the stationary reference frame, α μ and α ˜ μ are two contravariant 2 N × 2 N square matrices (corresponding with the matrix representations of Clifford algebras Cℓ1,2, Cℓ1,3, Cℓ1,4, …, Cℓ1,N (for N ≥ 2) and their generalizations [11,24,32], see also Appendix A) that by the isomorphic linear relations (65) are expressed in terms of two corresponding contravariant 2 N × 2 N matrices β μ and β μ , and S is a 2 N × 1 parametric column matrix. These matrices in (1 + 1), (1 + 2), (1 + 3), (1 + 4) and (1 + 5) space-time dimensions are given uniquely as follows, respectively:
For (1 + 1)-dimensional case we get:
β 0 = [ 0 0 0 1 ] , β 0 = [ 1 0 0 0 ] , β 1 = [ 0 1 0 0 ] , β 1 = [ 0 0 1 0 ] , S = [ ( u 0 v 1 u 1 v 0 ) w s ] ;
where u 0 , u 1 , v 0 , v 1 , w , s are arbitrary parameters.
For (1 + 2)-dimensional case we obtain where u 0 , u 1 , u 2 , v 0 , v 1 , v 2 , w , s are arbitrary parameters):
β 0 = [ 0 0 0 ( σ 0 + σ 1 ) ] , β 0 = [ σ 0 + σ 1 0 0 0 ] , β 1 = [ 0 σ 2 σ 2 0 ] , β 1 = [ 0 σ 3 σ 3 0 ] , β 2 = [ 0 σ 1 σ 0 0 ] , β 2 = [ 0 σ 0 σ 1 0 ] , S = [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 1 v 2 u 2 v 1 ) w s ] , σ 0 = [ 1 0 0 0 ] , σ 1 = [ 0 0 0 1 ] , σ 2 = [ 0 1 0 0 ] , σ 3 = [ 0 0 1 0 ]
For (1 + 3)-dimensional case we obtain:
β 0 = [ 0 0 0 ( γ 0 + γ 1 ) ] , β 0 = [ ( γ 0 + γ 1 ) 0 0 0 ] , β 1 = [ 0 γ 2 γ 3 0 ] , β 1 = [ 0 γ 3 γ 2 0 ] , β 2 = [ 0 γ 4 γ 5 0 ] , β 2 = [ 0 γ 5 γ 4 0 ] , β 3 = [ 0 γ 6 γ 7 0 ] , β 3 = [ 0 γ 7 γ 6 0 ] , S = [ ( u 0 v 1 u 1 v 0 ) w ( u 0 v 2 u 2 v 0 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 3 v 2 u 2 v 3 ) w ( u 1 v 3 u 3 v 1 ) w ( u 2 v 1 u 1 v 2 ) w s ] , γ 0 = [ 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 ] , γ 1 = [ 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 ] , γ 2 = [ 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 ] , γ 3 = [ 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 ] , γ 4 = [ 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 ] , γ 5 = [ 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 ] , γ 6 = [ 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 ] , γ 7 = [ 0 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 ] .
where u 0 , u 1 , u 2 , u 3 , v 0 , v 1 , v 2 , v 3 , w , s are arbitrary parameters. Moreover, the 4 × 4 matrices γ i (68) generate the Lorentz Lie algebra in (1 + 3) dimensions.
For (1 + 4)-dimensional case we have:
β 0 = [ 0 0 0 ( η 0 + η 1 ) ] , β 0 = [ ( η 0 + η 1 ) 0 0 0 ] , β 1 = [ 0 η 2 η 2 0 ] , β 1 = [ 0 η 3 η 3 0 ] , β 2 = [ 0 η 4 η 5 0 ] , β 2 = [ 0 η 5 η 4 0 ] , β 3 = [ 0 η 6 η 7 0 ] , β 3 = [ 0 η 7 η 6 0 ] , β 4 = [ 0 η 8 η 9 0 ] , β 4 = [ 0 η 9 η 8 0 ] , S = [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 4 v 0 u 0 v 4 ) w 0 0 0 0 ( u 3 v 4 u 4 v 3 ) w ( u 2 v 4 u 4 v 2 ) w ( u 1 v 4 u 4 v 1 ) w ( u 2 v 3 u 3 v 2 ) w ( u 1 v 3 u 3 v 1 ) w ( u 1 v 2 u 2 v 1 ) w s ] , η 0 = [ 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] ,   η 1 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 ] , η 2 = [ 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] , η 3 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 ] , η 4 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 ] , η 5 = [ 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] , η 6 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 ] , η 7 = [ 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] , η 8 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 ] , η 9 = [ 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] .
where u 0 , u 1 , u 2 , u 3 , u 4 , v 0 , v 1 , v 2 , v 3 , v 4 , w , s are arbitrary parameters. Furthermore, similar to the 4 × 4 γ i matrices in (68), the 8 × 8 matrices η i (69) generate the Lorentz Lie algebra in (1 + 4) dimensions.
For (1 + 5)-dimensional case the size of matrices β μ and β μ is 32 × 32 . S is also a 32 × 1 column matrix, given by:
S = [ S S ] ,   S = [ ( u 0 v 1 u 1 v 0 ) w ( u 0 v 2 u 2 v 0 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 0 v 4 u 4 v 0 ) w 0 0 0 ( u 0 v 5 u 5 v 0 ) w 0 0 0 0 0 0 0 ] ,   S = [ 0 0 0 ( u 5 v 4 u 4 v 5 ) w 0 ( u 3 v 5 u 5 v 3 ) w ( u 5 v 2 u 2 v 5 ) w ( u 1 v 5 u 5 v 1 ) w 0 ( u 4 v 3 u 3 v 4 ) w ( u 2 v 4 u 4 v 2 ) w ( u 4 v 1 u 1 v 4 ) w ( u 2 v 3 u 3 v 2 ) w ( u 3 v 1 u 1 v 3 ) w ( u 1 v 2 u 2 v 1 ) w s ] .
where u 0 , u 1 , u 2 , u 3 , u 4 , u 5 , v 0 , v 1 , v 2 , v 3 , v 4 , v 5 , w , s are arbitrary parameters.
Similar to the formulations (66)–(70), for the higher dimensional cases of invariant quadratic relation (52), the column matrix S and square matrices β μ and β μ (defining the square matrices α μ and α μ that correspond to the matrix representations of Clifford algebras and their generalization, see Section 3.3 and also Appendix A) are obtained using similar algebraic structures, where in (1 + N) space-time dimensions the size of square matrices β μ and β μ is 2 N × 2 N and the size of column matrix S is 2 N × 1 .

Remark 1. General Algebraic Formulation of the Column Matrix S Given bythe Matrix Equation (64)

In addition, as noted in Section 3.3, the matrix Equation (64) represents uniquely the equivalent form of the Lorentz invariant energy-momentum quadratic relation (52) (as the norm of the D-momentum), based on the axiomatic relations (17) and (17-1) and relevant general results obtained in Section 2.2 and Section 2.4 for homogeneous quadratic equations over the integral domain over ℤ. Therefore (as mentioned in Section 3.3), the general algebraic formulation of the entries of column matrices S obtaining in subsequent higher space-time dimensions, are similar to formulations of the obtained matrices S (66)–(70) corresponding, respectively, to the (1 + 0), (1 + 1), (1 + 2), (1 + 3), (1 + 4) and (1 + 5)-dimensional cases of Lorentz invariant matrix equation (64). Hence, the algebraic formulation of column matrix S in (1 + N) space-time dimensions would be generally defined as follows: the last entry of S is represented solely by the arbitrary parameter s , 2 N 1 entries are definitely zero (see below) and all the other 2 N 1 1 entries of S could be represented uniformly by the following unique algebraic formulation (expressing in terms of the arbitrary parameters: u 0 , u 1 , u 2 , u 3 , , u N , v 0 , v 1 , v 2 , v 3 , , v N , w ) given on the basis of a one-to-one correspondence between these (non-zero) entries of matrix S and the entries h μ ν (for μ > ν) of a 2 N + 1 × 2 N + 1 square matrix H [ h μ ν ] defined in (1 + N) dimensions, by:
h μ ν = ( u ν v μ u μ v ν ) w
where μ, ν = 0, 1, 2,…, N, and h μ ν = 0 for μ = ν.
Note that the algebraic form (70-1) is equivalent to form (45-2) which, as it has been noted in Section 2.4, generates a symmetric algebra Sym(V) on the vector space V , where ( u μ , v ν ) V [37].
Hence, as a basic algebraic property of the form (70-1), a natural unique isomorphism is defined between the underlying vector space V of the symmetric algebra Sym(V) (which is generated by algebraic form (70-1)) and the Weyl algebra W(V). Moreover, based on this isomorphism, the Weyl algebra W(V) could be defined as a (first) quantization of the symmetric algebra Sym(V), where the generators of the Weyl algebra W(V) would be represented by the corresponding (covariant) differential operators (such as i μ , as per quantum mechanics usage).
In Section 3.4, we use these general and basic algebraic properties of the column matrix S , in particular, in the procedure of quantization of the algebraic matrix Equation (64).
In addition to the above algebraic properties of the parametric entries of column matrix S , that are represented uniformly by the algebraic Formula (70-1), in terms of the arbitrary parameters: u 0 , u 1 , u 2 , u 3 , , u N , v 0 , v 1 , v 2 , v 3 , , v N , w , the following basic properties hold as well:
Displaying the column matrix S by two half-sized 2 N 1 × 1 column matrices S and S (containing respectively the upper and lower entries of S , similar to the Formula (57-2) and (62-2) representing the (1 + 5)-dimensional case of matrix S ) such that: S = [ S S ] , then we have:
(1)
The number of entries of the column matrix S that are zero, is exactly: ( 2 N 1 N ) , and the other N entries are represented solely either by the formulation: h μ 0 = ( u 0 v μ u μ v 0 ) w , or by its negative form, i.e.,: h μ 0 = h 0 μ = ( u μ v 0 u 0 v μ ) w , where μ = 1, 2,…, N, and h μ 0 denote the N entries (except the first entry h 00 that is zero) of the first column of square matrix H [ h μ ν ] (defined by the Formula (70-1));
(2)
The number of entries of the column matrix S that are zero, is exactly: ( 2 N 1 N ( N 1 ) 2 1 ) , and except the last entry (represented by arbitrary parameter s ), all the other ( N ( N 1 ) 2 ) entries are represented solely either by the formulation: h μ ν = ( u ν v μ u μ v ν ) w , or by its negative form, i.e.,: h μ ν = h ν μ = ( u μ v ν u ν v μ ) w , where μ > ν, μ, ν = 1, 2,…, N, and h μ ν denote the components of square matrix H [ h μ ν ] , and the last entry of column matrix S is also represented by the arbitrary parameter s .
(3)
If we exchange S and S in the column matrix S = [ S S ] , that could be shown by,
S ( C h ) = [ 0 I I 0 ] S = [ 0 I I 0 ] [ S S ] = [ S S ]
then based on the general formulation of matrix S (defined uniquely by Formula (66)–(70) for various space-time dimensions) and its algebraic properties (1) and (2) (mentioned above), it is concluded directly that the matrix Equation (64) given with the new column matrix S ( C h ) (70-2), i.e., equation: ( α μ p μ m 0 α ˜ μ k μ ) S ( C h ) = 0 , is which could be defined solely in (1 + 2) space-time dimensions for s = 0 , u 1 = 0 , v 1 = 0 , in (1 + 3) space-time dimensions for s = 0 , and in (1 + 4) space-time dimensions for s = 0 , u 1 = 0 , v 1 = 0 (which is reduced and be equivalent to the (1 + 3)-dimensional case of matrix Equation (64)). In (1 + 1) and (1 + 5) and higher space-time dimensions, the matrix equation ( α μ p μ m 0 α ˜ μ k μ ) S ( C h ) = 0 would be which are defined if and only if all the entries of column matrix S ( C h ) are zero. This means that the matrix Equation (64): ( α μ p μ m 0 α ˜ μ k μ ) S = 0 , is symmetric in the exchange of S and S (in the column matrix S = [ S S ] ), solely in (1 + 2)-dimensional space-time for s = 0 ,   u 1 = 0 ,   v 1 = 0 , and in (1 + 3)-dimensional space-time for s = 0 . In Section 3.12, this particular algebraic property of the column matrix S would be used for concluding a new crucial and essential issue in fundamental physics.
In the following Section, the natural isomorphism between the symmetric algebra Sym(V) (generated uniquely by the algebraic form (70-1)) and the Weyl algebra W(V), in addition to the general algebraic properties of column matrix S mentioned above (in Remark 1), would be used and applied directly in the procedure of first quantization of the Lorentz invariant system of linear (64).

3.4. A New Unique Mathematical Derivation of the Fundamental (Massive) Field Equations of Physics (Representing the Laws Governing the Fundamental Forces of Nature)

By first quantization (followed by a basic procedure of minimal coupling to space-time geometry) of the Lorentz invariant system of linear Equation (64) (representing uniquely the equivalent form of energy-momentum quadratic relation (52), see Section 3.3) expressed in terms of the Clifford algebraic matrices (65)–(70),..., two classes of general covariant field equations are derived uniquely as follows (given by (1 + N) space-time dimensions):
( i α μ μ m 0 ( R ) α ˜ μ k μ ) Ψ R = 0 ,
( i α μ D μ m 0 ( F ) α ˜ μ k μ ) Ψ F = 0
where i μ and i D μ are the general relativistic forms of energy-momentum quantum operator (where μ is the general covariant derivative, and D μ is gauge covariant derivative, for detail see the ordinary tensor formalisms of these equations, represented by Formula (78-1)–(79-3), in Section 3.5), m 0 ( R ) and m 0 ( F ) are the fields’ invariant masses, k μ = ( c / g 00 , 0 , , 0 ) is the general covariant velocity in stationary reference frame (that is a time-like covariant vector), α μ and α ˜ μ are two contravariant 2 N × 2 N square matrices (compatible with the matrix representations of certain Clifford algebras, see Section 3.3 and also Appendix A) defined by Formula (65)–(70) in Section 3.3. In the field Equation (72), ΨR is a column matrix as a (first) quantized form of the algebraic column matrix S (defined by relations (64)–(70-2)), determined and represented uniquely by Formula (73)–(77), …, in various space-time dimensions. The column matrix ΨR contains the components of field strength tensor R μ ν ρ σ (equivalent to the Riemann curvature tensor), and also the components of covariant quantity ϕ ρ σ ( G ) that defines the corresponding source current tensor by relation: J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) (which appears in the course of the derivation of field Equation (71), see Section 3.6 for details). In a similar manner, in the tensor field Equation (72), ΨF is also a column matrix as a (first) quantized form of the algebraic column matrix S (defined by relations (64)–(70-2)), determined and represented uniquely by Formula (73)–(77), …, in to various space-time dimensions. The column matrix ΨF contains both the components of tensor field F μ ν (defined as the gauge field strength tensor), and also the components of covariant quantity ϕ ( F ) that define the corresponding source current vector by relation: relation: J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) (which appears in the course of the derivation of field Equation (72), see Section 3.6 for details). Moreover, the general covariance formalism of the field Equations (71) and (72), would be also shown in Section 3.5.
In addition, in Section 3.11, based on a basic class of discrete symmetries for the field Equations (71) and (72), along with definite mathematical axiomatic formalism of the derivation of these equations, it is shown that these equations could be defined solely in (1 + 2) and (1 + 3) space-time dimensions. It is shown that (1 + 3) dimensional cases of these equations represent uniquely a new formalism of bispinor fields of spin-2 and spin-1 particles, respectively. It is also shown that the (1 + 2)-dimensional cases of these equations, represent asymptotically new massive forms of bispinor fields of spin-3/2 and spin-1/2 particles, respectively.
Moreover, in Section 3.12, based on the definite mathematical formalism of this axiomatic derivation approach, the basic assumption (3) in Section 3.1, along with the C, P and T symmetries (represented basically by their corresponding quantum matrix operators) of the fundamentally derived general covariant field Equations (71) and (72), it is concluded that the universe could be realized solely by the (1 + 2) and (1 + 3)-dimensional space-times (where this conclusion, in particular, is based on the T-symmetry). In Section 3.13 and Section 3.14, it is proved that ‘CPT’ is the only (unique) combination of C, P, and T symmetries that could be defined as a symmetry for interacting fields. In addition, in Section 3.14, on the basis of these discrete symmetries of the field equations (71) and (72), it is shown that only left-handed particle fields (along with their complementary right-handed fields) could be coupled to the corresponding (any) source currents.
Furthermore, in Section 3.15, it is argued that the metric of background curved space-time is diagonalized for the spin-1/2 fermion field equations (defined by the field Equation (110) as a generalized form of (1 + 2)-dimensional case of Equation (72)), where this property generates a certain set of additional symmetries corresponding uniquely to the SU(2)L⊗U(2)R symmetry group for spin-1/2 fermion fields (represented by two main groups of “1 + 3” generations, corresponding respectively to two subgroups of leptons and two subgroups of quarks), in addition to the SU(2)L⊗U(2)R and SU(3) gauge symmetry groups for spin-1 boson fields coupled to the spin-1/2 fermionic source currents. Moreover, based on these uniquely determined gauge symmetries, four new charge-less spin-1/2 fermions (represented by “ze , zn ; zu , zd”, where two right-handed charhe-less quarks zu and zd emerge specifically in two subgroups with anti-quarks such that: ( s ¯ ,   u ¯ ,   b ¯ , zu) and ( c ¯ ,   d ¯ ,   t ¯ , zd), and also three new massive spin-1 bosons (represented by W ˜ + , W ˜ , Z , where in particular Z is the complementary right-handed particle of ordinary Z boson), are predicted by this new mathematical axiomaticapproach.
As a particular result, in Section 3.6, based on the definite and unique formulation of the derived Maxwell’s equations (and also Yang-Mills equations, defined by the (1 + 3)-dimensional case of the field Equation (72), compatible with specific gauge symmetry groups as shown in Section 3.15, Formula (109)–(110-12)), it is also concluded that magnetic monopoles could not exist in nature.

3.5. Axiomatic Derivation of General Covariant Massive Field Equations (71) and (72)

First, it should be noted that via initial quantization (followed by a basic procedure of minimal coupling to space-time geometry) of the algebraic systems of linear Equation (64) (as a matrix equation given by the Clifford algebraic matrices (65)–(70), …, in various space-time dimensions), two categories of general covariant field Equations (with a definite matrix formalism compatible with the Clifford algebras and their generalizations, see Section 3.3 and also Appendix A) are derived solely, represented by the tensor Equations (71) and (72) in terms of two tensor fields R ρ σ μ ν and F μ ν , respectively. In fact, as it has been mentioned in Remark 1 (in Section 3.3), there is a natural isomorphism between the Weyl algebra and the symmetric algebra generated by the algebraic form (70-1) which represents the general formulation of the entries of algebraic column matrix S in the matrix Equation (64). In addition, the procedure of minimal coupling to space-time geometry would be simply defined as a procedure which, starting from a theory in flat space-time, substitutes all partial derivatives by corresponding covariant derivatives and the flat space-time metric by the curved space-time (pseudo-Riemannian) metric. Moreover, as mentioned in Remark 1 (in Section 3.3), on the basis of this natural isomorphism, the Weyl algebra could be also represented as a quantization of the symmetric algebra generated by the algebraic form (70-1). Hence, using this natural isomorphism, by first quantization (followed by a basic procedure of minimal coupling to space-time geometry) of matrix Equation (64), two definite classes of general covariant massive (tensor) field equations are determined uniquely, expressed in terms of two basic connection forms (denoted by two derivatives μ and D μ corresponding respectively to the diffeomorphism (or metric) invariance and gauge invariance), along with their corresponding curvature forms, denoted respectively by R ρ σ μ ν (as the gravitational field strength tensor, equivalent to Riemann curvature tensor) and F μ ν (as the gauge field strength tensor). This natural isomorphism could be represented by the following mappings (corresponding to the curvature forms R ρ σ μ ν and F μ ν , respectively):
( u μ v ν u ν v μ ) w ( μ ν ν μ ) ω R = R ( μ , ν ) ω R R ρ σ μ ν ,
( u μ v ν u ν v μ ) w ( D μ D ν D ν D μ ) ω F ( i g F ) F μ ν .
where R σ μ ν ρ = ( ν Γ σ μ ρ + Γ λ ν ρ Γ σ μ λ ) ( μ Γ σ ν ρ + Γ λ μ ρ Γ σ ν λ ) , F μ ν = D ν A μ D μ A ν , and g F , A μ are respectively the corresponding coupling constant and gauge field (that is defined generally as a Lie algebra-valued 1-form represented by a unique vector field [39]). Based on this natural unique isomorphism represented by the mappings (71-1) and (72-1), the column matrices ΨR and ΨF (in the expressions of field Equations (71) and (72), respectively) would be determined uniquely various dimensional space-times, represented by Formula (73)–(77), ….
In addition, as detailed in Section 3.4, the last entry of algebraic column matrix S in matrix Equation (64) (as it has been shown in the relations (64)–(70)), is represented by the arbitrary algebraic parameter s . In the course of the derivation of field Equations (71) and (72) (via the first quantization procedure mentioned above, and the mappings (71-1) and (72-1)), the arbitrary parameter s could be substituted solely by two covariant quantities ϕ ρ σ ( R ) and ϕ ( F ) that define the corresponding covariant source currents ϕ ρ σ ( R ) and J ν ( F ) (given bythe field Equations (71) and (72), respectively) by the conditional relations: J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) and J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) .
As another basic issue concerning the general covariance formulation of tensor field Equations (71) and (72), we should note that each of these Equations (as a system of equations) includes also an equation corresponding to the 2nd Bianchi identity, as follows, respectively:
( λ + i m 0 ( R ) k λ ) R ρ σ μ ν + ( μ + i m 0 ( R ) k μ ) R ρ σ ν λ + ( ν + i m 0 ( R ) k ν ) R ρ σ λ μ = 0 ,
( D λ + i m 0 ( F ) k λ ) F μ ν + ( D μ + i m 0 ( F ) k μ ) F ν λ + ( D ν + i m 0 ( F ) k ν ) F λ μ = 0
However, the tensor field R ρ σ μ ν as the Riemann curvature tensor, obeys the relation (71-2) tensor, if and only if a torsion tensor is defined in as: T τ μ ν = ( i m 0 ( R ) / 2 ) ( g τ μ k ν g τ ν k μ ) , and subsequently the relation (71-2) be equivalent to the 2nd Bianchi identity of the Riemann tensor. Consequently, the covariant derivative ν should be also defined with this torsion, that we may show it by ν . Moreover, as presented in Section 3.6, concerning the relation (72-2), we may also define a torsion field as: Z τ μ ν = ( i m 0 ( F ) / 2 ) ( g τ μ k ν g τ ν k μ ) , and write the relations (71-2) and (72-2) (representing the 2nd Bianchi identities) as follows:
λ R ρ σ μ ν + μ R ρ σ ν λ + ν R ρ σ λ μ = T λ μ τ R ρ σ τ ν + T μ ν τ R ρ σ τ λ + T ν λ τ R ρ σ τ μ ,
D λ F μ ν + D μ F ν λ + D ν F λ μ = 0
where the general relativistic form of gauge derivative D μ has been defined with torsion field Z τ μ ν . We use the derivatives μ and D μ in the ordinary tensor representations (i.e., the Formula (78-1)–(79-3)) of the field Equations (71) and (72) in Section 3.6. In addition, based on the formulations of torsions T τ μ ν and Z τ μ ν (that have appeared naturally in the course of derivation of the field Equations (71) and (72)) and the general properties of torsion tensors (in particular, the property that allows a torsion tensor to always be treated as an independent tensor field, or equivalently, as part of the space-time geometry [40,41,42]), it could be concluded directly that torsion field T τ μ ν generates the invariant mass of corresponding gravitational field, and torsion field Z τ μ ν generates the invariant mass of corresponding gauge field, respectively. Hence, based on our axiomatic derivation approach including the mappings (71-1) and (72-1) (mentioned above), the (1 + 1), (1 + 2), (1 + 3), (1 + 4), (1 + 5), …, dimensional cases of column matrices ΨR and ΨF in the specific expressions of general covariant massive (tensor) field Equations (71) and (72), are determined uniquely as follows, respectively; For (1 + 1)-dimensional space-time we have:
Ψ R = [ R ρ σ 10 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 ϕ ( F ) ] ,   J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ; ,
For (1 + 2)-dimensional space-time we obtain:
Ψ R = [ R ρ σ 10 R ρ σ 02 R ρ σ 21 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 F 02 F 21 ϕ ( F ) ] ,   J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ;
For (1 + 3)-dimensional space-time we have:
Ψ R = [ R ρ σ 10 R ρ σ 20 R ρ σ 30 0 R ρ σ 23 R ρ σ 31 R ρ σ 12 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 F 20 F 30 0 F 23 F 31 F 12 ϕ ( F ) ] ,   J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ;
For (1 + 4)-dimensional space-time we get:
Ψ R = [ R ρ σ 10 R ρ σ 02 R ρ σ 30 0 R ρ σ 04 0 0 0 0 R ρ σ 43 R ρ σ 42 R ρ σ 41 R ρ σ 32 R ρ σ 31 R ρ σ 21 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 F 02 F 30 0 F 04 0 0 0 0 F 43 F 42 F 41 F 32 F 31 F 21 ϕ ( F ) ] ,   J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ;
For (1 + 5)-dimensional space-time we obtain:
Ψ R = [ R ρ σ 10 R ρ σ 02 R ρ σ 30 0 R ρ σ 04 0 0 0 0 R ρ σ 43 R ρ σ 42 R ρ σ 41 R ρ σ 32 R ρ σ 31 R ρ σ 21 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 F 02 F 30 0 F 04 0 0 0 0 F 43 F 42 F 41 F 32 F 31 F 21 ϕ ( F ) ] ,   J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) ;   Ψ R = [ R ρ σ 10 R ρ σ 20 R ρ σ 30 0 R ρ σ 40 0 0 0 R ρ σ 50 0 0 0 0 0 0 0 0 0 0 R ρ σ 45 0 R ρ σ 53 R ρ σ 25 R ρ σ 51 0 R ρ σ 34 R ρ σ 42 R ρ σ 14 R ρ σ 32 R ρ σ 13 R ρ σ 21 ϕ ρ σ ( R ) ] , Ψ F = [ F 10 F 20 F 30 0 F 40 0 0 0 F 50 0 0 0 0 0 0 0 0 0 0 F 45 0 F 53 F 25 F 51 0 F 34 F 42 F 14 F 32 F 13 F 21 ϕ ( F ) ] , J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F )
where in the relations (73)–(77), J ν ρ σ ( R ) and J ν ( F ) are the source currents expressed, necessarily, in terms of the covariant quantities ϕ ρ σ ( G ) and ϕ ( F ) (as the initially given quantities), respectively. For higher-dimensional space-times, the column matrices Ψ R and Ψ F (with similar formulations) are determined uniquely as well.

3.6. General Covariant (Tensor) Representations of the Field Equations (71) and (72)

From the Field Equations (71) and (72) (derived uniquely with certain matrix formalisms compatible with the Clifford and Weyl algebras), the following general covariant field equations, with ordinary tensor formalisms, are obtained (but not vice versa), respectively:
λ R ρ σ μ ν + μ R ρ σ ν λ + ν R ρ σ λ μ = T λ μ τ R ρ σ τ ν + T μ ν τ R ρ σ τ λ + T ν λ τ R ρ σ τ μ ,
μ R ρ σ μ ν ( i m 0 ( R ) / ) k μ R ρ σ μ ν = J ρ σ ( R ) ν
R σ μ ν ρ = ( ν Γ σ μ ρ + Γ λ ν ρ Γ σ μ λ ) ( μ Γ σ ν ρ + Γ λ μ ρ Γ σ ν λ ) , J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , T τ μ ν = i m 0 ( R ) 2 ( g τ μ k ν g τ ν k μ ) .
and
D λ F μ ν + D μ F ν λ + D ν F λ μ = 0 ,
D μ F μ ν = J ν ( F )
F μ ν = D ν A μ D μ A ν J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) , Z τ μ ν = i m 0 ( F ) 2 ( g τ μ k ν g τ ν k μ ) .
where Γ σ μ ρ is the affine connection: Γ σ μ ρ = Γ ¯ σ μ ρ K σ μ ρ , Γ ¯ σ μ ρ is the Christoffel symbol (or the torsion-free connection), K σ μ ρ is the contorsion tensor defined by: K ρ σ μ = ( i m 0 ( R ) / 2 ) g ρ μ k σ (that is anti-symmetric in the first and last indices), T ρ σ μ is the torsion given by: T ρ σ μ = K ρ μ σ K ρ σ μ (that generates the invariant mass of the gravitational field), k μ = ( c / g 00 , 0 , , 0 ) (where we supposed c = 1 ) is the covariant velocity of particle (or the static observer) in the stationary reference frame, and A μ is the gauge potential vector field. Moreover, in general covariant field Equations (79-1)–(79-3), the covariant derivative D μ has been defined specifically with the torsion field Z τ μ ν (generating the invariant mass of gauge field strength tensor F μ ν ).
It should be emphasized again that the tensor field Equations (78-1)–(78-3) and (79-1)–(79-3) (which are obtained respectively from the original Equations (71) and (72), but not vice versa) show merely the general covariance formalism (including torsions fields T μ ν τ and Z μ ν τ ) of the axiomatically derived field Equations (71) and (72). The crucial issue here is that the original field Equations (71) and (72) could not be obtained from the tensor Equations (78-1)–(78-3) and (79-1)–(79-3). In fact, the tensor Equations (78-1)–(78-3) and (79-1)–(79-3) don’t represent completely the definite matrix formalism (compatible with certain Clifford and Weyl algebras) of the axiomatic field Equations (71) and (72). Hence, based on this mathematical axiomatic formalism and derivation approach of Equations (71) and (72) (presented in Section 3.3, Section 3.4 and Section 3.5), it is concluded that the fundamental force fields of physics cannot be described completely via the ordinary tensor representations of these fields (in the current standard classic and quantum relativistic field theoretic formalism of physics), such as the representations (78-1)–(78-3) and (79-1)–(79-3); and as shown in Section 3.3, Section 3.4 and Section 3.5, on the basis of this new mathematical formalism, all the fundamental force fields of physics could be represented (and described) solely by the axiomatically determined and formulated field Equations (71) and (72) with their definite covariant matrix formalisms (given and specified by Formula (65)–(70) for various space-time dimensions, compatible with certain Clifford and Weyl algebras).

3.7. Derivation of the Einstein field Equations

Along with the massive gravitational field Equations (78-1)–(78-3) (obtained uniquely from the originally derived field Equation (71)) that are expressed solely in terms of R μ ν ρ σ as the field strength tensor and also torsion’s depended terms, we also assume the following relation as basic definition for the Ricci tensor (where the Riemann curvature tensor and Ricci tensor obey the interchange symmetries: R μ ν ρ σ R ρ σ μ ν , R μ ν R ν μ , because of the torsion [43]):
( σ + i m 0 ( R ) k σ ) R μ ν ρ σ = ( ν + i m 0 ( R ) k ν ) R μ ρ ( μ + i m 0 ( R ) k μ ) R ν ρ
where the relation (78-4) particularly remains unchanged by the transformation:
R μ ν R μ ν + Λ g μ ν
(where as would be shown, Λ is equivalent to the cosmological constant). It should be noted that by taking Λ = 0 , from the 2nd Bianchi identity of the Riemann curvature tensor and relation (78-4) it could be shown that the Ricci tensor is also the contraction of the Riemann tensor, i.e., R μ ν = R μ ν σ σ (which is equivalent to the ordinary definition of the Ricci tensor). However, this ordinary definition for the Ricci tensor, necessarily, doesn’t imply the above transformation. In fact, in the following, we show that this basic transformation is necessary for having the cosmological constant in the gravitational field Equations (including the Einstein field equations which could be derived from the above equations and relations) expressed in terms of the Ricci and stress-energy tensors. As a direct result, a unique equivalent expression of gravitational field equations, in terms of the Ricci tensor R μ ν and stress-energy tensor T μ ν , could be also determined from the basic definition (78-4) (for Ricci curvature tensor, based on this axiomatic formalism), and field Equations (78-1)–(78-3), along with the following expression for current J ρ σ ν ( R ) (defined in terms of the stress-energy tensor T μ ν , T ( = T μ μ ) , and metric g μ ν , in D-dimensional space-time):
J ρ σ ν ( R ) = 8 π [ ( σ + i m 0 ( R ) k σ ) T ρ ν ( ρ + i m 0 ( R ) k ρ ) T σ ν ] + 8 π B [ ( σ + i m 0 ( R ) k σ ) T g ρ ν ( ρ + i m 0 ( R ) k ρ ) T g σ ν ] ,
where T μ ν T ν μ for m 0 ( R ) 0 , B = 0 for D = 1 , 2 , and B = 1 / ( D 2 ) for D 3 , the Einstein field Equations (as the massless case) are determined directly as follows:
R μ ν = 8 π ( T μ ν B T g μ ν ) Λ g μ ν

3.8. Showing that Magnetic Monopoles Cannot Exist in Nature

As a direct consequence of the uniquely derived general covariant field Equation (72) that are specified by the matrices (73)–(77) and (65)–(70) (or the general covariant field Equations (79-1)–(79-3) obtained from the original Equation (72)), which, in fact, represent the electromagnetic fields equivalent to a generalized massive form of the Maxwell’s Equations (as well as a generalized massive form of the Yang-Mills fields corresponding to certain gauge symmetry groups, see Section 3.15), it is concluded straightforwardly that magnetic monopoles could not exist in nature.

3.9. Gauge Invariance of (Massive) Tensor Field Equations(79-1)–(79-3) Corresponding to the Maxwell’s (and Yang-Mills) Equations in (1 + 3) Dimensions, and Dirac Equation in (1 + 2) Dimensions

The Lagrangian density specified for the tensor field F μ ν in the field Equations (79-1)–(79-3) is (supposing J ν ( F ) = 0 ) [39]:
L ( F ) = ( 1 / 4 g ) F μ ν F μ ν
where g is the metric’s determinant. Moreover, the trace part of torsion field Z τ μ ν in (79-3) is obtained as:
Z μ ν μ = Z ν = N ( i m 0 ( F ) / 2 ) k ν = N α k ν
where (1 + N) is the number of space-time dimensions and α = i m 0 ( F ) 2 . Now based on the definition of covariant vector k μ (as a time-like covariant vector), we simply get: ϕ : k ν = ν ϕ . This basic property, along with and Formula (81), imply the general covariant massive field Equations (79-1)–(79-3) (formulated originally with the torsion field (79-3) generating the invariant mass m 0 ( F ) of field F μ ν ), and the corresponding Lagrangian density (80), be invariant under the U(1) Abelian gauge group [17,39,44,45,46,47,48]. However, in Section 3.15, we show that assuming the spin-1/2 fermion fields (describing generally by the field Equation (110-9) compatible with specific gauge symmetry group (110-12), as shown in Section 3.15, Formula (110-9)–(110-12)) and their compositions as the source currents of the (1 + 3)-dimensional cases of general covariant massive field Equation (72) (describing the spin-1 boson field), then this field equation would be invariant under two types of gauge symmetry groups, including: SU(2)L⊗U(2)R and SU(3), corresponding with a group of seven bosons and a groups of eight bosons (as shown in Section 3.15, Formula (114-4)–(114-9)).

3.10. Identifying a New Particular Massive Gauge Boson

According to Refs. [44,45,46,47,48], in agreement with the recent astronomical data, we can directly establish a lower bound for a constant quantity which is equivalent to the constant α = i m 0 ( F ) 2 (defined by the relation (80)) as: | α | 21 . Hence, a new massive particle (corresponding to the U(1) symmetry group) would be identified with the invariant mass:
m γ 4.90571 × 10 50 kg
that is generated by a coupling torsion field of the type (79-3) of the background curved space-time. In addition, it should be noted that, in general, based on the covariant massive field Equations (71) and (72) derived by our axiomatic approach (or field Equations (78-1)–(78-3) and (79-1)–(79-3) obtained from (71) and (72)), the invariant masses of the elementary particles are generated by torsion fields of the types (78-3) (for spin-3/2 and spin-2 particles) and (79-3) (for spin-1/2 and spin-1 particles, see Section 3.15). Hence, this approach could be also applied for massive neutrinos concluding that their masses are generated by the coupling torsion fields (of the type (79-3)). Such massive particle fields coupled with the torsions (of the type (79-3)) of the background space-time geometry could be completely responsible for the mysteries of dark energy and dark matter [49,50].

3.11. Quantum Representations of C, P and T Symmetries of the Axiomatically Derived General Covariant Massive (Tensor) Field Equations (71) and (72)

As shown in Section 3.3, Section 3.4 and Section 3.5, the general covariant massive (tensor) field Equations (71) and (72) as the unique axiomatically determined Equations (representing the fundamental field of physics, as assumed in Section 3.1), are represented originally with definite matrix formalisms constructed from the combination of two specific matrix classes including the column matrices (73)–(77),… compatible with the Weyl algebras (based on the isomorphism (71-1)–(72-1)), and the square matrices (65)–(70), … that are compatible with the Clifford algebras and their generalizations; see Section 3.3, Section 3.4 and Section 3.5 and also Appendix A for detail).
In agreement with the principles of relativistic quantum theory [6], and also as another primary assumption in addition to the basic assumption (3) defined in Section 3.1, we basically represent the C, P and T symmetries of the source-free cases of by the following quantum matrix operators (with the same forms in both flat and curved space-time), respectively:
(Note: In Section 3.13, we show that only a certain simultaneous combination of the C, P and T transformations could be defined for the field Equations (71) and (72) with non-zero source currents.)
(1)- 
Parity Symmetry (P-Symmetry):
P ^ = γ P = [ I 0 0 I ]
where I is the identity matrix, and the size of matrix γ P in (1 + N)-dimensional space-time is 2 N × 2 N . The operator P ^ obeys the relations:
det ( P ^ ) = 1 , P ^ 2 = 1 , P ^ = P ^ 1 = P ^ = P ^ T
(2)- 
Time-Reversal Symmetry (T-Symmetry):
T ^ = T ^ 0 K ^ = i γ P γ C h K ^
where the operator K ^ denotes complex conjugation, the operator γ P defined by Formula (83) and the operator γ C h in (1 + 1) and (1 + 2) space-time dimensions, is given by:
γ C h = [ 0 I I 0 ] ,
and in (1 + 3) and higher space-time dimensions, γ C h is denoted by:
γ C h = [ 0 i I i I 0 ]
where the size of matrix γ C h in (1 + N)-dimensional space-time is: 2 N × 2 N . Moreover, in (1 + 1) and (1 + 2) space-time dimensions, the time reversal operator T ^ (84) and the Hermitian operator T ^ 0 = i γ P γ C h (specified in the Formula (84)) obey the relations:
T ^ 2 = 1 , T ^ 0 = T ^ 0 1 = T ^ 0 = T ^ 0 T ,
and in (1 + 3) and higher space-time dimensions, T ^ and T ^ 0 obey the relations:
T ^ 2 = 1 , T ^ 0 = T ^ 0 1 = T ^ 0 = T ^ 0 T
Concerning the time reversal symmetry, it should be noted that relations (84-3) are solely compatible with the fermionic fields and relations (84-4) are solely compatible with the bosonic fields. In addition, it should be noted that these basic quantum mechanical properties (i.e., the relations (84-3) and (84-4)) of the time reversal symmetry (84), are fully compatible with corresponding properties of the field tensors F μ ν and R ρ σ μ ν presented in Section 3.15, where the tensor field F μ ν (describing by general covariant field Equation (72)) represents (asymptotically) solely a massive bispinor field of spin-1/2 particles (as a new general covariant massive formulation of the Dirac equation) in (1 + 2) space-time dimensions, and also represents a massive bispinor field of spin-1—as new massive general covariant (matrix) formulations of both Maxwell and Yang-Mills field equations are compatible with specified gauge symmetry groups—in (1 + 3) space-time dimensions; and tensor field R ρ σ μ ν (describing by general covariant field Equation (71)) represents (asymptotically) solely a bispinor field of spin-3/2 particles (as a new massive general covariant form of the Rarita–Schwinger equation) in (1 + 2) space-time dimensions , and also represents a massive bispinor field spin-2 particles (equivalent to a generalized massive form of the Einstein equations) in (1 + 3) space-time dimensions.
(3)- 
Charge Conjugation Symmetry (C-Symmetry):
( Ψ R ) C = C ˜ Ψ R = i I K ^ Ψ R , ( Ψ F ) C = = C ˜ Ψ F = i I K ^ Ψ F
where C ˜ = i I K ^ , I is the identity matrix, the operator K ^ denotes complex conjugation, and the charge conjugation operator C ^ defined by: ( Ψ R ) C = C ^ ( Ψ ¯ R ) T , ( Ψ F ) C = C ^ ( Ψ ¯ F ) T . The charge conjugation operator C ^ obeys the following relations:
C ^ C ^ = 1 , C ^ = C ^ 1 = C ^ = C ^ T
As a basic additional issue, it is worth noting that the time-reversal operator (84) could be also expressed basically in terms of the parity matrix operator γ P (83), matrix operator C ˜ = i I K ^ given bythe definition of charge-conjugated transformation (85), and matrix γ C h defined by Formula (84-1) and (84-2), as follows:
γ P γ C h C ˜ = T ^
where we have: γ P γ C h = γ C h γ P .

Remark 2. The Main Properties of Matrix Operator γ C h (Defined by Unitary Matrices (84-1) and (84-2))

In this Section, the main properties of matrix operator γ C h (defined by unitary matrices (84-1) and (84-2)) have been presented. Each of the general covariant (tensor) field Equations (71) and (72) (including their source-free and non-source-free cases), as a system of differential equations, is symmetric and has the same spectrum as a result of multiplying by matrix γ C h . The multiplied column matrices Ψ R ( C h ) = γ C h Ψ R and Ψ F ( C h ) = γ C h Ψ F then obey the Equations (71) and (72), respectively, but with opposite sign in mass term such that: ( i α μ μ + m 0 ( R ) α ˜ μ k μ ) Ψ R ( C h ) = 0 , ( i α μ μ + m 0 ( F ) α ˜ μ k μ ) Ψ F ( C h ) = 0 .
As a general additional issue concerning the column matrices Ψ R ( C h ) = γ C h Ψ R and Ψ F ( C h ) = γ C h Ψ F , should be also added that the sign change of the mass terms introduced in the field Equations (71) and (72) is immaterial (the same property also holds for the ordinary formulation of Dirac equation, and so on [51]). In other words, the field Equation (71) of the form ( i α μ μ ± m 0 ( R ) α ˜ μ k μ ) Ψ R = 0 are equivalent, and similarly the field Equation (72) of the form ( i α μ μ ± m 0 ( F ) α ˜ μ k μ ) Ψ F = 0 would be equivalent as well. However, since the algebraic column matrix S in the matrix Equation (64) (derived and represented uniquely in terms of the matrices (66)–(70),… corresponding to various space-time dimensions), is not symmetric by multiplying by matrix γ C h (84-1) and (84-2) (except for (1 + 2) and (1 + 3)-dimensional cases of column matrix S , based on the definite algebraic properties of matrix S presented in Remark 1, (in Section 3.3), it is concluded that except the (1 + 2) and (1 + 3)-dimensional cases of the fundamental field Equations (71) and (72), these field equations could not be defined with the column matrices of the types Ψ R ( C h ) ( = γ C h Ψ R ) and Ψ F ( C h ) ( = γ C h Ψ F ) (if assuming that the column matrices Ψ R and Ψ F are defined with field Equations (71) and (72), i.e., they have the formulations similar to the formulations of originally derived column matrices (73)–(77),… corresponding to various space-time dimensions). This conclusion follows from this fact that the filed Equations (71) and (72) have been derived (and defined) uniquely from the matrix Equation (64) via the axiomatic derivation approach (including the first quantization procedure) presented in Section 3.4 and Section 3.5. In Section 3.12, using this property (i.e., multiplication of column matrices Ψ R and Ψ F defined in the fundamental field Equations (71) and (72), by matrix γ C h (84-1) and (84-2) from the left), this crucial and essential issue would be concluded directly that by assuming the time-reversal invariance of the general covariant filed Equations (71) and (72) (represented by the transformations T ^ Ψ R and T ^ Ψ F , where the quantum operator T ^ is given uniquely by Formula (84), i.e.,: T ^ = T ^ 0 K ^ = i γ P γ C h K ^ ), these fundamental field equations could be defined solely in (1 + 2) and (1 + 3) space-time dimensions (with the column matrices of the forms (96-1) and (98-2), respectively). Subsequently, in Section 3.13, it would be also shown that only a definite simultaneous combination of all the transformations C ^ , P ^ , T ^ and also matrix γ C h (given by quantum operators (83)–(87)) could be defined for the field Equations (71) and (72) with non-zero source currents. In addition, the matrix operator γ C h in (1 + 1) and (1 + 2) space-time dimensions obeys the relations:
( γ C h ) 2 = 1 , γ C h = ( γ C h ) 1 = ( γ C h ) = ( γ C h ) T ,
and in the (1 + 3) and higher dimensions obeys the following relations as well:
( γ C h ) 2 = 1 , γ C h = ( γ C h ) 1 = ( γ C h ) = ( γ C h ) T
Furthermore, in Section 3.14, the matrix γ C h would be also used basically for defining and representing the left-handed and right-handed components of the column field matrices Ψ R and Ψ F defined originally in the field Equations (71) and (72).

3.12. Showing That the Universe Could Be Realized Solely by the (1 + 2) and (1 + 3)-Dimensional Space-Times

The proof of this essential property of nature within the new mathematical axiomatic formalism presented in this article is mainly based on the T-symmetry (represented by quantum matrix operators (84)) of the fundamentally derived general covariant field Equations (71) and (72). As shown in Section 3.11, the source-free cases (as basic cases) of field Equations (71) and (72) are invariant under the time-reversal transformation defined by matrix operator (84). Moreover, in Section 3.13, it would be also shown that these field equations with non-zero source currents are solely invariant under the simultaneous transformations of all the C ^ , P ^ , and T ^ (83)–(85), multiplied by matrix γ C h (given by Formula (84-1) and (84-2)). Now, following the definite mathematical formalism of the axiomatic derivation approach of fundamental field Equations (71) and (72), assuming that any column matrix X R , or Y F , expressible in the tensor formulation of general covariant field Equation (71) or (72), is basically definable, if and only if, it could be also derived originally as a column matrix via the axiomatic derivation approach presented in Section 3.4 and Section 3.5.
On this basis, it could be shown that the time-reversal transformed forms of the column matrices Ψ R and Ψ F given bythe expressions of source-free cases of field Equations (71) and (72), are definable solely in (1 + 2) and (1 + 3)-dimensional space-times. Based on this result, along with the basic assumption (3) in Section 3.1, it is concluded directly that the universe could be realized solely with the (1 + 2) and (1 + 3)-dimensional space-times. We show this in the following in detail.
As noted, in fact, the above conclusion follows directly from the formulations of uniquely determined time-reversal transformed forms of column matrices Ψ R and Ψ F given by the expressions of source-free cases of field Equations (71) and (72). Denoting these column matrices by ϒ R = T ^ Ψ R and ϒ F = T ^ Ψ F , where the time-reversal operator (84) is defined by: T ^ = T ^ 0 K ^ = i γ P γ C h K ^ , they would be determined as follows in various dimensions:
For (1 + 1)-dimensional space-time we have:
ϒ R ( x , t ) = T ^ Ψ R ( x , t ) = T ^ 0 Ψ R ( x , t ) = [ 0 i R ρ σ 10 ( x , t ) ] ,   ϒ F ( x , t ) = T ^ Ψ F ( x , t ) = T ^ 0 Ψ F ( x , t ) = [ 0 i F 10 ( x , t ) ] ;
For (1 + 2)-dimensional space-time it is obtained:
ϒ R ( x , t ) = T ^ Ψ R ( x , t ) = T ^ 0 Ψ R ( x , t ) = [ i R ρ σ 21 ( x , t ) 0 i R ρ σ 10 ( x , t ) i R ρ σ 02 ( x , t ) ] ,   ϒ F ( x , t ) = T ^ Ψ F ( x , t ) = T ^ 0 Ψ F ( x , t ) = [ i F 21 ( x , t ) 0 i F 10 ( x , t ) i F 02 ( x , t ) ] ;
For (1 + 3)-dimensional space-time we get:
ϒ R ( x , t ) = T ^ Ψ R ( x , t ) = T ^ 0 Ψ R ( x , t ) = [ i R ρ σ 23 ( x , t ) i R ρ σ 31 ( x , t ) i R ρ σ 12 ( x , t ) 0 i R ρ σ 10 ( x , t ) i R ρ σ 20 ( x , t ) i R ρ σ 30 ( x , t ) 0 ] ,   ϒ F ( x , t ) = T ^ Ψ F ( x , t ) = T ^ 0 Ψ F ( x , t ) = [ i F 23 ( x , t ) i F 31 ( x , t ) i F 12 ( x , t ) 0 i F 10 ( x , t ) i F 20 ( x , t ) i F 30 ( x , t ) 0 ] ;
For (1 + 4)-dimensional space-time we have:
ϒ R ( x , t ) = T ^ Ψ R ( x , t ) = T ^ 0 Ψ R ( x , t ) = [ 0 i R ρ σ 43 ( x , t ) i R ρ σ 42 ( x , t ) i R ρ σ 41 ( x , t ) i R ρ σ 32 ( x , t ) i R ρ σ 31 ( x , t ) i R ρ σ 21 ( x , t ) 0 i R ρ σ 10 ( x , t ) i R ρ σ 02 ( x , t ) i R ρ σ 30 ( x , t ) 0 i R ρ σ 04 ( x , t ) 0 0 0 ] ,   ϒ F ( x , t ) = T ^ Ψ F ( x , t ) = T ^ 0 Ψ F ( x , t ) = [ 0 i F 43 ( x , t ) i F 42 ( x , t ) i F 41 ( x , t ) i F 32 ( x , t ) i F 31 ( x , t ) i F 21 ( x , t ) 0 i F 10 ( x , t ) i F 02 ( x , t ) i F 30 ( x , t ) 0 i F 04 ( x , t ) 0 0 0 ] ;
For (1 + 5)-dimensional space-time we obtain:
ϒ R ( x , t ) = T ^ Ψ R ( x , t ) = T ^ 0 Ψ R ( x , t ) = [ 0 0 0 i R ρ σ 45 ( x , t ) 0 i R ρ σ 53 ( x , t ) i R ρ σ 25 ( x , t ) i R ρ σ 51 ( x , t ) 0 i R ρ σ 34 ( x , t ) i R ρ σ 42 ( x , t ) i R ρ σ 14 ( x , t ) i R ρ σ 32 ( x , t ) i R ρ σ 13 ( x , t ) i R ρ σ 21 ( x , t ) 0 i R ρ σ 10 ( x , t ) i R ρ σ 20 ( x , t ) i R ρ σ 30 ( x , t ) 0 i R ρ σ 40 ( x , t ) 0 0 0 i R ρ σ 50 ( x , t ) 0 0 0 0 0 0 0 ] ,   ϒ F ( x , t ) = T ^ Ψ F ( x , t ) = T ^ 0 Ψ F ( x , t ) = [ 0 0 0 i F 45 ( x , t ) 0 i F 53 ( x , t ) i F 25 ( x , t ) i F 51 ( x , t ) 0 i F 34 ( x , t ) i F 42 ( x , t ) i F 14 ( x , t ) i F 32 ( x , t ) i F 13 ( x , t ) i F 21 ( x , t ) 0 i F 10 ( x , t ) i F 20 ( x , t ) i F 30 ( x , t ) 0 i F 40 ( x , t ) 0 0 0 i F 50 ( x , t ) 0 0 0 0 0 0 0 ] ;
Now based on the formulations of the derived time-reversal transformed column matrices ϒ R and ϒ F (88)–(92), although they could be expressed merely by the tensor formulations of field Equations (71) and (72), however, except the (1 + 2) and (1 + 3)-dimensional cases of these transformed column matrices, all the other cases cannot be derived originally as a column matrix via the axiomatic derivation approach presented in Section 3.4 and Section 3.5 (following the formulation of originally derived column matrices (73)–(77)). Below this conclusion (and subsequent remarkable results) is discussed in more detail.
In addition, it is also worth noting that on the basis of our derivation approach, since there are not the corresponding isomorphism (that could be represented by the unique mappings (71-1) and (72-1), in Section 3.5) between the entries of column matrices ϒ R and ϒ F (88), (91), (92),… and the entries (with the exactly same indices) of column matrix S (in the algebraic matrix Equation (64), where its last entry, i.e., arbitrary parameter “s” is zero compatible with the source-free cases of the field Equations (71) and (72)) that are given uniquely as follows in (1 + 1) and (1 + 4), (1 + 5),... and higher space-time dimensions, respectively, using the definitions (66)–(70),… (in Section 3.3), and also the algebraic properties of column matrix S (presented in Remark 1 in Section 3.3) representing in terms of two half-sized 2 N 1 × 1 column matrices S and S such that: S = [ S S ] (where u 0 , u 1 , u 2 , u 3 , , u N , v 0 , v 1 , v 2 , v 3 , , v N , w are arbitrary parameters):
For (1 + 1)-dimensional space-time we have:
S = [ S S ] = [ ( u 0 v 1 u 1 v 0 ) w 0 ] ;
For (1 + 4), (1 + 5)-dimensional space-times we get, respectively:
S = [ S S ] ,   S = [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 4 v 0 u 0 v 4 ) w 0 0 0 ] S = [ 0 ( u 3 v 4 u 4 v 3 ) w ( u 2 v 4 u 4 v 2 ) w ( u 1 v 4 u 4 v 1 ) w ( u 2 v 3 u 3 v 2 ) w ( u 1 v 3 u 3 v 1 ) w ( u 1 v 2 u 2 v 1 ) w 0 ] ;   S = [ S S ] ,   S = [ ( u 0 v 1 u 1 v 0 ) w ( u 0 v 2 u 2 v 0 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 0 v 4 u 4 v 0 ) w 0 0 0 ( u 0 v 5 u 5 v 0 ) w 0 0 0 0 0 0 0 ] ,   S = [ 0 0 0 ( u 5 v 4 u 4 v 5 ) w 0 ( u 3 v 5 u 5 v 3 ) w ( u 5 v 2 u 2 v 5 ) w ( u 1 v 5 u 5 v 1 ) w 0 ( u 4 v 3 u 3 v 4 ) w ( u 2 v 4 u 4 v 2 ) w ( u 4 v 1 u 1 v 4 ) w ( u 2 v 3 u 3 v 2 ) w ( u 3 v 1 u 1 v 3 ) w ( u 1 v 2 u 2 v 1 ) w 0 ] ;     ,
it would be directly concluded that in (1 + 1) and (1 + 4), (1 + 5),... and higher space-time dimensions, the column matrices ϒ R and ϒ F could not be defined as the column matrices in unique formulations of the axiomatically derived general covariant field Equations (71) and (72). In other words, for the (1 + 2)-dimensional cases of the transformed column matrices ϒ R and ϒ F (89), the corresponding isomorphism (represented uniquely by the mappings (71-1) and (72-1)) could be defined between the components of these matrices and the entries of column matrix S (67), for s = 0 (compatible with ϕ ρ σ ( R ) = 0 , ϕ ( F ) = 0 ), if and only if: i R ρ σ 02 ( x , t ) = 0 and i F 02 ( x , t ) = 0 . This could be shown as follows:
( S = [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 1 v 2 u 2 v 1 ) w 0 ] D e r i v a t i o n   P r o c e d u r e ( F i r s t   Q u a n t i z a t i o n ) ϒ R = [ i R ρ σ 21 ( x , t ) 0 i R ρ σ 10 ( x , t ) i R ρ σ 02 ( x , t ) ] ,   ϒ F = [ i F 21 ( x , t ) 0 i F 10 ( x , t ) i F 02 ( x , t ) ] ) u 2 v 0 u 0 v 2 = 0 , i R ρ σ 02 ( x , t ) = 0   , i F 02 ( x , t ) = 0 .
where for appeared parametric condition: u 2 v 0 u 0 v 2 = 0 , as it would be shown in Remark 3 (in Section 3.12), it could be supposed solely: u 2 = v 2 , v 0 = u 0 , implying conditions: R ρ σ 02 = 0 and F 02 = 0 , which could be assumed for the field strength tensors R ρ σ μ ν and F μ ν in (1 + 2)-dimensional space-time (without vanishing these tensor fields), based on their basic definitions given by Formula (71-1) and (72-1).
Hence, definite mathematical framework of our axiomatic derivation approach (presented in Section 3.4), in addition to the time-reversal invariance (represented by the quantum operator (84)) of source-free cases of general covariant field Equations (71) and (72), imply the (1 + 2)-dimensional case of column matrices ΨR and ΨF given by relations (74) (where we assumed ϕ ρ σ ( R ) = 0 , ϕ ( F ) = 0 ), could be given solely as follows, to be compatible with the above assumed conditions (i.e., being compatible with the mathematical framework of axiomatic derivation of field Equations (71) and (72), and also the time-reversal invariance defined by quantum operator (84)), and consequently, as the column matrices could be defined in the formulations of the fundamental tensor field Equations (71) and (72), respectively:
Ψ R = [ R ρ σ 10 0 R ρ σ 21 0 ] , Ψ F = [ F 10 0 F 21 0 ]
The formulations (96) that are represented the column matrices ΨR and ΨF in the field Equations (71) and (72) compatible with the above basic conditions, are also represented these matrices in the field Equations (71) and (72) with non-zero source currents compatible with two basic conditions (similar to above conditions) including a unique combination of the C, P and T symmetries (that have been represented by quantum operators (83)–(85)) for these cases of field Equations (71) and (72), and also the mathematical framework of axiomatic derivation of Equations (71) and (72). In fact, as it has been shown in Remark 2 (in Section 3.11), the field Equations (71) and (72) with non-zero source currents could have solely a certain combination (given by Formula (86) and (87)) of the C, P and T symmetries (that are represented by the operators (83), (84) and (85)). This unique combined symmetry in addition to the unique formulations (96) of source-free cases of column matrices ΨR and ΨF in (1 + 2)-dimensional space-time, implies these matrices could take solely the following forms to be defined in the formulations of the fundamental tensor field Equations (71) and (72) (with non-zero source currents):
Ψ R = [ R ρ σ 10 0 R ρ σ 21 ϕ ρ σ ( R ) ] ,   Ψ F = [ F 10 0 F 21 ϕ ( F ) ] ,   J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F )
In the same manner, concerning the (1 + 4)-dimensional cases of column matrices ϒ R and ϒ F (91), there would be a mapping between the entries of these matrices and entries (with the same indices) of algebraic column matrix S (69), where s = 0 (compatible with ϕ ρ σ ( R ) = 0 , ϕ ( F ) = 0 ), if and only if: i R ρ σ 10 ( x , t ) = 0 , i R ρ σ 41 ( x , t ) = 0 , i R ρ σ 31 ( x , t ) = 0 , i R ρ σ 21 ( x , t ) = 0 , i F 10 ( x , t ) = 0 , i F 41 ( x , t ) = 0 , i F 31 ( x , t ) = 0 , i F 21 ( x , t ) = 0 , i.e.,:
( S = [ ( u 0 v 1 u 1 v 0 ) w ( u 2 v 0 u 0 v 2 ) w ( u 0 v 3 u 3 v 0 ) w 0 ( u 4 v 0 u 0 v 4 ) w 0 0 0 0 ( u 3 v 4 u 4 v 3 ) w ( u 2 v 4 u 4 v 2 ) w ( u 1 v 4 u 4 v 1 ) w ( u 2 v 3 u 3 v 2 ) w ( u 1 v 3 u 3 v 1 ) w ( u 1 v 2 u 2 v 1 ) w 0 ] D e r i v a t i o n   P r o c e d u r e ( F i r s t   Q u a n t i z a t i o n ) ϒ R = [ 0 i R ρ σ 43 ( x , t ) i R ρ σ 42 ( x , t ) i R ρ σ 41 ( x , t ) i R ρ σ 32 ( x , t ) i R ρ σ 31 ( x , t ) i R ρ σ 21 ( x , t ) 0 i R ρ σ 10 ( x , t ) i R ρ σ 02 ( x , t ) i R ρ σ 30 ( x , t ) 0 i R ρ σ 04 ( x , t ) 0 0 0 ] ,   ϒ F = [ 0 i F 43 ( x , t ) i F 42 ( x , t ) i F 41 ( x , t ) i F 32 ( x , t ) i F 31 ( x , t ) i F 21 ( x , t ) 0 i F 10 ( x , t ) i F 02 ( x , t ) i F 30 ( x , t ) 0 i F 04 ( x , t ) 0 0 0 ] ) u 1 = 0 ,   v 1 = 0 , i R ρ σ 10 ( x , t ) = 0 ,   i F 10 ( x , t ) = 0 , i R ρ σ 21 ( x , t ) = 0 ,   i F 21 ( x , t ) = 0 , i R ρ σ 31 ( x , t ) = 0 ,   i F 31 ( x , t ) = 0 , i R ρ σ 41 ( x , t ) = 0 ,   i F 41 ( x , t ) = 0 .
This means that in (1 + 4) space-time dimensions, the mathematical framework of our axiomatic derivation approach (described in Section 3.4) in addition to the time reversal invariance (defined by the quantum operator (84)) of the source-free case of the derived general covariant fundamental field Equations (71) and (72) imply the column matrices ΨR and ΨF (76) (for ϕ ρ σ ( R ) = 0 , ϕ ( F ) = 0 ) could take solely the following forms (in general) to be defined in the formulations of the field Equations (71) and (72):
Ψ R = [ 0 R ρ σ 02 R ρ σ 30 0 R ρ σ 04 0 0 0 0 R ρ σ 43 R ρ σ 42 0 R ρ σ 32 0 0 0 ] ,   Ψ F = [ 0 F 02 F 30 0 F 04 0 0 0 0 F 43 F 42 0 F 32 0 0 0 ]
which are equivalent to the (1 + 3)-dimensional source-free cases of column matrices ΨR and ΨF (represented uniquely by Formula (75)). In addition, similar to the formulations (96-1), as it has been shown in Remark 2 (in Section 3.11), the field Equations (71) and (72) with non-zero source currents have a certain (and unique) combination of the C, P and T symmetries (that have been defined by the operators (83), (84) and (85)). This combined symmetry in addition to the forms (98), imply also the (1 + 4)-dimensional cases of column matrices ΨR and ΨF represented by Formula (76) could take solely the following forms (in general) to be defined in the formulations of fundamental field Equations (71) and (72):
Ψ R = [ 0 R ρ σ 02 R ρ σ 30 0 R ρ σ 04 0 0 0 0 R ρ σ 43 R ρ σ 42 0 R ρ σ 32 0 0 ϕ ρ σ ( R ) ] ,   Ψ F = [ 0 F 02 F 30 0 F 04 0 0 0 0 F 43 F 42 0 F 32 0 0 ϕ ( F ) ] , J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F )
Consequently, the (1 + 4)-dimensional cases of column matrices ΨR and ΨF that are originally given by formulations (76), are reduced to Formula (98-1) which are equivalent to the (1 + 3)-dimensional cases of these matrices (given originally by column matrices of the forms (75)), i.e.,:
Ψ R = [ R ρ σ 10 R ρ σ 20 R ρ σ 30 0 R ρ σ 23 R ρ σ 31 R ρ σ 12 ϕ ρ σ ( R ) ] ,   Ψ F = [ F 10 F 20 F 30 0 F 23 F 31 F 12 ϕ ( F ) ] ,   J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( D ν + i m 0 ( F ) k ν ) ϕ ( F ) .
Moreover, as it would be also noted in Section 3.15, it is noteworthy to add that the tensor field R ρ σ μ ν in column matrix ΨR (98-2) (expressing the general representation of column matrices definable in the formulation of (1 + 3)-dimensional case of general covariant field Equation (71)), in fact, equivalently represents a massive bispinor field of spin-2 particles in (1 + 3) space-time dimensions (which could be identified as a definite generalized massive matrix formulation of the Einstein gravitational field, as it has been also shown in Section 3.6), and the tensor field F μ ν in the column matrix ΨF (98-2) (expressing the general representation of column matrices definable in the formulation of (1 + 3)-dimensional case of general covariant field Equation (72)), in fact, equivalently represents a massive bispinor field of spin-1 particles in (1 + 3) space-time dimensions (which could be identified as definite generalized massive formulation of the Maxwell electromagnetic field, as it has been also shown in Section 3.6 and Section 3.9; and also Yang-Mills fields compatible with specific gauge groups, as it would be shown in Section 3.15).
Summing up, in this Section (Section 3.12) we showed that the axiomatic approach of derivation of the field Equations (71) and (72) (described in Section 3.1, Section 3.3 and Section 3.4) in addition to their time reversal invariance (represented basically by the quantum operator (84)), imply these fundamentally derived equations could be solely defined in (1 + 2) and (1 + 3) space-time dimensions. “Hence, based on the later conclusion and also the basic assumption (3) (defined in Section 3.1), we may conclude directly that the universe could be realized solely with the (1 + 2) and (1 + 3)-dimensional space-times, and cannot have more than four space-time dimensions.
Based on the axiomatic arguments and relevant results presented and obtained in this Section, in the following Sections we consider solely the (1 + 2) and (1 + 3)-dimensional cases of general covariant field Equations (71) and (72) that are defined solely with the column matrices of the forms (96-1) and (98-2).
Equivalent (Asymptotically) Representations of the Bispinor Fields of Spin-3/2 and Spin-1/2 Particles, Respectively, by General Covariant Field Equations (71) and (72) (Formulated Solely with Column Matrices of the Types (96-1)) in (1 + 2) Space-Time Dimensions
It is noteworthy that according to the Ref. [52] and also based on the basic properties of the Riemann curvature tensor R ρ σ μ ν in (1 + 2) space-time dimensions [53] (in particular the identity: R ρ σ μ ν = ε ρ σ α ε μ ν β G α β , where G α β is the Einstein tensor and), it would be concluded that R ρ σ μ ν which is defined by (1 + 2)-dimensional case of the general covariant massive field Equation (71) (which could be defined solely with a column matrix of the type ΨR (96-1)), represents asymptotically a general covariant bispinor field of spin-3/2 particles (that would be asymptotically equivalent to the Rarita–Schwinger equation). In a similar manner, according to the Ref. [52], and also following the basic properties of field strength tensor F μ ν in (1 + 2) space-time dimensions (that as a rank two anti-symmetric with three independent components holding, in particular, the identities: F μ ν = ε μ ν α T α , T α = ( 1 / 2 ) ε α μ ν F μ ν , showing that F μ ν could be represented equivalently by a vector T α with three independent components as well) it would be concluded that F μ ν which is defined by (1 + 2)-dimensional case of the general covariant massive (tensor) field Equation (72) (which could be defined solely with a column matrix of the type ΨF (96-1)), represents asymptotically a general covariant bispinor field of spin-1/2 particles (that would be asymptotically equivalent to the Dirac equation [52]). Furthermore, as shown in Section 3.15, the general covariant field Equations (72) (representing asymptotically the spin-1/2 fermion fields) is also compatible with the SU(2)L⊗U(2)R symmetry group (representing “1 + 3” generations for both lepton and quark fields including a new charge-less fermion).

Remark 3. Equivalent Representations of the Bispinor Fields of Spin-2 and Spin-1 Particles, Respectively, by General Covariant Field Equations (71) and (72) (Defined Solely with Column Matrices of the Types (98-2)) in (1 + 3) Space-Time Dimensions

It should be also noted that according to the Refs. [5,6,7,51,54,55,56], the basic properties of the Riemann curvature tensor including the relevant results presented in Section 3.6, it would be concluded that the field strength tensor R ρ σ μ ν (i.e., the Riemann tensor) the in (1 + 3) space-time dimensions by general covariant massive (tensor) field Equation (71) (formulated solely with a column matrix of the type ΨR (98-2)), represents a general covariant bispinor field of spin-2 particles (as a generalized massive formulation of the Einstein gravitational field equation). In a similar manner, according to the Refs. [5,6,7,51,54,55,56], the field strength tensor F μ ν which is defined in (1 + 3) space-time dimensions by the general covariant massive (tensor) field Equation (72) (formulated solely by a column matrix of the type ΨF (98-2)), and represents a general covariant bispinor field of spin-1 particles (representing new generalized massive formulations of the Maxwell’s equations, and also Yang-Mills field equations). Furthermore, as it would be shown in Section 3.15, the general covariant field Equation (72) (representing the spin-1 boson fields coupling to the spin-1/2 fermionic currents) is also compatible with the SU(2)L⊗U(2)R and SU(3) symmetry groups.
Moreover, based on these determined gauge symmetries for the derived fermion and boson field equations, four new charge-less spin-1/2 fermions (represented by “ze , zn ; zu , zd”, where two right-handed charhe-less quarks zu and zd emerge specifically in two subgroups with anti-quarks such that: ( s ¯ ,   u ¯ ,   b ¯ , zu) and ( c ¯ ,   d ¯ ,   t ¯ , zd)), and also three new massive spin-1 bosons (represented by W ˜ + , W ˜ , Z , where in particular Z is the complementary right-handed particle of ordinary Z boson), are predicted uniquely by this new mathematical axiomatic approach.

3.13. Showing that Only a Definite Simultaneous Combination of the Quantum Mechanical Transformations C ^ , P ^ , T ^ and γ C h . (Given Uniquely by the Matrix Operators (83)–(87)) Could Be Defined for the General Covariant Massive (Tensor) Field Equations (71) and (72) with Non-Zero Source Currents

As it has been shown in Remark 1 (in Section 3.11) and Section 3.12, since the algebraic column matrix S in the matrix Equation (64) (derived and represented uniquely in terms of the matrices (66)–(70), … corresponding to various space-time dimensions), is not symmetric by multiplying by matrix γ C h (84-1) and (84-2) (except for (1 + 2) and (1 + 3)-dimensional cases of column matrix S , based on the definite algebraic properties of matrix S presented in Remark 1 in Section 3.3), it is concluded that except the (1 + 2) and (1 + 3)-dimensional cases of the fundamental field Equations (71) and (72), these field equations could not be defined with column matrices of the type Ψ R ( C h ) ( = γ C h Ψ R ) and Ψ F ( C h ) ( = γ C h Ψ F ) (if assuming that the column matrices Ψ R and Ψ F are defined with field Equations (71) and (72), i.e., they have the formulations similar to the formulations of originally derived column matrices (73)–(77), … corresponding to various space-time dimensions). This conclusion follows from this fact that the filed Equations (71) and (72) have been derived (and defined) uniquely from the matrix Equation (64) via the axiomatic derivation approach (including the first quantization procedure) presented in Section 3.4 and Section 3.5. As it has been shown in Section 3.12, using this property (i.e., multiplication of column matrices Ψ R and Ψ F , defined in the unique expressions of fundamental field Equations (71) and (72), by matrix γ C h from the left), this crucial and essential issue is concluded directly that by assuming the time-reversal invariance of the general covariant filed Equations (71) and (72) (represented by the transformations T ^ Ψ R and T ^ Ψ F , where the quantum operator T ^ is given uniquely by Formula (84), i.e.,: T ^ = T ^ 0 K ^ = i γ P γ C h K ^ ), these fundamental field equations could be defined solely in (1 + 2) and (1 + 3) space-time dimensions (with the column matrices of the forms (96-1) and (98-2), respectively).
Hence, the definite mathematical formalism of the axiomatic approach of derivation of fundamental field Equations (71) and (72), along with the C, P and T symmetries (represented by the quantum matrix operators (83)–(87), in Section 3.11) of source-free cases (as basic cases) of these equations, in fact, imply these equations with non-zero source currents, would be invariant solely under the simultaneous combination of all the transformations C ^ , P ^ , and T ^ (83)–(85), multiplied by matrix γ C h (defined by Formula (84-1) and (84-2)). This unique combined transformation could be expressed uniquely as follows, respectively, for the particle fields (represented by column matrices Ψ R ( r , t ) , Ψ F ( r , t ) ) and their corresponding antiparticle fields (represented by column matrices Ψ R ( r , t ) , Ψ F ( r , t ) given solely with reversed signs of the temporal and spatial coordinates):
Z ^ C O M B Ψ R ( r , t ) = γ C h T ^ P ^ C ^   Ψ R ( r , t ) , Z ^ C O M B Ψ F ( r , t ) = γ C h T ^ P ^ C ^   Ψ F ( r , t ) ;
Z ˜ ^ C O M B Ψ R ( r , t ) = γ C h T ^ P ^ C ^   Ψ R ( r , t ) , Z ˜ ^ C O M B Ψ F ( r , t ) = γ C h T ^ P ^ C ^   Ψ F ( r , t ) .
The unique combined form of transformation Z ^ C O M B (99) (and also Z ˜ ^ C O M B (100), where Z ˜ ^ C O M B = Z ^ C O M B ) is based on the following two basic issues:
Firstly, it follows from the definite formulations of uniquely determined column matrices (73)–(77), … (corresponding to various space-time dimensions, however, as noted above, based on the arguments presented in Section 3.12, the only definable column matrices in the formulations of field Equations (71) and (72), are of the types ΨR and ΨF represented by Formula (96-1) and (98-2), in (1 + 2) and (1 + 3) space-time dimensions, respectively), where the source currents J ρ σ ν ( R ) and J ν ( F ) should be expressible by these conditional relations (in terms of the arbitrary covariant quantities ϕ ρ σ ( R ) and ϕ ( F ) ), respectively: J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( ν + i m 0 ( F ) k ν ) ϕ ( F ) . In other words, the unique formulation of derived combined symmetries Z C O M B and Z ˜ C O M B represented by the quantum operators Z ^ C O M B (99) and Z ˜ ^ C O M B (100), in particular, is a direct consequent of the above conditional expressions for source currents J ρ σ ν ( R ) and J ν ( F ) . As noted in Section 3.5, these relations appear as necessary conditions in the course of the axiomatic derivation of general covariant field Equations (71) and (72). In fact, in the field Equations (71) and (72) the uniquely derived column matrices ΨR and ΨF (73)–(77),…, not only contain all the components of tensor fields R ρ σ μ ν and F μ ν , but also contain the components of arbitrary covariant quantities ϕ ρ σ ( R ) and ϕ ( F ) (as the initially given quantities) which define the source currents J ρ σ ν ( R ) and J ν ( F ) by the above expressions, respectively, i.e.,: J ρ σ ν ( R ) = ( ν + i m 0 ( R ) k ν ) ϕ ρ σ ( R ) , J ν ( F ) = ( ν + i m 0 ( F ) k ν ) ϕ ( F ) . Now based on these conditional expressions in addition to this natural and basic circumstance that the source currents J ρ σ ν ( R ) and J ν ( F ) should be also transferred respectively as a rank three tensor and a vector, under the parity, time-reversal and charge conjugation transformations (defined by Formula (83)–(85)) of the field Equations (71) and (72), it would be concluded directly that the transformations (99) and (100) are the only simultaneous combinations of transformations C ^ , P ^ , T ^ (also including the matrix γ C h , necessarily, as it would be shown in the following paragraph), which could be defined for the field Equations (71) and (72) with “non-zero” source currents.
Secondly, appearing the matrix operator γ C h in simultaneous combinations γ C h T ^ P ^ C ^ and γ C h T ^ P ^ C ^ in the combined transformations (99) and (100), follows simply from the basic arguments presented in Section 3.12. In fact, in these uniquely determined combinations, the simultaneous multiplication by matrix γ C h (from the left) is a necessary condition for that the transformed column matrices: Z ^ C O M B Ψ R ( r , t ) , Z ˜ ^ C O M B Ψ R ( r , t ) , Z ^ C O M B Ψ F ( r , t ) , Z ˜ ^ C O M B Ψ F ( r , t ) (given bythe transformations (99) and (100)) could be also defined in the field Equations (71) and (72), based on the formulations of column matrices (96-1) and (98-2), as mentioned in Section 3.12 (however, it is worth notingthat this argument is not merely limited to the definability of column matrices of the types (96-1) and (98-2), and it could be also represented on the basis of unique formulations of all the originally derived column matrices (73)–(77), … corresponding to various space-time dimensions).
In the next Section, we show how the ‘CPT’ theorem in addition to the unique formulations of the combined transformations (99) and (100) (representing the only definable transformation forms, including C, P and T quantum mechanical transformations, for the field Equations (71) and (72) with “non-zero” source currents), imply only the left-handed particle fields (along with their complementary right-handed fields) could be coupled to the corresponding (any) source currents.

3.14. Showing that Only the Left-Handed Particle (along with Their Complementary Right-Handed Antiparticle) Fields Could Be Coupled to the Corresponding Source Currents

On the basis of the ‘CPT’ theorem [6,7], it would be concluded directly that the unique combined forms of transformations (99) and (100) (representing the only combination of C ^ , P ^ , and T ^ transformations multiplied by matrix γ C h , that could be defined as a symmetry for general covariant field Equations (71) and (72) with non-zero source currents), should be equivalent only to simultaneous combination of C ^ , P ^ , and T ^ transformations (that have been defined uniquely by Formula (83)–(85)). Moreover, based on the ‘CPT’ theorem, the simultaneous combination of transformations C ^ , P ^ , and T ^ should: “interchange the particle field and its corresponding antiparticle field; inverts the spatial coordinates r r ; reverse the spin of all particle fields; leave the direction of the momentum invariant; and, therefore, should interchange the left-handed and right-handed components of both particle field and its corresponding antiparticle field”. Hence, we should have:
Z ^ C O M B [ Ψ R ( r ) ] ( L e f t ) = γ C h T ^ P ^ C ^ [ Ψ R ( r ) ] ( L e f t ) = γ C h [ Ψ R ( r ) ] ( R i g h t ) = [ Ψ R ( r ) ] ( R i g h t ) , Z ^ C O M B [ Ψ F ( r ) ] ( L e f t ) = γ C h T ^ P ^ C ^ [ Ψ F ( r ) ] ( L e f t ) = γ C h [ Ψ F ( r ) ] ( R i g h t ) = [ Ψ F ( r ) ] ( R i g h t ) ;
Z ^ C O M B [ Ψ R ( r ) ] ( R i g h t ) = γ C h T ^ P ^ C ^ [ Ψ R ( r ) ] ( R i g h t ) = γ C h [ Ψ R ( r ) ] ( L e f t ) = [ Ψ R ( r ) ] ( L e f t ) , Z ^ C O M B [ Ψ F ( r ) ] ( R i g h t ) = γ C h T ^ P ^ C ^ [ Ψ F ( r ) ] ( R i g h t ) = γ C h [ Ψ F ( r ) ] ( L e f t ) = [ Ψ F ( r ) ] ( L e f t ) ;
Z ˜ ^ C O M B [ Ψ R ( r ) ] ( L e f t ) = γ C h T ^ P ^ C ^ [ Ψ R ( r ) ] ( L e f t ) = γ C h [ Ψ R ( r ) ] ( R i g h t ) = [ Ψ R ( r ) ] ( R i g h t ) , Z ˜ ^ C O M B [ Ψ F ( r ) ] ( L e f t ) = γ C h T ^ P ^ C ^ [ Ψ F ( r ) ] ( L e f t ) = γ C h [ Ψ F ( r ) ] ( R i g h t ) = [ Ψ F ( r ) ] ( R i g h t ) ;
Z ˜ ^ C O M B [ Ψ R ( r ) ] ( R i g h t ) = γ C h T ^ P ^ C ^ [ Ψ R ( r ) ] ( R i g h t ) = γ C h [ Ψ R ( r ) ] ( L e f t ) = [ Ψ R ( r ) ] ( L e f t ) , Z ˜ ^ C O M B [ Ψ F ( r ) ] ( R i g h t ) = γ C h T ^ P ^ C ^ [ Ψ F ( r ) ] ( R i g h t ) = γ C h [ Ψ F ( r ) ] ( L e f t ) = [ Ψ F ( r ) ] ( L e f t ) .
where the column matrices Ψ R ( r ) and Ψ F ( r ) represent the particle field and Ψ R ( r ) and Ψ F ( r ) denote the transformed forms of column matrices of Ψ R ( r ) and Ψ F ( r ) , respectively, under the simultaneous combination of transformations C ^ , P ^ , and T ^ (83)–(85). Furthermore, in agreement and based on the definitions and properties of quantum operators C ^ , P ^ , T ^ and matrix γ C h given by Formula (83)–(87), the left-handed and right-handed components of column matrices of the types (96-1) and (98-2) (representing the unique formulations of column matrices that could be defined in the field Equations (71) and (72), as mentioned in Section 3.12) are defined solely as follows for the column matrices Ψ R ( r ) , Ψ F ( r ) and also Ψ R ( r ) , Ψ F ( r ) (as the transformed forms of column matrices Ψ R ( r ) and Ψ F ( r ) under the C ^ P ^ T ^ transformation, respectively):
[ Ψ R ( r ) ] ( L e f t ) = 1 2 [ Ψ R ( r ) + γ C h Ψ R ( r ) ] , [ Ψ F ( r ) ] ( L e f t ) = 1 2 [ Ψ F ( r ) + γ C h Ψ F ( r ) ] , [ Ψ R ( r ) ] ( R i g h t ) = 1 2 [ Ψ R ( r ) γ C h Ψ R ( r ) ] , [ Ψ F ( r ) ] ( R i g h t ) = 1 2 [ Ψ F ( r ) γ C h Ψ F ( r ) ] ;
[ Ψ R ( r ) ] ( L e f t ) = 1 2 [ Ψ R ( r ) + γ C h Ψ R ( r ) ] , [ Ψ F ( r ) ] ( L e f t ) = 1 2 [ Ψ F ( r ) + γ C h Ψ F ( r ) ] , [ Ψ R ( r ) ] ( R i g h t ) = 1 2 [ Ψ R ( r ) γ C h Ψ R ( r ) ] , [ Ψ F ( r ) ] ( R i g h t ) = 1 2 [ Ψ F ( r ) γ C h Ψ F ( r ) ] .
where we have:
Ψ R ( r ) = [ Ψ R ( r ) ] ( L e f t ) + [ Ψ R ( r ) ] ( R i g h t ) , Ψ F ( r ) = [ Ψ F ( r ) ] ( L e f t ) + [ Ψ F ( r ) ] ( R i g h t ) , Ψ R ( r ) = [ Ψ R ( r ) ] ( L e f t ) + [ Ψ R ( r ) ] ( R i g h t ) , Ψ F ( r ) = [ Ψ F ( r ) ] ( L e f t ) + [ Ψ F ( r ) ] ( R i g h t ) .
Now using the definitions (105) and (106) in the Formula (101)–(104), we obtain:
γ C h [ Ψ R ( r ) ] ( L e f t ) = 1 2 γ C h [ Ψ R ( r ) + γ C h Ψ R ( r ) ] = 1 2 [ γ C h Ψ R ( r ) + Ψ R ( r ) ] = [ Ψ R ( r ) ] ( L e f t ) , γ C h [ Ψ F ( r ) ] ( L e f t ) = 1 2 γ C h [ Ψ F ( r ) + γ C h Ψ F ( r ) ] = 1 2 [ γ C h Ψ F ( r ) + Ψ F ( r ) ] = [ Ψ F ( r ) ] ( L e f t ) ;
γ C h [ Ψ R ( r ) ] ( R i g h t ) = 1 2 γ C h [ Ψ R ( r ) γ C h Ψ R ( r ) ] = 1 2 [ γ C h Ψ R ( r ) Ψ R ( r ) ] [ Ψ R ( r ) ] ( R i g h t ) , γ C h [ Ψ F ( r ) ] ( R i g h t ) = 1 2 γ C h [ Ψ F ( r ) γ C h Ψ F ( r ) ] = 1 2 [ γ C h Ψ F ( r ) Ψ F ( r ) ] [ Ψ F ( r ) ] ( R i g h t ) ;
γ C h [ Ψ R ( r ) ] ( L e f t ) = 1 2 γ C h [ Ψ R ( r ) + γ C h Ψ R ( r ) ] = 1 2 [ γ C h Ψ R ( r ) Ψ R ( r ) ] [ Ψ R ( r ) ] ( L e f t ) , γ C h [ Ψ F ( r ) ] ( L e f t ) = 1 2 γ C h [ Ψ F ( r ) + γ C h Ψ F ( r ) ] = 1 2 [ γ C h Ψ F ( r ) Ψ F ( r ) ] [ Ψ F ( r ) ] ( L e f t ) ;
γ C h [ Ψ R ( r ) ] ( R i g h t ) = 1 2 γ C h [ Ψ R ( r ) γ C h Ψ R ( r ) ] = 1 2 [ γ C h Ψ R ( r ) + Ψ R ( r ) ] = [ Ψ R ( r ) ] ( R i g h t ) , γ C h [ Ψ F ( r ) ] ( R i g h t ) = 1 2 γ C h [ Ψ F ( r ) γ C h Ψ F ( r ) ] = 1 2 [ γ C h Ψ F ( r ) + Ψ F ( r ) ] = [ Ψ F ( r ) ] ( R i g h t ) .
Based on the relations (101-1)–(104-1), it would be concluded directly that only the left-handed components of particle fields represented by [ Ψ R ( r ) ] ( L e f t ) , [ Ψ F ( r ) ] ( L e f t ) , and the right-handed components of their corresponding antiparticle fields represented by [ Ψ R ( r ) ] ( R i g h t ) , [ Ψ F ( r ) ] ( R i g h t ) (as the transformed forms of column matrices [ Ψ R ( r ) ] ( L e f t ) and [ Ψ F ( r ) ] ( L e f t ) under the C ^   P ^   T ^ transformation, respectively), obey the transformations (101-1) and (104-1) (as the necessary conditions given respectively by relations (101) and (104)). On the other hand, the right-handed components of particle fields represented by [ Ψ R ( r ) ] ( R i g h t ) , [ Ψ F ( r ) ] ( R i g h t ) , and the left-handed components of their corresponding antiparticle fields represented by [ Ψ R ( r ) ] ( L e f t ) , [ Ψ F ( r ) ] ( L e f t ) (as the transformed forms of column matrices [ Ψ R ( r ) ] ( R i g h t ) and [ Ψ F ( r ) ] ( R i g h t ) under the C ^   P ^   T ^ transformation, respectively), obey the transformations (102-1) and (103-1) (as the necessary conditions given respectively by relations (102) and (103)). Hence (and also following the basic assumption (3) defined in Section 3.1), it is concluded directly that only the left-handed particle fields (along with their complementary right-handed fields) could be coupled to the corresponding (any) source currents. This means that only left-handed bosonic fields (along with their complementary right-handed fields) could be coupled to the corresponding fermionic source currents; which also means that only left-handed fermions (along with their complementary right-handed fermions) can participate in any interaction with the bosons (which consequently would be only left-handed bosons or their complementary right-handed bosons).

3.15. Showing the Gauge Invariance of Axiomatically Derived General Covariant (Tensor) Field Equation (72) in (1 + 2)-dimensional Space-Time (Definable with Column Matrices of the Type ΨF (96-1), Representing the spin-1/2 Fermion Fields) under the SU(2)L⊗U(2)R Symmetry Group, and also Invariance of this Equation in (1 + 3)-dimensional Space-Time (Definable Column Matrices of the Type ΨF (98-2), Representing the Spin-1 Boson Fields Coupled to the Fermionic Source Currents) under the SU(2)L⊗U(2)R and SU(3) Symmetry Group

One of the natural and basic properties of the (1 + 2)-dimensional space-time geometry is that the metric tensor can be “diagonalized” [57]. Using this basic property, the invariant energy-momentum quadratic relation (52) (in Section 3.1) would be expressed as follows:
g 00 ( p 0 ) 2 g 00 ( p 0 s t ) 2 + g 11 ( p 1 ) 2 + g 22 ( p 2 ) 2 = 0
that is equivalent to: g 00 ( p 0 ) 2 + g 11 ( p 1 ) 2 + g 22 ( p 2 ) 2 = ( m 0 c ) 2 , where (as defined in Section 3.1) m 0 and p μ are the particle’s rest mass and momentum (3-momentum), p μ s t = m 0 k μ , and k μ = ( k 0 , 0 , 0 ) = ( c / g 00 , 0 , 0 ) denotes the covariant form of the 3-velocity of particle in stationary reference frame. As it would be shown in the following, a crucial and essential property of the quadratic relation (108) is its invariance under a certain set of sign inversions of the components of particle’s momentum: ( p 0 , p 1 , p 2 ) , along with similar inversions for the components: ( p 0 s t , p 1 s t , p 2 s t ) , where p 0 s t = m 0 k 0 , p 1 s t = p 2 s t = 0 . This set includes seven different types of the sign inversions (in total), which could be represented simply by the following symmetric group of transformations (based on the formalism of the corresponding Lorentz symmetry group of invariant relation (108)), respectively:
( p 0 , p 0 s t , p 1 , p 2 ) ( p 0 , p 0 s t , p 1 , p 2 ) = ( p 0 ( 1 ) , p 0 s t ( 1 ) , p 1 ( 1 ) , p 2 ( 1 ) )
( p 0 , p 0 s t , p 1 , p 2 ) ( p 0 , p 0 s t , p 1 , p 2 ) = ( p 0 ( 2 ) , p 0 s t ( 2 ) , p 1 ( 2 ) , p 2 ( 2 ) )
( p 0 , p 0 s t , p 1 , p 2 ) ( p 0 , p 0 s t , p 1 , p 2 ) = ( p 0 ( 3 ) , p 0 s t ( 3 ) , p 1 ( 3 ) , p 2 ( 3 ) )
( p 0 , p 0 s t , p 1 , p 2 ) ( p 0 , p 0 s t , p 1 , p 2 ) = ( p 0 ( 4 ) , p 0 s t ( 4 ) , p 1 ( 4 ) , p 2 ( 4 ) )
( p 0 , p 0 s t , p 1 , p 2 ) ( p 0 , p 0 s t , p 1 , p 2 ) = ( p 0 ( 5 ) , p 0 s t ( 5 ) , p 1 ( 5 ) , p 2 ( 5 ) )
( p 0 , p 0 s t , p 1 , p 2 ) ( p 0 , p 0 s t , p 1 , p 2 ) = ( p 0 ( 6 ) , p 0 s t ( 6 ) , p 1 ( 6 ) , p 2 ( 6 ) )
( p 0 , p 0 s t , p 1 , p 2 ) ( p 0 , p 0 s t , p 1 , p 2 ) = ( p 0 ( 7 ) , p 0 s t ( 7 ) , p 1 ( 7 ) , p 2 ( 7 ) )
Moreover, although, following noncomplex- algebraic values of momentum’s components p μ ( p μ = p μ ) , the corresponding complex representations of transformations (108-1)–(108-7) is not a necessary issue in general, however, if the invariant relation (108) is represented formally by equivalent complex form:
g 00 ( p 0 p 0 ) g 00 ( p 0 s t p 0 s t ) + g 11 ( p 1 p 1 ) + g 22 ( p 2 p 2 ) = 0
then, along with the set seven real-valued transformations (108-1)–(108-7), this relation would be also invariant under these corresponding sets of complex transformations (for a = 1 , 2 , 3 , , 7 ):
( p 0 , p 0 s t , p 1 , p 2 ) ( ± i p 0 ( a ) , ± i p 0 s t ( a ) , ± i p 1 ( a ) , ± i p 2 ( a ) ) , ( p 0 , p 0 s t , p 1 , p 2 ) ( i p 0 ( a ) , i p 0 s t ( a ) , i p 1 ( a ) , i p 2 ( a ) )
As the next step, in the following, using the transformations (108-1)–(108-7) (along with their corresponding complex forms (108-9)), a certain set of seven simultaneous (different) general covariant field Equations (corresponding to a group of seven bispinor fields of spin-1/2 particles) would be determined as particular cases of the (1 + 2)-dimensional form of general covariant field Equation (72) (defined with a column matrix of the type (96-1)).
Based on the definite formulation of (1 + 2)-dimensional case of system of linear Equation (64) (formulated in terms of the matrices (67)), for the energy-momentum relation (108) (along with the transformations (108-1)–(108-7)), the following set of seven systems of linear Equations (with different parametric formalisms) is determined uniquely. The general parametric solution of each of these systems of linear equations, obeys also the quadratic relation (108) (representing a set of seven forms, with different parametric formulations, of the general parametric solutions of quadratic relation (108)). This set of the seven systems of linear equations could be represented uniformly by a matrix equation as follows:
( α μ p μ ( a ) m 0 ( a ) α ˜ μ k μ ) S ( a ) = 0
where a = 1 , 2 , 3 , , 7 , p μ s t ( a ) = m 0 ( a ) k μ , α μ and α ˜ μ are two contravariant 4 × 4 real matrices (compatible with matrix representations of the Clifford algebra Cℓ1,2) defined solely by Formula (65) and (67), and parametric column matrix S ( a ) is also given uniquely as follows (formulated on the basis of definite parametric formulation of column matrix S (67) in (1 + 2) space-time dimensions):
S ( a ) = [ ( u 0 ( a ) v 1 ( a ) u 1 ( a ) v 0 ( a ) ) w ( u 2 ( a ) v 0 ( a ) u 0 ( a ) v 2 ( a ) ) w ( u 1 ( a ) v 2 ( a ) u 2 ( a ) v 1 ( a ) ) w s ]
which includes seven cases with specific parametric formulations expressed respectively in terms of seven groups of independent arbitrary parameters: u 0 ( a ) , u 1 ( a ) , u 2 ( a ) , v 0 ( a ) , v 1 ( a ) , v 2 ( a ) , and two common arbitrary parameters s and w (i.e., having the same forms in all of the seven cases of column matrix S ( a ) ). In addition, concerning the specific parametric expression (109-1) of column matrix S ( a ) in the formulation of matrix Equation (109), it is necessary to add that this parametric expression has been determined specifically by assuming (as a basic assumption in addition to the systematic natural approach of formulating the matrix Equation (109), based on the definite formulation of axiomatically determined matrix Equation (64)) the minimum value for total number of the arbitrary parameters in all of seven cases of column matrix S ( a ) , which implies equivalently the minimum value for total number of the arbitrary parameters in all of seven simultaneous (different) cases of matrix Equation (109) (necessarily with seven independent parametric solutions representing a certain set of seven different equivalent forms of the general parametric solution of quadratic relation (108), based on the general conditions of basic definition of the systems of linear equations corresponding to homogeneous quadratic and higher degree equations, presented in Section 2, and Section 2.2, Section 2.3, Section 2.4 and Section 3.1 concerning the homogenous quadratic equations).
In the following, in the derivation of the corresponding field Equations (from matrix Equation (109)), we will also use the above particular algebraic property of parameters s which has been expressed commonly in the expressions of all of seven simultaneous cases of matrix Equation (109) (and as shown below, concerning the (1 + 3)-dimensional corresponding form of matrix Equation (109), which holds the similar property).
In addition, along with the transformations (108-1)–(108-7) and algebraic matrix Equation (109), using the corresponding complex transformations (108-9), we may also formally have the following equivalent matrix Equation (with the complex expression):
( i α μ p μ ( a ) i m 0 ( a ) α ˜ μ k μ ) S ( a ) = 0
where a = 1 , 2 , 3 , , 7 , and p μ s t ( a ) = m 0 ( a ) k μ . Although, based on the real value of momentum p μ ( α ) ( p μ ( a ) = p μ ( a ) ) , the complex expression of each of the seven cases of algebraic matrix Equation (109), definitely, is not a necessary issue at the present stage. However, since the corresponding momentum operator p ^ μ ( a ) has a complex value (where p ^ μ ( a ) p ^ μ ( a ) ), in the following, using this basic property of the momentum operator, we derive a certain set of seven different simultaneous general covariant field equations from the matrix Equations (109) and (109-2) (based on the general axiomatic approach of derivation of general covariant massive field Equation (72), presented in Section 3.4, Section 3.5, Section 3.6, Section 3.7, Section 3.8, Section 3.9 and Section 3.10 in addition to certain forms of quantum representations of the C, P and T symmetries of this field equation, presented in Section 3.11, Section 3.12, Section 3.13 and Section 3.14). Furthermore, in the following, it would be also shown that the uniform representation of this determined set of seven simultaneous field equations, describe a certain group of seven simultaneous bispinor fields of spin-1/2 particles (corresponding, respectively, to a new right-handed charge-less fermion in addition to three right-handed anti-fermions, along with their three complementary left-handed fermions).
Furthermore, concerning the gravitational field Equation (71), it should be noted that following from the fact that the general covariant field Equation (71) should describe, uniquely and uniformly, the background space-time geometry via a certain form of the Riemann curvature tensor (which should be determined from the tensor field Equation (71)), the matrix Equation (109) could not be used for the derivation of a set of simultaneous different spin-3/2 fermion fields in (1 + 2) dimensions (there would be the same condition for the field Equation (71) in higher-dimensional space-times).
Hence, based on the axiomatic approach of derivation of (1 + 2)-dimensional case of field Equation (72) (defined solely by a column matrix of the form (96-1) in (1 + 2) space-time dimensions, as shown in Section 3.12), from the matrix Equation (109) and (109-2) (defined solely by column matrix (109-1)), and also taking into account the momentum operator’s property: p ^ μ p ^ μ , the following group of seven simultaneous (different) general covariant field equations could be determined:
( i α μ D μ ( f ) m 0 ( f ) α ˜ μ k μ ) Ψ F ( f ) = 0
specifying by the following group of transformations (based on the corresponding group of transformations (108-1)–(108-7) and (108-9)), for f = 1 , 2 , 3 , , 7 , respectively
( D 0 ( 1 ) , m 0 ( 1 ) , D 1 ( 1 ) , D 2 ( 1 ) ) = ( D 0 , m 0 , D 1 , D 2 )
( D 0 ( 2 ) , m 0 ( 2 ) , D 1 ( 2 ) , D 2 ( 2 ) ) = ( D 0 , m 0 , D 1 , D 2 )
( D 0 ( 3 ) , m 0 ( 3 ) , D 1 ( 3 ) , D 2 ( 3 ) ) = ( D 0 , m 0 , D 1 , D 2 )
( D 0 ( 4 ) , m 0 ( 4 ) , D 1 ( 4 ) , D 2 ( 4 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 )
( D 0 ( 5 ) , m 0 ( 5 ) , D 1 ( 5 ) , D 2 ( 5 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 )
( D 0 ( 6 ) , m 0 ( 6 ) , D 1 ( 6 ) , D 2 ( 6 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 )
( D 0 ( 7 ) , m 0 ( 7 ) , D 1 ( 7 ) , D 2 ( 7 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 )
where the column matrix Ψ F ( f ) would be also given as follows (based on the definite formulation of column matrix Ψ F (96-1) in Section 3.12, expressing the general representation of column matrices definable in the formulation of (1 + 2)-dimensional case of general covariant field Equation (72)):
Ψ F ( f ) = [ F 10 ( f ) 0 F 21 ( f ) ϕ F ] ,   J ν ( f ) = ( D ν ( f ) + i m 0 ( f ) k ν ) ϕ F
where in all of the seven simultaneous cases of field Equation (110) defined respectively by the column matrices Ψ F ( f ) (110-8) (for f = 1 , 2 , 3 , , 7 ), the scalar quantity ϕ F (that as a given initial quantity, defines the source currents J ν ( f ) (110-8)), necessarily, has the same value, based on the definite parametric formulation of the algebraic column matrix (109-1) (in particular, the common form of the corresponding arbitrary parameter s in the expressions of all of the seven simultaneous cases of matrix Equation (109)).
Following the definite formulations of set of seven general covariant (massive) field Equations (110) (specified, respectively, by the group of seven transformations (110-1)–(110-7)), the set of these could be represented uniformly by the following general covariant field equation as well (defined solely in (1 + 2) space-time dimensions):
( i α μ D μ m 0 α ˜ μ k μ ) Ψ F = 0
where the column matrix Ψ F given by:
Ψ F = [ F 10 0 F 21 ϕ F ] ,   J μ = ( D μ + i m 0 k μ ) ϕ F
and the field strength tensor F μ ν , scalar ϕ F , along with the source current J μ are defined as follows:
F μ ν = f = 1 7 F μ ν ( f ) τ f , I 2 ( D μ + i m 0 k μ ) ϕ F = f = 1 7 ( D μ ( f ) + i m 0 ( f ) k μ ) ϕ F = f = 1 7 J μ ( f ) τ f = J μ .
where I is the 2 × 2 identity matrix, and τ f = l f 2 (for f = 1 , 2 , 3 , , 7 ) are a set of seven 2 × 2 complex matrices given by,
τ f = l f 2 , { l 1 , l 2 , l 3 , l 4 , l 5 , l 6 , l 7 } = { [ 0 i i 0 ] , [ 0 1 1 0 ] , [ i 0 0 i ] , [ 1 0 0 1 ] , [ 0 1 1 0 ] , [ 0 i i 0 ] , [ 1 0 0 1 ] }
which as would be shown in the following, represents uniformly a combined gauge symmetry group of the form: SU(2)L⊗U(2)R , where the sub-set of three matrices “ τ 1 , τ 2 , τ 3 ” corresponds to SU(2)L group, and subset of four matrices “ τ 4 , τ 5 , τ 6 , τ 7 ” corresponds to U(2)R group.
Now based on the matrix formulation of field strength tensor F μ ν (defined by the general covariant field Equation (110-9)), and on the basis of C, P and T symmetries of this field Equation (as a particular form of the (1 + 2)-dimensional case of field Equation (72)) that have been represented basically by their corresponding quantum operators (in Section 3.5, Section 3.6, Section 3.7, Section 3.8, Section 3.9, Section 3.10, Section 3.11, Section 3.12, Section 3.13 and Section 3.14), it would be concluded that the general covariant field Equation (110-9) describes uniformly a group of seven spin-1/2 fermion fields corresponding to, respectively: “three left-handed fermions (for f = 1 , 2 , 3 ), in addition to their three complementary right-handed anti-fermions (for f = 5 , 6 , 7 ), and also a new single charge-less right-handed spin-1/2 fermion (for f = 4 ).” Hence, following the basic algebraic properties of seven matrices τ f (110-12), and the gauge symmetry group of the type: SU(2)L⊗U(2)R generated by these matrices, the three matrices τ 1 , τ 2 , τ 3 (corresponding with SU(2)L) represent respectively “three left-handed fermions”, and four matrices τ 4 , τ 5 , τ 6 , τ 7 (corresponding with U(2)R) represent respectively: “a new single right-handed charge-less spin-1/2 fermion, and three right-handed spin-1/2 fermions as the complementary particles of the three left-handed spin-1/2 fermions represented by matrices τ 1 , τ 2 , τ 3 ”.
Furthermore, as noted previously and shown below, as a natural assumption, by assuming the seven types of spin-1/2 fermion fields that are described by general covariant field Equation (110-9), as the source currents of spin-1 boson fields (that will be represented by two determined unique groups describing respectively by general covariant field Equations (114-4) and (114-5)), it would be concluded that there should be, in total, four specific groups of seven spin-1/2 fermion fields (each) with certain properties, corresponding to “1 + 3” generations of four fermions, including two groups of four leptons each, and two groups of four quarks each. Moreover, based on this basic circumstances, two groups of leptons would be represented uniquely by: “[( ν μ , e, ν τ ), ( ν ¯ μ , e+, ν ¯ τ ze)] and [(μ, ν e , τ), (μ+, ν ¯ e , τ+, zn)], respectively, where each group includes a new single right-handed charge-less lepton, represented by: ze and zn”; and two groups of quarks would be also represented uniquely by: “[( s , u , b ), ( s ¯ ,   u ¯ ,   b ¯ , zu)] and [( c , d , t ), ( c ¯ ,   d ¯ ,   t ¯ , zd)], respectively, where similar to leptons, each group includes a new single right-handed charge-less quark, represented by: zu and zd”. In addition, emerging two right-handed charhe-less quarks zu and zd specifically in two subgroups with anti-quarks ( s ¯ ,   u ¯ ,   b ¯ , zu) and ( c ¯ ,   d ¯ ,   t ¯ , zd), could explain the baryon asymmetry, and subsequently, the asymmetry between matter and antimatter in the universe.
Assuming the spin-1/2 fermion fields describing by general covariant massive field Equations (110-9) (defined by column matrix (110-10) in (1 + 2) space-time dimensions with a digonalized metric), as the coupling source currents of spin-1 boson fields (describing generally by (1 + 3)-dimensional case of general covariant field Equation (72) formulated with a column matrix of the type ΨF (98-2)), it is concluded that the (1 + 3)-dimensional metric could be also diagonalized for corresponding spin-1 boson fields. This conclusion follows directly from the above assumption that the (1 + 3)-dimensional metric of spin-1 boson fields (coupled to the corresponding fermionic source currents) would be also partially diagonalized such that for μ, ν = 0, 1, 2 and μν: g μ ν = 0 , which subsequently impliy g 03 = g 13 = g 23 = 0 . Hence, the invariant energy-momentum relation (52) will be expressed as follows in (1 + 3)-dimensional space-time with diagonalized metric:
g 00 ( p 0 ) 2 g 00 ( p 0 s t ) 2 + g 11 ( p 1 ) 2 + g 22 ( p 2 ) 2 + g 33 ( p 3 ) 2 = 0
that is equivalent to: g 00 ( p 0 ) 2 + g 11 ( p 1 ) 2 + g 22 ( p 2 ) 2 + g 33 ( p 3 ) 2 = ( m 0 c ) 2 , where (similar to the (1 + 2)-dimensional case in Section 3.1) m 0 and p μ are the particle’s rest mass and momentum (4-momentum), p μ s t = m 0 k μ , and k μ = ( k 0 , 0 , 0 , 0 ) = ( c / g 00 , 0 , 0 , 0 ) denotes the covariant form of the 4-velocity of particle in stationary reference frame. Now similar to the transformations (108-1)–(108-7), as it would be shown in the following, a crucial property of the quadratic relation (111) would be also its invariance under two certain sets of sign inversions of the components of particle’s momentum: ( p 0 , p 1 , p 2 , p 3 ) , along with similar inversions for the components: p 0 s t , p 1 s t , p 2 s t , p 3 s t (as particular cases), where p 0 s t = m 0 k 0 , p 1 s t = p 2 s t = p 3 s t = 0 . The first set of these includes seven different odd types of the sign inversions (i.e., with odd inversions), and the second set includes eight different even types of the sign inversions (i.e., with even inversions), which could be represented simply by the following two symmetric groups of transformations (based on the formalism of the Lorentz symmetry group of invariant relation (111)), respectively:
The first group includes,
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 1 ) , p 0 s t ( 1 ) , p 1 ( 1 ) , p 2 ( 1 ) , p 3 ( 1 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 2 ) , p 0 s t ( 2 ) , p 1 ( 2 ) , p 2 ( 2 ) , p 3 ( 2 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 3 ) , p 0 s t ( 3 ) , p 1 ( 3 ) , p 2 ( 3 ) , p 3 ( 3 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 4 ) , p 0 s t ( 4 ) , p 1 ( 4 ) , p 2 ( 4 ) , p 3 ( 4 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 5 ) , p 0 s t ( 5 ) , p 1 ( 5 ) , p 2 ( 5 ) , p 3 ( 5 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 6 ) , p 0 s t ( 6 ) , p 1 ( 6 ) , p 2 ( 6 ) , p 3 ( 6 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 7 ) , p 0 s t ( 7 ) , p 1 ( 7 ) , p 2 ( 7 ) , p 3 ( 7 ) )
And the second group is given by, respectively:
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 8 ) , p 0 s t ( 8 ) , p 1 ( 8 ) , p 2 ( 8 ) , p 3 ( 8 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 9 ) , p 0 s t ( 9 ) , p 1 ( 9 ) , p 2 ( 9 ) , p 3 ( 9 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 10 ) , p 0 s t ( 10 ) , p 1 ( 10 ) , p 2 ( 10 ) , p 3 ( 10 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 11 ) , p 0 s t ( 11 ) , p 1 ( 11 ) , p 2 ( 11 ) , p 3 ( 11 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 12 ) , p 0 s t ( 12 ) , p 1 ( 12 ) , p 2 ( 12 ) , p 3 ( 12 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 13 ) , p 0 s t ( 13 ) , p 1 ( 13 ) , p 2 ( 13 ) , p 3 ( 13 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 14 ) , p 0 s t ( 14 ) , p 1 ( 14 ) , p 2 ( 14 ) , p 3 ( 14 ) )
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) = ( p 0 ( 15 ) , p 0 s t ( 15 ) , p 1 ( 15 ) , p 2 ( 15 ) , p 3 ( 15 ) )
where, similar to the transformations (108-9) (as equivalent complex representations of the determined group of transformations (108-1)–(108-7), in (1 + 2)-dimensional space-time), following noncomplex-algebraic values of momentum’s components p μ ( p μ = p μ ), the corresponding complex representations of transformations (111-1)–(111-15) is not a necessary issue in general, however, if the invariant relation (111) is represented formally by equivalent complex form:
g 00 ( p 0 p 0 ) g 00 ( p 0 s t p 0 s t ) + g 11 ( p 1 p 1 ) + g 22 ( p 2 p 2 ) + g 33 ( p 3 p 3 ) = 0
then, along with the set fifteen real-valued transformations (111-1)–(111-15), this relation would be also invariant under these corresponding sets of complex transformations (for b = 1 , 2 , 3 , , 15 ):
( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( ± i p 0 ( b ) , ± i p 0 s t ( b ) , ± i p 1 ( b ) , ± i p 2 ( b ) , ± i p 3 ( b ) ) , ( p 0 , p 0 s t , p 1 , p 2 , p 3 ) ( i p 0 ( b ) , i p 0 s t ( b ) , i p 1 ( b ) , i p 2 ( b ) , i p 3 ( b ) )
In addition, in the following, we show that using the transformations (111-1)–(111-15) (along with their corresponding complex forms (111-17)), a set of fifteen different general covariant field equations would be determined, including two certain groups of simultaneous field Equations (corresponding, respectively, to a group of seven bispinor fields and a group of eight bispinor fields of spin-1 particles) as the particular cases of the (1 + 3)-dimensional form of field Equation (72) (defined with a column matrix of the type (98-2)).
Similar to the set of seven algebraic matrix Equations (109) (determined uniquely as the algebraic equivalent matrix representation of the energy-momentum relation (108)), based on the definite formulation of the system of linear Equation (64) in (1 + 3) space-time dimensions (formulated in terms of the matrices (68)), for the energy-momentum relation (111) (along with the transformations (111-1)–(111-15)) the following two sets of systems of linear equations are also determined uniquely, including respectively a set of seven and a set of eight systems of Equations (with different parametric formalisms). The general parametric solution of each of these systems of linear equations, obeys also the quadratic relation (111) (representing a set of fifteen forms, with different parametric formulations, of the general parametric solutions of quadratic relation (111)). Each of these sets of the systems of linear equations could be represented uniformly by a matrix equation as follows, respectively:
( α μ p μ ( b 1 ) m 0 ( b 1 ) α ˜ μ k μ ) S ( b 1 ) = 0 ,
( α μ p μ ( b 2 ) m 0 ( b 2 ) α ˜ μ k μ ) S ( b 2 ) = 0
where b 1 = 1 , 2 , 3 , , 7 , b 2 = 8 , 9 , , 15 , p μ s t ( b 1 ) = m 0 ( b 1 ) k μ , p μ s t ( b 2 ) = m 0 ( b 2 ) k μ , α μ and α ˜ μ are two contravariant 8 × 8 real matrices (compatible with matrix representations of the Clifford algebra Cℓ1,3) defined solely by Formula (65) and (68), and parametric column matrices S ( b 1 ) and S ( b 2 ) are also given uniquely as follows by two distinct expressions (formulated on the basis of definite parametric formulation of column matrix S (68) in (1 + 3) space-time dimensions):
S ( b 1 ) = [ ( u 0 ( b 1 ) v 1 ( b 1 ) u 1 ( b 1 ) v 0 ( b 1 ) ) w ( u 0 ( b 1 ) v 2 ( b 1 ) u 2 ( b 1 ) v 0 ( b 1 ) ) w ( u 0 ( b 1 ) v 3 ( b 1 ) u 3 ( b 1 ) v 0 ( b 1 ) ) w 0 ( u 3 ( b 1 ) v 2 ( b 1 ) u 2 ( b 1 ) v 3 ( b 1 ) ) w ( u 1 ( b 1 ) v 3 ( b 1 ) u 3 ( b 1 ) v 1 ( b 1 ) ) w ( u 2 ( b 1 ) v 1 ( b 1 ) u 1 ( b 1 ) v 2 ( b 1 ) ) w s ] ,   S ( b 2 ) = [ ( u 0 ( b 2 ) v 1 ( b 2 ) u 1 ( b 2 ) v 0 ( b 2 ) ) w ( u 0 ( b 2 ) v 2 ( b 2 ) u 2 ( b 2 ) v 0 ( b 2 ) ) w ( u 0 ( b 2 ) v 3 ( b 2 ) u 3 ( b 2 ) v 0 ( b 2 ) ) w 0 ( u 3 ( b 2 ) v 2 ( b 2 ) u 2 ( b 2 ) v 3 ( b 2 ) ) w ( u 1 ( b 2 ) v 3 ( b 2 ) u 3 ( b 2 ) v 1 ( b 2 ) ) w ( u 2 ( b 2 ) v 1 ( b 2 ) u 1 ( b 2 ) v 2 ( b 2 ) ) w s ]
which column matrix S ( b 1 ) includes seven cases with specific parametric formulations expressed respectively in terms of seven groups of independent arbitrary parameters: u 0 ( b 1 ) , u 1 ( b 1 ) , u 2 ( b 1 ) , v 0 ( b 1 ) , v 1 ( b 1 ) , v 2 ( b 1 ) , and two common arbitrary parameters s and w (i.e., having the same forms in all of the seven cases of column matrix S ( b 1 ) ), and column matrix S ( b 2 ) also includes eight cases with specific parametric formulations expressed respectively in terms of eight groups of independent arbitrary parameters: u 0 ( b 2 ) , u 1 ( b 2 ) , u 2 ( b 2 ) , v 0 ( b 2 ) , v 1 ( b 2 ) , v 2 ( b 2 ) , and two common arbitrary parameters s and w (with the same forms in all of seven cases of the column matrix S ( b 2 ) ). In addition, similar to the column matrix S ( a ) represented soley by Formula (109-1), the specific parametric expressions (113) of column matrices S ( b 1 ) and S ( b 2 ) in the formulation of matrix Equations (112-1) and (112-2), have been determined specifically by assuming (as a basic assumption in addition to the systematic natural approach of formulating the matrix Equations (112-1) and (112-2), based on the definite formulation of axiomatically determined matrix Equation (64)) the minimum value for total number of arbitrary parameters in both column matrices S ( b 1 ) and S ( b 2 ) , which implies equivalently the minimum value for total number of arbitrary parameters in all of the fifteen simultaneous (different) cases of matrix Equations (112-1) and (112-2) (necessarily with fifteen independent parametric solutions representing totally a certain set of fifteen different equivalent forms of the general parametric solution of quadratic relation (111), based on the general conditions of basic definition of the systems of linear equations corresponding to homogeneous quadratic and higher degree equations, presented in Section 2, and Section 2.2, Section 2.3, Section 2.4 and Section 3.1 concerning the homogenous quadratic equations). In the following, similar to the fundamental general covariant field Equation (109-2), in the derivation of the corresponding field Equations (from matrix Equations (112-1) and (112-2), respectively), we will also use the above particular algebraic properties of parameters s and s which, respectively, have been expressed commonly in the expressions of all of seven simultaneous cases of matrix Equation (112-1), and in the expressions of all of eight cases of matrix Equation (112-1).
Moreover, similar to the invariant relation (108) and derived matrix Equation (109), along with the transformations (111-1)–(111-15) and algebraic matrix Equations (112-1) and (112-2), using the corresponding complex transformations (111-17), we may also formally have the following equivalent matrix Equations (with the complex expression), respectively (for b 1 = 1 , 2 , 3 , , 7 , b 2 = 8 , 9 , , 15 ):
( i α μ p μ ( b 1 ) i m 0 ( b 1 ) α ˜ μ k μ ) S ( b 1 ) = 0 ,
( i α μ p μ ( b 2 ) i m 0 ( b 2 ) α ˜ μ k μ ) S ( b 2 ) = 0
Where p μ s t ( b 1 ) = m 0 ( b 1 ) k μ , p μ s t ( b 2 ) = m 0 ( b 2 ) k μ . Similar to the matrix Equation (109-2), though based on the real value of momentum p μ ( p μ = p μ ) , the complex expression of each of the seven cases of algebraic matrix Equation (112-1), and also each of the eight cases of algebraic matrix Equation (112-2), definitely, is not a necessary issue at the present stage. However, since the corresponding momentum operator p ^ μ has a complex value (where p ^ μ p ^ μ ) , in the following, using this basic property of the operator p ^ μ , we derive, distinctly, two certain groups of the general covariant field equations, including a group of seven different simultaneous field equations from the matrix Equations (112-1) and (112-3), and a group of eight different simultaneous field equations from the matrix Equations (112-2) and (112-4) (based on the general axiomatic approach of derivation of general covariant massive field Equations (72) presented in Section 3.4, Section 3.5, Section 3.6, Section 3.7, Section 3.8, Section 3.9 and Section 3.10, and the quantum representations of C, P and T symmetries of this equation, presented in Section 3.10, Section 3.11, Section 3.12, Section 3.13 and Section 3.14). Furthermore, in the following, it is also shown that each of these determined two sets of seven and eight simultaneous field equations describe, respectively, a uniform group of seven spin-1boson fields (corresponding to two left-handed massive charged bosons, along with their two complementary right-handed bosons; a left-handed massive charge-less boson, along with its complementary right-handed boson; and a single right-handed massless and charge-less boson), and a uniform group of eight spin-1 boson field (corresponding to eight massless charged bosons).
Hence, similar to the (1 + 2)-dimensional general covariant field Equation (114), based on the axiomatic approach of derivation of the (1 + 3)-dimensional case of field Equation (72) (defined solely by a column matrix of the form (98-2) in (1 + 3) space-time dimensions, as shown in Section 3.12), from the matrix Equations (112-1), (112-3) and (112-2), (112-4) (defined solely by column matrices (113)) and taking into account this basic momentum operator’s property: p ^ μ p ^ μ , the following two unique groups of seven and eight simultaneous general covariant field equations are determined solely, respectively:
( i α μ D μ ( b 1 ) m 0 ( b 1 ) α ˜ μ k μ ) Φ Z ( b 1 ) = 0
( i α μ D μ ( b 2 ) m 0 ( b 2 ) α ˜ μ k μ ) Φ G ( b 2 ) = 0
specifying by the following two groups of transformations (based on their two corresponding groups of (sign) transformations (111-1)–(111-7), (111-8)–(111-15) and (111-17)), for b 1 = 1 , 2 , 3 , , 7 and b 2 = 1 , 2 , 3 , , 8 , respectively:
The first group includes,
( D 0 ( 1 ) , m 0 ( 1 ) , D 1 ( 1 ) , D 2 ( 1 ) , D 3 ( 1 ) ) = ( D 0 , m 0 , D 1 , D 2 , D 3 )
( D 0 ( 2 ) , m 0 ( 2 ) , D 1 ( 2 ) , D 2 ( 2 ) , D 3 ( 2 ) ) = ( D 0 , m 0 , D 1 , D 2 , D 3 )
( D 0 ( 3 ) , m 0 ( 3 ) , D 1 ( 3 ) , D 2 ( 3 ) , D 3 ( 3 ) ) = ( D 0 , m 0 , D 1 , D 2 , D 3 )
( D 0 ( 4 ) , m 0 ( 4 ) , D 1 ( 4 ) , D 2 ( 4 ) , D 3 ( 4 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
( D 0 ( 5 ) , m 0 ( 5 ) , D 1 ( 5 ) , D 2 ( 5 ) , D 3 ( 5 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
( D 0 ( 6 ) , m 0 ( 6 ) , D 1 ( 6 ) , D 2 ( 6 ) , D 3 ( 6 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
( D 0 ( 7 ) , m 0 ( 7 ) , D 1 ( 7 ) , D 2 ( 7 ) , D 3 ( 7 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
and the second group is given as follows, respectively:
( D 0 ( 1 ) , m 0 ( 1 ) , D 1 ( 1 ) , D 2 ( 1 ) , D 3 ( 1 ) ) = ( D 0 , m 0 , D 1 , D 2 , D 3 )
( D 0 ( 2 ) , m 0 ( 2 ) , D 1 ( 2 ) , D 2 ( 2 ) , D 3 ( 2 ) ) = ( D 0 , m 0 , D 1 , D 2 , D 3 )
( D 0 ( 3 ) , m 0 ( 3 ) , D 1 ( 3 ) , D 2 ( 3 ) , D 3 ( 3 ) ) = ( D 0 , m 0 , D 1 , D 2 , D 3 )
( D 0 ( 4 ) , m 0 ( 4 ) , D 1 ( 4 ) , D 2 ( 4 ) , D 3 ( 4 ) ) = ( D 0 , m 0 , D 1 , D 2 , D 3 )
( D 0 ( 5 ) , m 0 ( 5 ) , D 1 ( 5 ) , D 2 ( 5 ) , D 3 ( 5 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
( D 0 ( 6 ) , m 0 ( 6 ) , D 1 ( 6 ) , D 2 ( 6 ) , D 3 ( 6 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
( D 0 ( 7 ) , m 0 ( 7 ) , D 1 ( 7 ) , D 2 ( 7 ) , D 3 ( 7 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
( D 0 ( 8 ) , m 0 ( 8 ) , D 1 ( 8 ) , D 2 ( 8 ) , D 3 ( 8 ) ) = ( i D 0 , i m 0 , i D 1 , i D 2 , i D 3 )
where the column matrices Φ Z ( b 1 ) and Φ G ( b 1 ) are also given as follows, written on the basis of definite formulations of algebraic column matrices (113) and in addition to the unique formulation of column matrix (98-2) (expressing the general representation of column matrices definable in the formulation of (1 + 3)-dimensional case of general covariant field Equation (72)):
Φ Z ( b 1 ) = [ Z 10 ( b 1 ) Z 20 ( b 1 ) Z 30 ( b 1 ) 0 Z 23 ( b 1 ) Z 31 ( b 1 ) Z 31 ( b 1 ) ϕ Z ] ,   J μ ( b 1 ) = ( D μ ( b 1 ) + i m 0 ( b 1 ) k μ ) ϕ Z ,   Φ G ( b 2 ) = [ G 10 ( b 2 ) G 20 ( b 2 ) G 30 ( b 2 ) 0 G 23 ( b 2 ) G 31 ( b 2 ) G 31 ( b 2 ) ϕ G ] ,   J μ ( b 2 ) = ( D μ ( b 2 ) + i m 0 ( b 2 ) k μ ) ϕ G
where in all of the seven simultaneous (different) field Equations (112-1) formulated with column matrix Φ Z ( b 1 ) (for b 1 = 1 , 2 , 3 , , 7 ), and also in all of the eight simultaneous (different) field Equations (112-2) formulated with column matrix Φ G ( b 2 ) (for b 2 = 1 , 2 , 3 , , 8 ), the scalar quantity ϕ Z (as initially given quantity) defines a commonly set of seven source currents J μ ( b 1 ) , and scalar quantity ϕ G also defines commonly set of eight source currents J μ ( b 2 ) .
Following the definite formulations of a set of seven field Equation (114-1), and set of eight field Equation (114-2) specified, respectively, by the transformations (114-1-1)–(114-1-7) and (114-2-1)–(114-2-8), these two sets of the field equations could be represented uniformly by the following general covariant field equations as well (defined solely in (1 + 3) space-time dimensions), respectively:
( i α μ D μ m 0 α ˜ μ k μ ) Φ Z = 0
( i α μ D μ m 0 α ˜ μ k μ ) Φ G = 0
where the column matrices Φ Z and Φ G are given by:
Φ Z = [ Z 10 Z 20 Z 30 0 Z 23 Z 31 Z 31 ϕ Z ] ,   J μ = ( D μ + i m 0 k μ ) ϕ Z ,   Φ G = [ G 10 G 20 G 30 0 G 23 G 31 G 31 ϕ G ] ,   J μ = ( D μ + i m 0 k μ ) ϕ G
and the field strength tensors Z μ ν , G μ ν and scalars ϕ Z and ϕ G , along with the source currents J μ and J μ are defined as follows:
Z μ ν = b 1 = 1 7 Z μ ν ( b 1 ) τ b 1 ,   G μ ν = b 2 = 1 8 G μ ν ( b 2 ) λ b 2 , I 2 ( D μ + i m 0 k μ ) ϕ Z = b 1 = 1 7 ( D μ ( b 1 ) + i m 0 ( b 1 ) k μ ) ϕ Z = b 1 = 1 7 J μ ( b 1 ) τ b 1 = J μ , I 3 ( D μ + i m 0 k μ ) ϕ G = b 2 = 1 8 ( D μ ( b 2 ) + i m 0 ( b 2 ) k μ ) ϕ G = b 2 = 1 8 J μ ( b 2 ) λ b 2 = J μ .
where I 2 , I 3 are 2 × 2 and 3 × 3 identity matrices, and τ b 1 = l b 1 2 (for b 1 = 1 , 2 , 3 , , 7 ) are the following set of seven 2 × 2 complex matrices:
τ b 1 = l b 1 2 , { l 1 , l 2 , l 3 , l 4 , l 5 , l 6 , l 7 } = { [ 0 i i 0 ] , [ 0 1 1 0 ] , [ i 0 0 i ] , [ 1 0 0 1 ] , [ 0 1 1 0 ] , [ 0 i i 0 ] , [ 1 0 0 1 ] }
which are similar to the set of matrices (110-12), and uniformly represent a combined gauge symmetry group of the form: SU(2)L⊗U(2)R, where the subset of three matrices “ τ 1 , τ 2 , τ 3 ” corresponds to SU(2)L group, and subset of four matrices “ τ 4 , τ 5 , τ 6 , τ 7 ” corresponds to U(2)R group.
The matrices λ b 2 = ( 1 / 2 ) λ b 2 (for b 2 = 1 , 2 , 3 , , 8 ) are also the following set of eight 3 × 3 complex matrices equivalent to the Gell-Mann matrices (representing the SU(3) gauge symmetry group):
λ b 2 = 1 2 λ b 2 ,   { λ 1 , λ 2 , λ 3 , λ 4 , λ 5 , λ 6 , λ 7 , λ 8 } = { [ 0 1 0 1 0 0 0 0 0 ] , [ 0 i 0 i 0 0 0 0 0 ] , [ 1 0 0 0 1 0 0 0 0 ] , [ 0 0 1 0 0 0 1 0 0 ] , [ 0 0 i 0 0 0 i 0 0 ] , [ 0 0 0 0 0 1 0 1 0 ] , [ 0 0 0 0 0 i 0 i 0 ] , 1 3 [ 1 0 0 0 1 0 0 0 2 ] }
Now based on the definite matrix formalisms of the field strength tensors Z μ ν and G μ ν (114-7) (described respectively by general covariant massive field Equations (114-4) and (114-5)), and on the basis of C, P and T symmetries of these field Equations (as two particular forms of the (1 + 3)-dimensional case of tensor field Equation (72)), represented by their corresponding quantum operators (defined in Section 3.11, Section 3.12, Section 3.13 and Section 3.14), it could be concluded that the field Equation (114-4) describes uniformly a definite group of seven simultaneous bispinor fields of spin-1 particles (corresponding to seven matrices τ 1 , τ 2 , τ 3 , τ 4 , τ 5 , τ 6 , τ 7 (114-8)), including, respectively: “three left-handed massive bosons that could be denoted by W , W ˜ + , Z (represented respectively by three matrices τ 1 , τ 2 , τ 3 , corresponding with SU(2)L), a right-handed charge-less spin-1 boson and also three right-handed spin-1 (massive) bosons denoted by ψ , W + , W ˜ , Z (represented respectively by four matrices τ 4 , τ 5 , τ 6 , τ 7 , corresponding with U(2)R), as the complementary particles of three left-handed bosons represented by matrices τ 1 , τ 2 , τ 3 ”.
In addition, following the definite representations of these determined seven bosons, it could be concluded that four bosons ( ψ , W , W + , Z ) correspond to the known bosons including photon (determined as a right-handed charge-less boson, compatible with the positive-frequency corresponding to the right-handed circular polarization state of photon), and W , W + , Z bosons. Hence, particles W ˜ + , W ˜ , Z represent three new massive spin-1 bosons (where, in particular, Z is the complementary right-handed particle of ordinary Z boson), predicted uniquely by this new mathematical axiomatic approach. Furthermore, the field Equation (114-5) would also describe uniformly a definite group of eight spin-1 boson fields (corresponding respectively to eight matrices λ 1 , λ 2 , λ 3 , λ 4 , λ 5 , λ 6 , λ 7 , λ 8 , representing the SU(3) gauge group).
Furthermore, by assuming the group of seven spin-1/2 fermion fields (described by field Equations (110-9)–(110-12)) as the source currents of spin-1 boson fields, it would be also concluded that the group of seven uniform spin-1 boson fields Z μ ν ( b 1 ) represented by ( W , W ˜ + , Z ) and ( ψ , W + , W ˜ , Z ) (describing by general covariant field Equation (114-4)), and the group of eight uniform spin-1 boson fields G μ ν ( b 2 ) (describing by general covariant field Equation (114-5)), hold certain properties (including the electrical and color charges, so on) compatible with the known properties of ordinary bosons W , W + , Z and photon, and also eight gluon fields (with their known definite properties, including the color charges represented by ‘color octet’ [6,7]). In addition, based on the group representation of three additional new bosons that are predicted uniquely by this new mathematical axiomatic approach, denoted by: W ˜ + , W ˜ , Z , these new bosons could have properties similar to the ordinary bosons W , W + , Z ; where in particular new boson Z (as the complementary right-handed particle of ordinary Z boson), can mix with Z boson.
Moreover, as mentioned in Section 3.15, by assuming (as a basic natural assumption) the seven types of spin-1/2 fermion fields describing by general covariant field Equation (110-9), as the source currents of the uniquely determined two groups of seven and eight spin-1 boson fields (describing respectively by general covariant field Equations (114-4) and (114-5)), it would be concluded that there should be, in total, four specific groups of seven spin-1/2 fermion fields (each) with certain properties, corresponding to “1 + 3” generations of four fermions, including two groups of four leptons each, and two groups of four quarks each. Moreover, based on this basic circumstances, two groups of leptons would be represented uniquely by: “[( ν μ , e, ν τ ), ( ν ¯ μ , e+, ν ¯ τ ze)] and [(μ, ν e , τ), (μ+, ν ¯ e , τ+, zn)], respectively, where each group includes a new single right-handed charge-less lepton, represented by: ze and zn”; and two groups of quarks would be also represented uniquely by: “[( s , u , b ), ( s ¯ ,   u ¯ ,   b ¯ , zu)] and [( c , d , t ), ( c ¯ ,   d ¯ ,   t ¯ , zd)], respectively, where similar to leptons, each group includes a new single right-handed charge-less quark, represented by: zu and zd”. In addition, emerging two right-handed charhe-less quarks zu and zd specifically in two subgroups with anti-quarks ( s ¯ ,   u ¯ ,   b ¯ , zu) and ( c ¯ ,   d ¯ ,   t ¯ , zd), could explain the baryon asymmetry, and subsequently, the asymmetry between matter and antimatter in the universe.

4. Conclusions

The main results obtained in this article, are mainly, the outcomes of the new algebraic axiom (17) (along with the basic assumptions (2)–(3) defined in Section 3.1). This new axiom as a definite generalized form of the ordinary axiom of “no zero divisors” of integral domains (including the domain of integers), has been formulated soley in terms of square matrices (with integer entries, appeared as primary objects for representing the integer elements in their corresponding algebraic axiomatic formalism). In Section 3 of this article, as a new mathematical approach to origin of the laws of nature, using a new basic algebraic axiomatic (matrix) formalism based on the ring theory and Clifford algebras (presented in Section 2), “it is shown that certain mathematical forms of fundamental laws of nature, including laws governing the fundamental forces of nature (represented by a set of two definite classes of general covariant massive field equations, with new matrix formalisms), are derived uniquely from only a very few axioms”; where as a basic additional assumption (that is the assumption (2) in Section 3.1), in agreement with the rational Lorentz symmetry group, it has been also assumed that the components of relativistic energy-momentum (D-momentum) can only take the rational values. Concerning the basic assumption of rationality of relativistic energy-momentum, it is necessary to add (as mentioned in Section 3.1) that the rational Lorentz symmetry group is not only dense in the general form of Lorentz group, but also is compatible with the necessary conditions required basically for the formalism of a consistent relativistic quantum theory [15]. In essence, the main scheme of the new mathematical axiomatic approach to fundamental laws of nature presented in Section 3, is as follows. First in Section 3.1, based on the assumption of rationality of D-momentum, by linearization (along with a parameterization procedure) of the Lorentz invariant energy-momentum quadratic relation, a unique set of Lorentz invariant systems of homogeneous linear Equations (with matrix formalisms compatible with certain Clifford, and symmetric algebras) has been derived. Then in Section 3.4, by first quantization (followed by a basic procedure of minimal coupling to space-time geometry) of these determined systems of linear equations, a set of two classes of general covariant massive (tensor) field Equations (with matrix formalisms compatible with certain Clifford, and Weyl algebras) has been derived uniquely as well. Each class of the derived general covariant field equations also includes a definite form of torsion field appeared as generator of the corresponding field’ invariant mass. In addition, in Section 3.4, Section 3.5, Section 3.6, Section 3.7, Section 3.8, Section 3.9, Section 3.10 and Section 3.11, it has been shown that the (1 + 3)-dimensional cases of two classes of derived field equations represent a new general covariant massive formalism of bispinor fields of spin-2, and spin-1 particles, respectively. In fact, these uniquely determined bispinor fields represent a unique set of new generalized massive forms of the laws governing the fundamental forces of nature, including the Einstein (gravitational), Maxwell (electromagnetic) and Yang-Mills (nuclear) field equations. Moreover, it has been also shown that the (1 + 2)-dimensional cases of two classes of these field equations represent (asymptotically) a new general covariant massive formalism of bispinor fields of spin-3/2 and spin-1/2 particles, respectively, corresponding to the Dirac and Rarita–Schwinger equations.
As a particular consequence, in Section 3.6, it has been shown that a certain massive formalism of general relativity–with a definite form of torsion field appeared originally as the generator of gravitational field’s invariant mass–is obtained only by first quantization (followed by a basic procedure of minimal coupling to space-time geometry) of a certain set of special relativistic algebraic matrix equations. In Section 3.9, it has been also proved that Lagrangian densities specified for the originally derived new massive forms of the Maxwell, Yang-Mills and Dirac field equations, are also gauge invariant, where the invariant mass of each field is generated solely by the corresponding torsion field. In addition, in Section 3.10, in agreement with recent astronomical data, a new particular form of massive boson has been identified (corresponding to U(1) gauge group) with invariant mass: mγ ≈ 4.90571 × 10−50 kg, which is specially generated by a coupled torsion field of the background space-time geometry.
Moreover, in Section 3.12, based on the definite mathematical formalism of this new axiomatic approach, along with the C, P and T symmetries (represented basically by the corresponding quantum matrix operators) of uniquely derived two fundamental classes of general covariant field equations, it has been concluded that the universe could be realized solely with the (1 + 2) and (1 + 3)-dimensional space-times (where this conclusion, in particular, is based on the time-reversal symmetry). In Section 3.13 and Section 3.14, it has been proved that ‘CPT’ is the only (unique) combination of C, P, and T symmetries that could be defined as a symmetry for interacting fields. In addition, in Section 3.14, on the basis of these discrete symmetries of derived field equations, it has been also shown that only left-handed particle fields (along with their complementary right-handed fields) could be coupled to the corresponding (any) source currents. Furthermore, in Section 3.15, it has been shown that metric of the background space-time is diagonalized for the uniquely derived fermion field Equations (defined and expressed solely in (1 + 2)-dimensional space-time), where this property generates a certain set of additional symmetries corresponding uniquely to the SU(2)L⊗U(2)R symmetry group for spin-1/2 fermion fields (representing “1 + 3” generations of four fermions, including a group of eight leptons and a group of eight quarks), and also the SU(2)L⊗U(2)R and SU(3) gauge symmetry groups for spin-1 boson fields coupled to the spin-1/2 fermionic source currents. Hence, along with the known elementary particles, eight new elementary particles, including: four new charge-less right-handed spin-1/2 fermions (two leptons and two quarks, represented by “ze , zn and zu , zd”), a spin-3/2 fermion, and also three new spin-1 massive bosons (represented by W ˜ + , W ˜ , Z , where in particular, the new boson Z is complementary right-handed particle of ordinary Z boson), have been predicted uniquely by this fundamental axiomatic approach. As a particular result, in Section 3.6, based on the definite and unique formulation of the derived Maxwell’s Equations (and also determined Yang-Mills equations, represented uniquely with two specific forms of gauge symmetries), it has been also concluded generally that magnetic monopoles could not exist in nature.
The new results obtained in this article, which are connecting with a number of longstanding essential issues in science and philosophy, demonstrate the wide efficiency of a new fundamental algebraic-axiomatic formalism presented in Section 2 of this article.

Acknowledgments

Special thanks are extended to Vitaly L. Ginzburg (Russia), Dmitry V. Shirkov (Russia), Leonid A. Shelepin (Russia), Vladimir Ya. Fainberg (Russia), Wolfgang Rindler (USA), Roman W. Jackiw (USA), Roger Penrose (UK), Steven Weinberg (USA), Ezra T. Newman (USA), Graham Jameson (UK), Sergey A. Reshetnjak (Russia), Sir Michael Atiyah (UK) (who, in particular, kindly encouraged me to continue this work as a new unorthodox mathematical approach to fundamental physics), and many others for their support and valuable guidance during my studies and research.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

The matrix Equation (64) in Minkowski flat space-time (with metric signature (+ – ... –)) would be represented simply by:
( α μ p μ m 0 I ) S = 0
where I is the identity matrix, and column matrix S is defined uniquely by Formula (66)–(70), … in (1 + 1), (1 + 2), (1 + 3), (1 + 4), (1 + 5), … space-time dimensions. The general contravariant forms of real matrices α μ that generate the Clifford algebra Cℓ1,N (for N ≥ 2) in (1 + N)-dimensional space-time, are (as mentioned in Section 3.3), are expressed by Formula (66)–(70),… in various space-times dimensions. Moreover, following the axiomatic approach of derivation of matrix Equation (64), matrices α μ in Minkowski flat space-time also hold the Hermiticity and anti-Hermiticity properties such that: α 0 = ( α 0 ) (compatible with ( α 0 ) 2 = 1 ), and α μ = ( α μ ) (compatible with ( α μ ) 2 = 1 , for μ = 1, 2, 3, …).
These matrices in the (1 + 1), (1 + 2), (1 + 3) and (1 + 4)-dimensional Minkowski space-time (as special cases of their general contravariant forms (65)–(69),…), have the following representations, respectively:
For (1 + 1)-dimensional space-time we have:
α 0 = [ 1 0 0 1 ] ,   α 1 = [ 0 1 1 0 ]
For (1 + 2)-dimensional case we get:
α 0 = [ σ 0 + σ 1 0 0 ( σ 0 + σ 1 ) ] = [ 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ] , α 1 = [ 0 σ 2 σ 3 σ 2 + σ 3 0 ] = [ 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 ] , α 2 = [ 0 σ 1 + σ 0 σ 1 σ 0 0 ] = [ 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 ] .
In (1 + 3) dimensions, we have:
α 0 = [ γ 0 + γ 1 0 0 ( γ 0 + γ 1 ) ] = [ 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 ] , α 1 = [ 0 γ 2 γ 3 γ 2 γ 3 0 ] = [ 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 ] , α 2 = [ 0 γ 4 + γ 5 γ 4 + γ 5 0 ] = [ 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 ] , α 3 = [ 0 γ 6 γ 7 γ 6 γ 7 0 ] = [ 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 ] .
This article has been invited and presented at the following international conferences:
  • The 2016 SIAM International Conference on Mathematical Aspects of Materials Science, Philadelphia, PA, USA, 2016. (https://www.siam.org/meetings/ms16).
  • The 17th International Conference on Quantum Foundations: Quantum and Beyond, International Centre for Mathematical Modeling in Physics, Engineering and Cognitive Sciences (ICMM), Linnaeus University, Sweden, 2016. (https://lnu.se/en/qb/).
  • The 4th International Conference on New Frontiers in Physics, CERN Organized and Supported Conference (Europe), ICNFP-2015. (https://indico.cern.ch/e/icnfp2015/).
  • The 2016 International Conference on Algebraic Geometry and Mathematical Physics, University of Tromsø, Norway, 2016. (https://site.uit.no/).
  • The XXXVII Max Born International Symposium, International Conference on Non-commutative Geometry, Quantum Symmetries and Quantum Gravity (II), Wroclaw University, Poland, 2016. (http://ift.uni.wroc.pl/~mborn37/).
  • The 2016 GRavitational-wave Astronomy International Conference in Paris, Institute d’Astrophysique de Paris (IAP), The University of Paris VI—Sorbonne University, CNRS, LERU, EUA, Paris, France, ‘Supported by the European Union’s 7th Framework Program: FP7/PEOPLE-2011-CIG,’ 2016. (http://www.iap.fr/vie_scientifique/ateliers/GravitationalWave/2016/scripts/abstract.aspx).
  • The 22nd Internnational Australian Institute of Physics Congress (AIP), University of Queensland, Brisbane, Australia, 2016. (http://appc-aip2016.org.au).
  • The 21st International Conference on General Relativity and Gravitation, Columbia University, New York, NY, USA, 2016. (http://www.gr21.org).

References

  1. Oxford University. Oxford English Dictionary Law of Nature; Oxford University Press: Oxford, UK, 2015. [Google Scholar]
  2. Wigner, E. The Unreasonable Effectiveness of Mathematics in the Natural Sciences. Commun. Pure Appl. Math. 1960, 13, 1–14. [Google Scholar] [CrossRef]
  3. Swartz, N. The Concept of Physical Law; Cambridge University Press: New York, NY, USA, 1985. [Google Scholar]
  4. Barrow, J. Theories of Everything: The Quest for Ultimate Explanations, 2nd ed.; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  5. Joos, H. Zur Darstellungstheorie der inhomogenen Lorentzgruppe als Grundlage quantenmechanischer Kinematik. Fortschr. Phys. 1962, 10, 65–146. [Google Scholar] [CrossRef]
  6. Weinberg, S. The Quantum Theory of Fields; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  7. Pal, P.B. An Introductory Course of Particle Physics; CRC Press, Taylor & Francis Publications: Boca Raton, FL, USA, 2014. [Google Scholar]
  8. Ohlsson, T. Relativistic Quantum Physics: From Advanced Quantum Mechanics to Introductory Quantum Field Theory; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  9. Sexl, R.U.; Urbantke, H.K. Relativity, Groups Particles, Special Relativity and Relativistic Symmetry in Field and Particle Physics; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  10. Aitchison, I.J.R.; Hey, A.J.G. Gauge Theories in Particle Physics: From Relativistic Quantum Mechanics to QED, 3rd ed.; IOP Publications: Bristol, UK, 2003. [Google Scholar]
  11. Zahedi, R.A. Linearization Method in the Ring Theory. Bull. Lebedev Phys. Inst. 1997, 6, 21–27. [Google Scholar]
  12. Zahedi, R.A. On Applied Aspects of the Ring Theory. Bull. Lebedev Phys. Inst. 1997, 4, 15–20. [Google Scholar]
  13. Zahedi, R.A. On the Connection Between Methods of the Ring Theory and the Group Approach. Bull. Lebedev Phys. Inst. 1997, 8, 5–11. [Google Scholar]
  14. Zahedi, R.A. The Linearization Method Based on the Ring Theory and its Applications in Physics; Russian Academy of Sciences (RAS) Publications: Moscow, Russia, 1998; Available online: http://dlib.rsl.ru/viewer/01000049410 (accessed on 12 September 2017).
  15. Hill, E.L. Relativistic Theory of Discrete Momentum Space and Discrete Space-Time. Phys. Rev. 1955, 100, 1780–1783. [Google Scholar] [CrossRef]
  16. Dickson, L.E. History of the Theory of Numbers, Volume II: Diophantine Analysis; Dover Publications: Mineola, NY, USA, 2005. [Google Scholar]
  17. Guidry, M.W. Gauge Field Theories: An Introduction with Applications; John Wiley & Sons: Hoboken, NJ, USA, 1991. [Google Scholar]
  18. Lucas, J.R. Gödel’s Incompleteness Theorems; Oxford University Press: Oxford, UK, 1999. [Google Scholar]
  19. Kaye, R. Models of Peano Arithmetic; Oxford University Press: Oxford, UK, 1991. [Google Scholar]
  20. Mendelson, E. Introduction to Mathematical Logic, 4th ed.; Chapman & Hall: London, UK, 1997. [Google Scholar]
  21. Smullyan, R. Conceptual Roots of Mathematics; Routledge: Abingdon-on-Thames, UK, 1991. [Google Scholar]
  22. Durbin, J.R. Modern Algebra: An Introduction, 3rd ed.; John Wiley and Sons: Hoboken, NJ, USA, 1993. [Google Scholar]
  23. Blyth, T.S.; Robertson, E.F. Groups, Rings and Fields: Algebra through Practice, Book 3; Cambridge University Press: Cambridge, UK, 1985. [Google Scholar]
  24. Tesser, S.B. Generalized Clifford Algebras and Their Representations. In Clifford Algebras and Their Applications in Mathematical Physics; Micali, A., Boudet, R., Helmstetter, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; pp. 133–141. [Google Scholar]
  25. Byrnes, J. (Ed.) Computational Noncommutative Algebra and Applications; Springer: Berlin, Germany, 2004. [Google Scholar]
  26. Ramakrishnan, A. Generalized Clifford Algebra and its applications—A new approach to internal quantum numbers. In Proceedings of the Conference on Clifford Algebra, its Generalization and Applications, Matscience (IMSc), Madras, India, 30 January–1 February 1971; Volume 20, pp. 87–96. [Google Scholar]
  27. Morris, O. On a Generalized Clifford Algebra. Q. J. Math. 1967, 18, 7–12. [Google Scholar] [CrossRef]
  28. Morris, O. On a Generalized Clifford Algebra II. Q. J. Math. 1968, 19, 289–299. [Google Scholar] [CrossRef]
  29. Pappacena, C.J. Matrix pencils and a generalized Clifford algebra. Linear Algebra its Appl. 2000, 313, 1–20. [Google Scholar] [CrossRef]
  30. Chapman, A.; Kuo, J. On the generalized Clifford algebra of a monic polynomial. Linear Algebra Appl. 2015, 471, 184–202. [Google Scholar] [CrossRef]
  31. Coskun, E.; Kulkarni, R.S.; Mustopa, Y. On representations of Clifford algebras of Ternary cubic forms. In New Trends in Noncommutative Algebra; American Mathematical Society: Providence, RI, USA, 2012; pp. 91–99. [Google Scholar]
  32. Ablamowicz, R.; Parra, J.; Lounesto, P. (Eds.) Clifford Algebras with Numeric and Symbolic Computations; Springer: Berlin/Heidelberg, Gemany, 2012. [Google Scholar]
  33. Smith, T. Decomposition of generalized Clifford algebras. Q. J. Math. 1991, 42, 105–112. [Google Scholar] [CrossRef]
  34. Zahedi, R. On the Logical Origin of the Laws Governing the Fundamental Forces of Nature: A New Axiomatic Matrix Approach; Preprint; Hokkaido University: Sapporo, Japan, January 2015; Available online: https://eprints.lib.hokudai.ac.jp/dspace/handle/2115/59279 (accessed on 5 September 2017).
  35. Mordell, L.J. Diophantine Equations, Pure and Applied Mathematics; Academic Press: Cambridge, MA, USA, 1969. [Google Scholar]
  36. Ireland, K.; Rosen, M. A Classical Introduction to Modern Number Theory, 2nd ed.; Springer: New York, NY, USA, 1998. [Google Scholar]
  37. Bourbaki, N. Algebra; Springer: Berlin, Germany, 1988. [Google Scholar]
  38. Forshaw, J.R.; Smith, A.G. Dynamics and Relativity; Manchester Physics Series; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  39. Frampton, P. Gauge Field Theories, 3rd ed.; Wiley-VCH: Weinheim, Germany, 2008. [Google Scholar]
  40. Hehl, F.W. Four lectures in Poincare gauge field theory. In Cosmology and Gravitation; Bergmann, P.G., de Sabbata, V., Eds.; Plenum Press: New York, NY, USA, 1980; pp. 5–61. [Google Scholar]
  41. Aldrovandi, R.; Pereira, J.G. An Introduction to Geometrical Physics; World Scientific: Singapore, 1995. [Google Scholar]
  42. De Andrade, V.C.; Pereira, J.G. Gravitational Lorentz Force and the Description of the Gravitational Interaction. Phys. Rev. D 1997, 56, 4689–4695. [Google Scholar] [CrossRef]
  43. Wald, R.M. General Relativity; University of Chicago Press: Chicago, IL, USA, 2010. [Google Scholar]
  44. Fresneda, R.; Baldiotti, M.C.; Pereira, T.S. Maxwell Field with Torsion. Braz. J. Phys. 2015, 45, 353–358. [Google Scholar] [CrossRef]
  45. Gasperini, M. Constraints on unified theories from the experimental tests of the equivalence principle. Gen. Relat. Gravit. 1984, 16, 1031–1037. [Google Scholar] [CrossRef]
  46. Nobili, A.M.; Lucchesi, D.M.; Crosta, M.T.; Shao, M.; Turyshev, S.G.; Peron, R.; Catastini, G.; Anselmi, A.; Zavattini, G. On the universality of free fall, the equivalence principle, and the gravitational redshift. Am. J. Phys. 2013, 81, 527. [Google Scholar] [CrossRef]
  47. Hehl, F.W.; Obukhov, Y.N. How Does the Electromagnetic Field Couple to Gravity, in Particular to Metric, Nonmetricity, Torsion, and Curvature? In Gyros, Clocks, Interferometers...: Testing Relativistic Gravity in Space; Lecture Notes in Physics; Springer: Berlin/Heidelberg, Germany, 2001; Volume 562, pp. 479–504. [Google Scholar]
  48. Baldiotti, M.C.; Elias, W.S.; Molina, C.; Pereira, T.S. Thermodynamics of quantum photon spheres. Phys. Rev. D 2014, 90, 104025. [Google Scholar] [CrossRef]
  49. Kouwn, S.; Oh, P.; Park, C.-G. Massive Photon and Dark Energy. Phys. Rev. D 2016, 93, 083012. [Google Scholar] [CrossRef]
  50. Chandler, D.L. Seeing the Dark: New MIT-Led Experiment Could Finally Shed Light on the Mysteries of Dark Matter. MIT News, 2013. Available online: https://news.mit.edu/2013/dark-matter-1025, https://dmtpc.mit.edu/DarkLight/DarkLightProposal_PAC39.pdf (accessed on 5 September 2017).
  51. Pilkuhn, H. Relativistic Quantum Mechanics; Springer: Berlin/Heidelberg, Germany, 2003; p. 156. [Google Scholar]
  52. Elton, D.; Vassiliev, D. The Dirac Equation without Spinors; Operator Theory: Advances and Applications; Springer: Berlin/Heidelberg, Germany, 1999. [Google Scholar]
  53. Deser, S.; Jackiw, R.; t’Hooft, G. Three-Dimensional Einstein Gravity: Dynamics of Flat Space. Ann. Phys. 1984, 152, 220–235. [Google Scholar] [CrossRef]
  54. Illge, R. Massive Fields of Arbitrary Spin in Curved Space-Times. Commun. Math. Phys. 1993, 158, 433–457. [Google Scholar] [CrossRef]
  55. Weinberg, S. Feynman rules for any spin. Phys. Rev. 1964, 133, B1318–B1332. [Google Scholar] [CrossRef]
  56. Weinberg, S. Feynman rules for any spin. II. Massless particles. Phys. Rev. 1964, 134, B882–B896. [Google Scholar] [CrossRef]
  57. Rindler, W. Relativity: Special, General, and Cosmological; Oxford University Press: New York, NY, USA, 2006. [Google Scholar]
1
2
Besides, we may argue that our presented axiomatic matrix approach (for a direct derivation and formulating the fundamental laws of nature uniquely) is not subject to the Gödel’s incompleteness theorems [18]. As in our axiomatic approach, firstly, we have basically changed (i.e., replaced and generalized) one of the main Peano axioms (when these axioms algebraically are augmented with the operations of addition and multiplication [19,20,21]) for integers, which is the algebraic axiom of “no zero divisors”. Secondly, based on our approach, one of the axiomatic properties of integers (i.e., axiom of “no zero divisors”) could be accomplished solely by the arbitrary square matrices (with integer components). This axiomatic reformulation of algebraic properties of integers thoroughly has been presented in Section 2 of this article.
3
Some of recent citations to this paper include: Gudder, S. Reconditioning in Discrete Quantum Field Theory. Int. J. Theor. Phys. 2017, 1–14; McKenzie, A. The Level I Multiverse is not the same as the Level III Multiverse. arXiv 2017, arXiv:1702.02019; McKenzie, A. Some Remarks on the Mathematical Structure of the Multiverse. PhilSci-Archive, University of Pittsburgh Publication, Pittsburgh, USA, 2016. Paper ID: 11921, URL: http://philsci-archive.pitt.edu/11921/, arXiv 2016, arXiv:1602.04247.

Share and Cite

MDPI and ACS Style

Zahedi, R. A Unique Mathematical Derivation of the Fundamental Laws of Nature Based on a New Algebraic-Axiomatic (Matrix) Approach. Universe 2017, 3, 67. https://doi.org/10.3390/universe3040067

AMA Style

Zahedi R. A Unique Mathematical Derivation of the Fundamental Laws of Nature Based on a New Algebraic-Axiomatic (Matrix) Approach. Universe. 2017; 3(4):67. https://doi.org/10.3390/universe3040067

Chicago/Turabian Style

Zahedi, Ramin. 2017. "A Unique Mathematical Derivation of the Fundamental Laws of Nature Based on a New Algebraic-Axiomatic (Matrix) Approach" Universe 3, no. 4: 67. https://doi.org/10.3390/universe3040067

APA Style

Zahedi, R. (2017). A Unique Mathematical Derivation of the Fundamental Laws of Nature Based on a New Algebraic-Axiomatic (Matrix) Approach. Universe, 3(4), 67. https://doi.org/10.3390/universe3040067

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop