_{1}

^{*}

The Standard Model of Particle Physics treats four fields—the gravitational, electromagnetic, weak and strong fields. These fields are assumed to converge to a single field at the big bang, but the theory has failed to produce this convergence. Our theory proposes
* one *primordial field and analyzes the evolution of this field. The key assumption is that
*only* the primordial field exists—if any change is to occur, it must be based upon self-interaction, as there is nothing other than the field itself to interact with. This can be formalized as the
*Principle*
*of *
*Self-interaction* and the consequences explored. I show that this leads to the linearized Einstein field equations and discuss the key ontological implications of the theory.

Newton, in his analysis of gravity, concluded that it made no sense to imagine a truly empty vacuum with nothing in it; something must be there, to transmit force from one place to another. This “something” is the gravitational field. Later Einstein reached the same conclusion [

The concept of field evolved from Faraday to Maxwell to Hertz [^{st} field [

The current status of gravitational field theory is marked by two questions. The Standard Model poses that electromagnetic, gravitomagnetic, weak and strong fields converge at the big bang, yet this convergence fails without supersymmetry, which has effectively been eliminated at the LHC [

The plan of this paper is as follows:

Section 1, the Introduction, traces the history of fields in physics and discusses the failure of convergence to a primordial field.

Section 2, the Principle of Self-interaction, introduces the concept of primordial field as the original entity that existed or came into existence “in the beginning”. As nothing else existed, any physical interaction could only be with the primordial field itself. Based on this formulation of physics, I propose a “change operator”, acting on the field and represent the change as the field acting on itself. This yields the self-interaction equation. I then solve this equation for a scalar aspect of the field and interpret this scalar as time, in which case the solution has a frequency property. I then consider a vector field aspect and formulate a vector equation.

Section 3 treats “contact with Newton’s equation” as a means of linking the symbolic formulation of self-interaction to physical reality. We find that Newton’s equation of gravity can be derived from the self-interaction equation, and compare the primordial self-interaction to the interaction between separate entities.

Section 4 introduces Hestenes’ Geometric Calculus as the most appropriate mathematical formalism for physics, based on the fact that every geometric algebra entity has both an algebraic interpretation and a geometric interpretation. In addition, Geometric Calculus deals with multi-vectors composed of different types of entities. For example, we found a scalar solution and a vector solution to the self-interaction equation; therefore we combine both solutions into a multi-vector representation. The geometric product of two vectors is introduced as the fundamental operation and the dual operator, i, is also introduced. The self-interaction equation is redefined in the Geometric Calculus formalism.

Section 5 expands the redefined self-interaction equation in terms of the fundamental constituents introduced in section 4, and then the expansion of the equation is regrouped in terms of “like terms” appearing on each side of the equation. This process yields four equations which are presented in terms of the fundamental constituents.

Section 6 interprets the four equations, derived in section 5, in terms of the associated physics; the result is a set of equations known as the “weak field equations” of general relativity.

Section 7 derives the gravitomagnetic wave equations and discusses the recent detection of gravitomagnetic waves.

Section 8 discusses momentum in the primordial field and describes how terms that would seem to cancel mathematically can exist if separated in space. The addition of these terms completes the linearized field equations first derived by Heaviside and later derived from Einstein’s general relativistic nonlinear field equations.

Section 9 discusses the fact that it has been impossible to successfully apply adjunct linear “Lorenz condition” stipulation to Einstein’s “generally covariant” formulation of general relativity. The Geometric Calculus is used to formulate the gravitomagnetic gauge field equations and to derive the Lorenz gauge condition.

Section 10 presents the key finding of the theory of gravity based on the Principle of Self-interaction: the fact that “field strength” does not appear in this theory. That distinguishes our linear gravitomagnetic field equations from the equivalent field equations derived from general relativity. Relativists assume that these equations apply only for “weak” gravitational fields, whereas our theory is “strength-independent”, and applies for all strengths.

Section 11 discusses the fact that Einstein’s nonlinear field equations can be derived iteratively from the linearized field equations we have derived from our theory of the self-interacting primordial field.

Section 12 presents conclusions and reiterates the “strength-independent” nature of the self-interaction solution, enabling a physical interpretation of the Kasner metric solution of Einstein’s equations. Other applications will be treated in future papers. We summarize by emphasizing that self-interaction theory of gravity presents a novel reinterpretation of gravity that addresses current confusions associated with the success of the post-Newtonian approach to relativity.

Perhaps the simplest assumption upon which to base a universe is that the universe either existed, or came into existence, as a primordial entity. That is, “in the beginning” this primordial entity, and nothing else, existed. If nothing else existed, there was nothing to interact with the primordial entity except itself. To be specific we call this entity a physical field.

Our physics tools are generally designed to relate changes in one physical entity to another entity. For example the basic equation

∇ f = s (1)

relates change (represented by operator ∇ , undefined) in a field (represented by f, undefined) to a source s. If our primordial field is represented byf, and change is represented by mathematical operator ∇ , then s does not exist apart from f. And change must have occurred if the primordial field evolved to the current state of our universe. Today, changes occur when things “interact with each other”, but, if nothing else existed to interact with the primordial field, any interaction could only be the field interacting with itself, and we denote this by ff. So change based on self-interaction is described by the equation

∇ f = f f (2)

This fundamental Principle of Self-interaction describes changes in the primordial field f as it interacts with itself.

In the following, we assume minimal knowledge of physics, while at the same time we assume knowledge of logic and mathematics, as necessary. As we cannot model the universe based only on two symbols ∇ and f, we assume some aspect or property p of the field and ask how the field changes with respect to this aspect. In other words ∇ = ∂ / ∂ p where f = f ( p ) . The self-interaction equation becomes:

∂ ∂ p f ( p ) = f ( p ) f ( p ) . (3)

This fundamental equation has solution, f = − p − 1 .

∂ ∂ p ( − p − 1 ) = ( − p − 1 ) ( − p − 1 ) ⇒ − ( − p − 2 ) = + p − 2 (4)

This implies a scalar property or aspect of the primordial universe. The most fundamental scalar in physics is almost certainly time t, therefore we initially identify p = t and f ( t ) = − 1 / t . The field appears to have a frequency aspect and we postpone interpretation of the − sign. Parameter t is cosmological time, the same time everywhere in the universe. It represents a distance (duration in time) from a beginning to the present state. To find a 3-space vector solution to ∇ f = f f we promote ∇ and f to ∇ → and f → for parameter r → = { x , y , z } . The self-interaction equation becomes

∇ → f → = f → f → . (5)

This combination is not well defined in vector calculus, so we expand the meaning of the change operation to project change onto the field, and use the inner product or dot product to represent self-interaction:

∇ → ⋅ f → = f → ⋅ f → . (6)

Physical fields have energy density proportional to the square of the fields; ρ E = f → ⋅ f → = f 2 . For unity speed of light ( c = 1 ) we have mass equivalent density ρ m = ρ E and we find ∇ → ⋅ f → = ρ E .

Einstein’s general relativity field equations are unphysical differential geometry equations unless, and until, they make contact with Newton’s equation; we apply the same criterion to the primordial field equation. Specifically, we write

∇ → ⋅ G → = G → ⋅ G → . (7)

If we assume that f → is gravitational field G → , we recall that, unlike the electromagnetic field, the gravitational field energy is negative, as it is necessary to add energy to free a body captured in a gravity well. Thus energy density ρ E is negative and (for c = 1 ) ρ m = ρ E and the self-interaction equation becomes

∇ → ⋅ G → = − ρ m . (8)

If Newton’s gravitational constant g = 1 , this is seen to be Newton’s equation of the gravitational field. Therefore our primordial field is tentatively identified as the gravitational field. Although it may be assumed that gravity is sourced by mass, Calabi [

G → ( r → ) = 1 r → ≡ r → r 2 such that G → ⋅ G → = 1 r 2 = ρ E (9)

We compare this with Newton’s force law

F → = − g m m ′ r 2 = − g m m ′ r → r 3 ≡ m ′ G → N (10)

describing the force of gravity on test mass m ′ a distance r from the source of the field, m. Let us use unit test mass m ′ = 1 and keep g = 1 and set G → N = F → / m ′ = ( m / r 3 ) r → = ρ m r → . Although we have suppressed the display of m ′ and g, these constants enter into any dimensional check, so that if Newton’s gravitational constant has units | g | = l 3 / m t 2 we find | G → N | = | g ρ r | = | l 3 m t 2 m l 3 l 1 | = l t 2 , which correctly has dimensions of acceleration.

In the following, we set Newton’s gravitational constant g = 1 . This scalar dimensional constant is always present and has dimensions | g | = l 3 / m t 2 whether we display it or suppress its display. Similarly, we let test mass m ′ = 1 retaining the result shown in Equation (10): F → = m ′ G → N . Thus:

Newtonian: G → N = g m r → r 3 (11)

Primordial: G → = g r → r 2 (12)

If the forces are identical these two equations imply m = r . We choose the origin of the gravitational system ( 0 , 0 , 0 ) to be the center-of-mass and consider the test mass m ′ to be located at r → = ( x , y , z ) ≡ ( r , θ , ϕ ) . In Newtonian problems, the mass m is located at the origin and m ′ at r → .

We see that the gravitational field derived from Newton’s force law has dependence G → N ~ r → / r 3 while the solution of the primordial field self-interaction Equation (9) has dependence G → ~ r → / r 2 . How can we explain this difference?

The difference between the gravitational acceleration due to fixed mass m, distance r from the test mass m ′ , and a primordial field G → tested at the same point (see

field G → at r → . For the primordial field, assumed spherically symmetric, Birkhoff’s Shell theorem [

m = r 3 ( ρ ) = r 3 ( G → ⋅ G → ) = r 3 ( 1 r 2 ) = r . (13)

This is the result implied by Equations (11) and (12). Thus the Newtonian force on the test mass at r → , due to mass m inside the sphere, is

G → N = F → m ′ = g m r → r 3 ⇒ g r r → r 3 = g r → r 2 = G ( r ) , (14)

and self-interaction of primordial field at r → is correctly given by the self-interaction equation. In order to derive this result, we modified equation ∇ → f → = f → f → by specializing in the inner product represented by “ ⋅ ”. We next generalize our tool-set to solve the unmodified equation.

Our goal is to create a physical model, or a physical theory of reality, based on minimal knowledge of physics at the time of creation and on the most effective mathematics. As our world has both logical relationships and shapes, the only mathematical field in which every mathematical term has both an algebraic and a geometric interpretation is Hestenes’ geometric calculus [

∫ M d m x ∂ F = ∮ ∂ m d m − 1 x F (15)

This theorem is compatible with and contains Gauss’s theorem, Stokes theorem, Green’s theorem, and the Cauchy integral formula, in coordinate-free formalism [

u → v → = u → ⋅ v → + u → ∧ v → . (16)

The geometric product of two vectors yields a multi-vector consisting of inner product u → ⋅ v → , which is a scalar projecting one vector into the other, and outer product u → ∧ v → , which is a bivector, a directed area representing the rotation of u → into v → . The scalar product is identical to the vector dot product. The bivector can be related to the vector cross product as seen in

u → ∧ v → = i u → × v → (17)

Bivector u → ∧ v → is a directed area with no defined shape, while u → × v → is an axial vector, which is not included as a vector in vector analysis, but which is required to represent magnetic fields. The term i is the dual operator that transforms cross-product into wedge, as shown, or vice versa: u → × v → = − i u → ∧ v → . An axial vector can be envisioned as a vector cross product, but its rotational aspect can be represented as a bivector. The negative sign associated with − t − 1 is interpreted to mean left-handed circulation of the local field with rotational frequency ~ t − 1 .

In geometric calculus, as in vector calculus, the derivative operator ∇ → is viewed as a vector. Therefore the geometric product of ∇ → with field f → is as follows

∇ → f → = ∇ → ⋅ f → + ∇ → ∧ f → . gradient =divergence +curl(18)

This relation gradient = divergence + curl is not true in any other mathematical formalism.

Scalar derivative ∂ t operating on scalar function f ( t ) yields f ( t ) = − t − 1 for self-interaction equation ∂ t f ( t ) = f ( t ) f ( t ) , while a primordial field with aspects of distance in time and space leads to a directional field f → . If these aspects are separable with respect to time and space, we resolve our field into two primary subfields, G → ( r ) and C → ( t ) and express the primordial field

f ¯ = G → + i C → (19)

where G → ( r ) is a vector and C ↔ ( t ) is a bivector formed by the dual operator i operating on the C → ( t ) field vector. The nature of G → ~ r → − 1 and the nature of C → ~ − t − 1 , where the ~ symbol implies proportionality; scalar constants may be required to match experimental measurement of G → ( r ) and C ↔ ( t ) . The overbar denotes a multivector. If ∇ → is the vector derivative with respect to space, and ∂ t is the scalar derivative with respect to time, then for changes in space and time, we generalize change operator ∇ ¯ to include both derivatives, and field f ¯ to include both subfields:

∇ ¯ f ¯ = f ¯ f ¯ ⇒ ( ∇ → + ∂ t ) ( G → + i C → ) = ( G → + i C → ) ( G → + i C → ) . (20)

( ∇ → + ∂ t ) ( G → + i C → ) = ( G → + i C → ) ( G → + i C → ) (21a)

First we multiply out all terms, noting that the dual operator i commutes with all vectors.

∇ → G → + ∂ t G → + i ∇ → C → + i ∂ t C → = G → G → + i G → C → + i C → G → − C → C → (21b)

Next we expand the geometric products on both sides and then group like terms.

∇ → ⋅ G → + i ∇ → × G → + ∂ t G → + i ∇ → ⋅ C → − ∇ → × C → + i ∂ t C → = G → ⋅ G → + i G → × G → + i G → ⋅ C → − G → × C → + i C → ⋅ G → − C → × G → − C → ⋅ C → − i C → × C → (21c)

Observing that a curl of a vector with itself is identically zero we delete terms G → × G → and C → × C → . We also note that G → × C → + C → × G → = 0 . The remaining terms should be grouped by like terms.

∇ → ⋅ G → + i ∇ → × G → + ∂ t G → + i ∇ → ⋅ C → − ∇ → × C → + i ∂ t C → = G → ⋅ G → + i G → ⋅ C → + i C → ⋅ G → − C → ⋅ C → (21d)

First we group scalars; next scalars multiplied by the dual operator, then vector terms and finally vectors multiplied by i. This expansion of the self-interaction equation yields four equations:

∇ → ⋅ G → = G → ⋅ G → − C → ⋅ C → (22a)

i ∇ → ⋅ C → = i 2 G → ⋅ C → (22b)

∂ t G → − ∇ → × C → = 0 (22c)

i ∇ → × G → + i ∂ t C → = 0 (22d)

These equations derive from the self-interaction of the primordial field according to our Self-interaction Principle, based on the simplest assumptions. They are quite explicit, yet to proceed further we need to make use of what more we know of physical reality. For example, we know that physical fields are real and have energy. Ohanian and Ruffini state: [

When we apply our knowledge to the first scalar Equation (22a), we interpret G → as the gravitational field and G → ⋅ G → as the self-energy density of the gravitational field. Ignoring the complication of C → ⋅ C → we have ∇ → ⋅ G → = ρ . Calling on additional information, we know that gravitational energy is negative since we must apply positive energy to a body captured in a gravitational field in order for it to escape the field, hence the mass density ρ ~ − G → ⋅ G → and we obtain:

∇ → ⋅ G → = − ρ Newton’s gravitational equation (23a)

This implies that C → ⋅ C → has positive energy so that equivalent mass density − C → ⋅ C → in Equation (22a) contributes correctly to Newton’s equation. It has recently been shown that rotational energy in molecules is equivalent to mass [

Our next interpretation also relies on analogy with the electromagnetic field, where E → ⋅ B → = 0 , since the fields are orthogonal to each other. Obviously, our C-field is the gravitomagnetic field, sometimes called cogravitation [

∇ → ⋅ C → = 0 . (23b)

Again analogous to ∇ → ⋅ B → = 0 , this implies that no gravitomagnetic “pole” exists and also that the gravitomagnetic field can be derived from a gauge equation C → = ∇ → × A → since ∇ → ⋅ ∇ → × A ¯ ≡ 0 .

We obtained terms − C → × G → and − G → × C → . When we add these and note that G → × C → = − C → × G → we initially assume that these terms cancel. We will revisit these terms later. Until then we have:

∇ → × C → = ∂ t G → (23c)

By now it’s obvious from the electromagnetic analogy that the last equation is

∇ → × G → = − ∂ t C → (23d)

Grouping these for convenience we obtain:

∇ → ⋅ G → = − ρ , ∇ → ⋅ C → = 0 , ∇ → × G → = − ∂ t C → ∇ → × C → = + ∂ t G → (24)

These equations, derived from the Self-interaction Principle, were derived by Oliver Heaviside in 1893, and later from Einstein’s relativistic field equations, as the “weak field equations”.

The Maxwell-like field equations invite the following procedure, based on the vector identity:

∇ → × ( ∇ → × V → ) = ∇ → ( ∇ → ⋅ V → ) − ∇ 2 V → , (25)

where the last term can be written − ∇ → ⋅ ( ∇ → V → ) . The first term on the right vanishes at all times for V → = C → since ∇ → ⋅ C → ≡ 0 . For no mass density ρ = 0 , and for minimal field density ρ ≈ 0 , we have ∇ → ⋅ G → ≈ 0 , leaving the relation ∇ → × ( ∇ → × V → ) = − ∇ 2 V → . Substitute first G → and then C → into this identity.

∇ → × ( ∇ → × G → ) = − ∇ 2 G → ⇒ ∇ → × ( − ∂ C → ∂ t ) = − ∂ ∂ t ( ∇ → × C → ) ⇒ − ∂ 2 G → ∂ t 2 (26a)

∇ → × ( ∇ → × C → ) = − ∇ 2 C → ⇒ ∇ → × ( + ∂ G → ∂ t ) = + ∂ ∂ t ( ∇ → × G → ) ⇒ − ∂ 2 C → ∂ t 2 (26b)

Summarizing, we have obtained the wave equations

− ∇ 2 G → + ∂ 2 G → ∂ t 2 = 0 and − ∇ 2 C → + ∂ 2 C → ∂ t 2 = 0 (27)

Dimensional analysis indicates that a velocity-squared term is needed, so we assume v = 1 and include the symbolic speed in the equation.

− ∇ 2 G → + 1 v 2 ∂ 2 G → ∂ t 2 = 0 , − ∇ 2 C → + 1 v 2 ∂ 2 C → ∂ t 2 = 0 (28)

The 2017 [

c = 1 ε 0 μ 0 ⇒ [ ε ( g ) μ ( g ) ] − 1 / 2 = [ ( − 1 4 π g ) ε ( − 4 π g c 2 ) μ ] − 1 / 2 = ( c − 2 ) − 1 / 2 = c (29)

Since it is now known that velocity v in Equations (28) is equal to the speed of light we observe that the Principle of Self-interaction predicts gravitational waves.

We now revisit the two terms appearing in the expansion of the Self-interaction Equation (21c) which we deleted based on G → × C → + C → × G → = 0 . We did so due to antisymmetry G → × C → = − C → × G → . However since E → × B → is the momentum energy density of the electromagnetic field; we interpret − C → × G → as the momentum energy density of the gravitomagnetic field. The energy density has equivalent mass density, and momentum density implies that energy is moving with velocity v → . Hence we rewrite − C → × G → as − ρ v → and obtain the complete set of Heaviside equations:

∇ → ⋅ G → = − ρ , ∇ → ⋅ C → = 0 , ∇ → × G → = − ∂ t C → ∇ → × C → = − ρ v → + ∂ t G → (30)

This requires some physical explanation. First, we note that G → × C → cancels C → × G → if they are the same vectors. Based on the significance of the gravitational

waves detected from in-spiraling neutron stars, we re-examine this interpretation. Obviously, observed from Earth, the momentum expressed by G → × C → is G → ( r → ) × C → ( t ) as symbolized in

But imagine that the earth is positioned at the far side of the inspiraling stars. In this case, the momentum vector observed would be G → ( − r → ) × C → ( t ) and would have the opposite apparent circulation, changing the sign of the momentum term. The position of Earth is arbitrary, and physics tells us that the inspiraling stars produced gravitational radiation in both directions, so we conclude that, instead of canceling, both momentum terms exist, and that is why we include the ρ v → term in Equation (30). Finally, we generalized from the velocity v = c of the gravitational radiation to encompass subluminal momentum density with v < c . Physical reasoning causes us to restore the terms that we originally canceled for mathematical reasons.

We anticipate problems with energy-momentum tensors in general relativity; unambiguous gauge fields simply cannot be defined, therefore we will reformulate our field Equation (21) by defining ρ = G → ⋅ G → − C → ⋅ C → and p → = ρ v → ~ G → × C → where mass density current p → (momentum density) is the analog of the electromagnetic charge density current j → . The source multivector becomes p ¯ = ρ + p → ≡ ρ ( 1 + v → ) and the reformulated multivector field equation becomes:

( ∇ → + ∂ t ) f ¯ = p ¯ . (31)

As we have seen, this structure supports the wave operator:

∂ μ ∂ μ = ( ∇ 2 − ∂ t 2 ) . (32)

The following, modeled after Arthur [

“There is no way to successfully apply adjunct linear “Lorenz condition” stipulation to the Einstein equation when it is presented in the customary “generally covariant” form:

G μ ν = − ( 8 π g / c 2 ) T μ ν . (33)

Let us multiply Equation (31) by operator ( ∇ → − ∂ t ) to obtain the sourceless wave equation on the left-hand side and the term ( ∇ → − ∂ t ) ( ρ + p → ) on the right,

and recall that the scalar and vector potentials of electromagnetism give rise to scalar wave equations with ρ and j → as multivector source j ¯ = ρ + j → . Arthur suggests that we find a single wave equation relating the four momentum A ¯ to j ¯ . Our analogy replaces current density j → by momentum density p → to obtain

( ∇ 2 − ∂ t 2 ) A ¯ = p ¯ (34)

which, using Equation (31), we rewrite as

( ∇ → + ∂ t ) ( ∇ → − ∂ t ) A ¯ = ( ∇ → + ∂ t ) f ¯ (35)

This allows us to factor ( ∇ → + ∂ t ) from each side to obtain

f ¯ = ( ∇ → − ∂ t ) A ¯ + f ¯ ′ (36)

where f ¯ ′ is any solution of the homogeneous (source free) Equation (28); ( ∇ → + ∂ t ) f ¯ ′ = 0 . For simplicity, we choose f ¯ ′ = 0 , and recall that ( ∇ 2 − ∂ t 2 ) is a scalar. Therefore, if we multiply A ¯ by a scalar in Equation (34), we see that, since the multivector p ¯ is a scalar plus a vector, and A ¯ must have the same form as p ¯ , we must write a multivector

A ¯ = − ϕ + A → (37)

which is analogous to the electromagnetic gauge field four-vector. Next expand Equation (36):

f ¯ = ( ∇ → − ∂ t ) A ¯ = ( ∇ → − ∂ t ) ( − ϕ + A → ) (38)

G → + i C → = − ∇ → ϕ + ∇ → ⋅ A → + ∇ → ∧ A → + ∂ t ϕ − ∂ t A →

As always, each type of term must satisfy the equation separately; so, matching scalar, vector, and bivector terms, we obtain the following equations:

0 = ∂ t ϕ + ∇ → ⋅ A →

G → = − ∇ → ϕ − ∂ t A → (39)

i C → = ∇ → ∧ A → ⇔ C → = ∇ → × A →

Since the field f ¯ has no scalar terms, we set the scalar terms to zero. From Newton’s theory, we have ϕ ~ m / r , thus the gravitational field agrees with Newton plus a gauge term. The last Equation (39c) follows from (23b) ∇ → ⋅ C → = 0 . A dimensional analysis performed almost anywhere along the way will suggest that the gauge field A → has dimensions of velocity A → ~ v → . In analogy with electromagnetic theory, the product of charge q with gauge field A → yields electromagnetic momentum p → q = q A → , therefore our analogous product of mass with gauge field yields gravitomagnetic gauge field momentum p → m = m A → ≡ m v → . Since this is the momentum in gravitomagnetism, we assume our assignment A → ⇒ v → is correct. We check this by examining term − ∂ t A → in Equation (39), where we find that the term represents an acceleration ∂ t v → that is dimensionally compatible with gravitational acceleration G → .

Finally, scalar Equation (39a) is dimensionally correct, since ∇ → ⋅ A → ~ ∂ ∂ x ∂ x ∂ t ~ ∂ ∂ t .

But most significantly, this equation, in electromagnetic field theory is the Lorenz condition, and thus we consider this the gravitomagnetic Lorenz condition that is missing in general relativity. Great effort has been expended in relativity to establish an adjunct stipulation of the Lorenz condition ∂ μ A μ = 0 . Kauffmann substitutes this into the “poster child” of gauge imposition in electromagnetic theory, ∂ μ ∂ μ A ν − ∂ ν ∂ μ A μ = j ν , and simplifies the relation to

∂ μ ∂ μ A ν = j ν . (40)

This “stipulation” is equivalent to Equation (34), suggesting that our approach has been correct. In summary, the much desired gravitomagnetic Lorenz gauge condition that is still missing from general relativity is obtained rather directly from the Principle of Self-interaction.

The key factor concerning our derivation of the Heaviside-Einstein equations (30) is the fact that “field strength” never enters the equation. Our description of the field as “primordial” implies field strengths associated with the big bang, in strong contrast to the century old perception of weak field approximation. The Self-Interaction Principle replaces the “weak field approximation” with the “all-field equations”—the equations hold for all finite strengths of the gravitational field. That is the key lesson to be learned from this theory.

Will observes that “most of our understanding of gravitational radiation has come from approximations to Einstein’s equations.” And Padmanabhan analyzes a Lagrangian for the two body problem in the post-Newtonian approximation and finds that the perihelion precession per orbit “miraculously matches with the corresponding expression for a test body in the Schwarzschild metric. No simple reason for this conclusion is known and it is an issue worth thinking about”.

The self-interaction theory contradicts this prevailing view by applying to any gravitational field, regardless of strength. It addresses Will’s statement that “we have no good understanding of why this approximation to general relativity should be so effective.” Will reviews binary pulsars and inspiralling compact binaries, including black holes, and the surprising fact that the approximate calculations agree with those of numerical relativity for very strong fields. He notes no obvious reason to expect weak field equations to work for inspiralling black holes, but they do.

The Self-Interaction Principle provides the reason. It is “strength-independent”; it never makes assumptions about “weak field” approximation.

The mass-energy density approach of field theory is equivalent to Einstein’s metric-based theory of curved space-time. Although Einstein’s metric-based general relativity is by far the most familiar theory of gravity, a number of approaches have formulated gravity as a gauge theory. In 1954, Gupta constructed a theory in which the “source” couples to the massless spin-2 field h μ ν as the energy-momentum tensor, including the energy momentum of the h μ ν field itself. The coupling induces a cubic term in the Lagrangian, resulting in a corresponding cubic term T 3 μ ν in the energy momentum tensor, which is then included in the source [

The most significant aspect of the gauge approach is that, per Feynman [

Although linear equations are transformable into non-linear via iterative analysis, the nonlinear equations did not just “snap into being”. As Padmanabhan noted, it is necessary to know beforehand that the final field equations have to match with those in Einstein’s theory, in order to introduce the extra assumptions to obtain it, but these extra assumptions are essentially equivalent to the result we’re attempting to derive! So Einstein’s equations did not just fall out of an analysis of gravity; or even the Equivalence Principle; many issues still have yet to be resolved. One simply cannot obtain an expression for the energy-momentum tensor for the spin-2 field that is unique and gauge invariant; instead one can obtain a large class of non-unique theories. This has resulted in numerous “Einstein-like” theories, a number of which are finally being eliminated by the real gravitational wave data that is increasingly being detected.

Ohanian and Ruffini observe that “almost all of the result that had been the subject of experimental investigation can be described by the linear approximation …the deflection of light,the time delay of light,gravitational time dilation,gravitational lensing,and gravitational radiation emerge from the linear approximation.” Recently Will [

Poisson and Will [

In contrast, the Self-Interaction Principle makes no field strength assumption, other than the implicit assumption that the strength of the field at the big bang is included in the theory.

Despite equivalence of the linear formalism to Einstein’s nonlinear form, the “weak field approximation” terminology has misled physicists to believe that real gravitation is described by space-time curvature corresponding to the nonlinear formalism, although Feynman, Padmanabhan, Weinberg, and others insist that curved spacetime is not a necessary conception of gravity.

Familiarity with the Schwarzschild and Kerr metrics convince many that the “proper” theory of gravity is general relativity. Nevertheless, these metric solutions are static; they represent a geometric solution that does not evolve over time, given a fixed mass, M. For a dynamic space-time, an exact metric solution to Einstein’s field equations has existed for over 90 years, yet its interpretation has been “obscure and questionable”. This Kasner metric has recently been interpreted [

We have worked from the assumption that our universe evolved from a single primordial field. The corresponding Principle of Self-interaction produces the known gravitomagnetic field equations; however, in contrast with the “weak field approximation” assumption, a self-interacting field remains self-interacting—it does not become non-self-interacting due simply to a physicist’s deleting nonlinear terms for ease of solution. Therefore the most significant aspect of the derivation of the gravitational field equations from the self-interaction principle is that there is absolutely no mention of field strength. Derivation from the Principle of Self-interaction instead yields all-strength-field-equations of the gravitational field.

The space-time curvature aspect of gravity is based on an approximate equivalence principle that holds only at a mathematical point. On the other hand, there is another assumed equivalence that some find remarkable; Ohanian and Ruffini state:

“That the exact nonlinear equations are implied by the linear equations… is a remarkable feature of Einstein’s theory.”

We restate this as follows:

That Einstein’s exact nonlinear equations are implied by our linear equations… is a remarkable feature of the Principle of Self-interaction.

This offers a new ontological understanding of physical reality.

The author declares no conflicts of interest regarding the publication of this paper.

Klingman, E.E. (2021) The Primordial Principle of Self-Interaction. Journal of Modern Physics, 12, 65-81. https://doi.org/10.4236/jmp.2021.122007