The Earth-ionosphere cavity supports global electromagnetic standing waves driven by lightning — the Schumann resonances, with fundamental modes near 7.8, 14.3 and 20.8 Hz. Because a hypothetical photon rest mass would shift these eigenfrequencies relative to Maxwell electrodynamics, the observed resonance spectrum yields an upper bound on the photon mass. Kroll (1971) pioneered this idea with perfect-conductor shells; Sucu and Tekincay (2019) reformulated it with a massive spin-1 (Zitterbewegung/Duffin-Kemmer-Petiau-family) wave equation, obtaining mγ<1.3×10−50m_\gamma < 1.3\times10^{-50} kg, but with phenomenological cavity quality factors, a thin-shell geometry, and a single-mode fit. The open problem is to place this bound on a first-principles footing: derive the cavity eigenfrequencies and quality factors as functions of the photon mass from realistic ionospheric conductivity — including laterally non-uniform, time-varying structure and boundary impedance — and extract the bound from multi-mode observational data with error propagation. Malta and Helayël-Neto (2022) completed the spherically symmetric version of this program in Proca theory (mγ≤2.5×10−14m_\gamma \le 2.5\times10^{-14} eV/c2^2 conservative); day-night asymmetry, geomagnetic-field effects, source-distribution modeling in amplitude spectra, and extension to planetary cavities remain open.

Background

Schumann resonances (SR) are the global electromagnetic standing waves of the cavity formed by the Earth's surface and the lower ionosphere, excited by worldwide lightning. The resonant modes are extremely low frequency (ELF, here ~3 Hz to ~30 Hz) transverse-magnetic (TM) waves — modes whose magnetic field has no radial component — observed near 7.8, 14.3, 20.8 and 27 Hz. The ionospheric wall is a poor conductor at these frequencies, so each mode has a finite quality factor QQ (the ratio of resonance frequency to linewidth, Q≈4Q\approx 4–66 for the first modes), set by how ELF waves leak into and dissipate in the ionosphere. Classical modeling of SR with realistic, height-dependent atmospheric conductivity profiles is a mature subject (e.g., Mushtak and Williams 2002; Greifinger, Mushtak and Williams 2007).

A massive photon is described by the de Broglie-Proca equations — Maxwell's equations augmented by a mass term, equivalent at the level of wave propagation to replacing the dispersion relation ω=ck\omega = ck by ω2=c2k2+mγ2c4/ℏ2\omega^2 = c^2k^2 + m_\gamma^2c^4/\hbar^2, and by a length scale λγ=ℏ/(mγc)\lambda_\gamma = \hbar/(m_\gamma c) beyond which Coulomb fields are exponentially screened. Because λγ\lambda_\gamma for any accessible mass is enormous, small-mass tests use the largest available systems: laboratory limits are around 10−1410^{-14} eV/c2^2 (Williams, Faller and Hill 1971), fast radio bursts (FRBs — millisecond-long, apparently extragalactic radio transients) constrain mγ≲10−14m_\gamma \lesssim 10^{-14} eV/c2^2 (Bonetti et al. 2017; Shao and Zhang 2017), and solar-wind magnetohydrodynamics (plasma modeling of the solar wind's magnetic field) reaches roughly 10−1810^{-18} eV/c2^2 (reviewed by Goldhaber and Nieto 2010). Kroll (1971) noticed that the Earth-ionosphere cavity, whose TM eigenfrequencies are of order 10 Hz, could be sensitive to mγ∼10−14m_\gamma \sim 10^{-14} eV/c2^2, and derived the bound mγ<4.2×10−49m_\gamma < 4.2\times10^{-49} kg by solving the massive scalar (Klein-Gordon) equation between perfectly conducting spherical shells with a mass-sensitivity coefficient and a nominal 70 km ionosphere height.

Sucu and Tekincay (2019) re-examined this with a spin-1 relativistic wave equation derived from Barut's Zitterbewegung model — a first-order equation for massive or massless spin-1 bosons that coincides with the spin-1 Duffin-Kemmer-Petiau (DKP) equation and reproduces Maxwell theory in the massless limit — solved in an annular cavity with poorly conducting walls. They recovered the classical SR frequencies in the massless case and, for the massive case, the bound m0<1.3×10−50m_0 < 1.3\times10^{-50} kg from the experimentally observed first mode (ν1=7.8\nu_1 = 7.8 Hz) and the first-mode quality factor. Their cavity damping is encoded in a phenomenological 'regularization parameter' ϵj\epsilon_j through Qj=ln⁡(1+fj/ϵj)Q_j = \ln(1 + f_j/\epsilon_j), with ϵj\epsilon_j fitted to reproduce observed frequencies rather than derived from ionospheric electrodynamics; the geometry uses a fixed ionosphere radius bb with the thin-shell approximation (b−a)/a≪1(b-a)/a \ll 1 (aa the Earth radius); and the bound uses a single mode. These are the specific approximations the open problem asks to remove.

Subsequently, Malta and Helayël-Neto (2022) carried out the analogous program in standard de Broglie-Proca theory: they solved the massive-vector wave equation with realistic spherically stratified atmospheric conductivity profiles (two numerical aeronomical profiles plus Cole's analytic profile — a closed-form fit to measured atmospheric conductivity — all exhibiting the sharp conductivity rise, or 'knee', near 60 km that effectively bounds the cavity), computed the eigenfrequencies and quality factors of the first three modes as functions of mγm_\gamma, and compared with observed values and their variabilities, obtaining mγ≤1.4×10−14m_\gamma \le 1.4\times10^{-14} eV/c2^2 (optimal) and ≤2.5×10−14\le 2.5\times10^{-14} eV/c2^2 (conservative), about ten times tighter than Kroll's estimate. Notably, the nominal Sucu-Tekincay limit (1.3×10−501.3\times10^{-50} kg ≈7.3×10−15\approx 7.3\times10^{-15} eV/c2^2) is numerically smaller, but it rests on a fitted phenomenological QQ rather than a controlled conductivity model. What remains unmodeled in every treatment to date: the day-night asymmetry of the ionosphere and the geomagnetic field (requiring two- or three-dimensional numerics such as finite-difference time-domain (FDTD) analysis, which steps the wave equations forward in time on a discrete spatial grid), the inclusion of mγm_\gamma in theoretical amplitude spectra built from realistic global lightning source distributions, and the extension to other planets, where Jupiter's roughly eleven-times-larger cavity could push sensitivity toward mγ≲10−15m_\gamma \lesssim 10^{-15} eV/c2^2.

Problem Statement

Establish a robust upper bound on the photon rest mass from Schumann resonances of the Earth-ionosphere cavity on a first-principles footing. The spherically symmetric, static-conductivity instance of this program has already been completed in Proca theory (Malta and Helayël-Neto 2022, mγ≤2.5×10−14m_\gamma \le 2.5\times10^{-14} eV/c2^2 conservative); the open problem is the corresponding analysis with a realistic laterally inhomogeneous, time-varying ionosphere. Specifically: derive the cavity's TM eigenfrequencies fℓf_\ell and quality factors QℓQ_\ell (ℓ=1,2,3,…\ell = 1, 2, 3, \ldots) as functions of the photon mass mγm_\gamma from massive-vector electrodynamics (Proca or equivalent spin-1 formulation) with a physically realistic ionosphere — the required core being a height-dependent, laterally inhomogeneous and time-varying (day-night) atmospheric/ionospheric conductivity profile, with inhomogeneous and frequency-dependent boundary impedance and geomagnetic-field structure desirable but not required, rather than perfect-conductor shells, fixed-radius thin-shell geometry (b−a)/a≪1(b-a)/a\ll1, or phenomenologically fitted damping parameters — and, by comparing these massive-model predictions with multi-mode observational resonance data including error propagation for both observational variability and modeling uncertainty, extract the tightest defensible limit on mγm_\gamma. If the realistic analysis does not improve on the published conservative bound, reporting the limit the analysis itself supports, together with a quantitative demonstration that the added realism leaves that bound essentially unchanged, also resolves the task. The same analysis could also quantify how residual geophysical variability (diurnal and seasonal ionospheric changes, global lightning source distributions) limits the attainable sensitivity, and, as an optional planetary extension of the same framework, assess what bound larger cavities such as Jupiter's could provide given remotely sensed Schumann spectra.

References

  1. Y. Sucu, C. Tekincay: Photon in the Earth-ionosphere cavity: Schumann resonances, Astrophysics and Space Science 364, 56 (2019). DOI 10.1007/s10509-019-3547-7, arXiv:1803.10685. https://doi.org/10.1007/s10509-019-3547-7
  2. P. C. Malta, J. A. Helayël-Neto: Constraining the photon mass via Schumann resonances, Physical Review D 106, 116014 (2022). DOI 10.1103/PhysRevD.106.116014, arXiv:2208.00527. https://doi.org/10.1103/PhysRevD.106.116014
  3. N. M. Kroll: Concentric Spherical Cavities and Limits on the Photon Rest Mass, Physical Review Letters 27, 340 (1971). DOI 10.1103/PhysRevLett.27.340. https://doi.org/10.1103/PhysRevLett.27.340
  4. E. R. Williams, J. E. Faller, H. A. Hill: New Experimental Test of Coulomb's Law: A Laboratory Upper Limit on the Photon Rest Mass, Physical Review Letters 26, 721 (1971). DOI 10.1103/PhysRevLett.26.721. https://doi.org/10.1103/PhysRevLett.26.721
  5. A. S. Goldhaber, M. M. Nieto: Photon and graviton mass limits, Reviews of Modern Physics 82, 939 (2010). DOI 10.1103/RevModPhys.82.939. https://doi.org/10.1103/RevModPhys.82.939
  6. L. Shao, B. Zhang: Bayesian framework to constrain the photon mass with a catalog of fast radio bursts, Physical Review D 95, 123010 (2017). DOI 10.1103/PhysRevD.95.123010. https://doi.org/10.1103/PhysRevD.95.123010
  7. L. Bonetti, J. Ellis, N. E. Mavromatos, A. S. Sakharov, E. K. Sarkisyan-Grinbaum, A. D. A. M. Spallicci: FRB 121102 casts new light on the photon mass, Physics Letters B 768, 326 (2017). DOI 10.1016/j.physletb.2017.03.014. https://doi.org/10.1016/j.physletb.2017.03.014
  8. V. C. Mushtak, E. R. Williams: ELF propagation parameters for uniform models of the Earth-ionosphere waveguide, Journal of Atmospheric and Solar-Terrestrial Physics 64, 1989 (2002). DOI 10.1016/S1364-6826(02)00222-5. https://doi.org/10.1016/S1364-6826(02)00222-5
  9. P. S. Greifinger, V. C. Mushtak, E. R. Williams: On modeling the lower characteristic ELF altitude from aeronomical data, Radio Science 42, RS2S12 (2007). DOI 10.1029/2006RS003500. https://doi.org/10.1029/2006RS003500