Skip to content

Phase 2: measuring gravity on the ground and in tests of the equivalence principle #3

Closed

Description

@dan-garden

Report section

This phase gives an absolute gravimeter its tide correction, gives the atom sensor its ground mode, and gives gravitics an Eötvös calculator with real limits.

  • Item 8: metrology::{GM_MOON, EARTH_ROTATION} (defect D9)
    • Formula or algorithm: Constants with named conventions
    • Primary source: IERS Conventions 2010 (TN36)
    • Test value to reproduce: The existing derived-constant test; longman keeps its own cgs constants so its vector holds exactly
    • Joins: every crate
  • Item 9: gravimetry::tides::longman(lat, lon, h, t) → (moon, sun, total)
    • Formula or algorithm: Longman's closed form from mean orbital elements, times the gravimetric factor 1 + h₂ − (3/2)k₂
    • Primary source: Longman 1959
    • Test value to reproduce: 40.7914°, 282.1414°E, 370 m, 2015-04-23 00:00 UT → moon 0.0324029651226, sun −0.0288682178454, total 0.00353474727722 mGal (LongmanTide tests); reproduced to < 10⁻⁹ mGal
    • Joins: gravimetry, metrology
  • Item 10: gravimetry::tides::{delta, gamma}
    • Formula or algorithm: δₙ = 1 + (2/n)hₙ − ((n+1)/n)kₙ; γₙ = 1 + kₙ − hₙ
    • Primary source: Mäkinen 2021
    • Test value to reproduce: Conventional δ = 1.16 and 1 + k − h ≈ 0.69 (Mäkinen); δ₂(0.612, 0.303) = 1.1575 with Longman's own Love numbers
    • Joins: gravimetry ↔ gravitational::tides (k₂)
  • Item 11: gravimetry::tides::permanent(lat, h), reusing geopotential::TideSystem
    • Formula or algorithm: W_T = A(r/a)²(sin²φ − 1/3) with A = −2.9166 m²s⁻²; station gravity and geoid offsets
    • Primary source: Mäkinen 2021, Eqs. 18–20
    • Test value to reproduce: g_T(φ) = −30.49 + 90.95 sin²φ + 0.31 sin⁴φ µGal; spherical form −30.485 µGal at the equator; mean-minus-zero-tide geoid 9.9 cm at the equator
    • Joins: gravimetry ↔ geopotential
  • Item 12: Diffraction order in gravimetry::atom (defect D2), with the gradiometer added to examples/buried_tunnel.rs, which already uses Stray's geometry
    • Formula or algorithm: Δφ = n_d k_eff ΔΓ L T², with n_d = 1 on the ground and 2 in microgravity
    • Primary source: Stray et al. 2022
    • Test value to reproduce: 150 E over 1 m at T = 85 ms → 17.5 mrad (computed 17.45); projection limit ≈ 44 E/√Hz (computed 45.9); 466 E/√Hz reaches 20 E "within 10 min" (averaging_time gives 543 s)
    • Joins: gravimetry ↔ gravitational::detectors
  • Item 13: gravitics::equivalence::{eta, combine, coupling_bound, sphere_form_factor} (identifiers say eta, to avoid the unit gravimetry::local::EOTVOS)
    • Formula or algorithm: η = 2(a_A − a_B)/(a_A + a_B); δ = ∓(α/(G_N u²))ε²(Q/A_r)_⊕Δ(Q/A_r); φ(x) = 3(x cosh x − sinh x)/x³
    • Primary source: Touboul et al. 2022; Fayet 2025
    • Test value to reproduce: MICROSCOPE η = (−1.5 ± 2.3 ± 1.5)×10⁻¹⁵ (quadrature 2.746×10⁻¹⁵); |g_{B−L}| < 1.1×10⁻²⁵ for spin 1 (computed 1.067×10⁻²⁵) and < 1.3×10⁻²⁵ for spin 0 (1.283×10⁻²⁵); α/(G u²) = 1.2536×10³⁶
    • Joins: gravitics::modified (yukawa_factor, boson_mass_ev), metrology
  • Item 14: Test vectors for geo5d (outside physica): existing fractional_frequency_shift, geodetic_precession and shapiro_delay; the planned perihelion_advance
    • Formula or algorithm: —
    • Primary source: Bothwell et al. 2022; Zheng et al. 2023; Will 2014; GRAVITY 2020; González Hernández et al. 2020
    • Test value to reproduce: JILA −9.8(2.3)×10⁻²⁰ per mm against g/c² = 1.090×10⁻¹⁹; Wisconsin (−12.4 ± 0.7 ± 2.5)×10⁻¹⁹ per cm against −10.9; Cassini geometry 262.4 µs round trip (twice shapiro_delay); lunar de Sitter 1.919″ per century ("about 2″"); S2 Schwarzschild precession 12.21′ per orbit ("12′"); solar redshift 633.14 m/s against 638 ± 6
    • Joins: geo5d ↔ geopotential

Done when each item's test reproduces its published value or identity, the module is in the crate's README table, and the check scripts still pass. Values in the verify-before-coding issue are not used as tests until their source is read.


From the cross-domain frontiers report. Test values come from the check scripts in notes/; every one passed on 2026-10-25.

Activity

  1. dan-garden commented on Oct 25, 2026

    @dan-garden
    CollaboratorAuthor

    Phase 2 is in #16, with item 14 in FusionVault/geo5d#1. These are the findings beyond the report.

    Tides (items 9–11)

    • LongmanTide's vector reproduces to 3×10⁻¹¹ mGal. longman takes a Julian epoch, which holds the time to about 4 µs, and that sets the residual.
    • Averaged over the 18.61-year nodal cycle, Longman's tide is δ₂ times the permanent tide. The match with Mäkinen's −g_T is 0.1% at the equator and poles (0.19 µGal worst, at 45°). This checks the time-domain tide against the spectral M0S0 average independently; it is the non-zero mean Honkasalo (1964) pointed out.
    • Two quirks in Longman as coded:
      • The degree-2 factor 1.1575 also scales the degree-3 lunar term, whose own δ₃ is 1.07.
      • His cgs masses make the lunar tide 0.09% and the solar 0.2% larger than today's GMs.
      • Both are kept so the vector holds, and the docs say so.
    • The permanent tide is exact at any height in Cartesian form, W_T = (A/3a²)(2z² − p²).
      • It matches Mäkinen's Eqs. 19–21 to their printed digits.
      • Eq. 19's height factor 1 + 2h/a is good to 4×10⁻⁶ at 1 km.
      • A = −2.9166 follows from IERS H₀ = −0.31460 m via his Eq. 16 (−2.91663).
    • Defect, now documented: a mean-tide C̄20 gives the wrong gravity. TideSystem::MeanTide puts W_T into the Earth's own field, falling as r⁻³ instead of growing as r². That is right only on r = R ("there is no mean-tide C20", Mäkinen). Its gravity is +45.7 µGal at the equator against −30.49, a factor −3/2, and a test pins it.
    • Small offsets in the C̄20 conventions.
      • The EGM direct term −3.11080×10⁻⁸/√5 differs from IERS's A₀H₀ = −1.39141×10⁻⁸ by 1.7×10⁻⁴.
      • TN36 quotes the −4.1736×10⁻⁹ offset for EGM2008. Its own Eq. 6.14 with k₂₀ = 0.30190 gives −4.2007×10⁻⁹.

    Gradiometer (item 12)

    • The tunnel signal comes out close to Stray's. Over the example's tunnel the 1 m-baseline difference is 152 E and 17.7 mrad, against Stray's 150 E and 17.5 mrad.
    • Over a line mass the baseline difference is exactly the point gradient at the geometric-mean depth √(z₁z₂).
    • Averaging time: at the 466 E/√Hz reached outdoors, a 10:1 detection takes 937 s per station; at the projection limit it takes 9 s.

    Equivalence (item 13)

    • Fayet doesn't print the Earth's charges per mass, but Eq. 7's coefficients give them:
      • (B − L)/A_r = 0.5138 and L/A_r = 0.4870, which sum to 1.0008 nucleons per u.
      • B's coefficient has only three digits and Δ(B/A_r) only two, so its value of 1.010 is good to 1%.
    • All six bounds reproduce. The bounds for B, L and B − L at both spins match his Eqs. 24–25. For B − L they match note 40's unrounded 1.07 and 1.28×10⁻²⁵, which differ by exactly √(6.5/4.5).
    • The form factor matches Fayet's tables, through the existing functions. φ(x) matches Table IV to its printed digits. Eq. 80's five shells, with yukawa_factor and boson_mass_ev, reproduce Table V's F and G at 710 km for m = 10⁻¹⁴ to 2×10⁻¹³ eV.
    • The closed-form φ(x) fails near zero. It returns −4.96 at x = 10⁻⁸ and loses 8 digits at 10⁻⁴, so it is summed as a series below x = 1.
    • The registry isn't ordered by uncertainty. MICROSCOPE's Pt–Pt control has a smaller combined error (2.55×10⁻¹⁵) than Ti–Pt (2.75×10⁻¹⁵).

    geo5d (item 14)

    • Mercury: JPL's J2000 elements give 42.980″ per century against Will's 42.98″ (Eq. 65).
    • S2: with GRAVITY 2020's own Table E.1 it is 12.16′ against their 12.1′. The 12.21′ above used later GRAVITY elements.
    • Cassini: the exact round trip is 262.42 µs, 0.5% under Will's Eq. 63 shorthand (263.84).
    • The solar redshift paper has a sign slip. González Hernández et al. call the Earth-distance correction "+4 cm s⁻¹", but their 633.10 m/s is 633.35 − 0.04 − 0.21. The correction lowers the value.
  2. added a commit that references this issue on Oct 25, 2026
  3. dan-garden commented on Oct 25, 2026

    @dan-garden
    CollaboratorAuthor

    Released 2026-10-25 as compatible z-bumps: geopotential 0.3.5, gravimetry 0.7.1, gravitics 0.5.5 and physica 0.3.4 (tagged <crate>-vX.Y.Z), plus geo5d 0.5.2 (v0.5.2) for item 14. All are confirmed on crates.io.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    enhancementNew feature or request

    Type

    No type

    Projects

    No projects

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions