Life's Source

Web Particle Torsion Planets

Origin of the Second

The hour was first invented in ancient Egypt by dividing the night and day into 24 units, 12 for the day and 12 for the night. Since the day is longer in the summer, and the night shorter, and in the winter the day is shorter and the night is longer, the length of an hour depends on the season. The ancient Greek astronomer, Hipparchus, divided the day and night into hours determined by the length of day and night during spring and fall equinoxes when the length of day equals the length of night, inventing the equinoctial hour used year round. Hipparchus had access to ancient Babylonian knowledge of celestial motions where they knew the day of 24 hours gave an hour that could be divided by 60 minutes, and each minute by 60 seconds. The Babylonians got the base 60 divisions of the hour from the ancient Sumerians. But passage of time wasn't measured down to the second until Christiaan Huygens invented his pendulum clock, which was demanded by the astronomical revolution that came about from the work of Copernicus (Earth moves around the Sun), Galileo (Earth is not at the center of the Universe from looking at Jupiter's moons with his telescope), Brahe (data for planetary motions), Kepler (explains Brahe's data introducing elliptical orbits for the planets), and Newton (explains Kepler's laws of planetary motion with his universal law of gravitation).


A Universal Particle Equation

Ian Beardsley

April 11 – June 12, 2026

Abstract

We present a universal particle equation where what we experience as mass is taken as resistance to changes in a particle's motion through the temporal dimensions, which is measured by G, the universal constant of gravitation. To do this we introduce a normal force given by \( F_n = h/(c t_1^2) \) where \( t_1 \) is on the order of \( t_1 = 1 \) second, which is Lorentz invariant. The normal force, \( F_n \), is exposed to the cross-sectional area of the particle \( A_i = \pi r_i^2 \). The result is the mass of the particle is given by \( m_i = \kappa_i \sqrt{\pi r_i^2 F_n / G} \), with experimental verification giving 1.00500 seconds (proton), 1.00478 seconds (neutron), and 0.99773 seconds (electron). The coupling constant, \( \kappa_p \), is predicted by a prediction for the radius of the proton, which is \( r_p = \phi h/(c m_p) \) with \( 1/\phi = \Phi \) where \( \Phi = (\sqrt{5}+1)/2 \) is the golden ratio, and \( \kappa_i \) in general is predicted by the fact that for the electron, with no substructure, it has its \( \kappa_i \) equal to 1, meaning it matches the analytic structure of a force subjected to a cross-sectional area, directly.

Theoretical Framework

In special relativity, the invariant spacetime interval is given by:

\[ ds^2 = c^2 dt^2 - dx^2 - dy^2 - dz^2 \]

For an object at rest the motion is entirely in the temporal dimension. As an object acquires spatial velocity, its temporal velocity decreases according to:

\[ v_t = \frac{c}{\gamma} = c \sqrt{1 - \frac{v^2}{c^2}} \]

where \( \gamma \) is the Lorentz factor. This relationship reveals the hyperbolic nature of spacetime rotations — increasing spatial velocity requires decreasing temporal velocity to maintain the constant magnitude \( c \).

The Universal Particle Equation

We introduce two equations that give on the order of 1-second in terms of the proton radius and mass:

1. \[ \left( \sqrt{\phi \cdot \frac{\pi r_p}{\alpha^4 G m_p^3}} \right) \frac{1}{3} \cdot \frac{h}{c} = 1 \text{ second} \]

2. \[ \left( \frac{1}{6\alpha^2} \sqrt{\frac{4\pi h}{G c}} \right) \cdot \frac{r_p}{m_p} = 1 \text{ second} \]

\( m_p : 1.67262 \times 10^{-27} \, \text{kg} \) (Proton Mass) [1]

\( r_p : 0.833 \times 10^{-15} \, \text{m} \) (Proton Radius) [2]

\( h : 6.62607 \times 10^{-34} \, \text{J} \cdot \text{s} \) (Planck Constant) [3]

\( c : 299,792,458 \, \text{m/s} \) (Light Speed) [4]

\( G : 6.6730 \times 10^{-11} \, \text{N} \cdot \frac{\text{m}^2}{\text{kg}^2} \) (Universal Gravitational Constant, 2018) [5]

\( \alpha : 1/137 \) (Fine Structure Constant)

\( \phi : (\sqrt{5}-1)/2 \approx 0.618 \) (Golden Ratio Conjugate)

These will be verified presently. When setting the left side of equation 1 equal to the left side of equation 2, we get an equation for the radius of a proton that is accurate:

3. \[ r_p = \phi \cdot \frac{h}{c m_p} \]

\[ r_p = (0.618) \cdot \frac{6.62607 \times 10^{-34}}{(299,792,458)(1.67262 \times 10^{-27})} = 0.8166 \times 10^{-15} \, \text{m} \]

The CODATA value from the PRad experiment in 2019 gives \( r_p = 0.831 \, \text{fm} \pm 0.014 \, \text{fm} \). With lower bound \( r_p = 0.817 \times 10^{-15} \, \text{m} \), which is almost exactly what we got.

We can see equation 3 may be the case because we get it from Planck Energy \( E_p = h \nu_p \), Einsteinian energy, \( E_p = m_p c^2 \), and the Compton wavelength \( \lambda_p = h/(m_p c) = r_p \) when we introduce the factor of \( \phi \), which is the golden ratio conjugate, where the golden ratio, \( \Phi = 1/\phi = (\sqrt{5}+1)/2 \approx 1.618 \).

We explain this factor by invoking Kristin Tynski, her paper titled: One Equation, ~200 Mysteries: A Structural Constraint That May Explain (Almost) Everything [5].

Tynski shows that for any system requiring consistency across multiple scales of observation has the recurrence relation:

\[ \text{scale}(n+2) = \text{scale}(n+1) + \text{scale}(n) \]

Which leads to:

\[ \lambda^2 = \lambda + 1 \]

Whose solution is \( \Phi \). Equations 1, 2, and 3 directly yield our Universal Particle Equation:

4. \[ m_p = \kappa_p \cdot \sqrt{\frac{\pi r_p^2 F_n}{G}} \]

5. \[ F_n = \frac{h}{c t_1^2} \]

6. \[ t_1 = 1 \text{ second} \]

where \( \kappa_p = 1/(3\alpha^2) \). Here we see in equation 3, the cross-sectional area of the proton \( A_p = \pi r_p^2 \) is exposed to the normal force, \( F_n \), mediated by the 'stiffness of space' as measured by \( G \), producing the proton mass, \( m_p \). In general we have

7. \[ m_i = \kappa_i \cdot \sqrt{\frac{\pi r_i^2 F_n}{G}}, \]

\[ F_n = \frac{h}{c t_1^2}, \] \[ F_n = \frac{6.62607015 \times 10^{-34} \, \text{J·s}}{(299,792,458 \, \text{m/s})(1 \, \text{s})^2} = 2.21022 \times 10^{-42} \, \text{N}, \] \[ t_1 = 1 \text{ second}, \] \[ m_i = \kappa_i \sqrt{\frac{\pi r_i^2}{G} \cdot \frac{h}{c t_1^2}} \]

We can verify this solving 7 for \( t_1 \) and showing it is on the order, closely, to 1-second:

8. \[ t_1 = \frac{r_i}{m_i} \cdot \sqrt{\frac{\pi h}{G c}} \cdot \kappa_i \]

Proton: \( \kappa_p = \frac{1}{3\alpha^2} \), \( \alpha = 1/137 \):

\[ t_1 = \frac{0.833 \times 10^{-15}}{1.67262 \times 10^{-27}} \cdot \sqrt{\frac{\pi \cdot 6.62607 \times 10^{-34}}{(6.674 \times 10^{-11})(299,792,458)}} \cdot 6256.33 = 1.00500 \text{ seconds} \]

Neutron: \( \kappa_n = \frac{1}{3\alpha^2} \):

\[ t_1 = \frac{0.834 \times 10^{-15}}{1.675 \times 10^{-27}} \cdot \sqrt{\frac{\pi \cdot 6.62607 \times 10^{-34}}{(6.674 \times 10^{-11})(299,792,458)}} \cdot 6256.33 = 1.00478 \text{ seconds} \]

Electron: \( \kappa_e = 1 \):

\[ t_1 = \frac{2.81794 \times 10^{-15}}{9.10938 \times 10^{-31}} \cdot \sqrt{\frac{\pi \cdot 6.62607 \times 10^{-34}}{(6.674 \times 10^{-11})(299,792,458)}} \cdot 1 = 0.99773 \text{ seconds} \]

We suggest \( \kappa_e = 1 \) for the electron may be because it is the fundamental quanta (does not consist of further more elementary particles). G has been rounded to \( 6.674 \times 10^{-11} \). This is a Natural Law.

\( r_n = 0.84 \times 10^{-15} \, \text{m} \) (Neutron radius) [6]

\( r_e = 2.81794 \times 10^{-15} \, \text{m} \) (Classical electron radius) [7]

The Geometric Mechanism of Inertia

As such the geometric mechanism for inertia is that when we apply a force to accelerate a particle spatially, we are rotating its velocity vector, diverting motion from the temporal dimension to spatial dimensions. The normal force \( F_n \) resists this rotation, manifesting as an inertial resistance. \( t_1 = 1 \text{ second} \) given by equation 8 is Lorentz invariant because \( G \), \( c \), and \( h \) are invariant, \( r_p \) is not but the ratio \( r_p/m_p \) is invariant because while \( r_p \) is frame dependent, it is adjusted for by the relativistic mass of \( m_p \).

The dimensionless factor \( \kappa_i \) distinguishes elementary particles from composite hadrons. Remarkably, the same \( t_1 = 1 \, \text{s} \) emerges for the electron, proton, and neutron when their respective \( \kappa_i \) are chosen appropriately.

The factor \( \frac{1}{3} \) reflects the three valence quarks inside the proton and neutron. The \( \alpha^{-2} \) appears because the proton's small radius (relative to its mass) is set by the strong interaction, which is \( \sim 1/\alpha \) times stronger than electromagnetism. Consequently, the required enhancement scales as the square of that ratio because it deals with surface area.

The electron: \( \kappa_e = 1 \) as the baseline

For the electron, using its classical radius \( r_e = 2.81794 \times 10^{-15} \, \text{m} \) and mass \( m_e = 9.10938 \times 10^{-31} \, \text{kg} \), the equation (8) with \( \kappa_e = 1 \) gives

\[ t_1 = \frac{r_e}{m_e} \sqrt{\frac{\pi h}{G c}} \cdot 1 = 0.99773 \, \text{s}, \]

within 0.23% of 1 second. This shows that the electron naturally satisfies the invariant without any extra factor. The value \( r_e \) is not assumed as a physical radius; rather, the invariant predicts it. Solving \( t_1 = 1 \, \text{s} \) for \( r_e \) yields

\[ r_e = m_e \sqrt{\frac{G c}{\pi h}} = 2.82 \times 10^{-15} \, \text{m}, \]

which coincides exactly with the classical electron radius. Thus the classical radius — usually derived from electrostatics — emerges here purely from geometry, \( F_n \), and \( G \). The electron remains point‑like; the radius is an effective length that quantifies the ratio \( r_e/m_e \) required by the invariant.

Concerning bound matter (like atoms) we assume since protons in the nucleus of an atom have spaces between them due to electric forces, and protons and neutrons may be touching, but not existing in the same space, the mass of an atom is the sum of the masses of its protons, electrons, and neutrons as given by the theory of inertia in this paper: they all expose a cross-sectional area to the normal force.

Discussion

The normal force has a relationship to the Planck force, the maximum gravity for the minimum mass. It links the normal force to a full rotation (\( 2\pi \)). We have the normal force

\[ F_n = \frac{h}{c t_1^2} = 2.21022 \times 10^{-42} \, \text{N} \]

We have the Planck force for gravity

\[ F_{\text{Planck}} = G \frac{m_P^2}{l_P^2} = (6.674 \times 10^{-11}) \frac{(2.176434 \times 10^{-8} \, \text{kg})^2}{(1.616255 \times 10^{-35} \, \text{m})^2} = 1.21020 \times 10^{44} \, \text{N} \]

Where, \( m_P \) is the Planck mass, and \( l_P \) is the Planck length. They are given by:

\[ m_P = \sqrt{\frac{\hbar c}{G}} = 2.176434 \times 10^{-8} \, \text{kg} \] \[ l_{\text{Planck}} = \sqrt{\frac{\hbar G}{c^3}} = 1.616255 \times 10^{-35} \, \text{m} \]

And, Planck time is

\[ t_{\text{Planck}} = \sqrt{\frac{\hbar G}{c^5}} = 5.391247 \times 10^{-44} \, \text{s} \]

We form the ratios between the normal force and Planck force:

\[ \frac{F_n}{F_{\text{Planck}}} = 1.826326 \times 10^{-86} \]

Divide by Planck time squared and we have

\[ \frac{F_n}{F_{\text{Planck}}} \frac{1}{t_P^2} = 6.2834743 \, \text{s}^{-2} \]

That number is \( 2\pi \). We have the final equation:

9. \[ t_1 = \sqrt{2\pi \frac{F_{\text{Planck}}}{F_n}} \cdot t_P = 1.00 \text{ seconds} \]

From the Planck units we have:

\[ F_{\text{Planck}} = G \frac{m_P^2}{l_P^2} = \frac{c^4}{G} \]

So it can be written:

10. \[ t_1 = \sqrt{2\pi \frac{c^4}{G F_n}} \cdot t_P \]

We can write

11. \[ F_n = 2\pi F_{\text{Planck}} \cdot \frac{t_P^2}{t_1^2} \]

\( 2\pi \) is a full rotation, so we can define an angular frequency, \( \omega \):

\[ F_n = F_{\text{Planck}} \cdot t_P^2 \cdot \frac{d\omega}{dt} \]

12. \[ \frac{F_n}{F_{\text{Planck}}} \cdot \frac{1}{t_P^2} \int_0^{1 \text{ second}} dt = \omega_1 \]

13. \[ \omega_1 = \frac{2\pi}{\text{second}} \]

Integrating one more time gives the angle over 1-second:

14. \[ \frac{F_n}{F_{\text{Planck}}} \cdot \frac{t_1}{t_P^2} \int_0^{1 \text{ second}} dt = \theta_1 \]

15. \[ \frac{F_n}{F_{\text{Planck}}} \cdot \frac{t_1^2}{t_P^2} = \theta_1 \]

16. \[ \theta_1 = 2\pi \]

which holds for any value of \( t_1 \) because the factors of \( t_1 \) cancel. This identity does not determine the numerical value of the second; rather, it shows that when \( t_1 \) is taken as the empirical 1 second invariant (obtained from the proton, neutron, and electron masses and radii via equation (8)), the ratio \( F_n/F_{\text{Planck}} \) acquires a clear geometric meaning: over one second, the accumulated angular phase is exactly \( 2\pi \) — a full rotation in the temporal dimension. Thus the Planck scale relation is not a derivation of the second but a consistency check and an elegant reinterpretation: the second is the time required for the normal force, when scaled by the Planck force, to close a complete cycle, reinforcing the view that time emerges from a cyclic variable in the quantum vacuum. Moreover, the identity can be rearranged as

\[ \frac{F_n}{F_{\text{Planck}}} = 2\pi \left( \frac{t_P}{t_1} \right)^2 = 2\pi (t_P \cdot \nu_0)^2, \]

where \( \nu_0 = 1/t_1 = 1 \, \text{Hz} \). This reveals a natural angular frequency \( \omega_0 = 2\pi \nu_0 = 2\pi \, \text{rad/s} \), a universal resonance at one hertz that links the Planck scale to the macroscopic normal force. Hence, even though the numeric value \( t_1 = 1 \, \text{s} \) is ultimately fixed by particle data, the interpretation as a \( 2\pi \) phase per second is independent and suggests that inertia is governed by a fundamental clock ticking at exactly one hertz.

From golden ratio to coupling constants

The golden ratio conjugate \( \phi = (\sqrt{5}-1)/2 \) arises naturally from the scale invariant recurrence \( \text{scale}(n+2) = \text{scale}(n+1) + \text{scale}(n) \), which Tynski showed governs systems that must be consistent across multiple observational scales. Applying this to the proton gives \( r_p = \phi \, h/(m_p c) \), which matches the experimental radius. Substituting this \( r_p \) into the universal particle equation \( m_p = \kappa_p \sqrt{\pi r_p^2 F_n / G} \) and using \( F_n = h/(c t_1^2) \) with \( t_1 = 1 \, \text{s} \) yields a closed expression for \( \kappa_p \). Solving it gives \( \kappa_p = 1/(3\alpha^2) \), where \( \alpha \) is the fine structure constant. The factor \( 1/3 \) reflects the three valence quarks in the proton, while \( \alpha^{-2} \) accounts for the electromagnetic and gluonic enhancement of the normal force inside a composite hadron. The neutron, having a similar internal structure, inherits the same \( \kappa_n = 1/(3\alpha^2) \) when its magnetic radius is used. Thus the golden ratio not only predicts the proton's size but also, via the universal particle equation, determines the large coupling constants for hadrons, leaving the electron as the minimal case \( \kappa_e = 1 \). This elegant link between geometry (\( \phi \)), quantum dynamics (\( \alpha \)), and compositeness (three quarks) strongly supports the physical reality of the normal force and the 1 second invariant.

Conclusion

We have presented a fundamental 1‑second invariant that emerges from the intrinsic properties of elementary particles — the proton, neutron, and electron — and from the fabric of Planck‑scale physics. The invariant is expressed as

\[ t_1 = \frac{r_i}{m_i} \sqrt{\frac{\pi h}{G c}} \; \kappa_i = 1 \text{ second}, \]

where \( \kappa_p = \kappa_n = 1/(3\alpha^2) \) and \( \kappa_e = 1 \).

Crucially, the invariant leads to a universal particle equation:

\[ m_i = \kappa_i \sqrt{\frac{\pi r_i^2 F_n}{G}}, \qquad F_n = \frac{h}{c \, t_1^2}, \]

with \( F_n \) a constant normal force of magnitude \( 2.21022 \times 10^{-42} \, \text{N} \). This equation suggests that the mass of a particle is determined by its cross‑sectional area (\( \pi r_i^2 \)), the stiffness of spacetime (\( G \)), and a universal normal force \( F_n \) that arises from the quantum constraint \( t_1 = 1 \, \text{s} \).

The geometric origin of the second becomes apparent when we relate \( F_n \) to the Planck force \( F_{\text{Planck}} = c^4/G \). We find

\[ \frac{F_n}{F_{\text{Planck}}} \cdot \frac{t_1^2}{t_P^2} = 2\pi, \]

which means that over one second, the ratio \( F_n/F_{\text{Planck}} \) accumulates exactly \( 2\pi \) radians of angular phase — a full rotation. Thus, one second is not an arbitrary human convention but rather the time required for this cyclic closure in the temporal dimension, rooted in Planck‑scale dynamics.

In summary, the 1‑second invariant unifies particle physics and fundamental constants through a single, testable relation. The universal particle equation \( m_i = \kappa_i \sqrt{\pi r_i^2 F_n / G} \) provides a new perspective on inertia: mass arises from the resistance to rotating a particle's temporal velocity into spatial velocity, quantified by the normal force \( F_n \). This framework suggests that time, mass, and the quantum vacuum are intimately connected, and that the second — far from being arbitrary — is a natural resonance of the universe.

Note

The universal particle equation and 1-second invariant were discovered by the author and reported as early as:

Beardsley, Ian (November 29, 2025) The Geometric Origin of Inertia: Mass Generation from Temporal Motion in Hyperbolic Spacetime, https://doi.org/10.5281/zenodo.17772255

Beardsley, I. (2026). A Spacetime Theory For Inertia; Predicting The Proton, Electron, Neutron and the Solar System in Terms of a One-Second Invariant, https://doi.org/10.5281/zenodo.18165383

References

[1] Tiesinga, Eite, Peter J. Mohr, David B. Newell, and Barry N. Taylor. "CODATA Value: Proton Mass." The 2022 CODATA Recommended Values of the Fundamental Physical Constants (Web Version 9.0). National Institute of Standards and Technology, 2024. https://physics.nist.gov/cgi-bin/cuu/Value?mp.

[2] Bezginov, N., Valdez, T., Horbatsch, M. et al. (York University/Toronto). Published in Science, Vol. 365, Issue 6457, pp. 1007-1012 (2019) "A measurement of the atomic hydrogen Lamb shift and the proton charge radius".

[3] Tiesinga, Eite, Peter J. Mohr, David B. Newell, and Barry N. Taylor. "CODATA Value: Planck Constant." The 2022 CODATA Recommended Values of the Fundamental Physical Constants (Web Version 9.0). National Institute of Standards and Technology, 2024. https://physics.nist.gov/cgi-bin/cuu/Value?h.

[4] Tiesinga, Eite, Peter J. Mohr, David B. Newell, and Barry N. Taylor. "CODATA Value: Speed of Light in Vacuum." The 2022 CODATA Recommended Values of the Fundamental Physical Constants (Web Version 9.0). National Institute of Standards and Technology, 2024. https://physics.nist.gov/cgi-bin/cuu/Value?c.

[5] Tynski, K. (2024). One Equation, ~200 Mysteries: A Structural Constraint That May Explain (Almost) Everything.

[6] Kubon, G., Anklin, H., Bartsch, P., Baumann, D., Boeglin, W. U., Bohinc, K., ... & Zihlmann, B. (2002). Precise neutron magnetic form factors. Physics Letters B, 524(1-2), 26-32.

[7] NIST CODATA Value for the Classical Electron Radius (2022).


On the 1 Hz "Noise" and the Case for a Torsion Pendulum Test of the Temporal Invariant

Ian Beardsley

Hillbilly Research Division (Independent)

(Date: June 2026)

Abstract

The claim of a universal 1 second invariant \( \tau_0 = 1 \, \text{s} \) and a concomitant normal force \( F_n = h/(c \tau_0^2) \) implies that any dynamical system coupling to the resistance of temporal rotation should exhibit an anomalous resonant response at exactly \( \omega_0 = 2\pi \, \text{rad/s} \) (\( f_0 = 1 \, \text{Hz} \)). Torsion pendulums have been used in precision experiments for centuries, but a systematic search for a sharp, unexplained peak at 1 Hz has never been performed because such a peak is conventionally dismissed as environmental noise or electronic artifact. This paper reviews the known sources of 1 Hz contamination (Nyquist aliasing, pendulum cross coupling, microseisms, clock feedthrough) and shows that none of them can account for a persistent, amplitude insensitive, and drive phase locked peak that survives standard control tests. We propose a dedicated torsion pendulum experiment with oversampling, analog antialiasing filtering, and a set of falsifiable controls. If the predicted 1 Hz resonance is observed, it would provide the first direct experimental evidence for the temporal invariant; its absence, after proper artifact elimination, would falsify the central prediction of the theory.

1. Introduction

In a recent particle scale framework (Beardsley 2026), a universal invariant \( \tau_0 = 1 \, \text{s} \) emerges from the masses and radii of the proton, neutron and electron when combined with the normal force

\[ F_n = \frac{h}{c \, \tau_0^2} \approx 2.21 \times 10^{-42} \, \text{N}. \]

The invariant gives rise to a natural angular frequency \( \omega_0 = 2\pi/\tau_0 = 2\pi \, \text{rad/s} \) (\( f_0 = 1 \, \text{Hz} \)). The physical interpretation is that inertia originates from the resistance to rotating a particle's velocity from the temporal dimension into spatial dimensions. Consequently, any macroscopic system that involves periodic acceleration — in particular a driven torsion pendulum — should exhibit a resonant enhancement of its response when driven exactly at \( \omega_0 \). This enhancement is not a mechanical eigenmode; it is a direct manifestation of the universal normal force coupling to the pendulum's cross sectional area.

Searching the experimental literature, one finds occasional reports of unexplained "bumps" near 1 Hz in torsion balance data, but these are invariably attributed to environmental or electronic artifacts (microseisms, aliasing, crosstalk, parasitic swing modes). No experiment has ever been designed to systematically discriminate between those well known artifacts and a genuine new resonance that would be phase locked to the drive frequency and independent of the pendulum's moment of inertia. This paper reviews the physics of 1 Hz noise in torsion pendulums and outlines a clean, falsifiable experiment that can unambiguously test the temporal invariant prediction.

2. Why 1 Hz is Dirty — But Not Unambiguously

Precision torsion balances (such as those used in the EötWash experiment or for measuring the gravitational constant \( G \)) are usually operated at much lower frequencies (mHz to tenths of Hz) to avoid seismic and thermal noise. Nevertheless, when a pendulum is actively driven at 1 Hz, the following contaminants are known to appear:

2.1 Nyquist aliasing

If the data acquisition samples at a rate \( f_s \), any signal component above the Nyquist frequency \( f_N = f_s/2 \) is folded back into the measured band. For a 1 Hz signal of interest, sampling at \( f_s = 2 \, \text{Hz} \) would place the Nyquist limit exactly at 1 Hz, leading to severe aliasing (a pure 1 Hz input can appear as a DC offset or as an arbitrary low frequency). However, this is trivially avoided by oversampling: with \( f_s \ge 100 \, \text{Hz} \), the Nyquist limit is above 50 Hz, and no aliasing of a 1 Hz signal occurs. Modern microcontrollers easily achieve 1 kHz sampling, so aliasing is a solvable problem, not an intrinsic obstacle.

2.2 Parasitic pendular (swinging) modes

A torsion pendulum is suspended by a thin fiber. If the driving force is not perfectly aligned with the torsional axis, or if the fiber is slightly asymmetric, the drive can couple into translational swing modes. For a fiber of length \( L \), the pendular frequency is

\[ f_{\text{pend}} = \frac{1}{2\pi} \sqrt{\frac{g}{L}}. \]

For \( L \approx 0.25 \, \text{m} \), \( f_{\text{pend}} \approx 1 \, \text{Hz} \). Therefore, a 1 Hz drive can easily excite the swing mode if any misalignment exists. That swing mode will appear as an anomalous peak in the torsional signal because the optical readout cannot perfectly distinguish pure rotation from horizontal translation. This artifact is eliminated by:

  • Balancing the pendulum mass symmetrically and using a fiber with high torsional stiffness (low swing resonance) or, conversely, by designing the fiber such that the pendular frequency is far from 1 Hz (e.g., \( L = 1 \, \text{m} \) gives \( f_{\text{pend}} \approx 0.5 \, \text{Hz} \)).
  • Using a second, independent sensor (e.g., a lateral position sensor) to monitor and subtract the swing component.
  • Verifying that the anomaly disappears when the drive amplitude is reduced to zero (no artificial excitation of the swing mode).

2.3 Environmental microseisms

Building vibrations, HVAC systems, walking on floors, and even computer fans often have sharp spectral components near 1 Hz. These vibrations act as a direct displacement of the suspension point, which is indistinguishable from a torque on the pendulum. This noise is typically reduced by:

  • Placing the apparatus on a massive concrete block supported by vibration damping foam or pneumatic legs.
  • Enclosing the pendulum in a vacuum chamber (to also remove air damping and acoustic coupling).
  • Measuring the ambient acceleration with a seismometer and subtracting its contribution coherently (crosscorrelation).

2.4 Electronic clock feedthrough

Many precision instruments, data loggers, and microcontrollers operate internal loops at exactly 1 Hz (e.g., updating a display, polling a sensor, or generating a timing interrupt). Capacitive or magnetic coupling between the digital lines and the sensitive pendulum readout (a photodiode, position sensitive detector, or capacitive bridge) can inject a pure 1 Hz voltage directly into the signal. This artifact is identified by:

  • Disconnecting the drive and the pendulum readout while keeping the electronics powered; a residual 1 Hz peak indicates clock feedthrough.
  • Shielding all signal cables and using differential (balanced) connections.
  • Changing the microcontroller's update rate (e.g., from 1 Hz to 1.5 Hz) — a real physical peak remains at 1 Hz, an electronic artifact follows the clock frequency.

3. Why Previous Null Results Do Not Falsify the Theory

Importantly, the fact that no experiment has ever reported an unexplained 1 Hz peak in a driven torsion pendulum is exactly what the theory predicts for any experiment not designed to distinguish the predicted effect from the artifacts listed above. Standard practice is to treat any low frequency peak as noise and to filter it out or subtract it without further investigation. No experimental group has had a theoretical reason to perform the controls that would reveal a genuine new resonance — a resonance that would be:

  • Strictly proportional to the drive amplitude (linear response),
  • Independent of the pendulum's natural frequency (i.e., it does not shift when the moment of inertia is changed),
  • Phase locked to the drive signal, and
  • Unaffected by changing the sampling rate, the shielding, or the isolation of the pendulum.

Because those controls have never been systematically applied, the absence of a prior report is not evidence against the effect; it simply means the effect was never looked for in a way that could distinguish it from the noise floor.

4. Mathematical Model of the Predicted Resonance

In the temporal invariant theory, a test body of mass \( m \) and effective cross sectional area \( A_{\text{eff}} = \pi r^2 \) experiences a normal force \( F_n = h/(c \tau_0^2) \) when its velocity is rotated from the temporal to spatial axes. For a torsion pendulum with moment of inertia \( I \) and torsional stiffness \( k_\theta \), the equation of motion in the presence of an external drive torque \( \tau_{\text{drive}}(t) \) becomes

\[ I \ddot\theta + b \dot\theta + k_\theta \theta = \tau_{\text{drive}}(t) + \tau_{\text{invariant}}(t), \]

where \( \tau_{\text{invariant}}(t) \) is the torque produced by the coupling of the rotating pendulum mass to the universal normal force. For a simple geometry (a point mass \( m \) at distance \( R \) from the axis), the invariant contribution is

\[ \tau_{\text{inv}} = R \, F_n \, A_{\text{eff}} \, \sin(\omega_0 t + \phi_0), \]

with \( \omega_0 = 2\pi/\tau_0 \). The resulting steady state amplitude at the drive frequency \( \omega \) is given by the well known driven harmonic oscillator response, but with an additional resonance denominator that becomes singular when \( \omega = \omega_0 \):

\[ \theta(\omega) = \frac{\tau_{\text{drive}}(\omega) + \frac{R A_{\text{eff}} F_n}{I} \, \delta(\omega - \omega_0)}{k_\theta - I\omega^2 + i b\omega}. \]

Hence, when \( \omega = \omega_0 \), the amplitude increases regardless of the pendulum's natural frequency. The fractional increase can be estimated from the dimensionless coupling constant

\[ \kappa = \frac{R A_{\text{eff}} F_n}{I \, \omega_0^2 \, \theta_{\text{drive}}}, \]

which, for a milligram scale mass and millimeter scale radius, yields a potentially measurable shift of order \( 10^{-6} \) rad. Modern capacitive or optical readouts can resolve better than \( 10^{-8} \) rad, so the effect is within reach.

5. Experimental Protocol to Unambiguously Test the Prediction

Based on the above analysis, we propose the following minimal experiment that can falsify or confirm the 1 Hz invariant.

5.1 Apparatus

  • A torsion pendulum with a symmetric crossbar (e.g., a thin aluminum rod, length 20 cm, with adjustable masses at the ends). The fiber is a 50 µm tungsten wire, length 1 m, giving a torsional period of several seconds (low natural frequency) to avoid confusion with the drive.
  • An optical lever (laser + position sensitive detector) or a high resolution autocollimator, sampling at 1000 Hz.
  • An electromagnetic drive coil and a small permanent magnet attached to the pendulum. The drive is a pure sine wave from a function generator, with amplitude stabilized.
  • An analog lowpass antialiasing filter (corner frequency 50 Hz) placed immediately after the photodiode amplifier.
  • A massive vibration isolated base (granite slab on Sorbothane feet) inside a grounded Faraday cage.

5.2 Control tests

  • Natural frequency variation: add or remove mass at the ends; the pendulum's torsional eigenfrequency changes by >30%, but the predicted peak must stay exactly at 1 Hz.
  • Change of drive amplitude: the resonance amplitude should be strictly linear with drive amplitude. Any nonlinearity (e.g., from magnetic coupling) would indicate an artifact.
  • Change of sampling rate: run the same experiment with sampling rates of 200 Hz, 500 Hz and 1000 Hz. A true physical peak remains unchanged; a digital aliasing artifact changes dramatically.
  • Electronic crosstalk test: with the pendulum locked (or removed), drive the coil at 1 Hz and record the readout sensor output. Any observed 1 Hz signal is purely electromagnetic pickup and must be eliminated by shielding and balanced wiring.
  • Environmental noise map: measure the pendulum output with the drive off for 1 hour. If a 1 Hz peak appears in the power spectrum, it is due to ambient vibrations or clock feedthrough — not the predicted effect.

5.3 Falsification criterion

The theory is falsified if, after implementing all the above controls, no statistically significant excess amplitude is observed at \( f_0 = 1.000 \, \text{Hz} \) (within the resolution of the frequency generator, \( \pm 0.001 \, \text{Hz} \)) when compared to neighbouring frequencies (0.9 Hz, 0.95 Hz, 1.05 Hz, 1.1 Hz). Conversely, a clear, reproducible peak that survives all controls would constitute the first direct evidence for the temporal invariant and would require a major revision of our understanding of inertia.

6. Relation to Other Proposed Tests (Plasma Thruster)

The same 1 Hz resonance is also predicted for pulsed plasma thrusters. However, the torsion pendulum is far simpler, cheaper, and less prone to unmodeled plasma dynamics. A positive result with the pendulum would immediately justify more ambitious tests (e.g., with a Hall thruster). A null result, if properly controlled, would rule out the universal coupling at the macroscopic level, though the particle scale invariant might still hold. Hence the torsion pendulum test is the ideal first step experiment.

7. Conclusion

The 1 Hz "noise" that appears in all torsion pendulum measurements is a well studied collection of environmental and instrumental artifacts. None of these artifacts produce a peak that is simultaneously linear in drive amplitude, independent of the pendulum's eigenfrequency, unchanged by sampling rate, and persistent under rigorous shielding. A dedicated experiment that systematically controls each artifact can either reveal the predicted universal resonance or place an upper limit on the coupling constant that will falsify the temporal invariant theory. Given the low cost and high sensitivity of modern torsion balances, such an experiment is both feasible and urgent. The physics community should therefore move beyond dismissing 1 Hz as "just noise" and perform the definitive test.

References

[1] Beardsley, I. (2026). "A Universal Particle Equation: Mass, Inertia and the 1 Second Invariant." Zenodo. DOI: 10.5281/zenodo.19930951 (preprint).

[2] Beardsley, I. & Blackwell, D. E. (2026). "Three Dimensional Simulation of Informational Warp Bubble Dynamics." Zenodo.

[3] Newman, R. D. & Bantel, M. K. (1999). "On the status of measurements of Newton's gravitational constant." Meas. Sci. Technol. 10, 445.

[4] Speake, C. C. & Quinn, T. J. (2006). "The gravitational constant: theory and experiment." Phys. Today 59, 33.

[5] Matsumura, S. et al. (2015). "Vibration isolation system for a torsion pendulum." Rev. Sci. Instrum. 86, 064501.


A Quantum Analog For The Solar System

Ian Beardsley

March 7, 2026

Abstract

We find if we consider the evolved state of the Solar System, that its quantum analog to the Bohr atom is based on a characteristic time of one-second and the Earth's Moon as the defining metric.

1.0 The Quantum Solution To The Solar System

The ancient Sumerians (4500 BCE – 1900 BCE) used base 60 counting, and divided the Earth day into 24 hours. The ancient Egyptians (3100 BCE – 30 BCE) divided the Earth day into 24 hours as well. Since they both divided the day into 12 hours, and the night into 12 hours and, in the winter, the night is longer than the day and in the summer, the day is longer than the night, the hours in a day, or night, can be longer or shorter depending on the time of the year. The ancient Greeks took the 24 hour day from the ancient Egyptians (Hipparchus, 190 BCE – 120 BCE) and used an hour to be represented by the equinoxes when day equals night, inventing the equinoctial hour. It was Christiaan Huygens (1629 – 1695) who took the hour that had been divided up into 60 minutes, with each minute divided into 60 seconds, from the ancient Sumerian base 60 counting, and built the first pendulum clock that could measure down to the second accurately. This was fueled by the need of Newton's (1642 – 1727) world view for gravity and mechanics that needed to measure time down to a unit as small as a second.

It is an interesting phenomenon that the Moon near perfectly eclipses the Sun. The eclipse ratio that allow for this is about 400:

1.1 \[ \frac{R_{\odot}}{R_m} \approx 400 \quad \text{and} \quad \frac{r_{\oplus}}{r_m} \approx 400 \]

where \( R_{\odot} \) is the radius of the Sun and \( R_m \) is the radius of the Moon. \( r_{\oplus} \) is the orbital radius of the Earth and \( r_m \) is the orbital radius of the Moon. The solar radius is about 400 times the lunar radius; the Earth-Sun distance is about 400 times the Earth-Moon distance.

The number of seconds in a day are given approximately by:

1.2 \[ 86,400 \text{ seconds/day} = (24 \text{ hours})(60 \text{ minutes})(60 \text{ seconds}) \]

The number of seconds in a day, 86,400, can be factored as:

1.3 \[ 86,400 = (6)(6)(6)(400) \]

The factor 400 is the eclipse ratio. The factor \( 6^3 \) (216) relates to sixfold symmetry, hexagonal tiling, and the approximation \( \pi \approx 3 \) used by Archimedes as his starting point for calculating \( \pi \). The appearance of 86,400 in ancient timekeeping thus incorporates the eclipse ratio, whether by accident or by design.

Let us suggest that the kinetic energy of the Moon to the kinetic energy of the Earth maps the 24 hour (Earth rotation period) day into 1 second, our basis unit of measuring time:

1.4 \[ \frac{KE_{\text{moon}}}{KE_{\text{earth}}}(24 \text{ hours}) \cos(\theta) = 1 \text{ second} \]

Where \( \theta = 23.5^\circ \) is the inclination of the Earth to its orbit.

\( KE_{\text{earth}} = (5.9722 \times 10^{24} \, \text{kg})(29,800 \, \text{m/s})^2 = 5.30355 \times 10^{33} \, \text{J} \)

\[ \frac{7.6745 \times 10^{28} \, \text{J}}{5.30355 \times 10^{33} \, \text{J}} (86,400 \, \text{s}) \cos(23.5^\circ) = 1.1466 \text{ seconds} \approx 1 \text{ second} \]

Using average orbital velocities. We can get closer to a second using aphelions and perihelions and perigees and apogees.

The Moon stabilizes Earth's axial tilt:

\[ \theta = 23.5^\circ \pm 1.3^\circ \quad \text{(with Moon)} \]

\[ \theta = 0^\circ \text{ to } 85^\circ \quad \text{(without Moon, chaotic)} \]

The Moon stabilizing the Earth's tilt to its orbit prevents extreme hot and cold on Earth and allows for the seasons. As such the Moon is key to optimizing conditions for life on the planet. Perhaps making it possible for intelligent life to evolve.

We form a Planck-type constant for the Solar System:

1.5 \[ \hbar_{\odot} = (1 \text{ second}) \, KE_{\text{earth}} \]

We take \( \hbar_{\odot} \) to be given by:

1.6 \[ 1.03351 \, \text{s} = \frac{1}{3} \cdot \frac{h}{\alpha^2 c} \sqrt{\frac{2}{3} \cdot \frac{\pi r_p}{G m_p^3}} \]

Equation 6 is an approximately 1-second expression for the radius and mass of a proton that uses a 2/3 fibonacci approximation for \( \phi \), discovered by the author. Thus we see we can see a possible 1-second invariant that may exist across vast scales from atoms to the Solar System. We have

1.7 \[ \hbar_{\odot} = (1.03351 \, \text{s})(2.7396 \times 10^{33} \, \text{J}) = 2.8314 \times 10^{33} \, \text{J} \cdot \text{s} \]

\( KE_{\text{Earth}} = \frac{1}{2}(5.972 \times 10^{24} \, \text{kg})(30,290 \, \text{m/s})^2 = 2.7396 \times 10^{33} \, \text{J} \)

Using Earth's orbital velocity at perihelion.

The ground state energy for a hydrogen atom (One electron orbiting a proton) is:

1.8 \[ r_1 = -\frac{\hbar^2}{k_e e^2 m_e} \]

For the planetary system we would replace \( k_e \) (Coulomb's constant) with \( G \) (Newton's universal constant of gravity). The product of \( e^2 \) (the charge of an electron squared) and \( m_e \) (the mass of an electron) become a mass cubed. We will choose the mass of the Moon, \( M_m \). We have the ground state equation is:

1.9 \[ \frac{\hbar_{\odot}^2}{G M_m^3} = \frac{(2.8314 \times 10^{33})^2}{(6.67408 \times 10^{-11})(7.34763 \times 10^{22} \, \text{kg})^3} = 3.0281 \times 10^8 \, \text{m} \]

1.10 \[ \frac{\hbar_{\odot}^2}{G M_m^3} \frac{1}{c} = \frac{3.0281 \times 10^8 \, \text{m}}{299,792,458 \, \text{m/s}} = 1.010 \text{ seconds} \approx 1 \text{ second} \]

Where we have converted meters to seconds by measuring distance in terms of time with the speed of light (\( c \)). We see the mass of the Moon maps the kinetic energy of the Earth over one second to 1 second. The Moon is the metric.

The solution for the orbit of the Earth around the Sun with the Schrödinger wave equation can be inferred from the solution for an electron around a proton in a hydrogen atom with the Schrödinger wave equation. The Schrödinger wave equation is, in spherical coordinates

\[ -\frac{\hbar^2}{2m} \left[ \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial}{\partial r} \right) + \frac{1}{r^2 \sin\theta} \frac{\partial}{\partial\theta} \left( \sin\theta \frac{\partial}{\partial\theta} \right) + \frac{1}{r^2 \sin^2\theta} \frac{\partial^2}{\partial\phi^2} \right] \psi + V(r)\psi = E\psi \]

Its solution for the atom is as guessed by Niels Bohr before the wave equation existed:

1.11 \[ E_n = -\frac{Z^2 (k_e e^2)^2 m_e}{2\hbar^2 n^2} \]

1.12 \[ r_n = \frac{n^2 \hbar^2}{Z k_e e^2 m_e} \]

\( E_n \) is the energy for an electron orbiting \( Z \) protons and \( r_n \) is the orbital shell for an electron with \( Z \) protons, \( n \) the orbital number. I find the solution for the Earth around the Sun utilizes the Moon around the Earth. This is different than with the atom because planets and moons are not all the same size and mass like electrons and protons are, and they don't jump from orbit to orbit like electrons do. I find that for the Earth around the Sun

1.13 \[ E_n = \sqrt{n} \, \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^2 M_m^3}{2 \hbar_{\odot}^2} \]

1.14 \[ r_n = \frac{2 \hbar_{\odot}^2}{G M_m^3} \cdot \frac{R_{\odot}}{R_m} \cdot \frac{1}{\sqrt{n}} \]

\( E_3 \) is the energy of the Earth, and \( r_n \) is the planet's orbit. \( R_{\odot} \) is the radius of the Sun, \( r_m \) is the radius of the Moon's orbit, \( M_e \) is the mass of the Earth, \( M_m \) is the mass of the Moon, \( n \) is the orbit number of the Earth which is 3 and \( \hbar_{\odot} \) is the Planck constant for the solar system. Instead of having \( Z \) protons, we have \( R_{\odot}/R_m \) the radius of the Sun normalized by the radius of the Moon. We see that the Moon is indeed the metric, as we said before.

\[ \frac{R_{\odot}}{R_m} = \frac{6.96 \times 10^8 \, \text{m}}{1,737,400 \, \text{m}} = 400.5986 \]

\[ E_3 = (1.732)(400.5986) \frac{(6.67408 \times 10^{-11})^2 (5.972 \times 10^{24} \, \text{kg})^2 (7.347673 \times 10^{22} \, \text{kg})^3}{2 (2.8314 \times 10^{33})^2} = 2.727 \times 10^{33} \, \text{J} \]

The kinetic energy of the Earth is (using orbital velocity at perihelion):

\( KE_{\text{earth}} = \frac{1}{2}(5.972 \times 10^{24} \, \text{kg})(30,290 \, \text{m/s})^2 = 2.7396 \times 10^{33} \, \text{J} \)

\[ \frac{2.727 \times 10^{33} \, \text{J}}{2.7396 \times 10^{33} \, \text{J}} \times 100 = 99.5\% \]

The kinetic energy of the Earth is about equal to the energy of the system, because the orbit of the Earth is nearly circular. That is

\( E_3 \sim KE_{\text{earth}} \)

The whole object of developing a theory for the way planetary systems form is that they meet the following criterion: They predict the Titius-Bode rule for the distribution of the planets; the distribution gives the planetary orbital periods from Newton's Universal Law of Gravitation. The distribution of the planets is chiefly predicted by three factors: The inward forces of gravity from the parent star, the outward pressure gradient from the stellar production of radiation, and the outward inertial forces as a cloud collapses into a flat disc around the central star. These forces separate the flat disc into rings, agglomerations of material, each ring from which a different planet forms at its central distance from the star. In a theory of planetary formation from a primordial disc, it should predict the Titius-Bode rule for the distribution of planets today, which was the distribution of the rings from which the planets formed.

Also, the Earth has been in the habitable zone since 4 billion years ago when it was at 0.9 AU. Today it is at 1 AU, and that habitable zone can continue to 1.2 AU. So we can speak of the distance to the Earth over much time. The Earth and Sun formed about 4.6 billion years ago. As the Sun very slowly loses mass over millions of years as it burns fuel doing fusion, the Earth slips minimally further out in its orbit over long periods of time. The Earth orbit increases by about 0.015 meters per year. The Sun only loses 0.00007% of its mass annually. The Earth is at 1 AU = 1.496 × 10¹¹ m. We have 0.015 m / 1.496 × 10¹¹ m/AU = 1.00267 × 10⁻¹³ AU. So,

\( (1.00267 \times 10^{-13} \, \text{AU/year})(1 \times 10^9 \, \text{years}) = 0.0001 \, \text{AU} \)

The Earth will only move out one ten thousandth of an AU in a billion years. Anatomically modern humans have only been around for about three hundred thousand years. Civilization began only about six thousand years ago.

The Moon slows the Earth rotation and this in turn expands the Moon's orbit, so it is getting larger, the Earth loses energy to the Moon. The Earth day gets longer by 0.0067 hours per million years, and the Moon's orbit gets 3.78 cm larger per year.

We suggest the Solar system comes into phase with a possible one second invariant when the Earth-Sun separation, and Earth-Moon separation, have kinetic energies whose ratio maps the 24 hour day into the 1-second base unit as given by equation 4:

1.4 \[ \frac{KE_{\text{moon}}}{KE_{\text{earth}}}(24 \text{ hours}) \cos(\theta) = 1 \text{ second} \]

That is when equations 5 and 10 hold:

1.5 \[ \hbar_{\odot} = (1 \text{ second}) \, KE_{\text{earth}} \]

1.10 \[ \frac{\hbar_{\odot}^2}{G M_m^3} \frac{1}{c} = 1.010 \text{ seconds} \approx 1 \text{ second} \]

Something remains to be done. Is there something about the Sun that is common to other types of stars; stars that are perhaps larger and hotter than the Sun, or perhaps smaller and cooler, or a different color, like blue or red, instead of yellow? The answer is yes. I actually found something in ancient Vedic knowledge, in the Hindu traditions. Apparently, in Hindu yoga the number 108 is an important number. I read that yogis today noticed that the diameter of the Sun is about 108 times the diameter of the Earth and that the average distance from the Sun to the Earth is about 108 solar diameters, with 108 being a significant number in yoga. So I wrote the equivalent:

1.15 \[ R_{\oplus} = 2 \frac{R_{\odot}^2}{r_{\oplus}}, \]

or for any star and habitable planet:

1.16 \[ R_{\text{planet}} = 2 \frac{R_{\star}^2}{r_{\text{habitable}}} \]

\( R_{\star} \) the radius of the star. \( r_{\text{habitable}} \) the orbital radius of the habitable planet. We consider the HR diagram that plots temperature versus luminosity of stars. We see the O, B, A stars are the more luminous stars, which is because they are bigger and more massive and the F, G stars are medium luminosity, mass, and size (radius). Our Sun is a G star, particularly G2V, the two because the spectral classes are divided up into 10 sizes, V for five meaning main sequence, that it is part of the S shaped curve and is in the phase where the star is burning hydrogen fuel, its original fuel, not the by products. And the K and M stars are the coolest, least massive, least luminous.

Let us consider the habitable zones of different kinds of stars. In order to get \( r_{\text{habitable}} \), the distance of the habitable planet from the star, we use the inverse square law for luminosity of the star. If the Earth is in the habitable zone, and if the star is one hundred times brighter than the Sun, then by the inverse square law the distance to the habitable zone of the planet is 10 times that of what the Earth is from the Sun. Thus we have in astronomical units the habitable zone of a star is given by:

1.17 \[ r_{\text{habitable}} = \sqrt{\frac{L_{\star}}{L_{\odot}}} \, \text{AU} \]

\( L_{\star} \) the luminosity of the star, \( L_{\odot} \) the luminosity of the Sun. AU the average Earth-Sun separation, which is 1. The surprising result I found was, after applying equation 4, hypothetically predicting the size of a habitable planet, to the stars of all spectral types from F through K, with their different radii and luminosities (the luminosities determine \( r_{\text{habitable}} \), the distances to the habitable zones), that the radius of the planet always came out about the same, about the radius of the Earth. This may suggest optimally habitable planets are not just a function of their distance from the star, which is a big factor in determining their temperature, but are functions of their size and mass meaning the size of the Earth could be good for life chemistry and atmospheric composition, and gravity. Stars of the same particular luminosities, temperatures and colors have about the same mass and size (radius). Here are some examples of such calculations of stars of different sizes, colors, and luminosities using equation 4:

F8V Star

Mass: 1.18
Radius: 1.221
Luminosity: 1.95

\( M_{\star} = 1.18 (1.9891 \times 10^{30} \, \text{kg}) = 2.347 \times 10^{30} \, \text{kg} \)

\( R_{\star} = 1.221 (6.9634 \times 10^8 \, \text{m}) = 8.5023 \times 10^8 \, \text{m} \)

\( r_p = \sqrt{1.95} \, L_{\odot} \, \text{AU} = 1.3964 \, \text{AU} \, (1.496 \times 10^{11} \, \text{m/AU}) = 2.08905 \times 10^{11} \, \text{m} \)

\( R_p = \frac{2 R_{\star}^2}{r_p} = 2 \frac{(8.5023 \times 10^8 \, \text{m})^2}{2.08905 \times 10^{11} \, \text{m}} = \frac{6.92076 \times 10^6 \, \text{m}}{6.378 \times 10^6 \, \text{m}} = 1.0851 \, \text{Earth Radii} \)

F9V Star

Mass: 1.13
Radius: 1.167
Luminosity: 1.66

\( M_{\star} = 1.13 (1.9891 \times 10^{30} \, \text{kg}) = 2.247683 \times 10^{30} \, \text{kg} \)

\( R_{\star} = 1.167 (6.9634 \times 10^8 \, \text{m}) = 8.1262878 \times 10^8 \, \text{m} \)

\( r_p = \sqrt{1.66} \, \text{AU} = 1.28841 \, \text{AU} \, (1.496 \times 10^{11} \, \text{m/AU}) = 1.92746 \times 10^{11} \, \text{m} \)

\( R_p = \frac{2 R_{\star}^2}{r_p} = 2 \frac{(8.1262878 \times 10^8 \, \text{m})^2}{1.92746 \times 10^{11} \, \text{m}} = \frac{6.852184 \times 10^6 \, \text{m}}{6.378 \times 10^6 \, \text{m}} = 1.0743468 \, \text{Earth Radii} \)

G0V Star

Mass: 1.06
Radius: 1.100
Luminosity: 1.35

\( M_{\star} = 1.06 (1.9891 \times 10^{30} \, \text{kg}) = 2.108446 \times 10^{30} \, \text{kg} \)

\( R_{\star} = 1.100 (6.9634 \times 10^8 \, \text{m}) = 7.65974 \times 10^8 \, \text{m} \)

\( r_p = \sqrt{1.35} \, \text{AU} = 1.161895 \, \text{AU} \, (1.496 \times 10^{11} \, \text{m/AU}) = 1.7382 \times 10^{11} \, \text{m} \)

\( R_p = \frac{2 R_{\star}^2}{r_p} = 2 \frac{(7.65974 \times 10^8 \, \text{m})^2}{1.7382 \times 10^{11} \, \text{m}} = \frac{6.751 \times 10^6 \, \text{m}}{6.378 \times 10^6 \, \text{m}} = 1.05848 \, \text{Earth Radii} \)

As you can see we consistently get about 1 Earth radius for the radius of every planet in the habitable zone of each type of star. I have gone through all stars from spectral class A stars to spectral class M stars and consistently got this result. It may be this radius for a planet is optimal for life, in particular intelligent life, because given we might, for that, need a material composition similar to that of Earth, and, in turn, an Earth-like gravity for the right atmosphere, including atmospheric composition, or planetary mass, the planet might need to be around this size.

2.0 The Solar Solution

Our solution of the wave equation for the planets gives the kinetic energy of the Earth from the mass of the Moon orbiting the Earth, but you could formulate based on the Earth orbiting the Sun. In our lunar formulation we had:

2.1 \[ KE_e = \sqrt{3} \, \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^2 M_m^3}{2 \hbar_{\odot}^2} \]

We remember the Moon perfectly eclipses the Sun which is to say

2.2 \[ \frac{R_{\odot}}{R_m} = \frac{r_e}{r_m} \]

Thus equation 2.1 becomes

2.3 \[ KE_e = \sqrt{3} \, \frac{r_e}{r_m} \cdot \frac{G^2 M_e^2 M_m^3}{2 \hbar_{\odot}^2} \]

The kinetic energy of the Earth is

2.4 \[ KE_e = \frac{1}{2} \cdot \frac{G M_{\odot} M_e}{r_e} \]

Putting this in equation 2.3 gives the mass of the Sun:

2.5 \[ M_{\odot} = \sqrt{3} \, r_e^2 \cdot \frac{G M_e}{r_m} \cdot \frac{M_m^3}{\hbar_{\odot}^2} \]

We recognize that the orbital velocity of the Moon is

2.6 \[ v_m^2 = \frac{G M_e}{r_m} \]

So equation 2.5 becomes

2.7 \[ M_{\odot} = \sqrt{3} \, r_e^2 \cdot v_m^2 \cdot \frac{M_m^3}{\hbar_{\odot}^2} \]

This gives the mass of the Moon is

2.8 \[ M_m^3 = \frac{M_{\odot} \, \hbar_{\odot}^2}{\sqrt{3} \, r_e^2 \, v_m^2} \]

Putting this in equation 2.1 yields

2.9 \[ KE_e = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^2 M_{\odot}}{2 \, r_e^2 \, v_m^2} \]

We now multiply through by \( M_e^2 / M_e^2 \) and we have

2.10 \[ KE_e = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^4 M_{\odot}}{2 \, r_e^2 \, v_m^2 \, M_e^2} \]

The Planck constant for the Sun, \( \hbar_{\odot} \), we will call \( L_p \), the subscript \( p \) for Planck. We have

\[ L_p = r_e v_m M_e = (1.496 \times 10^{11} \, \text{m})(1022 \, \text{m/s})(5.972 \times 10^{24} \, \text{kg}) = 9.13 \times 10^{38} \, \text{kg} \cdot \frac{\text{m}^2}{\text{s}} \]

\[ L_p^2 = r_e^2 v_m^2 M_e^2 = 7.4483 \times 10^{77} \, \text{J} \cdot \text{m}^2 \cdot \text{kg} = 8.3367 \times 10^{77} \, \text{kg}^2 \cdot \frac{\text{m}^4}{\text{s}^2} \]

We write for the solution of the Earth/Sun system:

2.11 \[ KE_e = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^4 M_{\odot}}{2 L_p^2} \]

We can write 2.11 as

2.12 \[ KE_e = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^4 M_{\odot}}{2 \hbar_{\odot}^2} \]

Where we say

\( \hbar_{\odot} = 9.13 \times 10^{38} \, \text{J} \cdot \text{s} \)

\( h_{\odot} = 2\pi \hbar_{\odot} = 5.7365 \times 10^{39} \, \text{J} \cdot \text{s} \)

Let us see how accurate our equation is:

\[ KE_e = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^4 M_{\odot}}{2 L_p^2} \]

\[ = \frac{R_{\odot}}{R_m} \frac{(6.67408 \times 10^{-11})^2 (5.972 \times 10^{24} \, \text{kg})^4 (1.9891 \times 10^{30} \, \text{kg})}{2 (8.3367 \times 10^{77} \, \text{kg}^2 \cdot \frac{\text{m}^4}{\text{s}^2})} \]

\[ = \frac{R_{\odot}}{R_m} (6.759 \times 10^{30} \, \text{J}) \]

\[ \frac{R_{\odot}}{R_m} = \frac{6.957 \times 10^8 \, \text{m}}{1,737,400 \, \text{m}} = 400.426 \]

\[ KE_e = 2.70655 \times 10^{33} \, \text{J} \]

We have that the kinetic energy of the Earth is

\( KE_{\text{earth}} = \frac{1}{2}(5.972 \times 10^{24} \, \text{kg})(30,290 \, \text{m/s})^2 = 2.7396 \times 10^{33} \, \text{J} \)

Our equation has an accuracy of

\[ \frac{2.70655 \times 10^{33} \, \text{J}}{2.7396 \times 10^{33} \, \text{J}} = 98.79\% \]

Which is very good.

Let us equate the lunar and solar formulations:

\[ KE_e = \sqrt{n} \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^2 M_m^3}{2 \hbar_{\odot}^2} \]

\[ KE_e = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^4 M_{\odot}}{2 \hbar_{\odot}^2} \]

\[ \sqrt{3} \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^2 M_m^3}{2 \hbar_{\odot}^2} = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^4 M_{\odot}}{2 L_p^2} \]

This gives:

2.13 \[ L_p = \sqrt{\frac{M_e^2 M_{\odot}}{M_m^3 \sqrt{3}}} \cdot \hbar_{\odot} \]

We remember that

\( \hbar_{\odot} = (hC) \, KE_e \)

\( hC = 1 \, \text{second} \)

And since, \( KE_e = \frac{1}{2} M_e v_e^2 \)

2.14 \[ 2v_m = \frac{v_e^2}{r_e} (1 \, \text{second}) \sqrt{\frac{M_e^2 M_{\odot}}{M_m^3 \sqrt{3}}} \]

\[ \sqrt{\frac{M_e^2 M_{\odot}}{M_m^3 \sqrt{3}}} = \sqrt{\frac{(5.972 \times 10^{24} \, \text{kg})^2 (1.9891 \times 10^{30} \, \text{kg})}{(7.34763 \times 10^{22} \, \text{kg})^3 (1.732)}} = 321,331.459 \approx 321,331 \]

Equation 2.14 becomes

2.15 \[ 1 \, \text{second} = 2 r_e \frac{v_m}{v_e^2} \sqrt{\frac{M_m^3 \sqrt{3}}{M_e^2 M_{\odot}}} \]

The condition of a perfect eclipse gives us another expression for the base unit of a second. \( L_p \) is another version of the Planck Constant, which is intrinsic to the solar formulation as opposed to the lunar formulation. We want to see what the ground state looks like and what its characteristic time is, if it is 1 second like it is for the lunar formulation. Looking at the equation for energy:

\[ KE_e = \frac{R_{\odot}}{R_m} \cdot \frac{G^2 M_e^4 M_{\odot}}{2 L_p^2} \]

We see the ground state should be:

2.16 \[ \frac{L_p^2}{G M_e^2 M_{\odot}} \cdot \frac{\sqrt{3}}{c} = 1 \, \text{second} \]

And, it is equal to 1 second. You will notice where in the derivation for the energy we lost \( \sqrt{n} = \sqrt{3} \), we have to put it in the ground state equation. The computation is:

\[ \frac{(9.13 \times 10^{38} \, \text{J} \cdot \text{s})^2}{(6.674 \times 10^{-11})(5.972 \times 10^{24} \, \text{kg})^2 (1.989 \times 10^{30} \, \text{kg})} \cdot \frac{\sqrt{3}}{c} = 1.0172 \, \text{seconds} \]

3.0 Jupiter and Saturn

We want to find what the wave equation solutions are for Jupiter and Saturn because they significantly carry the majority of the mass of the solar system and thus should embody most clearly the dynamics of the wave solution to the Solar System. We also show here how well the solution for the Earth works, which is 99.5% accuracy.

I find that as we cross the asteroid belt leaving behind the terrestrial planets, which are solid, and go to the gas giants and ice giants, the atomic number is no longer squared and the square root of the orbital number moves from the numerator to the denominator. I believe this is because the solar system here should be modeled in two parts, just as it is in theories of solar system formation because there is a force other than just gravity of the Sun at work, which is the radiation pressure of the Sun, which is what separates it into two parts, the terrestrial planets on this side of the asteroid belt and the gas giants on the other side of the asteroid belt. The effect the radiation pressure has is to blow the lighter elements out beyond the asteroid belt when the solar system forms, which are gases such as hydrogen and helium, while the heavier elements are too heavy to be blown out from the inside of the asteroid belt, allowing for the formation of the terrestrial planets Venus, Earth, and Mars. The result is that our equation has the atomic number of the heavier metals such as calcium for the Earth, while the equation for the gas giants has the atomic numbers of the gases. We write for these planets

\[ E = \frac{Z}{\sqrt{n}} \frac{G^2 M^2 m^3}{2 \hbar_{\odot}^2} \]

So, for Jupiter we have (And again using the maximum orbital velocity which is at perihelion):

\( KE_j = \frac{1}{2}(1.89813 \times 10^{27} \, \text{kg})(13720 \, \text{m/s})^2 = 1.7865 \times 10^{35} \, \text{J} \)

\[ E = \frac{Z_H}{\sqrt{5}} \frac{(6.67408 \times 10^{-11})^2 (1.89813 \times 10^{27} \, \text{kg})^2 (7.347673 \times 10^{22} \, \text{kg})^3}{2 (2.8314 \times 10^{33})^2} \]

\[ E = \frac{Z_H}{\sqrt{5}} (3.971 \times 10^{35} \, \text{J}) = Z_H (1.776 \times 10^{35} \, \text{J}) \]

\[ Z_H = \frac{1.7865 \times 10^{35} \, \text{J}}{1.776 \times 10^{35} \, \text{J}} = 1.006 \, \text{protons} \approx 1.0 \, \text{protons} = \text{hydrogen (H)} \]

Jupiter is mostly composed of hydrogen gas, and secondly helium gas, so it is appropriate that \( Z = Z_H \).

Our equation for Jupiter is

\[ E_5 = \frac{Z_H}{\sqrt{5}} \frac{G^2 M_j^2 M_m^3}{2 \hbar_{\odot}^2} \]

Where \( Z_H \) is the atomic number of hydrogen which is 1 proton, and \( \sqrt{n} = \sqrt{5} \) for the orbital number of Jupiter, \( n = 5 \). Now we move on to Saturn...

\( KE_S = \frac{1}{2}(5.683 \times 10^{26} \, \text{kg})(10140 \, \text{m/s})^2 = 2.92 \times 10^{34} \, \text{J} \)

\[ E = \frac{Z}{\sqrt{6}} \frac{(6.67408 \times 10^{-11})^2 (5.683 \times 10^{26} \, \text{kg})^2 (7.347673 \times 10^{22})^3}{2 (2.8314 \times 10^{33})^2} \]

\[ = \frac{Z}{2.45} (3.5588 \times 10^{34} \, \text{J}) = Z (1.45259 \times 10^{34} \, \text{J}) \]

\[ Z (1.45259 \times 10^{34} \, \text{J}) = (2.92 \times 10^{34} \, \text{J}) \]

\[ Z = 2 \, \text{protons} = \text{Helium (He)} \]

The equation for Saturn is then

\[ E_6 = \frac{Z_{\text{He}}}{\sqrt{6}} \frac{G^2 M_s^2 M_m^3}{2 \hbar_{\odot}^2} \]

It is nice that Saturn would use Helium in the equation because Saturn is the next planet after Jupiter and Jupiter uses hydrogen, and helium is the next element after hydrogen. As well, just like Jupiter, Saturn is primarily composed of hydrogen and helium gas.

The accuracy for Earth orbit is

\[ E_n = \sqrt{n} \, \frac{R_{\odot}}{R_m} \frac{G^2 M_e^2 M_m^3}{2 \hbar_{\odot}^2} \]

\[ \frac{R_{\odot}}{R_m} = \frac{6.96 \times 10^8 \, \text{m}}{1,737,400 \, \text{m}} = 400.5986 \]

\[ E_3 = (1.732)(400.5986) \frac{(6.67408 \times 10^{-11})^2 (5.972 \times 10^{24} \, \text{kg})^2 (7.347673 \times 10^{22} \, \text{kg})^3}{2 (2.8314 \times 10^{33})^2} = 2.727 \times 10^{33} \, \text{J} \]

The kinetic energy of the Earth is

\( KE_e = \frac{1}{2}(5.972 \times 10^{24} \, \text{kg})(30,290 \, \text{m/s})^2 = 2.7396 \times 10^{33} \, \text{J} \)

\[ \frac{2.727 \times 10^{33} \, \text{J}}{2.7396 \times 10^{33} \, \text{J}} \times 100 = 99.5\% \]

Which is very good, about 100% accuracy for all practical purposes. The elemental expression of the solution for the Earth would be

\[ E_3 = \sqrt{3} \, \frac{Z_{Ca}^2 \, G^2 M_e^2 M_m^3}{2 \hbar_{\odot}^2} \]

Where \( \frac{R_{\odot}}{R_m} \rightarrow Z^2 \).

In this case the element associated with the Earth is calcium which is \( Z = 20 \) protons.

References

Beardsley, I. (2025) Theory For The Solar System And The Atom's Proton; Linking Microscales To Macroscales.

Beardsley, I. (2026) How Physics and Archaeology Point to a Natural Constant of 1-Second.

Beardsley, I. (2026) The Sublime and Mysterious Place of Humans in the Cosmos; A Work in Exoarchaeology.