Home Energy Nuclear Electricity Climate Change Lighting Control Contacts Links



By Charles Rhodes, P. Eng., Ph.D.

This web page provides insight as to how energy and forces arise from electric, magnetic and gravitational fields.

An important concept in physics is that the total energy of a particle is finite but the particle's energy distribution in space stretches out to infinity. Hence a high vacuum contains energy due to electric and magnetic fields originating from distant particles.

The concept that electric and magnetic vector fields contain energy leads to well known equations for apparent electric, magnetic and gravitational forces as well as to a better understanding of physical phenomena.

However, the reader is cautioned that General Relativity indicates that gravity is actually a distortion of space-time rather than the imaginary field used for most orbital calculations.

Each isolated particle has corresponding surrounding electric, magnetic and gravitational vector fields that extend out to infinity but which contain only a finite amount of distributed energy.

When multiple particles co-exist their vector fields overlap everywhere in the universe. At every point in the universe there is a net electric field vector, a net magnetic field vector and a net gravitational field vector. The field overlap causes a change in the field energy content of every element of volume. The law of conservation of energy simultaneously causes an equal but opposite change in the kinetic energy content of the same element of volume and/or the emission/absorption of a photon. The apparent force between particle i and all the other particles in a cluster is really just the change in total field energy with respect to a change in the position of particle i in the cluster.

Assume that at every point in space the electric, magnetic and gravitational field vectors from different particles linearly add to give position dependent net electric, net magnetic and net gravitational field vectors.

Then the net electric field vector Fe at position X is given by:
Fe = Sum of all Fei
where Fei denotes the electric field vector at X due to the ith particle.

Then the net magnetic field vector Fm at position X is given by:
Fm = Sum of all Fmi
where Fmi denotes the magnetic field vector at X due to motion of the ith particle.

Then the net gravitational field vector Fg at position X is given by:
Fg = Sum of all Fgi
where Fgi denotes the gravitational field vector at X due to the energy of the ith particle.

Assume that the gravitational, electric and magnetic field vectors of a particle are all mutually mathematically orthogonal. In effect this assumption means that field unit vectors have six dimensions. In addition to the normal X, Y, Z cartesian co-ordinate unit vectors there are imaginary co-ordinate jX, jY, jZ unit vectors, where j = (-1)^0.5. The additional imaginary unit vectors allow gravitational field energy to be negative.

Assume that for each field type the rest energy density of the field is proportional to the net field vector squared.

Then the element of total field energy dEo contained in an element of volume dV = (dX dY dZ) at X can be expressed as:
dEo = [(Ke / 2) |Fe|^2 + (Km / 2) |Fm|^2 - (Kg / 2) |Fg|^2] dV
Fe = net electric field vector at X
Fm = net magnetic field vector at X
Fg = net gravitational field vector at X
Ke = positive natural constant
(Ke = Epsilon = electric permittivity of free space)
Km = positive natural constant
(Km = Mu = magnetic permeability of free space)
Kg = positive natural constant = 1 / (4 Pi G)
G = Newton gravitational constant

Note that when the gravitational fields of normal matter particles add the gravitational field energy density becomes more negative. When a field energy at position (X - Xo) becomes more negative the kinetic energy at position (X - Xo) instantaneously becomes more positive by the same amount so that the total energy at X is unchanged.

It is hypothesized that when the gravitational fields of normal matter and antimatter particles overlap the gravitational field energy becomes more positive. It is further hypothesized that when the gravitational fields of two antimatter particles overlap the gravitational field energy becomes more negative. These two hypotheses, although consistent with observational data, have not been experimentally confirmed.

Note that the above formulation of energy density implicitly assumes that the "dark matter" density is zero and that the energy content of a universe with no energy and no charge is zero. It remains to be seen whether "dark matter" really exists or whether it is a product of limited understanding of the behavior of galaxies. The dopler shift of the arms of rotating spiral galaxies suggests that the galaxies may have more mass than is indicated by the arm component luminosities.

Each particle i can be characterized as having a charge Qi, a magnetic moment Mi, a radius Ri and a rest energy Eoi. The rest energy Eoi includes both non-field energy and the energy content of the particle's electric, magnetic and gravitational fields.

The field energies are integrals over spacial volume of the field energy densities.

Most real systems involve multiple particles. A real system has a Centre of Momentum (CM), which serves as a system position reference point.

The law of conservation of energy requires that for any fully isolated system the total system energy with respect to an inertial observer is constant. Hence, if progressive overlap of vector fields causes a change in field energy the law of conservation of energy requires a corresponding opposite sign change in kinetic energy and photon energy to keep the total energy constant.

If two particles forming an isolated system interact the individual particle energies can remain unchanged or an element of energy dE can be transferred from one particle to the other or a third particle can be formed but the total system energy remains unchanged. Note that a process involving creation of a third particle from an interaction between two particles is usually not reversible because such reversal requires a three body interaction. The probability of a random three body interaction is usually extremely small compared to the probability of a random two body interaction.

A common example is two bodies interacting and liberating a photon in circumstances where the probability of photon capture is very small. Thus conservation of energy in combination with photon emission leads to formation of assemblies of particles mutually bound in potential wells and determines the direction of evolution of many processes in the local universe. eg Condensation of water vapor to form liquid water.

A nearly isolated cluster of particles can absorb or emit photons, which increase or decrease the total cluster energy. Hence the total cluster energy will gradually change until the rate of energy absorption equals the rate of energy emission. The Earth, in its orbit in space around the sun, is an example of a nearly isolated cluster of particles.

Vector fields from different particles add linearly. Hence, overlap of vector fields from multiple particles changes the net local vector field strength linearly. However, the local field energy density is proportional to the net local vector field strength squared. Hence the local field energy density changes nonlinearly as the vector field overlap changes.

A change in the total cluster field energy dEo with respect to a change in a position dXi of a particular particle is expressed as a force on that particle. If the velocity of a particle affects its potential energy (eg an electrically charged particle moving in a magnetic field which creates a secondary magnetic field) then there is an additional dynamic force component.

A force on a particle due to a potential energy gradient causes a corresponding change in energy motion. Energy motion is also known as momentum.

Let Xo be an observers position vector. The force Fi on particle i causes a change in kinetic energy dEki during a change in particle position vector:
d(Xi - Xo).
dEki = Fi * d((Xi - Xo))

Electric and magnetic vector fields have real unit vectors.

Overlap of electric vector fields from charges of opposite sign causes net far field vector cancellation and hence a reduction in total positive field energy, leading to an increase in kinetic energy which appears to be the result of an attractive force. Similarly, overlap of electric vector fields from charges of the same sign causes net far field vector addition and hence an increase in total positive field energy, leading to a reduction in kinetic energy which appears to be the result of a repulsive force.

Magnetic fields, when viewed as originating from small electric current loops, also diminish in proportion to (1 / |X - Xi|^2) but are orientation dependent with respect to sign. As with electric charges, opposite signs lead to an attractive force whereas same signs lead to a repulsive force.

Gravitational fields from normal matter have unit vectors that are proportional to (-1)^0.5, which when squared causes negative potential energy / unit volume. Overlap of gravitational vector fields from mutual proximity of normal matter masses causes far field vector addition and hence an increase in magnitude of the negative field energy. Hence as normal masses approach each other they gain positive kinetic energy to keep the total energy constant, leading to the appearance of an attractive force.

Hence forces are simply a result of the law of conservation of energy.

Note that if anti-matter produces a gravitational field with a unit vector of:
(- i) = [- (-1)^0.5]
as expected from its negative energy, then the resulting gravitational field energy density around antimatter will be negative like ordinary matter. Hence anti matter may be able to form gravitationally bound planets with other antimatter particles.

Pair production lifts energy from below the field free vacuum state to above the field free vacuum state. Hence an anti-matter particle is in effect an energy hole with negative energy. The gravitational field related to this hole should reverse direction. When an electron-positron pair annihilate each other the change in rest mass energy is twice the electron rest mass energy. Thus the energy of the positron is negative with respect to the field free vacuum reference.

The form of the vector field equations suggests that overlap of the gravitational vector field from ordinary matter with the gravitational vector field from antimatter results in the far field vector cancellation which makes the total potential energy less negative and hence causes a force that is repulsive. Hence normal matter and antimatter will gravitationally repel each other. Hence we do not expect to find any free antimatter in our solar system.

The deduced repulsive gravitational force between normal matter and antimatter may explain certain aspects of the intergalactic expansion of the universe that have been observed by astronomers. Two galaxies, one composed of normal matter and the other composed of anti-matter will likely repel each other. Some parties attribute this repulsion to "dark energy". The gravitational interaction between matter and antimatter is presently impossible to resolve in the laboratory because the electric and magnetic forces affecting single particles are many orders of magnitude larger than the gravitational force.

An isolated particle has a finite constant electric charge Q. This charge is the same regardless of the distance of the observer. The electric field E is a result of the contained charge. The surface area at radius R from the nominal position of a particle is (4 Pi R^2). Hence the electric field E at radius R is given by:
E = Q / (4 Pi Epsilon R^2)
Epsilon = natural constant

For an electric field:
Ke = (Epsilon)
Epsilon = permittivity of free space
= 8.85 X 10^-12 coulomb^2 newton^-1 m^-2

Epsilon is one of a handfull of independent natural constants that can only be determined by experimental measurement.

Consider a particle with charge Qi with radius Ri. For R > Ri the electric field Fe around isolated charged particle i is given by:
Fe = Qi / (4 Pi R^2 Epsilon)
Pi = 3.14159
R = radius from the particle center

The corresponding electric field energy density is given by:
(Ke / 2) Fe^2 = (Epsilon / 2)[Qi / (4 Pi R^2 Epsilon)]^2

Then the electric field energy Ee surrounding the charge Qi is given by:
Ee = Integral from Ri to infinity of:
[Epsilon / 2][Qi / 4 Pi R^2 Epsilon]^2] 4 Pi R^2 dR

= (Qi^2 / 8 Pi Epsilon Ri)
This Ee is the electric field rest energy associated with an isolated charge of radius Ri.

Thus an isolated charge Qi has an electric field rest energy Eei of:
Eei = (Qi^2 / 8 Pi Epsilon Ri)

Similarly an isolated charge Qj has an electric field rest energy Eej of:
Eej = (Qj^2 / 8 Pi Epsilon Rj)

If the two charges are both within radius Ro the total electric field rest energy Eet is given by:
Eet = ((Qi + Qj)^2 / 8 Pi Epsilon Ro)
= ((Qi^2 + Qj^2 + 2 Qi Qj) / (8 Pi Epsilon Ro)

Hence the change in total electric field potential energy required to bring two isolated charges together is:
Eet - Eei - Eej
= (2 Qi Qj) / (8 Pi Epsilon Ro)

Differentiating this expression with respect to Ro gives:
F = d(Et - Ei - Ej) / dRo
= - (1 / 4 Pi Epsilon)( Qi Qj / Ro^2)

Recall that the electrostatic force Fe attracting two charges Qi and Qj separated by distance Ro is given by:
Fe = - (1 / 4 Pi Epsilon)(Qi Qj / Ro^2)

Hence we have shown that the electric force is simply the change in electric field energy with respect to position that results from overlap of electric fields.

For a magnetic field the constant Km is given by: Km = (1 / Mu)
Mu = permiability of free space
= 4 Pi X 10^-7 webers / amp-m
= 4 Pi X 10^-7 T-m / amp

1 Tesla (T) = 1 weber / m^-2

Consider a toroidal solenoid. The magnetic field inside the solenoid is given by:
B = (Mu N I) / L
B = magnetic field strength N = number of turns
I = current through each turn
L = average magnetic path length

The magnetic field volume within the toroidal solenoid is:
where A is the cross sectional area of the magnetic field.

The solenoid self inductance is:
N B A / I
= (N A / I) (Mu N I / L)
= (Mu N^2 A) / L

The magnetic field energy Em is:
Em = (inductance) I^2 / 2
= (Mu N^2 I^2 A) / 2 L

The magnetic field energy density is:
(Mu N^2 I^2 A) / (2 L (L A))
= (Mu N^2 I^2) / 2 L^2
= (Mu / 2)(N I / L)^2
=(1 / 2 Mu) B^2
= (Mu / 2) H^2
= (Km / 2) H^2

Hence the magnetic field energy density is given by:
(Km / 2) Fm^2 = (Mu / 2) H^2
Fm = H = (B / Mu)

Thus we have shown the relationship between toroidal solenoid parameters and magnetic field energy density.

In an electrical contactor the closing force results from reducing the length of a magnetic circuit of approximately uniform cross-section. The stored magnetic energy is:
Em = (Mu N^2 I^2 A) / 2 L

The contactor closing force is given by:
dEm / dL = - (Mu N^2 I^2 A) / 2 L^2

Every element of volume has associated with it a gravitational vector field component that contains gravitational field energy. The gravitational field vector is mathematically orthogonal to both electric and magnetic field vectors. Like the electric and magnetic fields the gravitational field has a net vector flux per unit area that when squared is proportional to gravitational field energy density.

However, because gravitational fields contain energy and gravitational fields are themselves the result of energy the resulting forces are not exactly proportional to (1 / R^2).

Assume that each particle is surrounded by a radial gravitational vector field.

For gravity the total external radial vector flux from a particle is proportional to the contained energy.

The total gravitational vector flux due to contained energy Ec is given by:
Flux = + j Ec / Kg C^2
j = (-1)^0.5
Kg = natural constant
C = speed of light
Ec = total energy contained in sphere with radius R.

The surface area of a sphere of radius R is:
(4 Pi R^2).
Hence, for a single isolated particle the external local vector field strength diminishes approximately in proportion to:
1 / (R^2).

The gravitational field flux per unit area at radius R due to energy E contained inside radius R is given by:
Gravitational Field Flux / area = Fg
= j Ec / (C^2 Kg 4 Pi R^2)

Ec = M C^2
so that M = mass contained within radius R.

Assume that the gravitational field energy density is proportional to the square of the vector field flux per unit area. Then the gravitational field energy density at R is given by:
(Gravitational Field Energy / unit volume) = (Kg / 2) (Gravitational Field Flux / unit area)^2

= (Kg / 2) Fg^2
= (Kg / 2) [j Ec / (C^2 Kg 4 Pi R^2)]^2
= (Kg / 2) [j M / Kg 4 Pi R^2]^2
= - [1 / (32 Kg)] [(M^2) / (Pi^2 R^4)]

Note that the gravitational field energy density is negative.

Let R = Ro at the surface of a solid sphere of mass Mo. The gravitational field energy between radius R = Ro and R = infinity is:
Integral from R = Ro to R = infinity of:
- [1 / (32 Kg)] [(M^2) / (Pi^2 R^4)] 4 Pi R^2 dR
= Integral from R = Ro to R = infinity of:
- [1 / (8 Kg Pi)] [(M^2) / (R^2)] dR

Mi = mass inside a sphere of radius Ri where Ri > Ro
Mo = same mass inside radius Ro

However: Mi = Mo + Integral from R = Ro to R = Ri of (dM / dR) dR.

Mass density = (energy density / C^2). Then:
dM / dR = (Field energy density) 4 Pi R^2 / C^2
= - [1 / (32 Kg)] [(M^2) / (Pi^2 R^4)] [4 Pi R^2 / C^2]
= - [1 / (8 Kg)] [(M^2) / (Pi R^2 C^2)]

[dM / M^2] = - [1 / (8 Kg Pi C^2)] [dR / R^2]

Integrate from R = Ro, M = Mo to R = Ri, M = Mi giving:
{[- 1 / Mi] - [- 1 / Mo]} = - [1 / (8 Kg Pi C^2)]{[- 1 / Ri] - [- 1 / Ro]}
{[1 / Mi] - [1 / Mo]} = [1 / (8 Kg Pi C^2)][(1 / Ro) - (1 / Ri)]
{(Mo - Mi) / Mo Mi} = [1 / (8 Kg Pi C^2)] [(1 / Ro) - (1 / Ri)]
(Mo - Mi) = [(Mo Mi) / (8 Kg Pi C^2)][(1 / Ro) - (1 / Ri)]
(Eo - Ei) / C^2 = [(Mo Mi) / (8 Kg Pi C^2)][(1 / Ro) - (1 / Ri)]
(Eo - Ei) = [(Mo Mi) / (8 Kg Pi)] [(1 / Ro) - (1 / Ri)]

Kg = 1 / (G 4 Pi)

(Eo - Ei) = [(G Mo Mi) / (2)][(1 / Ro) - (1 / Ri)]
(Ei - Eo) = - [(G Mo Mi) / 2][(1 / Ro) - (1 / Ri)]
which is the gravitational field energy between R = Ri and R = Ro.

In order to evaluate the full gravitational field energy we let Ri go to infinity. Then:
(Ei - Eo) = - [(G Mo Mi) / 2 Ro]

Remember that in this formula Mo and Mi refer to the same object except that Mi includes its gravity field energy which is negative. Hence (Ei - Eo) is the gravity field energy formed by bringing bits of matter from R = infinity to R = Ro. Note that there is simultaneous creation of an equal amount of positive kinetic energy.

Mi = Mo + (Ei - Eo) / C^2

(Ei - Eo) = - [(G Mo Mo) / 2 Ro] - [(G Mo (Mi - Mo) / 2 Ro]
= - [(G Mo Mo) / 2 Ro] - [(G Mo (Ei - Eo) / 2 C^2 Ro]

(Ei - Eo)[1 + (G Mo / 2 C^2 Ro)] = - [(G Mo Mo) / 2 Ro]
(Ei - Eo) = - [(G Mo^2) / 2 Ro] / [1 + (G Mo / 2 C^2 Ro)]

Gravitational force
= (change in gravitational energy with respect to a change in Ro)
= d(Ei - Eo) / dRo = + {[(G Mo^2) / 2 Ro^2] / [1 + (G Mo / 2 C^2 Ro)]}
- {[(G Mo^2) / 2 Ro](G Mo / 2 C^2 Ro^2) / [1 + (G Mo / 2 C^2 Ro)]^2}
= + {[(G Mo^2) / 2 Ro^2] / [1 + (G Mo / 2 C^2 Ro)]}
- {[(G^2 Mo^3) / 4 C^2 Ro^3] / [1 + (G Mo / 2 C^2 Ro)]^2}
=+ {[(G Mo^2) / 2 Ro^2] / [1 + (G Mo / 2 C^2 Ro)]}
- {[G Mo / 2 C^2 Ro][(G Mo^2) / 2 Ro^2] / [1 + (G Mo / 2 C^2 Ro)]^2}
= {[(G Mo^2) / 2 Ro^2] / [1 + (G Mo / 2 C^2 Ro)]} {1 - ([G Mo / 2 C^2 Ro] / [1 + (G Mo / 2 C^2 Ro)])}

[(G Mo) / (2 C^2 Ro)] << 1
this gravitational force expression simplifies to:
d(Ei - Eo) / dRo ~ [(G Mo^2) / 2 Ro^2] [1 - (G Mo / C^2 Ro)]

Note that the gravitational force is not exactly proportional to 1 / R^2 due to the energy contained within the gravitational field itself. However, this effect is very small and is difficult to detect with nearly circular planetary orbits. Precise measurement of the advance of the perihelion of the elliptical orbit of Mercury was required to experimentally observe this effect.

Newtons gravitational force equation is:
Force = G Ma Mb / R^2

In Newtons method there is no concept of a field.

According to Newton the potential energy released in bringing mass dM from infinity to Ro is:
Integral from R = infinity to R = Ro of: dE = (G M dM) dR / R^2 = - G M dM / Ro

Then the potential energy released in building up M at Ro from M = 0 to M = Mo is:
DeltaE = - G M^2 / 2 Ro

Thus according to Newton moving mass Mo from distributed bits at infinity to a spherical shell at R = Ro causes a change in potential energy :
DeltaE = - G Mo^2 / 2 Ro

Thus the field method and Newton's method give similar results if:
(G Mo / C^2 Ro) << 1
or if
(G Mo / Ro) << C^2

Note that (2 G Mo / Ro)^0.5 is the theoretical escape velocity from a Newtonian system. This parameter becomes significant when the escape velocity is a significant fraction of the speed of light C as is the case in the proximity of a large black hole.

The escape velocity is the minimum radial velocity that an object must have to escape from a gravitational potential well. In Newtonian gravitation:
Gravitational acceleration is:
dV / dT = - G Mo / R^2
dV = - (G Mo) dT / R^2
= - (G Mo) (dT / dR) dR / R^2
= - (G Mo) (1 / V) dR / R^2

V dV = - (G Mo) dR / R^2

Integrating from V = Vo, R = Ro to V = Vf, R = Rf gives:
(Vf^2 / 2) - (Vo^2 / 2) = (G Mo) [(1 / Rf) - (1 / Ro)]

Escape is prevented if Vf = 0 at Rf = infinity
- (Vo^2 / 2) = (G Mo) [ - (1 / Ro)]
Vo^2 = (2 G Mo / Ro)

Thus a particle cannot escape from a gravitational potential well if its initial velocity Vo is:
Vo < (2 G Mo / Ro)^0.5

Hence if:
(2 G Mo / Ro) > C^2
then nothing can escape and there is a black hole.

Consider the value of (G Mo / C^2 Ro) near the surface of the sun:
G = 6.67384 X 10^-11 m^3 kg^-1 s^-2
Ms = 1.98892 X 10^30 kg
Rs = 6.955 X 10^8 m
C = 3 X 10^8 m / s

Substituting in the above formula gives:
[G Mo / (Ro C^2)]
= [6.67384 X 10^-11 m^3 kg^-1 s^-2 X 1.98892 X 10^30 kg] / [6.955 X 10^8 m X 9 X 10^16 m^2 / s^2]
= .212 X 10^-5
= 2.12 X 10^-6

Under favorable circumstances this error fraction of two parts per million in apparent gravitational force can be detected by precise astronomical observations.

This web page last updated April 22, 2018.

Home Energy Nuclear Electricity Climate Change Lighting Control Contacts Links