Quantum gravity (Planck)

Emulating gravitational and atomic orbits via n-body discrete rotating particle-particle orbital pairs at the Planck scale

This orbit model uses a geometrical approach instead of the dimensioned physical constants (G, h, c) to emulate gravitational and atomic orbits via an n-body network of rotating particle-to-particle orbital pairs. Each particle in the 'universe' is connected to every other particle by a discrete circular rotating orbital (representing a unit of momentum) with the particles at each orbital pole, thereby forming a universe-wide network of particle-to-particle orbital pairs [1].

Although dimensionless (the model does not use formulas that require dimensioned constants), orbit period and radius can be measured in Planck units, and so this approach is applicable to modelling gravitational orbitals at the Planck scale. In dimensioned terms, per unit of Planck time tp each orbital pair rotates by 1 unit of Planck length lp at velocity c (v = lp/tp) in 4-axis hypersphere co-ordinates, all orbitals are then summed and averaged to give the new particle co-ordinates. The process then repeats, gravitational orbits between macro-bodies in 3-D space emerging as the averaged (over time) sum of these individual rotating orbitals. The particles (of the orbiting bodies) are connected with each other directly via these orbitals, information regarding the macro orbiting bodies or associated barycenter(s) is not required.

In a simple orbit where a small mass (such as an electron or satellite) is orbiting a larger mass (nucleus, planet ...), the radius of orbit follows the formula

To emulate gravitational and atomic orbit(al)s, the (inverse) fine-structure constant alpha is used. Although the orbit itself occurs in the hypersphere at the speed of light, this orbital constant returns period and velocity for a 2-D plane (of 3-D space).

Atomic orbitals are a geometrical subset of these gravitational orbitals, and so likewise can be interpreted as physical units of momentum rather than regions of probability where the particle can be found. They are treated as analogous in physical characteristics to the photon, with electron transition between orbitals occurring when the orbital itself absorbs (or ejects) a photon. The electron has no direct role in the transition.

Although simply the sum of rotating circles, the equations that emerge to describe simple orbits resemble the Bohr model.

By simply selecting the start co-ordinates on a 2-D plane for each point (unit of mass) accordingly, we can 'design' the required orbits. No other parameters are required.

8-body (8 mass points, 28 orbitals), the resulting orbit is a function of the start positions of each point

Particle-particle orbitals Edit

Particles are assigned `point' (representing a discrete unit of mass) co-ordinates within a 4-axis hyper-sphere `universe'. Every particle is then connected to every other particle by a circular orbital. The hyper-sphere expands in incremental steps (FOR age = 1 TO ... representing a discrete unit of time), this expansion causes the orbitals to rotate in steps (representing a discrete unit of length) correspondingly.

In a simulation, for each value of incrementing variable age (the simulation clock-rate), all (n-body) orbitals rotate by 1 step (driven by the hypersphere expansion) and the co-ordinates calculated. These are then summed and averaged giving new particle co-ordinates. As this occurs for every particle before the next increment to age, the process can be updated in 'real time' on a serial processor. As the orbitals are circular, the barycenter for each orbital is its center, the particles at each orbital 'pole'.

As these steps are discrete, they may be measured in terms of the Planck units, each point as a unit of Planck mass, with the orbitals rotating 1 unit of Planck length lp per 1 unit of Planck time tp (i.e.: 1 increment to age) at velocity c = lp/tp in hyper-sphere coordinates.

Although orbital and so particle motion occurs at c, the hyper-sphere expansion [2] is equidistant and so `invisible' to the observer. Instead observers will register these 4-axis orbits (in hyper-sphere co-ordinates) as a circular motion on a 2-D plane (in 3-D space). An apparent time dilation effect emerges as a consequence.

N-body orbit simulation Edit

A simulation program [3] comprising n-body rotating orbitals is described. Each particle in the simulation is assigned initial (x, y) 2-D point co-ordinates (representing 3-D space), forming orbital pairs that rotate around each other on a 2-D plane according to an angle β as defined by the orbital pair radius.


The rationale for angle β is introduced in the section on atomic orbitals (below).

The total distance travelled, 1 unit of length per increment to age (1 time unit) is given in (x, y, z) co-ordinates, where the (z) axis represents the hypersphere expansion axis.

As the simulation treats each (particle-particle) orbital independently (independent of all other orbitals), no information regarding the points (other than their initial start coordinates) is required by the simulation.

Fixed radius, variable mass Edit

For the following simple orbits, 1 point is assigned as the orbiting point, the remaining points forming the 'central' mass. The only distinction being that the central mass points are assigned (x, y) co-ordinates relatively close to each other, and the orbiting point is assigned (x, y) co-ordinates distant from the central points (this becomes the orbital radius). The simulation however treats all points equally, the center points also orbiting each other according to their orbital radius (in the simulation there is no mass greater than 1 unit).

After every orbital has rotated 1 length unit anti-clockwise, the new co-ordinates for each rotation per point are then averaged and summed, the process then repeats. After 1 complete orbit (return to the start position by the orbiting point), the period tsim (as the number of increments to the simulation clock age) and the (x, y) plane orbit length lsim are noted.


1. i; number of 'physical' center points in the orbit (the center mass).

2. j = i*x + 1; number of virtual center points (to reduce computation time, i*x virtual points are added to increase center mass up to j = jmax

3. jmax; maximum number of mass points per orbital radius

4. x, y; start co-ordinates for each point (2-D plane).

5. rα; a radius constant, here rα = sqrt(2α) = 16.55512; where alpha = inverse fine structure constant = 137.035 999 084 (CODATA 2018).

Symmetrical 4 body orbit; (3 center mass points, 1 orbiting point, 6 orbital pairs). Note that all points orbit each other.

Example: i = 81, j = jmax = 32*81+1 = 2593 (3321 orbitals)


tsim = 58430803.8405

lsim = 3528109.1197


This gives the following equations for the 2-D plane

  ratio of mass to maximum mass per orbital radius
 , orbital radius.
 , barycenter
 , orbiting point velocity
 , orbited point(s) velocity
 , orbiting point period
 , distance travelled

Fixed mass, variable radius Edit

To model a 1kg satellite to earth orbit will require earth mass/Planck mass = 0.2744 x1033 points and 1kg/Planck mass = 45940510 points. We can reduce calculation by using only relative mass and then use the dimensionless ng to assign the start parameters. For example, from the standard gravitational parameters, the earth to moon mass ratio approximates 81:1.


There is 1 orbiting point (distant point) and 81 central points (points in close vicinity)


To calculate ng

  = 384400km
  = 0.00887m (Schwarzschild radius)
  = 17783.25

This gives

  = 1445178.5

Converting from dimensionless numbers to SI Planck units using lp and c;


We can use the actual radius and period to translate between values.

  = 27.322 days

The above assumes a circular orbit, to form an elliptical orbital we can use unaligned orbitals (PE vs KE).

Gravitational coupling constant Edit

In the above, particles were assigned a mass as a theoretical unit of Planck mass (a point). Conventionally, the Gravitational coupling constant αG characterizes the gravitational attraction between a given pair of elementary particles in terms of a particle (i.e.: electron) mass to Planck mass ratio;


For the purposes of this simulation, particles are treated as an oscillation between an electric wave-state (duration particle frequency) and a mass point-state (duration 1 unit of Planck time). The above value αG then represents the probability that 2 electrons will be in the mass point-state at any unit of Planck time (wave-particle duality at the Planck level represented by an electric-wave to mass-point oscillation [4]).

As mass is not treated as a constant property of the particle, measured particle mass becomes the averaged frequency of discrete point mass at the Planck level. If 2 dice are thrown simultaneously and a win is 2 'sixes', then approximately every 36 throws (frequency) of the dice will result in a win. The inverse of αG is the frequency of occurence of the mass point-state between the 2 electrons. As 1 second is 1042 units of Planck time, this occurs about once a minute. Gravity now has a similar magnitude to the strong force (at this, the Planck level), albeit this interaction occurs seldom (the Planck level), and so when averaged over time (the macro level), gravity appears weak.

If particles oscillate between an electric wave-state to Planck-mass (for 1 unit of Planck-time) point-state, then at any discrete unit of Planck time a number of particles in the universe will simultaneously be in the mass point-state. For example a 1kg satellite orbits the earth, for any given unit of time, satellite (B) will have   particles in the point-state. The earth (A) will have   particles in the point-state, and so the number of orbital links (the gravitational coupling constant) between the earth and the satellite will sum to the number of orbitals;



1. 1kg satellite at a synchronous orbit radius

  (earth as the center mass)
  (Schwarzschild radius)
  (synchronous orbit)

2. The energy required to lift a 1 kg satellite into geosynchronous orbit is the difference between the energy of each of the 2 orbits (geosynchronous and earth).

  (energy per orbital)
  (number of orbitals)

3. The orbital angular momentum of the planets derived from the angular momentum of the respective orbital pairs.


The orbital angular momentum of the planets;

mercury = .9153 x1039  
venus    = .1844 x1041  
earth    = .2662 x1041 
mars     = .3530 x1040 
jupiter   = .1929 x1044   
pluto   = .365 x1039

Orbital angular momentum combined with orbit velocity cancels ng giving an orbit constant. Adding momentum to an orbit will therefore result in a greater distance of separation and a corresponding reduction in orbit velocity accordingly.


32 mass points (496 orbitals) begin with random co-ordinates, after 232 steps they have clumped to form 1 large mass and 2 orbiting masses.

Freely moving points Edit

The simulation calculates each point as if freely moving in space, and so is useful with 'dust' clouds where the freedom of movement is not restricted (i.e.: in the above example, the earth particles do not follow gravitational orbits around each other). When measuring the orbit of a single point around a larger mass, after each complete orbit we can note that the orbit period and radius reduces (as a function of center mass and start radius distance).

In this animation, 32 mass points begin with random co-ordinates (the only input parameter here are the start (x, y) coordinates of each point; i, j, r ... are not preset). We then fast-forward 232 steps to see that the points have now clumped to form 1 larger mass and 2 orbiting masses. The larger center mass is then zoomed in on to show the component points are still orbiting each other, there are still 32 freely orbiting points, only the proximity between them has changed.

8-body circular orbit plus 1-body with opposing orbitals 1:2

PE vs. KE (opposing orbitals) Edit

Gravitational potential and kinetic energy are measures of alignment of the orbitals. In the above examples, all orbitals rotate in the same direction = kinetic energy. If all orbitals are unaligned the object will appear to 'fall' = potential energy.

In this example, for comparison, onto an 8-body orbit (blue circle orbiting the center mass green circle), is imposed a single point (yellow dot) with a ratio of 1 orbital (anti-clockwise around the center mass) to 2 orbitals (clockwise around the center mass) giving an elliptical orbit.

The change in orbit velocity (acceleration towards the center and deceleration from the center) derives automatically from the change in the orbital radius.

The orbital drift (as determined where the blue and yellow meet) is due to these orbiting points rotating around each other.

Precession Edit

semi-minor axis:  

semi-major axis:  

radius of curvature : 


arc secs per 100 years (drift):

  = 365.25 days

drift =  

Mercury (eccentricity = 0.205630)
T = 87.9691 days
a = 57909050 km (n = 378.2734) 
b = 56671523 km (l = 374.2096)
drift = 42.98
Venus (eccentricity = 0.006772) 
T = 224.701 days
a = 108208000 km (n = 517.085) 
b = 108205519 km (l = 517.079)
drift = 8.6247
Earth (eccentricity = 0.0167)
T = 365.25 days
a = 149598000 km (n = 607.989) 
b = 149577138 km (l = 607.946)
drift = 3.8388
Mars (eccentricity = 0.0934)
T = 686.980 days
a = 227939366 km (n = 750.485) 
b = 226942967 km (l = 748.843)
drift = 1.351

Illustration of B's cylindrical orbit relative to A's time-line axis

Hyper-sphere orbit Edit

An expanding hyper-sphere forms the scaffolding of the `universe'. The hyper-sphere expands in uniform incremental steps (the simulation clock-rate) as the origin of the speed of light, and so (hyper-sphere) time and velocity are constants. Particles are pulled along by this expansion, the expansion as the origin of motion, and so all objects, including orbiting objects, travel at, and only at, the speed of light in these hyper-sphere co-ordinates [5]. Time becomes time-line.

While B (satellite) has a circular orbit period on a 2-axis plane (the horizontal axis representing 3-D space) around A (planet), it also follows a cylindrical orbit (from B1 to B11) around the A time-line (vertical expansion) axis (td) in hyper-sphere co-ordinates. A is moving with the universe expansion (along the time-line axis) at (v = c), but is stationary in 3-D space (v = 0). B is orbiting A at (v = c), but the time-line axis motion is equivalent (and so `invisible') to both A and B, as a result the orbital period and velocity measures will be defined in terms of 3-D space co-ordinates by observers on A and B. In dimensionless terms;


For object B


For object A


Planck force Edit






Atomic orbitals Edit

The atomic orbital in this model[6] is a specific case of the gravitational orbital. It is also treated as a distinct unit of rotational momentum (rather than a region of probability), and has physical similarities to the photon (during orbital transition, it is the orbital radius which absorbs/ejects the photon thereby lengthening or shortening, the electron itself has a passive role). The orbital dimensions are a function of the fine structure constant alpha and the electron-nucleus wavelength. The gravitational orbital simulation can be applied to atomic orbitals, the difference being the angle of rotation β (into this angle is included an extra alpha term).

electron (blue dot) moves 1 step anti-clockwise along the alpha orbital circumference

Base orbital Edit

The basic alpha orbital is a circular orbital with a radius 2α. In contrast to the gravitational orbital, the atomic orbital has an additional   component (note: in this atomic orbital model the n variable is analogous to the principal quantum number).


The value for CODATA 2018 (inverse) alpha α = 137.035999084, thus the alpha orbital (where n = 1), will require about tref = 2π4α2 ~ 471964 steps to make 1 complete rotation. To reduce computational time, the electron-nucleus frequency component f of the orbital radius is added later. For a simple Hydrogen atom, f is taken from the sum of the wavelengths of the electron + proton (λH = λe + λp).

The frequency component f determines the duration of the alpha orbital step, each step corresponding to 1 frequency cycle, if we use this λH example, then the frequency is about fH = 0.11954 x1023 tp (units of Planck time), thus for the duration of 1 step;

time required = 0.11954 x1023 tp
distance travelled = wavelength/(2α n)

however for computational reasons the simulation calculates only the alpha orbital component.

orbital transition in steps

Orbital transition Edit

The photon is divided into segments  .


Photon absorption:

A photon will be absorbed by the orbital radius such that after each transition step the orbital radius will be reduced by rincr (as rincr has a minus value).


Conversely, by ejecting a photon the orbital radius will increase.

However an incoming photon is actually 2 photons as per the Rydberg formula.


As   is minus, the '(+ )' photon will decrease per transition step the orbital radius in   increments as above, however the '- (+ )', because of the Rydberg minus term, will conversely increase the orbital radius. And so for the duration of the (+ ) photon wavelength, the orbital radius does not change as the 2 photons cancel;


However the (+ ) has the longer wavelength, and so after the '(+ )' photon has been absorbed, and for the remaining duration of the (+ ) wavelength, the orbital radius will be extended (  = +0.000581).


As with gravitational orbitals, the velocity of rotation derives from the radius of the orbital (from  ) and so adjusts as the orbital radius changes. The simulation calculates the velocity per step to determine the final transition velocity component.

We may also note that of the base orbital radius r, although r = 2α, it is a construct of these   segments, as we transition from r = 0 to r = 2α in   increments, we require   steps.

Relativistic orbital Edit

If we include the frequency component f in the orbital radius (which is not practical for a simulation where f ~ 1023), we find that the electron travels 1 Planck length per unit of Planck time in hypersphere coordinates, as with gravitational orbitals.

H atom Edit

H ( ) = 2466061413.187GHz is the most precisely measured of the spectra and can be used as a reference.


Including the hypersphere relativistic term, the transition from an ninitial to the higher nfinal orbit by photon absorption;


H atom orbital transitions from n1-n2, n2-n3, n3-n1 via 2 photon capture (photons expand (contract) the orbital radius, the electron continues its anti-clockwise rotation.

In the animation (right), the electron is moving in an anti-clockwise direction, it is the orbital radius that is extending (or contracting) in   steps, the electron has a passive role. Note that we are calculating only for a 2-D plane (whereas the electron travels in 3-D space), furthermore there are no assumptions made regarding complex geometries and so we can only suggest models. For example, we can compensate (and so improve accuracy in line with experimental values) by truncating the base orbital period from   to  . Here we set   = 471961.8898 for the H atom ni=1 to nf=2 transition using the (CODATA 2018 [7]) wavelengths for the electron and proton. There are other approachs, the results are similiar.

electron; 2.42631023867e-12m

proton; 1.32140985539e-15m

to = 471961.8898 gives H(1s-2s) 2466061413MHz

Optimal values for these transitions;

  = 471960.1758 gives H(1s-3s) 2922743278MHz

  = 471959.8404 gives H(1s-4s) 3082581563MHz

We can extrapolate to get approximate values

  = 471959.715 for H(1s-5s) 3156563436MHz

  = 471959.626 for H(1s-6s) 3196751153MHz

The theorectical value for ionization of the H electron;

  = 471959.242135 for 3288086800MHz

One possible explanation could involve the Doppler effect as each ejected photon carries away different 'quantities' of momentum with it. If we set tref = 2π2α2α as the reference = 1.0 then the shift per photon would correspond to

H(1s-2s) = 1.002628

H(1s-3s) = 1.002980

H(1s-4s) = 1.003098

Ionization = 1.003297

In summary, photon #1 strikes and deletes a section of the orbital radius whilst photon #2 adds to the orbital radius (torbital), as a consequence the orbital radius does not change (the electron orbit continues as per normal). After photon #1 has been absorbed, the radius will begin expanding until photon #2 is also completely absorbed (ttransition). The electron is now in the new orbit(al).

Thus the quantum nature of the orbitals also derives from the photon wavelengths, this suggests that the photon wavelengths are not independent of the electron wavelength.

He dissociation of electron 1 (red) then electron 2 (blue) transition n=1 to n=2

Helium Edit

Higher atoms have more complex geometries but the principal remains the same. In this He animation, both orbiting electrons occupy the same orbital radius. As electron 1 (red) is ionized (absorbing momentum), electron 2 (blue) simultaneously drops to a lower orbital radius (emitting momentum), thus electron 2 subsidizes the ionization of electron 1.

After dissociation of the red electron, the blue electron then transitions from its new n= 1 to an n= 2 orbital. The axes measure only the alpha orbital component.

gravitational orbit with mass and distance ratio equivalent to diatomic H

Diatomic H Edit

Diatomic Hydrogen radius = 37pm. The H Bohr radius was set above at 2α * (λe + λp) = 105.89pm. To simulate as a 'gravitational' orbit using only an anti-clockwise rotation with no allowance for charge, we set;

electrons; mass = 1 point, start co-ordinates (-99, 0) and (99, 0)

protons; mass = 1836 points, (0, 37) and (0, -37)

orbit center = (0, 0)

number of point to point orbitals = 6747301

... thereby setting the distance from each electron to each proton = 105.89 respectively and electron to electron at 2*99.46. The H2 ionization energy (15.426eV) is 1.1344x greater than for the H atom (13.59844eV). Likewise combining the 2 electron-electron radius (105.89/99.46 = 1.06727) gives 1.13454.

If we reduce proton-proton separation, the protons act as a single center mass and the electrons follow a circular orbit. By increasing the proton-proton separation, the electron orbits increase proportionately. This separation distance (74pm) gives the optimal solution for a gravitational orbit.

External links Edit

References Edit

  1. Macleod, Malcolm J.; "3. Emulating gravity via n-body rotating particle-particle orbital pairs at the Planck scale". RG. Feb 2011. doi:10.13140/RG.2.2.11496.93445/15. 
  2. Macleod, Malcolm; "2. Programming cosmic microwave background for Planck unit Simulation Hypothesis modeling". RG. 26 March 2020. doi:10.13140/RG.2.2.31308.16004/7. 
  3. Macleod, Malcolm J.; "3. Quantum gravity via n-body particle-particle orbitals in a Planck scale simulation". RG. Feb 2011. doi:10.13140/RG.2.2.11496.93445/15. 
  4. Macleod, M.J. "Programming Planck units from a mathematical electron; a Simulation Hypothesis". Eur. Phys. J. Plus 113: 278. 22 March 2018. doi:10.1140/epjp/i2018-12094-x. 
  5. Macleod, Malcolm; "1. Programming relativity for Planck unit Simulation Hypothesis modeling". RG. 26 March 2020. doi:10.13140/RG.2.2.18574.00326/3. 
  6. Macleod, Malcolm J.; "4. Atomic orbitals and fine structure constant alpha in Planck scale simulations". RG. Feb 2011. doi:10.13140/RG.2.2.23106.71367/5. 
  7. https://physics.nist.gov/cgi-bin/cuu/Value?ecomwl