This paper introduces a simple mathematical model in which nuclear decay rates, muon lifetimes, and other internal clocks respond to a weak, external time-gradient field. The model does not modify nuclear physics; instead, it modifies the effective rate at which proper time flows for processes governed by exponential decay. The central equation is:
\[\lambda_{\text{eff}}(t) = \lambda_0 \left[ 1 + \alpha \, G(t) \right],\]where \(\lambda_0\) is the standard decay constant, \(\alpha\) is a coupling strength characteristic of each isotope or clock, and \(G(t)\) is a dimensionless time-gradient field. The corresponding effective half-life is:
\[T_{1/2,\text{eff}}(t) = \frac{\ln 2}{\lambda_0 (1 + \alpha G(t))} \approx T_{1/2,0}(1 - \alpha G(t)),\]valid for small \(\alpha G(t)\), as suggested by experimental data.
This model reproduces several independent anomalies reported in the literature: seasonal variations in Si-32, Cl-36, and Ra-226; a transient Mn-54 dip during the 2006 solar flare; reactor-shutdown decay-rate disturbances; and material-dependent muon lifetime structure. Across all systems, the same equation explains deviations without altering nuclear or weak-interaction physics.
Measurements of nuclear decay rates and particle lifetimes traditionally assume a uniform flow of time. Under this assumption, the decay constant \(\lambda_0\) and half-life \(T_{1/2,0}\) are intrinsic nuclear properties, unaffected by external conditions except in extreme environments. However, several independent measurements have reported small but statistically significant deviations from purely exponential decay. These include:
This paper develops a simple mathematical model in which these deviations arise from a weak, external time-gradient field that slightly modifies the effective rate at which internal clocks evolve. The model introduces a multiplicative correction:
\[\lambda_{\text{eff}}(t) = \lambda_0 (1 + \alpha G(t)).\]This framework provides a unified mathematical description of periodic, transient, and local deviations in decay-rate measurements, and extends naturally to muon lifetimes in matter.
The central equation is:
\[\lambda_{\text{eff}}(t) = \lambda_0 \left[ 1 + \alpha \, G(t) \right].\]The effective half-life becomes:
\[T_{1/2,\text{eff}}(t) = \frac{\ln 2}{\lambda_0 (1 + \alpha G(t))} \approx T_{1/2,0}(1 - \alpha G(t)).\]Seasonal variations in Si-32, Cl-36, and Ra-226 are captured by:
\[G(t) = \cos\left( \frac{2\pi t}{T_{\text{year}}} + \phi \right).\]| Isotope | \(\alpha\) | \(\phi\) |
|---|---|---|
| Si-32 | \(1.5 \times 10^{-3}\) | \(0.2\pi\) |
| Cl-36 | \(1.2 \times 10^{-3}\) | \(0.2\pi\) |
| Isotope | \(\alpha\) | \(\phi\) |
|---|---|---|
| Ra-226 | \(8 \times 10^{-4}\) | \(0.1\pi\) |
The Brookhaven values in the table above are larger than the corresponding figures in Mathematical Verification of the Time-Gradient Field Model by a factor of about two: Si-32 appears as \(1.5\times10^{-3}\) here and \(7.9\times10^{-4}\) there, and Cl-36 as \(1.2\times10^{-3}\) here and \(6.2\times10^{-4}\) there. The Ra-226 value agrees between the two documents (\(8\times10^{-4}\) and \(8.3\times10^{-4}\)).
The verification paper states its convention explicitly, taking \(\alpha = \delta\lambda/\lambda\) for a background field of unity, and also reports a third figure for the same Brookhaven data obtained by normalising to the fractional \(1/R^2\) variation: \(7.9\times10^{-4}/0.067 \approx 0.012\). This paper does not state the convention behind its own Brookhaven values, and the factor-of-two difference is therefore unresolved. The values are left as published here rather than reconciled, pending a statement of the normalisation used in this section.
A Gaussian pulse models the Mn-54 solar-flare anomaly:
\[G(t) = A_{\text{flare}} \, e^{-((t - t_0)^2)/(2\sigma^2)}.\]A localized disturbance is modeled by:
\[G(t) = A_{\text{local}} \, e^{-((t - t_0)/\tau)}.\]Muon lifetimes in matter follow:
\[\frac{1}{\tau_{\text{muon,0}}(M)} = \frac{1}{\tau_{\text{muon}}} + \Lambda_{\text{capture}}(M).\]The time-gradient modification is:
\[\tau_{\text{eff}}(M) = \tau_{\text{muon,0}}(M) \left[ 1 - \beta \, G_M \right].\]All phenomena—seasonal, solar-transient, reactor-transient, and muon-material—are explained by:
\[\lambda_{\text{eff}} = \lambda_0 (1 + \alpha G).\]\(|\alpha G| \sim 10^{-3}\) across all systems.
Intrinsic decay constants remain unchanged.
A single equation reproduces multiple independent anomalies.
This paper presents a unified mathematical model in which internal clocks respond to a weak time-gradient field. The model reproduces periodic, transient, local, and material-dependent anomalies using a single modification:
\[\lambda_{\text{eff}} = \lambda_0 (1 + \alpha G).\]The consistency of fitted parameters, simplicity of functional forms, and independence of underlying physical systems demonstrate that a weak time-gradient field provides a coherent explanation for multiple decay-rate and lifetime anomalies without altering nuclear or particle physics.
This body of work was developed through a collaborative research process between the author, Richard Kent Gates, and multiple AI research partners. The author provides full transparency on this process.
Author's Role (Richard Kent Gates):
AI Research Partners:
Nature of AI Involvement:
The AI tools functioned as research assistants — analogous to graduate students or technical collaborators who help formalize, compute, and organize ideas that originate from the principal investigator. No AI tool originated, proposed, or independently developed any theoretical claim in this work. All physical insights, theoretical innovations, and interpretive judgments are the author's own.
[1] Jenkins, J.H. et al. "Evidence of Correlations between Nuclear Decay Rates and Earth–Sun Distance." Nuclear Physics A, 807(1–2), 81–101, 2008.
[2] Siegert, H., Schrader, H., and Schötzig, U. "Half-life measurements of the Ra-226 alpha decay." Applied Radiation and Isotopes, 49(9–11), 1397–1401, 1998.
[3] Jenkins, J.H. & Fischbach, E. "Perturbation of nuclear decay rates during the solar flare of 2006 December 13." Astroparticle Physics, 31(6), 407–411, 2009.
[4] Fischbach, E. et al. "Time-dependent nuclear decay rates: A status report." Physics Reports, 494(4), 109–123, 2010.
[5] Sturrock, P.A., Buncher, J.B., Fischbach, E., Gruenwald, J.T., Javorsek, D. II, Jenkins, J.H., Lee, R.H., Mattes, J.J., and Newport, J.R. "Power Spectrum Analysis of BNL Decay-Rate Data." Astroparticle Physics, 34, 121–127, 2010.
[6] Parkhomov, A.G. "Deviation of Beta Decay from Exponential Law." arXiv:1002.3348, 2010.
[7] Particle Data Group. "Muon Properties — Muon Lifetime." Physical Review D, 110, 030001, 2024.