InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Physics » "Selected Topics in Applications of Quantum Mechanics", book edited by Mohammad Reza Pahlavani, ISBN 978-953-51-2126-8, Published: May 13, 2015 under CC BY 3.0 license. © The Author(s).

Chapter 2

Photons and Signals in the Age of Information

By Cynthia Kolb Whitney
DOI: 10.5772/59067

Article top


Illustration of waveform development.
Figure 1. Illustration of waveform development.
A typical barred spiral galaxy: “A Barred Spiral Galaxy ngc 1300 hubble photo”.
Figure 2. A typical barred spiral galaxy: “A Barred Spiral Galaxy ngc 1300 hubble photo”.
The skeleton for a barred spiral disc galaxy. Originally computed for [20].
Figure 3. The skeleton for a barred spiral disc galaxy. Originally computed for [20].
Log energy loss and gain rates vs. log system radius.
Figure 4. Log energy loss and gain rates vs. log system radius.
More and more balance points below the Hydrogen ground state.
Figure 5. More and more balance points below the Hydrogen ground state.
Positronium solutions.
Figure 6. Positronium solutions.
Same-charge solutions.
Figure 7. Same-charge solutions.

Photons and Signals in the Age of Information

Cynthia Kolb Whitney1

1. Introduction

The history of Physics contains some inexplicable mysteries, and one of them is this: back in the early 20th century, Einstein was working on the idea of the ‘Photon’ and on the idea of the ‘Signal’ at essentially the same time [1,2], but he did not relate them to each other. They seem to have arisen totally separate in his mind, and they led to totally separate subsequent developments. The Photon played a central role in the development of Quantum Mechanics (QM), and the Signal played the central role in the development of Special Relativity Theory (SRT).

QM and SRT are now the two great pillars of early 20th century Physics, but they seem to be in conflict over the issue of communication. In QM, Schrödinger’s Equation is basically the Fourier transform (just a restatement in terms of different variables) of a statement from Classical Mechanics: Total energy = kinetic energy + potential energy. This statement has no signal propagation speed involved in it. So QM appears to allow instantaneous communication over arbitrary distances. But in SRT, Einstein’s Second Postulate limits all communication to light speed, c.

Since QM and SRT conflict so dramatically on the issue of communication, at least one of them must, in some sense, be wrong. So we are left unsure about what to believe, or take as a foundation for future development.

Can we find out anything decisive from experiments or observations?

The foundation for QM is usually called the Quantum Hypothesis, and not the Quantum Postulate, because the granularity aspect of QM is experimentally testable in various ways. So far, the granularity aspect of QM never seems to fail. This fact stands in favor of QM. But many people still do not believe in the instantaneous communication aspect of QM that is manifest in Schrödinger’s Equation.

The foundation for SRT can be called the Light-Speed Postulate, but not the Light-Speed Hypothesis, because light speed really is not experimentally testable, since a test would involve at least two different spacetime points, and the correlation of data from two different spacetime points would involve the Light-Speed Postulate itself. Despite numerous claims to the contrary, SRT has not been tested in a way that actually could have falsified it, and only very indirect testing appears even feasible for SRT.

Einstein’s General Relativity Theory (GRT) flows from SRT, and GRT is testable, at least observationally, although not experimentally. But the new hypotheses that it offers for observational test are few in number, and great in technical difficulty. So even indirect testing of SRT through GRT does not look very promising.

An altogether different approach therefore seems needed: instead of demanding experiments or observations, we should be reviewing the founding ideas themselves, in light of new insights gathered in the intervening century. Evidently, at least one of the founding ideas, and possibly both of them, need to be updated, or else retired and replaced. This paper aims to identify possible update(s)/replacement(s) that may help.

Some of the needed insights come from engineering practice, rather than from theoretical physics. In the mid 20th century there was a flowering of Information Theory (IT), first in connection with wartime code breaking and code making, and then in connection with the post-war communication industry. All of that development led in turn to our modern computation industry, and our current ‘Age of Information’.

IT uses the concept of Entropy, taken from classical Thermodynamics, and with the application of a minus sign, provides a quantitative mathematical measure for Information. This measure can be used in support of all sorts of engineering concept analyses and design decisions, etc. A convenient reference about the IT concepts and their general applications is Leon Brillouin’s wonderful little book Science and Information Theory [3]. Flores Gallegos [4] discusses some particular applications in QM.

Viewed from our vantage point here in the early 21st century, IT actually provides a clear disqualifier for Einstein’s Second Postulate. The problem is this: the Second Postulate is based on the behavior of a classical infinite plane wave, and an infinite plane wave cannot convey any information whatsoever!

The reason for this perhaps startling assertion is that an infinite plane wave is to electromagnetic communication what a steady hum is to auditory communication: background at best. There is no music in the monotonous hum, and there is no message in the infinite plane wave.

Information requires structure: amplitude modulation, or frequency modulation, or on-off switching. An infinite plane wave does not have any such structure. Because of this deficit, we certainly need a new Signal model as the foundation for an updated SRT.

Possibly we also need a better Photon model at the foundation of QM. We do not actually have a detailed and universally accepted model for photons in QM. Attention has focused more on material systems, which are said to change state, and in so doing, emit or absorb photons that carry packets of energy and angular momentum. As for electric and magnetic fields, what we have is the notion from Quantum Electrodynamics of a ‘virtual photon’. This terminology reveals the desirability of having some unified, photon-like, approach for both Coulomb-Ampère and radiation fields, but it does not actually provide details.

Without detail to argue against, we are apparently free to develop a Photon model de novo, in a way that serves not only as the Photon model for QM, but also provides the more realistic Signal model needed for an updated SRT.

The new Photon/Signal model need not involve yet another new Postulate. Remember what Euclid taught the world through his Geometry: use no more Postulates than absolutely necessary. The reason is that unnecessary Postulates can conflict with other Postulates already in place, and so lead to Paradoxes.

In SRT, we do indeed have many Paradoxes, involving rods, clocks, trains, lightening strikes, snakes, barns, twins, and so on, and on. A prime suspect for their root cause is the unnecessary Second Postulate. In QM, Schrödinger’s Equation seemingly came full-formed from heaven, and to that extent was also a Postulate, and indeed one that conflicted with Einstein’s Second Postulate. The mysterious feel of quantum duality, for example, may suggest a possible QM Paradox yet to be fully articulated. If so, the new Photon/Signal model can offer a candidate approach to solve the problem.

The Photon/Signal model is just very familiar, old-fashioned mathematics: 1) partial differential equations (Maxwell’s first order coupled field equations), 2) their family of solutions (starting with Gaussian pulses, and generating by differentiations the higher-and-higher order Hermite polynomials multiplying the Gaussian pulses), and 3) boundary conditions (no backflow of energy behind the source, no overflow of energy beyond the receiver). This formulation is enough to determine the particular solution that fits any particular problem.

One general rule in Science is this: when a venerable and well-tested mathematical approach exists, try it first, before abandoning it in favor of a new approach. The irony is that Einstein chose not to apply usual the mathematical approach, and instead to introduce his additional Postulate. But he did it so long ago that, despite the long list of Paradoxes generated by SRT, his Second Postulate has itself become ‘venerable’, and therefore nearly impossible to unseat.

But with SRT updated with the new Photon/Signal model, based entirely on the old-fashioned mathematical approach, QM no longer needs to conflict with SRT. Atoms can have solution states in which energy loss by radiation is countered by an energy gain mechanism newly identified with the updated SRT. It is no longer necessary to postulate the Schrödinger equation just to prevent atomic death by loss of orbit energy to far-field radiation.

With n new starting point for QM, there can be a new development for QM. Some type of Quantum Gravity has long been sought, but with GRT being founded in SRT, with c-speed-only communication, and with QM being founded in Classical Physics, with instantaneous communication, that goal has been hard to reach. But the combination of an updated SRT and traditional Statistical Mechanics (SM) offers a way forward.

The new protocol for treating gravity is this: 1) First notice that gravitational attraction formally resembles a statistical residue from magnetic interactions between elements of charge-neutral matter that are carrying electrical currents, ‘current elements’ for short. Current elements were well described by Ampère before Maxwell ever came along. 2) Then apply SM to pairs of current elements.

QM also connects to modern technology problems. Chemistry is wonderfully rich application area for QM. But current-day Quantum Chemistry (QC) is not very easy to use, mainly because of heavy computation loads associated with numerical integrations. The new photon/signal model leads to a different approach for QC, one that is algebraic, rather than integral, in character. This paper includes some recent results from the application of Algebraic Chemistry (AC) to research on the chemistry of water; namely, the form of water known as ‘EZ water’, because it excludes positive ions.

Finally, QM may connect to Elementary Particle Physics in a manner yet to be fully developed. Just as the myriad compounds in Chemistry arises from not-very-many chemical elements, some significant part of the myriad of currently understood ‘elementary’ particles may arise from just two of them: the electron and the positron.

The suggestion for this idea lies in Chemistry’s Periodic Table (PT). The way that the electron spin states fill up with increasing nuclear charge suggests the existence of not only electron pairs that have opposing spins, but also electron rings that have aligned spins. Electron rings involve up to three, five, even seven, electrons. Basically, electrons within atoms form, not only opposite-spin couples, but also same-spin teams. Since electron rings apparently do occur in atoms, they may well also occur separated from atomic nuclei, and therefore looking like exotic elementary particles. And the same may be said of positrons. Thus we have a rich array of possibilities yet to explore.

2. The photon / signal model

The first part of the Photon/Signal model consists of the governing partial differential equations. These are Maxwell’s four first-order coupled field equations. Jackson [5] gives Maxwell’s equations in modern notation and Gaussian units as:


Here B is magnetic field and E is electric field. The constant 1/c=ε0μ0, where ε0 is electric permittivity and μ0 is magnetic permeability. In free space, D=ε0E, H=B/μ0, and charge density ρ and current density J are zero. Free space will be the case of interest henceforth in this paper.

The second part of the Photon/Signal model is the family of suitable finite-energy solutions. This is developed as follows. Let the word ‘pulse’ be the short description for a field profile that is rounded on top and sloping down on the sides, and fading gradually to zero; for example, a Gaussian function.

The mechanism for the waveform development is that the spatial derivatives applied in Maxwell’s equations change the original Gaussian pulse into successively longer wavelets consisting of successively higher order Hermite polynomials multiplying the original Gaussian. As a result, the energy in the wavelet gets more and more spread out along the propagation path.

Observe that waveform development is inexorable, just like ever-growing entropy in Thermodynamics. This is, I believe, where Entropy really enters into Physics. That is, Maxwell’s first order coupled field equations are what give to Physics its obvious Arrow of Time. In citing electromagnetism as the cause for irreversibility, this idea follows Bentwich [6]. (It does not, however, give any hope for reversing anything.)

Let the spatial variable argument for the pulse be x. Let the direction of the initial pulse be y. Figure 1 illustrates this Gaussian pulse, along with a snapshot showing how it evolves over one complete cycle through Maxwell’s first-order coupled field equations. Series 1 is the input pulse, and Series 2 is the waveform developed from it. A more complicated Figure was given in Whitney [7]. This simplified Figure 1 focuses on just the one issue: waveform development. The wavelet is shown bold because it has not been given sufficient attention before.


Figure 1.

Illustration of waveform development.

Now, to solve the stated propagation problem, one can pose a primary pair of pulses in E and B to guarantee travel, and, for more realism, add a second such pair, offset a quarter cycle in time and perpendicular in space, to model circular polarization, like a real physical photon exhibits.

The third part of the Photon/Signal model is the pair of propagation boundary conditions: no backflow of energy behind the source, and no overflow of energy beyond the receiver. To guarantee these boundary conditions, one can demand E=0 at these boundaries.

One can fulfill the required zero E fields by matching the signal leaving the source toward the receiver with 1) another fictitious signal going from the source in the opposite direction, so as to make the E field zero at the source, and 2) another fictitious signal approaching the receiver from the opposite direction, so as to make the E field zero at the receiver. One can even continue this boundary-fixing process, to clean up each tiny new departure from zero E at a boundary that each additional fictitious signal creates at the end of the path opposite to the end that it is meant to correct.

Carried to an infinite sum of corrections upon corrections, this is certainly a very complicated picture about electromagnetic fields. How then can we extract from it some simple statement about propagation speed? We have all been trained to just say c. “But relative to what?” we might ask. If we knew the source and receiver were stationary relative to each other, we could be sure that c=crelative to receiver, consistent with Einstein’s Second Postulate. Otherwise, we would have to say more.

If we would focus attention to moments when the bulk of the energy is very near the receiver, we could be fairly confident that c= crelative to receiver, again consistent with Einstein’s Second Postulate. But if we would focus attention to moments when the bulk of the energy is still very near the source, moments when the receiver is nothing more than a distant phantom, we would be hard pressed to argue against the proposition that c= crelative to source. (This would be consistent with the 1908 Ritz Proposal [8], which was much investigated in the early to mid 20th century as a candidate alternative to Einstein’s Second Postulate, but was ultimately rejected.)

If we would need to characterize an entire propagation scenario, we would have to go even further, and consider all moments along the way. We would have neither the Einstein Second Postulate, nor the Ritz Proposal, but rather something else more complicated. It certainly must conform to Einstein’s Second Postulate late in the scenario, when the bulk of the energy is near the receiver. But early in the scenario, when the bulk of the energy is still near the source, it must conform to Ritz’s Proposal. And in between, it must represent in some mathematically appropriate way an idea not previously considered: a transition from one reference for c to the other.

Here is one way to formulate the problem. Let variable x represent distance along the propagation path, and let variable t represent time into the propagation process. At any point x,t there are fields E (x,t) and B (x,t) with magnitudes E(x,t) and B(x,t), and from them a local energy density:


One way to characterize the propagation is with path integrals. [9-13] Consider the following ratio of two path integrals:


This ratio provides a function that begins at the source, and ends at the receiver, and at the temporal midpoint of the scenario, gives equal weight to both the source and the receiver. This behavior captures the proposed transition of the reference to which the light speed c is assumed relative.

3. Using the changing reference for c changes the results

The concept of the changing reference for c to be relative to appears relevant for an important scenario that was considered even before Einstein arrived on the scene. The scenario involves the potentials and fields created by rapidly moving sources, and it was addressed starting in the late nineteenth century and very early twentieth centuries. Researchers then made the same Assumption that Einstein later made his Second Postulate in founding SRT, but they were not attentive enough to see that there indeed was an Assumption, and to call it out as his Second Postulate. So Einstein is to be commended for calling attention to this Assumption.

The sources generally cited for this early (1898 to 1901) problem are A. Liénard [14] and E. Wiechert [15]. Although they worked at about the same time, they worked separately. They got the same results, as did all contemporary and subsequent investigators, because all persons working from then up until now have used the same input Assumption; namely, that the speed of light is always c with respect to the receiver of the light.

The Liénard- Wiechert results are given in [5], and in every other standard EM book. They are displayed in [7], and that short passage is quoted again here, for review and subsequent further discussion:

“The standard scalar and vector potentials are:


“where κ=1 n⋅β, β is source velocity normalized by light speed c, and n=R/R (a unit vector), and R=rsource (t-R/c)-rreceiver (t) (an implicit definition for the terminology ‘retarded’).

“The LW fields obtained from those potentials are then:


“The 1/R fields are radiation fields, and they make a Poynting vector (energy flow per unit area per unit time) that lies along nretarded:

P=c4πEradiative×Bradiative   =c4πEradiative×(nretarded×Eradiative)=c4π(Eradiative)2nretarded   .

“The 1/R2 fields are Coulomb-Ampère fields, and the Coulomb field


“does not lie along nretarded as one might initially expect; instead, it lies along (n-β)retarded. Assume that β does not change much over the total field propagation time, in which case (n-β)retarded is virtually indistinguishable from npresent.”

Thus the Coulomb attraction/repulsion and the radiation Poynting vector have distinctly different directions. This result does not look right physically. It looks as though advance information is being provided on one, but not the other, of two information channels.

Now consider the same problem using the new and more nuanced definition for the light speed reference. Observe that a line that connects the source position at the temporal midpoint of the scenario to the receiver position at the temporal midpoint of the scenario defines both the distance and the direction that the energy must travel in order to achieve the proposed transmission of the signal from the source to the receiver. This specification implies that we need potentials and fields to be, not retarded, but half-retarded. Now the potentials become:

Φ(r,t)=e[1/κR] halfretardedandA(r,t)=e[β/κR] halfretarded.

The fields become:


The Poynting vector becomes:

P=c4πEradiative×Bradiative   =c4πEradiative×(nhalfretarded×Eradiative)=c4π(Eradiative)2nhalfretarded.

The direction of the Coulomb field becomes:


This direction is the same as the direction of the radiation Poynting vector. That is, the Coulomb field and the Poynting vector are now reconciled to the same direction, instead of conflicting with each other.

4. The corrected force direction means the hydrogen atom can survive classically

Many authors have expressed the opinion that really explaining the Hydrogen atom requires some presently-unknown short-range repulsive force between the electron and the proton; see, for example, Lokajicek, et al [16]. But given the results just presented, no mysterious new repulsive force is needed. With the direction of the Coulomb field being nhalf-retarded, there is a tiny tangential component of Coulomb force aligned with the orbit velocity. So there is a torque on the atom, and the torque pumps energy into the atom, and that process can work to balance the energy loss due to radiation.

That is to say: having a more nearly correct model for potentials and fields created by rapidly moving charges makes it possible to explain the immortality of the Hydrogen atom without first postulating the immortality of the Hydrogen atom; i.e., postulating Schrödinger’s equation.

That is to say: we need not postulate Schrödinger’s equation; w can instead just carry out the old-fashioned math.

In my 2012 and 2013 Intech papers, I listed just the pertinent results. Here is more detail and derivation:

Let the masses of the electron and the proton be me and mp. Note that me<<mp, but mp is not infinite.

Let the orbit radii of the electron and the proton be re and rp. Note that rp<<re, but rp is not zero.

Let the charges on the electron and the proton be e and +e.

The magnitude of the nominally attractive force within the atom is F=e2/(re+rp)2.

Let the orbit frequency be Ω. The orbit speed of the electron is ve=reΩ and that of the proton is vp=rpΩ.

The magnitude of the tiny tangential force on the electron is Fe=Fvp/2c=FrpΩ/2c.

The magnitude of the tiny tangential force on the proton is Fp=Fve/2c=FreΩ/2c.

The magnitudes of the torques on the electron and on the proton are Te=reFe and Tp=rpFp, both equal to FrerpΩ/2c.

The total torque on the electron-proton system is Ttotal=Te+Tp=2T=FrerpΩ2/c.

The torque power delivered to the system is Ptorque=TΩ=FrerpΩ2/c.

The squared orbit frequency is determined from either F=mereΩ2 or F=mprpΩ2.

The more convenient of the two options is Ω2=F/mere. With that expression, the approximation rere+rp yields:


This is a reasonably simple expression. But more important than its simplicity is its very existence. The existence of any such expression means that there exists an energy gain mechanism to balance against the known energy loss mechanism, i.e. radiation. This situation provides a chance for balance, allowing the Hydrogen atom to avoid death by energy loss to radiation.

The Ptorque is actually quite large, and so it changes the whole emphasis of worry concerning the Hydrogen atom. The question becomes, not why does the Hydrogen atom not radiate and collapse to death, but rather why does the Hydrogen atom not torque itself up and expand way beyond its known size?

The fact is: there exists much, much more radiation than was previously worried about, and it is enough to produce the proper balance between radiation and torque.

In [17], I said that the extra radiation arises from finite signal propagation speed, which results in circular motion of the center of mass of the Hydrogen atom, which in turn produces Thomas rotation, and thereby scales up by a factor of 2 the overall rotation rate generating the radiation, which increases the radiation power by a factor of 24.

But what should one say about that center-of-mass circular motion? In Newtonian physics, where signal propagation speed is infinite, there is no such thing. In Maxwell physics, the emphasis is on fields, and the responses of individual charges, but not as much on the responses of whole charge systems, such as atoms. So the issue doesn’t come up there. In Einstein’s relativity physics, the emphasis is often on the observers of events more than on the events themselves. System center-of-mass circulation seems not to come up, although Thomas rotation does.

In [17], I noted that Thomas rotation is generally believed to be a result of the properties of Lorentz transformations, and hence of SRT. That is the belief because one can think of Lorentz transformations, not only in the usual, passive sense, as conversion from an observer in one inertial coordinate frame to another observer in another inertial coordinate frame, but also in the active sense, as the application of a ‘boost’ in velocity, the result of an acceleration, the result of a physical force. A series of non-co-linear boosts does indeed produce Thomas rotation.

But in [17] I also remarked that Thomas rotation does arise, not just from Lorentz transformations, but also from Galilean Transformations. That fact can be demonstrated in detail as follows:

For simplicity, let all motion be in the x,y plane. The scenario begins at time coordinate ct0 with one of the particles, say the electron, at rest at spatial coordinates x0,y0. Let an attraction from another particle act in the x direction. Let an increment of velocity ΔVx =ΔV be imposed, and let an increment of time Δt elapse. The coordinates of the electron then become:


Now let an attraction from another particle act in the y direction. Let an increment of velocity ΔVy =ΔV be imposed, and let another increment of time Δt elapse. The coordinates of the electron then become:


Observe that, if the Galilean velocity boosts had been applied in the opposite order, then the ending coordinates of the electron would have been:


Observe that: the squared incremental length changes have the same magnitude either way: (2ΔVΔt)2+(ΔVΔt)2 (ΔVΔt)2+(2ΔVΔt)2=5(ΔVΔt)2. That fact means the two possible sequences of Galilean boost applications differ only by a rotation. That means each one individually contains a rotation equal to half that total angle difference. This is the Thomas rotation.

Let me now go further, and assert that Thomas rotation will arise from any kind of velocity transformation – Lorentz, or Galilean, or any other new kind that may not have a name yet. Thomas rotation is a property of actual reality, not of any particular mathematical model for reality.

With the Thomas rotation included, the total radiation from the atomic system is:

Ptotal radiated=24_2e23c3ae2=(25_e6/me2)/3c3(re+rp)4.

The value of the separation re+rp for which Ptotal radiated=Ptorque is:


In the traditional approach to QM, re+rp=h2/4π2μe2, where μ is the reduced mass, defined by μ1=me1+mp1, and very nearly equal to me, and h is Planck’s constant, 6.626176×1034 Joule-sec, a fundamental constant given by Nature.

The present analysis does not require Planck’s constant as an input. Instead, it provides an estimate of Planck’s constant as an output:

h=4π2μe2(re+rp)4π2mee232mpe2/3me2c2  πe2c128mp/3me6.77×1034Joule-sec   .;

and, for convenience, also an estimate of the often-seen reduced Planck’s constant:


These estimates can be improved by taking due account of the fact that the sines of small angles are not exactly equal to the angles themselves, and the cosines of small angles are not exactly equal to unity. The sine corrections are third order in angle, while the cosine corrections are second order in angle, which is more significant. Including those corrections reduces the radiation power slightly, and so reduces the solution re+rp slightly, and so reduces the estimates of h and slightly – a step in the right direction.

Observe that, in making h an output from, rather than an input to, a theory, the present work follows both Enders [18] and Ralston [19].

5. EM interactions between neutral atoms: A candidate model for gravity?

Our current best understanding of gravity comes from GRT, which is founded in, and developed from, SRT. The parameter c from SRT appears in GRT, and reveals the lineage. Like electromagnetic signals, gravitational signals must have the finite propagation speed c. So might gravitational signals actually be electromagnetic signals? If we choose to update SRT with a more realistic signal model in place of the number c, does that require a similar update for GRT? This Section explores such questions.

First let it be noted that the Einstein gravitational field equations, like Maxwell’s coupled field equations, are a description at the microscopic scale, but the observable phenomena exist at an extremely macroscopic scale.

Figure 2 illustrates a fairly typical barred spiral galaxy. What this image seems to suggest is: there are mass concentrations at the two armpits of this galaxy. Maybe they are mega stars, or black holes. Maybe they orbit another mass concentration at the center of the galaxy, or maybe they just orbit each other. In any event, the two mass concentrations together create a rather structured field of gravitational potential, into which millions, or billions, of smaller stars are entrained, or temporarily detained, as they orbit the galaxy.

Figure 3 shows the skeleton of a potential field created by two super-massive bodies orbiting at half the signal propagation speed. The lines mark minima in gravitational potential as a function of angle around the galaxy. Observe that this skeleton approximately matches this galaxy image.


Figure 2.

A typical barred spiral galaxy: “A Barred Spiral Galaxy ngc 1300 hubble photo”.


Figure 3.

The skeleton for a barred spiral disc galaxy. Originally computed for [20].

Note 1: To persist over time, the orbit speed of the driving two-body system in Figs. 2 and 3 must be such that the rates of energy loss by gravitational radiation and energy gain by torquing balance each other. Like the Hydrogen atom, the galactic-size two-body problem can be solved this way.

Note 2: Since most of the captive stars in Fig 2 orbit at lesser speed than the driving two bodies, and their potential pattern, with its skeleton, Fig. 3, the individual stars in the outer reaches do not keep up with the rotating spiral potential pattern. An individual star sees a recurring ‘density wave’ of neighbor stars first approaching, then receding. Such density waves have long been known, but not well explained.

Note 3: Fig. 3 is constructed using brute-force calculation of potential for lots and lots of location points, and then using numerical search over angles for local minimum values, and then fitting a function to the minima. Note the slightly sinuous bar between the two driving bodies. Even this detail is suggested in Fig. 2.

Now, in order to model gravity in terms of electromagnetic interactions, we need an expression of electromagnetic interaction that is appropriate for neutral atoms. The best expression appears to be one that was known even before Maxwell. André Marie Ampère already had a well-developed theory about forces between what he called ‘current elements’. This term referred to charge-neutral material increments in electrical circuits. In modern times, P. Graneau wrote extensively about Ampère’s theory and experiments; see for example Graneau [21]

Ampère’s theory works perfectly well for ordinary closed circuits, as well as for incomplete broken circuits, such as may exist momentarily in transient situations, like explosive rupture of circuits. Ampère’s theory ought not be forgotten solely on the basis that more modern theory also works perfectly well for closed and stable electrical circuits. Indeed, in some technological applications involving transient situations like ruptures, Ampère’s theory explains more than the modern theory does.

One reason why Ampère’s formulation can sometimes be more powerful than a formulation based on Maxwell’s equations is that Ampère’s formulation can describe a multi-participant scenario straight away, whereas Maxwell’s equations require iteration through successive steps involving both Maxwell’s equations, and the Lorentz force law for the force F acting on a charge q moving with velocity v:

The iteration goes as follows: first, the input charge positions and motions generate E and B fields; second, the Lorentz force law tells how each charge q responds to the fields from the other charges; another step through Maxwell’s equations tells how all the fields change, and so on.

One particular scenario illustrates the difference between the approaches especially well. Consider a current-carrying wire. It is tedious to use the Maxwell-Lorentz-Maxwell-Lorentz iterative approach to arrive at the understanding that the moving electrons favor the surface of the wire, or even the exterior neighborhood near the wire, leaving the interior of the wire depleted of electrons, and therefore in a state of internal repulsion between the remaining positive nuclei. The Ampère’s formulation skips over all this detail, and describes the resulting consequence: the wire experiences internal longitudinal force, and in fact might even rupture. If it does rupture, one can easily tell that the event was not due to ordinary resistive heating and melting, since the fragments are found to be neither hot to touch nor melted in appearance.

The Ampère approach looks promising for gravity problems because any gravity problem is definitely a multi-participant scenario. And for the same reason, the following analysis also invokes ideas from modern Statistical Mechanics.

Ampère’s force formula can be written:


The indices m and n identify two interacting currents. The im and in are current magnitudes. The Δm and Δn are magnitudes of tiny directed length increments Δm and Δn through which the currents flow. The products of currents and directed length increments, im Δm and in Δm, are the current elements. The rm,n is the length of the vector separation rm,n between the current elements. The α, β, and γ are angles with respect to the connecting line between the two current elements, and with respect to each other. Current element im Δm is at angle α from the connecting line, and current element in Δn is at angle β from the connecting line. The γ is the angle between the two planes defined by the connecting line and each of the two current elements, as if the distance rm,n did not separate them. The value ranges are all full circle: 0<α<2π,

One can get a feel for the general behavior of Ampère’s force formula by considering the angle factor 3cosαcosβ2cosγ for a few special cases:

  1. Current elements side-by-side and parallel, as in parallel wires. Both current elements are perpendicular to the connecting line, so α and β are π/2 and cosα and cosβ are zero. But γ is zero, and cosγ=1, so the angle factor evaluates to 2. The force ΔFm,n is then negative. The current elements attract each other. If they reside in parallel wires, the wires attract each other. This you know from experience is true. In a plasma, instead of in a solid wire, it is called the ‘pinch effect’.

  2. Current elements side-by-side, but anti-parallel. This case is just opposite to Case 1 above: now γ=π and cosγ=1. The current elements repel each other. If they reside in a circuit, that circuit likes to straighten out any kinks and enclose more area. This you may know from experience is true.

  3. Current elements end-to-end, as in an electrical circuit. All three angles are zero, all three cosines are unity, and the angle factor evaluates to +1, so the force ΔFm,n is positive. The current elements repel each other.

The 1/(rm,n)2 aspect of the Ampère force law is just like Newton’s law for gravity. Ampère designed his law that way, because, in his time, the greatest prior achievement in Science was Newton’s conquest of gravity. Now we wish to return the favor, and exploit the Ampère Force Law to understand something novel about gravity.

What makes the Ampère current element so potentially appropriate for application to gravity? First of all, it is charge-neutral, like the masses in a gravity scenario. Secondly, its electrons are moving, and although its nucleus is moving too, that motion is not anywhere near as fast. So at all times and all places where matter exists, at the microscopic level a net electron current flows.

The concept that current elements generate forces that can attract or repel each other suggest that pairs of current elements – or pairs of atoms – can be regarded as a system that can have positive or negative total energy. The kinetic part of the energy may be disregarded, since the current elements may be essentially static, but the potential part of the energy is worth paying attention to.

The main novel feature that gravity presents is that we usually have, not two current elements, but huge numbers of atoms, and each atom must have some relationship with all other atoms. The complexity of the situation naturally conjures up ideas from Statistical Mechanics. Here, ideas from Statistical Mechanics are applied to gravity described in terms of Ampère forces between atoms that are viewed as current elements. Some atom-to-atom relationships are momentarily attractive, and some relationships are momentarily repulsive, and all relationships must vary over time. We can look at the population of atom pairs as a whole, and think of it as a statistical ensemble, in which every condition of attraction/repulsion is represented somewhere.

Every area of physics that has statistical ensembles has Gaussian probability functions. In Classical Thermodynamics, a Gaussian probability density function for a random variable, such as a component of a particle momentum vector, implies maximum entropy, subject to a prescribed value for the standard deviation of that random variable. In Quantum Mechanics, a Gaussian probability density function (the squared wave function amplitude) is associated with minimum uncertainty, meaning minimum product of standard deviations in Fourier conjugate variables, like position and momentum.

Sometimes it is not immediately obvious that a problem involves the equivalent of a Gaussian function, because the Gaussian itself involves a squared variable, such as x2 or p2. The squared variable is proportional to some energy E. So one sees probability density functions expressed in the form exp(E/<E>)/<E> where <E> is the average value of energy E, usually something like kT, where k is Boltzmann’s consistent and T is absolute temperature.

In the case of gravity, that energy E of interest is gravitational potential energy. It can have both positive and negative values. The central concept in Statistical Mechanics is that lower-energy states are populated more richly than higher-energy states are. This concept means that any two atoms, viewed as current elements, will be with respect to each other in a state of negative potential energy more often than in a state of positive potential energy. So they will, on average, attract each other more than repel each other. Therefore, <E> is negative. Since temperature cannot be negative, this is something novel.

To deal with negative <E> we really need a probability density function with a Gaussian factor of the form exp[E2/2<E2>], where <E2> is another parameter, positive, but also not related to temperature.

A few examples can illustrate how to find the parameters.

    Let the two energies be Emax and Emin=|Emax|. With the attractive–force, negative-energy state dominating the scenario, <E> is negative, <E>=|<E>|. The two Boltzmann factors are:

    exp(+Emax/|<E>|)for the state with negative energy,Emin=Emax


    exp(Emax/|<E>|)for the state with positive energy, Emax.

    The average energy Eavg must satisfy the definition:


    Simple trial calculations and numerical search of the results works well enough to solve the problem at hand. The solution is approximately:


    (There is also, of course, a positive, and presently irrelevant, solution of the same magnitude.)

      For this case we have we have line integration in place of point evaluation. Eavg becomes:


      Because the denominator is simpler, begin with that. It is:

      |Eavg|sinh(E/|Eavg|)|Emax0=[|Eavg|sinh(Emax/|Eavg|)]                                             =+|Eavg|sinh(Emax/|Eavg|)   .

      This is a positive number.

      The numerator is more complicated, but it can be evaluated using integration by parts:

         Emax0UdV=UV|Emax0Emax0VdU   .

      where U=E and dV=sinh(E/|Eavg|)dE, so that V=|Eavg|cosh(E/|Eavg|). The first term in the numerator evaluation is:

      UV|Emax0=U(Emax)V(Emax)=[Emax|<E>|cosh(Emax/|<E>|)]   =Emax|<E>|cosh(Emax/|<E>|)   .

      The second term in the numerator evaluation is:

      Emax0VdU=Emax0|<E>|cosh(E/|<E>|)dE=|Eavg|2sinh(E/|<E>|)|Emax0    = |<E>|2sinh(Emax/|<E>|)=|<E>|2sinh(Emax/|<E>|)   .

      The sought numerator divided by denominator for <E> is then:

      <E>=Emax|<E>|cosh(Emax/|<E>|)|Eavg|2sinh(Emax/|<E>|)|Eavg|sinh(Emax/|Eavg|)        =Emaxtanh(Emax/|<E>|)+|<E>|   .

      Since the sought <E> is negative, equal to |<E>|, we have:


      Again, numerical search is a practical approach for finding a solution. We find:


      Observe that, as should be expected, this solution is significantly smaller in magnitude than was the solution with only two energy values, ±Emax, to balance between, which came in at 0.833Emax.

        The continuous Gaussian profile,


        extends to infinity in both directions. This attribute is inappropriate for the problem at hand, which definitely possesses limits ±Emax beyond which the modeling problem does not extend. Therefore, let us turn to discrete approximations for a Gaussian. These are based on the binomial expansion for an arbitrary (a+b)n. The binomial coefficients are familiar to many people from Pascal’s famous triangle:

          1  1     1   1     2     1 1     3    3    1etc.

        The numbers in Pascal’s triangle are constructed with addition of neighboring numbers above. This is easy for small n, but small n means crude, and we need refined. So we want large n. So we need a formula involving multiplication instead of addition. That would be:

        1, n, n(n1)1×2,  n(n1)(n2)3!,  , (n)!(n)!(nn)!, n(n1)(n2)3!, n(n1)1×2,n,1.

        Observe that the binomial coefficients are symmetric around the middle of the list, like a Gaussian function is symmetric around zero argument. If n is an even number, the number of binomial coefficients is 2n1, odd, and the middle one, the maximum one, is n!/[(n/2)!]2. If n is an odd number, the number of binomial coefficients is 2n, an even number, and the middle two numbers, the maximum two, are both n!/{[(n+1)/2]![(n1)/2]!}.

        For our modeling problem, let n be an odd number. Let the binomial coefficients be represented as Bb with b=n to b=1. Let them be associated with equal energy increments ΔE=Emax/n starting from Emax and covering the range to zero. Associate the minimum binomial coefficient Bn with the increment starting with Emax, and the maximal binomial coefficient with the increment ending with E=0, and associate the other coefficients with the increments between those limits. The problem to solve is:


        Numerical investigations done to date suggest that the solution comes at approximately <E>=Emax/2n. Here the n for this discrete model is analogous to the standard deviation σ for the corresponding continuous Gaussian.

        The problem of modeling gravity therefore reduces to the problem of determining what value of n should be used. Here is the most pertinent fact: compared to anything electromagnetic, gravity is extremely weak. Consider two Hydrogen atoms at a given separation distance. Let us compare the gravitation force with the maximum Ampère force between them.

        The gravitational attraction is proportional to G(mp)2, where G is the universal gravitation constant, about 6.6×1011 Newton× meter2 per kilogram2, and mp is the mass of the proton, about 1.66×1027kg, so (mp)2 is about 2.76×1054 kg2. Overall,


        Clearly, the maximum Ampère force between atoms viewed as current elements is generously larger than the gravitational force between atoms viewed as charge-neutral masses – by about 32 orders of magnitude!

        But the typical Ampère force is nowhere near as big as the maximum Ampère force, due to the fact that it depends on three angles, any one of which can spoil it. Occurrence of the maximum Ampère force is very rare indeed. And occurrence of the minimum (negative) Ampère force is equally rare. Only the piddling near-zero Ampère forces are common, and even then, each tiny attractive force is mostly cancelled with the tiny repulsive force of equal magnitude but less frequent occurrence. All we have is a tiny residue of attractive force, due to the nature of Boltzmann factors and Statistical Mechanics.

        Of course, being tiny does not mean being insignificant. Like the tiny residue that is microwave background radiation, the tiny residue that is gravity is a possible key to understanding something about the Universe in which we live. Here is an example problem: at present, we know the actual particle radius of the electron is something extremely tiny, but we do not know what its numerical value is. There exist a number of length-dimensioned quantities associated with the electron, all called ‘radius’, but distinguished by specific names and numerical values. MacGregor [22] lists seven of them. Most are on the order of 1013 cm, although one is much smaller, and is presently only upper-bounded at <1016 cm.

        The radius attributed to the electron can have a role in the gravity problem. The ratio of an atomic radius to the electron radius can imply a candidate level of discretization for the binomial approximation to the Gaussian factor involved in the gravity problem. For an atomic radius, let us consider the first orbit radius of Hydrogen, rH1=0.529×108cm. For the electron radius, let us consider two of the possibilities from MacGregor [22].

        One of the electron radii is called classical. This one captures the Coulomb energy equivalence

        14πε0e2/re classical=mec2,

        which implies

        re classical=14πε0(e2/mec2)=2.82×1013cm.

        The ratio ρ=rH1/re classical is then:


        which implies


        The square root of this number would then be the dimensionless n for the discretization:


        This is a number already famous in Physics, but in a context other than gravity. It is the inverse of the so-called ‘fine structure constant’ α, defined as:

        This number plays a role in spectroscopy, where spectral lines occur in families, closely spaced but clearly distinguishable. There, the explanation comes from QM; clearly, another manifestation of natural discretization.

        But if the classical radius of the electron were used in the gravity problem, the ratio of the average Ampère force magnitude to the maximum Ampère force magnitude would be approximately


        This ratio is not appropriately small, so this is not the right discretization level for the gravity problem.

        Another one of the electron radii given by MacGregor [22] is called actual. It characterizes results of scattering experiments, and is the one presently only upper-bounded, at re actual <1016 cm. No one knows how much smaller it could eventually turn out to be. So how much smaller would it have to be, in order to account for the extreme weakness of gravity? Gravity requires 2/n1032, or.


        That in turn requires:

        re actual=rH1/ρ=0.529×108cm/4×10641073cm.

        At present, such a value for re actual certainly looks impossible to test with any kind of measurement. It is smaller than anything we yet know about any elementary particle. But that circumstance may be a good thing, because an extremely small electron makes it easier to understand what data from Chemistry reflects, discussed next. And an extremely small electron, along with a correspondingly small positron, helps explain aspects Elementary Particle Physics, discussed after that.

        6. Algebraic chemistry and EZ water

        Prof. Gerald Pollack of U. Washington wrote the authoritative book [23] about the physical phenomenon called ‘EZ water’. The EZ is short for Exclusion Zone, with the word ‘exclusion’ referring to a surface phenomenon that expels positive hydronium ions.

        Prof. Pollack gave a talk about EZ Water at the 2013 meeting of the Natural Philosophy Alliance at College Park, MD, USA. All the phenomena he described were surprising; some were truly puzzling. EZ water apparently makes extended orderly arrays of hexagonal units. How can that behavior comport with our understanding that Nature maximizes entropy? Explanations then available were not at all quantitative. That fact suggested a real need for a more quantitative approach.

        I had recently written my book about Algebraic Chemistry (AC). [24] The name reflects the fact that the technique has no integrals or other complicated math operations that would demand capabilities beyond those of a hand calculator. The worst operation is square root. So the AC approach looked promising for quick application to EZ water.

        The fundamental idea behind AC is that all atoms share some similarities with Hydrogen atoms: 1) They have a nucleus that is similar to a proton, but scaled up to nuclear charge Z and nuclear mass M ; 2) They have a population of electrons that is not entirely unlike a single electron; i.e., an interacting community that is somewhat coherent, and somewhat like one big electron orbiting the nucleus; 3) It is possible for the electron count to be different from the nuclear charge. This last possibility is what characterizes ions, and thereby creates all of Chemistry.

        We begin with a clue: Eq. (17) indicates that the radius of the Hydrogen atom scales with the mass of the proton. This fact suggests that the base orbit energy of the Hydrogen atom scales with the inverse of proton mass. It further suggests that for element with nuclear charge Z and mass M, the base orbit energy may scale with Z/M. If so, then when first-order ionization potentials for all elements are scaled by the inverse factor, M/Z, then the scaled first-order ionization potentials (called IP1,Z) might fall into some pattern.

        We proceed with an observation: A pattern indeed emerges: the rise on every period in the Periodic Table is exactly the same factor, 7/2.

        We make a Hypothesis: All IP1,Z contain information valuable for all other elements: population generic information. Each IP1,Z contains a universal baseline contribution IP1,1 about interaction between the nucleus and the population of electrons as a whole. For all elements beyond Hydrogen, there is also a contribution ΔIP1,Z about interactions among the electrons.

        The ΔIP1,2 can be very significant. For Helium, ΔIP1,2 is huge, meaning that two electrons bond together very strongly. And for Lithium, ΔIP1,3 is negative, meaning that two electrons actively work together to try to exclude a third electron.

        Over the periods, there is obvious detail about the electron-electron interactions. Within each period, there are obvious sub-periods keyed to the nominal angular-momentum quantum number that is being filled. Plotted on a log scale, all sub-period rises are straight lines. The slopes all appear to be rational fractions. We can display these rational fractions in a Table, as was done in [24] and [25]:


        A non-traditional parameter N is included in the display because, for l>0, it is possible to write a simple formula for the fraction:


        Also, all periods in the Periodic Table have length 2N2.

        All this numerical regularity suggests that there really is a reliable pattern here, and we can reasonably seek to exploit it. Here is the first exploitation that suggests itself: Given first-order ionization potentials of many elements, we can estimate the additional energy required to remove a second electron from each, and then a third, and so on. This was first done in [24]. Formulae were given for each individual electron removal or addition, and evaluated for a large number of elements.

        One point that Ref. [24] emphasized was that the energy to remove a second electron, or a third, and so on, is not the same thing as the so- called ‘second-order ionization potential’, ‘third-order ionization potential’, and so on. Those energies are very large, which implies that those events are very violent: ripping two, or three, or more, electrons off an atom all at once. Those energies do exhibit a lot of numerical regularity, but that isn’t important for understanding typical lab-bench chemistry, which is all about gentle events that occur one-at-a-time.

        Ref [25] presented the equivalent summed formulae for removal of, or addition of, one, two, and three electrons, removed or added one-at-a-time. Basically, use of these formulae save the user some repetitive arithmetic that would be incurred using the formulae from Ref. [24].

        Development of More Formulae:

        For the present paper, the formulae from [25] are extended from the illustrative cases of removing, or adding, one, two, and three electrons, to the general case of removing, or adding, N electrons.

        Ref. [25] used symbols W and H to distinguish between energy increments associated with electron-nucleus interaction, and energy increments associated with electron-electron interactions. That distinction is analogous to the distinction between work and heat in thermodynamics: the work part is something a human can control, and the heat part is something that Nature simply does, regardless of what the human does.

        We had:

        Wremoving e1 from the neutral atom=IP1,1(Z/MZ).

        The ion produced has a little less attraction between the nucleus and the now reduced electron cloud. So removing another electron should take a little less work:

        Wremoving e1 & e2=IP1,1[Z+Z(Z1)]/MZ.

        And then:

        Wremoving e1, e2, & e3=IP1,1[Z+Z(Z1)+Z(Z2)]/MZ.

        This pattern generalizes to:

        Wremoving e1 through eN=IP1,1[i=1NZ(Z+1i)]/MZ.

        We also had:

        Hremoving e1 from the neutral atom=(ΔIP1,ZΔIP1,Z1)(Z/MZ).

        We inferred in [25] that:

        Hremoving e1, & e2=[ΔIP1,ZZΔIP1,Z2(Z2)]/MZ.


        Hremoving e1, e2, & e3=[ΔIP1,ZZΔIP1,Z3(Z3)]/MZ.

        This pattern generalizes to:

        Hremoving e1 through eN=[ΔIP1,ZZΔIP1,ZN(ZN)]/MZ.

        Finally, in [25] we had:

        (W+H)removing e1, & e2=IP1,1[Z+Z(Z1)]/MZ+[ΔIP1,ZZΔIP1,Z2(Z2)]/MZ   ,


        (W+H)removing e1, e2, & e3=IP1,1[Z+Z(Z1)+Z(Z2)]/MZ+[ΔIP1,ZZΔIP1,Z3(Z3)]/MZ   ,

        This pattern generalizes to:

        (W+H)removing e1 through eN=IP1,1[i=1NZ(Z+1i)]/MZ+[ΔIP1,ZZΔIP1,ZN(ZN)]/MZ   ,

        Now let us turn to adding electrons. First, use the formula for the energy for removing an electron from a neutral atom of element Z to describe instead removing an electron from the singly charged negative ion of element Z, which has Z+1 electrons to start with:

        (W+H)removing e1 from negative ion=[IP1,1Z(Z+1)+ΔIP1,Z+1(Z+1)ΔIP1,ZZ]/MZ.

        Reversing the direction of the operation:

        (W+H)adding e1 to neutral atom=[IP1,1Z(Z+1)+ΔIP1,Z+1(Z+1)ΔIP1,ZZ]/MZ.

        This means the work for adding one electron into the nuclear field is:

        Wadding e1 to neutral atom =IP1,1Z(Z+1)/MZ.

        And the heat for re-adjusting the electron population is:

        Hadding e1 to neutral atom=[ΔIP1,Z+1(Z+1)+ΔIP1,ZZ]/MZ.

        Now let us add a second electron. This will require additional work:

        Wadding e2 after e1=IP1,1Z(Z+2)/MZ.

        And it will cause another heat adjustment:

        Hadding e2 after e1=ΔIP1,Z+2(Z+2)/MZ+ΔIP1,Z+1(Z+1)/MZ.

        This means total energy involved in adding two electrons is:

        (W+H)adding e1 & e2=IP1,1[Z(Z+1)+Z(Z+2)]/MZ+[ΔIP1,ZZΔIP1,Z+2(Z+2)]/MZ   .

        Likewise, the total energy involved in adding three electrons is:

        (W+H)adding e1, e2 & e3=IP1,1[Z(Z+1)+Z(Z+2)+Z(Z+3)]/MZ+[ΔIP1,ZZΔIP1,Z+3(Z+3)]/MZ   .

        This patten generalizes to:

        (W+H)adding e1 through eN=IP1,1[i=1NZ(Z+i)]/MZ+[ΔIP1,ZZΔIP1,Z+N(Z+N)]/MZ   ,

        Numerical Data to Insert in Formulae:

        Numerical data for elements up to number 118 are given in [24]. The numerical analysis of EZ water requires at most the data for the first ten elements. Expressed in electron volts, eV, these numerical data are:

        Hydrogen: Z=1 , M=1.008 , IP1,1=14.250eV , ΔIP1,1=0eV .
        Helium: Z=2 , M=4.003 , IP1,2=49.875eV , ΔIP1,2=35.625eV .
        Lithium: Z=3 , M=6.941 , IP1,3=12.469eV , ΔIP1,3=1.781eV .
        Beryllium: Z=4 , M=9.012 , IP1,4=23.327eV , ΔIP1,4=9.077eV .
        Boron: Z=5 , M=10.811 , IP1,5=17.055eV , ΔIP1,5=2.805eV .
        Carbon: Z=6 , M=12.011 , IP1,6=21.570eV , ΔIP1,6=7.320eV .
        Nitrogen: Z=7 , M=14.007 , IP1,7=27.281eV , ΔIP1,7=13.031eV .
        Oxygen: Z=8 , M=15.999 , IP1,8=27.281eV , ΔIP1,8=13.031eV .
        Fluorine: Z=9 , M=18.998 , IP1,9=34.504eV , ΔIP1,9=20.254eV .
        Neon: Z=10 , M=20.180 , IP1,10=43.641eV , ΔIP1,10=29.391eV .

        Ordinary Water:

        Here are some example calculations concerning possible ionic configurations of ordinary, normal water.

        Most people would guess that water is 2H++O2. But let us evaluate that ionic configuration. The transition HH+ takes:

          IP1,1/M1=14.250/1.008=14.1369 eV.

        So 2H+ takes    2×14.1369=28.2738 eV.

        The transition OO2 takes:

        IP1,1[Z(Z+1)+Z(Z+2)]/MZ+[ΔIP1,ZZΔIP1,Z+2(Z+2)]/MZ=14.250[8×9+8×10]/15.999+[ΔIP1,8×8ΔIP1,10×10]/15.999=14.250[72+80]/15.999+[13.031×829.391×10]/15.999=14.250[8.4853+8.9443]/15.999+[104.248293.910]/15.999=14.250[17.4296]/15.999+[189.662]/15.999=15.524211.8546=27.3788 eV   

        So the ionic configuration 2H++O2 requires 28.273827.3788=0.8950 eV. This is a positive energy requirement, which implies that some external assistance is needed to create this ionic configuration. So normal water may not be 2H++O2 after all!

        Another possibility is readily at hand though. The ionic configuration for normal water could be 2H+O2+. The transition HH takes:

        IP1,1Z(Z+1)/MZΔIP1,Z+1(Z+1)/MZ+ΔIP1,ZZ/MZ=IP1,11×2/M1ΔIP1,2×2/M1+ΔIP1,1×1/M1=14.250×1.4142/1.00835.625×2/1.008+0×1/1.008   =19.992470.6845+0=90.6769 eV   .

        So 2H takes 2×(90.6769)=181.3538 eV. Notice that this is a huge negative energy. It reflects the fact that electrons really like to make pairs. Indeed, their propensity to do so motivated the invention of the so-called ‘spin’ quantum number. Without spin, many electrons in atoms would be violating the ‘Pauli Exclusion Principle’, which says that only one electron can be in any particular quantum state. Electron pairs are famous in Condensed Matter Physics too, under the name ‘Cooper pairs’.

        Proceeding now to the transition OO2+, that takes:

        IP1,1[Z+Z(Z1)]/MZ+[ΔIP1,ZZΔIP1,Z2(Z2)]/MZ=IP1,1[8+8×7]/M8+[ΔIP1,8×8ΔIP1,6×6]/M8=14.250[8+56]/15.999+[13.031×87.320×6]/15.999=  13.7907+3.7707=17.5614 eV   .

        Thus the creation of the water molecule in the ionic configuration 2H+O2+ demands altogether 181.3538+17.5614= 163.7924 eV. This energy is solidly negative, which means that ordinary water is overwhelmingly in this ionic configuration, 2H+O2+. This is normal water.

        However, in situations where more than one version of anything can exist, both generally do exist, in proportions determined by their so-called Boltzmann factors, exp(E/kT). Here E is energy, k is Boltzmann’s constant, and T is absolute temperature. Boltzmann factors are the result of entropy maximization at work. Because of non-zero Boltzmann factors, there will exist a tiny, tiny fraction of the first ionic configuration, 2 H++O2.

        This analysis of normal water shows how quantitative approaches can sometimes unseat long-standing, but never-justified, assumptions in Chemistry.

        About EZ Water:

        The following transition is generally thought to represent the creation of EZ water:

        3 normal water molecules1 EZ water ion + 1 hydronium ionThat is,3(H2O)®(H3O2)-+(H3O)+.

        Here parentheses are used to avoid implying anything about what charge the individual atoms within any ion or radical may carry. A full numerical analysis should consider all possible, or at least all plausible, ionic configurations of every molecule or radical involved.

        One possible ionic configuration for the EZ water ion (H3O2) is 3H++2O2. The 3H+ takes total energy 3×IP1,1/Mz= 3×14.250/1.008=42.4107 eV. The 2O2 takes total energy 2×(27.3788)= 54.7576 eV. So the ionic configuration 3H++2O2 altogether takes 42.410754.7576= 12.3469 eV. This is a negative energy, so this ionic configuration certainly can occur.

        But there is also another possibility for the EZ water ion (H3O2). It could have the ionic configuration 3H+2O+. An H takes 90.6769 eV, so 3H takes 3×(90.6769)= 272.0307 eV. An O+ takes energy:

        (IP1,1×8+ΔIP1,8×8ΔIP1,7×7)/M8=(14.250×8+13.031×813.031×7)/15.999    =(114.000+104.24891.217)/15.999=7.9399 eV   .

        So 2O+ takes 2×7.9399=15.8798 eV. Then the ionic configuration 3H+2O+ takes 272.0307+15.8798= 256.1509 eV. This energy is much more negative than that of the first candidate ionic configuration for EZ water, 3H++2O2. This fact means EZ water is nearly always in this second candidate ionic configuration, 3H+2O+.

        One possible ionic configuration for the hydronium ion (H3O)+ is 3H++O2. This, I believe, is what most people would guess. But from the study of regular water, we know the candidate 3H+ would take 3×14.250=42.4107 eV, and that the candidate O2 would take 27.3788 eV, so the candidate ionic configuration 3H++O2 for hydronium would take 42.410727.3788=15.0319 eV. This energy is positive, so this ionic configuration for the hydronium ion is not promising.

        However, as was the case with the EZ water ion, there is another possibility for the hydronium ion. It (H3O)+ could have the ionic configuration 3H+O4+. The 3H would take 3×(90.6769)= 272.0307 eV, and the O4+ would take:

        IP1,1[Z+Z(Z1)+Z(Z2)+Z(Z3)]/MZ+[ΔIP1,ZZΔIP1,Z4(Z4)]/MZ=14.250×[8+56+48+40]/15.999+[13.031×89.077×4]/15.999=14.250×28.7361/15.999+67.940/15.999   = 25.5947+4.2465=29.8411 eV   .

        So for the hydronium ion (H3O)+, the second candidate ionic configuration 3H+O4+ would take 272.0307+29.8411= 242.1895 eV. This very negative energy explains why the reaction product that accompanies EZ water is a hydronium ion, rather than a naked proton plus a normal water molecule, which would take 14.1369 163.7925=149.6556eV, which is not as negative.

        The EZ water ion and the hydronium ion together take 268.6204242.1895=510.8099 eV. Compare this energy to the energy taken by three normal water molecules: 3×(163.7925)=491.3775 eV. The EZ water ion with the hydronium ion has lower energy than the three normal water molecules. That means that Nature will take any opportunity to make EZ water ions and hydronium ions.

        It appears that what creates the opportunity is a surface, plus a little energy to separate the ions. Any material body provides some gravity to create a surface, and if there is also some small energy source, such as sunlight, to help separate ions, and if there is also some normal water, the situation will automatically create EZ water too. Even an icy comet might be able to create some EZ water.

        7. Microphysics

        Just as the myriad compounds in Chemistry arise from not-very-many chemical elements, some significant part of the m


        Figure 4.

        Log energy loss and gain rates vs. log system radius.

        Observe that radiation dominates for radii below the crossing point, whereas torquing dominates for radii above the crossing point. This means the balance between the two effects is unstable: a small excursion from balance in either direction causes more excursion in the same direction. This is interesting. It means that Hydrogen does not like to exist as an isolated atom. It wants to engage in chemical reactions. In the Universe at large, you will find Hydrogen in H2 molecules, or other molecules, or you will find Hydrogen plasma, consisting of naked protons and free electrons, but you will not find many isolated Hydrogen atoms.

        1. We can discover even more about the Hydrogen atom if we include the appropriate angle sine and cosine factors in the calculations. These factors are oscillatory. So negative numbers sometimes occur with the torquing curve, and they cannot be plotted on a logarithmic scale. However, we are mainly interested in the points of balance between torquing and radiation, and the radiation curve is always non-negative, so the torquing curve is non-negative too at the balance points.

        Figure 5 shows that we have not just one balance point, but many balance points. The balance points occur in close pairs, one stable and one unstable. The two pairs furthest left on the plot do not show as crossings because of the finite resolution of the plot, but they are certainly present.

        More solution pairs are to be found, off the Figure to the left. Indeed, the solution pairs continue indefinitely, into smaller and smaller system radii. So Hydrogen has an infinite family of ‘sub-states’. Mills discusses these in [27].

        The smaller and smaller radii of the balance points in Fig. 5 correspond to higher and higher orbit speeds. This idea conflicts with a prohibition imposed by SRT: no physical particle possessing mass is allowed to move at a speed matching or exceeding light speed c. What does this conflict mean? I believe it means the prohibition should be understood more precisely to say: no physical particle possessing mass can be perceived to move at a speed matching or exceeding light speed c, if we agree to process all received data in accord with Einstein’s Second Postulate. If we do not agree to that, then there is no particle speed limit.


        Figure 5.

        More and more balance points below the Hydrogen ground state.

        1. We can also study positronium (a system consisting of one electron and one positron). See Fig. 6. As with Hydrogen, the oscillatory angle factors create a family of solutions for this system too. Half of them are stable, and half are unstable, and uncountably many of them occur at small radii and high speeds well in excess of c. Due to the finite resolution of the plot, only one pair of solutions is clearly visible in Fig. 6. But two more, at smaller radius and higher speed, are also certainly present. We can characterize these, and all high-speed solutions, without even knowing exactly what the radiation curve is like - its exact amplitude, or its r4 dependence. The one low-speed solution is just v00.02×c. The many high-speed solutions have to occur in pairs just above and below orbit speeds of the form vn=v0+n×2πc, where n is an arbitrary positive integer.


        Figure 6.

        Positronium solutions.

        A parenthetical note applies for Fig. 6: for same-mass systems, the amplitude of the radiation curve should be less by a factor of 4 because there is no center-of-mass motion. This sort of numerical detail does not significantly affect where the solutions fall. That is determined almost entirely by the cosine factors that produce the deep dips that intersect the peaks in the curve for rate of energy gain due torquing.

        1. The oscillatory nature of the angle factors can turn a situation of seeming repulsion into a situation of actual attraction. This phenomenon of sign reversal due to signal delay is well known to engineers, who often deal with oscillating signals in feedback control systems For the present application, consider two electrons in a circular orbit, and suppose they move at speed πc. One electron launches its signal radially outward. By the time this electron has executed half an orbit, this signal has expanded a distance equal to the orbit diameter. By then, the two electrons have exchanged places. So the expanding signal first contacts the second electron at exactly the signal launch point. Then the two electrons complete their orbit. At the end, the second electron finally understands its signal: it is to move radially. But by now, the two electrons have changed places again, and for the second electron, the direction commanded is inward. That situation is equivalent to attraction.

        Given this mechanism for attraction, we can also study homogeneous systems: two electrons, or two positrons, for example. Again, there exist both stable and unstable solutions, and there are infinitely many of each, corresponding to orbit speeds of the form vn=πc+n×2πc for arbitrary positive integer n. See Fig. 7.


        Figure 7.

        Same-charge solutions.

        An additional parenthetical note applies for Fig. 7: for same-charge systems, the angular pattern of the radiation is quadrupole, rather than dipole, so the amplitude of the radiation energy loss curve should decline as r6 rather than r4. Again, this sort of numerical detail does not significantly affect where the solutions fall, which is determined almost entirely by the cosine factors that produce the deep dips that intersect the peaks in the curve for rate of energy gain due torquing.

        The stable solutions for two electrons bring to mind the situation that is so well known in Chemistry: electron pairs. They are everywhere in Chemistry. The most famous case occurs for Helium. Helium is a noble gas, and it reacts with other elements only under extreme duress. Helium has two electrons, and pulling one electron away is very costly: Helium has the highest ionization potential of any element. The message is: two electrons definitely do form a stable subsystem within an atom. The standard QM explanation for this invokes the concept of electron spin, with two possible values, ±/2, allowing two electrons in the same overall energy state. Electron pairing also occurs famously in Solid State Physics, under the name of Cooper Pairs.

        8. Conclusion

        In the present Chapter, the new concept applied is the more realistic signal model for use in an improved version of SRT. The realistic signal model is based on Information Theory,

        The new concept is implemented with very standard mathematics: differential equations, their family of solutions, and the particular problem boundary conditions. These mathematical ingredients for a proper signal model were all available in 1905, but they were not used in SRT. Why? I believe the fundamental reason is the history: Information Theory was not yet available, so no researcher at that time would have been likely to detect the inadequacy of the infinite plane wave as a signal model.

        The present paper has shown that there are rewards for instead using the realistic photon/signal model. They include more insight into Quantum Mechanics, and into Gravity Theory, and potentially into Elementary Particle Physics. These are all subjects to be studied much more fully in the future.

        Many textbook treatments of SRT devote a lot of space to Lorentz Transformations (LT’s). The present work has not mentioned LT’s at all. To this author, LT’s just seem to describe the wrongly informed opinions of different observers. So I don’t really want to focus on LT’s. But I have to mention them, because repairing SRT to take proper account of the concepts of IT casts doubt on Einstein’s SRT, and hence on LT’s. Therefore, I hereby relegate the unavoidable discussion of coordinate transformations, LT’s and others, to the following Appendix.

        9. Appendix

        The situation in the late nineteenth century included the following fact: Maxwell’s first order coupled field equations appeared not to be invariant under Galilean transformation of coordinates (GT’s). Phipps [28] has written extensively about this apparent conflict between Maxwell’s Electromagnetic Theory and Newton’s Mechanics. In the early twentieth century, SRT brought in LT’s, and the conflict seemed to be resolved: Maxwell’s electromagnetic theory was clearly invariant under LT’s. This fact was taken as evidence in favor of Einstein’s SRT over Newton’s Mechanics.

        But there is a puzzle left to resolve: Maxwell’s first order coupled field equations appear to qualify as tensor equations. Mathematicians had developed tensors in the first place to enable the articulation of mathematical statements that would be coordinate-free. So tensor equations are by definition invariant to all invertible coordinate transformations.

        So what had happened here? I believe two circumstances had collided to create a very bad situation. One circumstance was that Mathematics had such a long history of developing eternal truths: the focus had been on arithmetic, geometry, and trigonometry – all of them eternal in character. Even archeo-astronomy was largely about the eternal repetition of events, and not about temporal evolution of events. Eternal truths really need not have a time dimension. They can, however, have as many spatial dimensions as may be desired, and that became the focus for much of tensor analysis. The other circumstance was that time became a really significant variable with the advent of modern Physics: Kepler, Galileo, Newton, and Maxwell. And time is a really different kind of variable than space is. Maxwell was very well aware of the difference, as he developed his electromagnetic theory in terms of Hamilton’s quaternions. The modern equivalent of the quaternion tool is the set of four 2×2 complex Pauli spin matrices:


        The first one, the time-like one, is the identity matrix. The other three, the space-like ones, produce the identity matrix when squared. When two of them are cross-multiplied, they generate a factor of i=1 times the third one, corresponding to the vector cross product in three-dimensional space.

        The collision between circumstances came in the formulation of differential operators. People were familiar with the scalar chain rule, d/dt =/t+v/x, and did not realize that more information was needed in the context of vector and tensor applications with time as well as space dimensions.

        The 2×2 Pauli matrices are easy to appreciate visually, so I will also discuss transformations of coordinates in terms of 2×2 matrices also - but only real ones, not complex ones. Let s stand for any spatial coordinate. A general coordinate transformation involving ct and s has the form:

        For the familiar Lorentz transformation, A=B=v/c, where v is the speed of the new coordinate frame relative to the old one. The letter v is lower case to remind us that v<c ; i.e. v/c<1. We have:


        For the long discarded Galilean transformation, A=V/c and B=0. The letter V is upper case to remind us that V is not limited, and might exceed c. So we have:


        For all such general coordinate transformations, there also exists a complement transformation:


        Its purpose is to preserve inner products; for example:


        Observe that:

        [ct's'][ct"s"]=[cts]11AB[1AB1]11AB[1AB1][cts]=[cts]11AB[1AB00BA+1][cts](ct)2s2   .

        For Lorentz transformation, the complement transformation is the inverse, or equivalently, the reverse transformation:


        But for Galilean transformation, the complement transformation is:


        This is the inverse transformation, but not the reverse transformation. It is rather the transpose of the reverse transformation. It looks so very strange because, for more than a century now, only Lorentz transformations of velocity have been used in mainstream theoretical Physics, and transposition does not change them.

        The vital role for this strange new thing lies with the differential operators. The story is much like it was for the coordinates: there are two complementing transformations, and they involve, not only inversion/reversal, but also transposition. Let ct represent differentiation with respect to the time-like coordinate, and s represent differentiation with respect to the spatial variable. Let us demand invariance of inner products involving differential operators; for example, like:


        i.e., always two statements – not just one statement. This level of detail was missing from the scalar chain rule, and that omission caused people to believe that Maxwell’s equations could not be shown to be invariant under GT. And so they welcomed LT instead. This is not to say we should now revert to using GT again. Indeed, because of the half-retardation issue discussed in Sect. 3, the best transformation to use may involve, not V/c, but rather V/2c. This question needs detailed future study.

        The use of 2×2 matrices can make the detail needed in such future study very clear. However, many mathematicians tend to prefer tensor notation. But current-day tensor notation uses only two index positions, both on the right: down called ‘covariant’, up, called ‘contravariant’. To represent the transformations needed for Physics, it would be helpful, and maybe necessary, to add two more index locations, up and down on the left, to acknowledge transposition, and using words like ‘trans-covariant’ and ‘trans-contravariant’ to emphasize what putting indices in those positions means.


        This Chapter is dedicated to the memory of a most courageous researcher in theoretical and applied electrodynamics: Dr. Peter Graneau, 1921-2014. He encouraged me, and many other researchers, to give serious attention to the History behind Physics.


        1 - Einstein, A., On Thr Electrodynamics of Moving Bodies, Annelen der Physik 17, 891-921 (1905).
        2 - Einstein, A., Annelen der Physik On the Theory of Light Production and Light Absorption, Annelen der Physik 20, 199-206 (1906).
        3 - Brillouin, L., Science and Information Theory, Second Edition, Dover Publications, Inc., Mineola, NY, (2013).
        4 - Flores-Gallegos, N., Shannon Informational Entropies and Chemical Reactivity, Chapt. 29, pp. 683-722 in Advances in Quantum Mechanics, Intech (2013).
        5 - Jackson, J.D., Classical Electrodynmics, Second Edition, John Wiley & Sons, New York, NY. (1975).
        6 - Bentwich, J., The Theoretical Ramifications of the Computational Unified Field Theory, Chapt. 28, pp. 671-681, Advances in Quantum Mechanics, Intech (2013).
        7 - Whitney, C.K., Better Unification for Physics in General Through Quantum Mechanics in Particular, Chapter 7, pp. 127-160 in Theoretical Concepts of Quantum Mechanics, Ed. M.R. Pahlavani, InTech. (2012).
        8 - Ritz, W., Researches critiques sur l’electrodynamique generale, Ann. Chim. et Phys. 13, pp. 145-275, (1908).
        9 - Putz, M.V., Path Integrals for Electronic Densities, Reactivity Indices, and Localization Functions in Quantum Systems, International Journal of Molecular Sciences 10, pp. 4816-4940 (2009).
        10 - Bracken, P., Quantum Mechanics Entropy and a Quantum Version of the H-Theorem, Chapt. 21, pp. 469-488 in Theoretical Concepts of Quamtum Mechanics 469-488 (2012).
        11 - Bracken, P., The Schwinger Action Principle and its Applications to Quantum Mechanics, Chapt. 8, pp 159-182 in Advances in Quantum Mechanics, Intech (2013).
        12 - Spitnev, V., Generalized Path Intergral Technique: Nanoparticles Incident on a Slit Grating, Matter Wave Interference, Chapt 9., pp. 183-212 in Advances in Quantum Mechanics, Intech (2013).
        13 - Bulnes, F., Correction, Restoration and Re-Composition if Quantum Mechanical Fields of Particles by Path Integrals and Their Applications, 489-514 in Advances in Quantum Mechanics, Intech (2013).
        14 - Liénard, A., Champ Electrique et Magnétique produit par une Charge Electrique Concentrée en un Point et Animée d’un Movement Quelconque, L’Eclairage Electrique, XVI, pp. 5-14, 53-59, and 106-112 (1898).
        15 - Wiechert, E., Elektrodynamische Elementargesetze, Archives Néerlandesises des Sciences Exactes et Naturelles, série II, Tome IV, pp. 549 – 573 (1901).
        16 - Lokajicek, M., V. Kundrát, and J. Procházka, Schrödinger Equation and (Future) Quantum Physics, Chapter 6 in Advances in Quantum Mechanics, Ed. Paul Bracken, InTech. (2013).
        17 - Whitney, C.K., On the Dual Concepts of ‘Quantum State’ and ‘Quantum Process’, Chapter 17 in Advances in Quantum Mechanics, Ed. Paul Bracken, InTech. (2013).
        18 - Enders, P., Quantization as Selection Rather than Eigenvalue Problem, Chapt 23, pp. 543-564, Advances in Quantum Mechanics, Intech (2013).
        19 - Ralston, J.P., Emergent un-Qusntum Mechanics, Chapt. 19, pp. 437-475, in Advances in Quantum Mechanics, Intech (2013)
        20 - Whitney, C.K., Reasonable doubt: Cosmology’s gift for physics, Journal of Computational Methods in Science and Engineering 2013, 13, pp. 291-302 (2013).
        21 - Graneau, P., Ampère–Neumann Electrodynamics of Metals, Hadronic Press, Nonantum, MA, USA, pp. 7 – 22 (1987).
        22 - MacGregor, M., The Enigmatic Electron- A Doorway to Particle Masses, 2nd ed; El Mac Books, Santa Cruz, CA, USA; pp. 3 –11 (2013).
        23 - Pollack, G., The Fourth Phase of Water, Beyond Solid Liquid Vapor; Ebner & Sons, Seattle (2013).
        24 - Whitney, C.K., Algebraic Chemistry: Applications and Origins, Nova Science Publishers (2013).
        25 - Whitney, C.K., EZ Water: A New Quantitative Approach Applied – Many Results Obtained, Water 5, January 5, pp. 105-120 (2014).
        26 - Whitney, C.K., The Neutron: a Challenge for Post-Maxwell Physics, Hadronic Journal, 3, pp. 201 – 250 (2008).
        27 - Mills, R.L., The Nature of the Chemical Bond Revisited and an Alternative Maxwellian Approach, Physics Essays 17, 3, pp. 342 – 389 (2004).
        28 - Phipps, T.E., Heretical Verities: Mathematical Themes in Physical Description; classic non-fiction library, Urbana, IL, USA, pp. 102 – 108 (1986).