arxiv: 2605.13960 · v1 · submitted 2026-05-13 · 🌌 astro-ph.EP · astro-ph.IM· gr-qc

Recognition: 2 theorem links

· Lean Theorem

Gravitational-wave Tomography of the Moon: Constraining Lunar Structure with Calibrated Gravitational Waves

Han Yan , Jan Harms

Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:56 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.IMgr-qc

keywords gravitational waveslunar structureseismic tomographyelastic parametersnormal modesperturbation theoryMoon interiormodal amplitudes

0 comments

The pith

Calibrated gravitational waves reduce estimation errors of the Moon's elastic parameters by an order of magnitude through tomographic inversion of its seismic response.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a perturbative framework that treats the Moon as a resonant detector for mid-frequency gravitational waves whose strain is independently known from detector networks. It maps spherically symmetric changes in elastic moduli, density, and interface locations to shifts in normal-mode frequencies and amplitudes using first-order perturbation theory. This turns observed lunar seismic responses into an inverse problem for internal structure. The central result is that such calibrated observations shrink parameter uncertainties by roughly a factor of ten compared with uncalibrated seismic data alone. Readers should care because the approach offers a new geophysical tomography channel that complements existing lunar seismology.

Core claim

A first-principles perturbative framework maps spherically symmetric perturbations of elastic and density structure to measurable changes in GW-driven modal amplitudes. The formalism combines a normal-mode representation of the elastic response, first-order perturbation theory for eigenvalues and eigenfunctions, and a linearized observation model that links frequency and amplitude observables to bulk and shear moduli, density, and interface locations. Observations of calibrated GWs thereby reduce estimation errors of the Moon's elastic parameters by about an order of magnitude.

What carries the argument

Linearized observation model from normal-mode perturbation theory that converts changes in modal amplitudes and frequencies into constraints on perturbations of bulk modulus, shear modulus, density, and layer boundaries.

If this is right

Elastic-parameter uncertainties shrink by a factor of approximately ten when GW strain is calibrated independently.
Modal frequencies and amplitudes become direct observables for inverting spherically symmetric interior perturbations.
Interface depths and density profiles can be recovered jointly with bulk and shear moduli.
The same formalism applies to any spherically symmetric body whose seismic response to known external strains can be measured.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Future lunar landers with broadband seismometers could use passing gravitational waves as calibration signals to refine interior models without new hardware.
The method might extend to other solar-system bodies if their normal-mode responses to external strains become measurable.
Discrepancies between predicted and recorded responses could flag the breakdown of the spherical-symmetry or first-order assumptions.
Joint analysis with Apollo-era seismic data could provide an independent consistency check on existing lunar velocity models.

Load-bearing premise

The incoming gravitational-wave strain amplitude is known independently from networks of detectors, and the Moon's seismic response is captured accurately by first-order perturbation theory around a spherically symmetric background model.

What would settle it

Seismic recordings of a known gravitational-wave event would falsify the central claim if the observed modal amplitude shifts deviate from the model's predicted changes by more than the claimed order-of-magnitude error reduction.

Figures

Figures reproduced from arXiv: 2605.13960 by Han Yan, Jan Harms.

**Figure 2.** Figure 2: FIG. 2. Relative errors for [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

read the original abstract

The recent success of gravitational-wave (GW) astronomy together with renewed plans for lunar geophysical instrumentation has revived interest in using the Moon as a resonant detector for mid-frequency (mHz-Hz) GWs. In realistic observational scenarios, the GW strain amplitude is expected to be constrained independently by networks of GW detectors, which motivates an inverse, \emph{tomographic} question: to what extent can measurements of the Moon's seismic response to known GWs be used to infer its internal structure? In this work, we develop a first-principles, perturbative framework that maps spherically symmetric perturbations of the elastic and density structure to measurable changes in observables, especially GW-driven modal amplitudes of the Moon. The formalism combines (i) a normal-mode representation of the elastic response, (ii) first-order perturbation theory for eigenvalues and eigenfunctions, and (iii) a linearized observation model that links frequency and amplitude observables to model parameters (bulk and shear moduli, density, and interface locations) and their perturbations. We show that the estimation errors of the Moon's elastic parameters can be reduced by about an order of magnitude with observations of calibrated GWs.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper sketches a perturbative method to turn calibrated gravitational-wave strains into better constraints on the Moon's elastic structure, but the linear approximation may not be accurate enough for the Moon's actual internal contrasts.

read the letter

The main thing here is a new way to combine independent GW strain measurements with the Moon's seismic response to get tighter constraints on its internal elastic parameters. They set up a first-principles perturbative mapping from density and modulus perturbations to changes in normal mode frequencies and amplitudes driven by known GWs. The formalism combines normal-mode representation, first-order perturbation theory for eigenvalues and eigenfunctions, and a linearized observation model linking those to bulk and shear moduli, density, and interface locations. The central quantitative claim is that estimation errors drop by about an order of magnitude once the GW amplitude is known from external networks. That specific combination for lunar tomography is not something I've seen laid out before. What they do well is keep the setup clean and grounded in standard Rayleigh-Schrödinger perturbation on a spherically symmetric background, which makes the mapping from structure changes to observable shifts straightforward to follow. The motivation from upcoming lunar instruments and existing GW detector networks is also timely and practical. The soft spot is the lack of any check on whether first-order perturbation holds for realistic lunar contrasts. The Moon has 10-30% jumps in density and shear modulus at the crust-mantle and core-mantle boundaries, and the abstract gives no bound on perturbation size or synthetic tests showing when second-order terms stay below the claimed precision. Without those, the Fisher-matrix error reduction remains an unverified projection. The paper is aimed at people working at the intersection of GW astronomy and planetary geophysics who want new observational strategies for the Moon. A reader looking for fresh ideas on multi-messenger constraints would get value from the formalism even if the numbers need more work. I would send it to peer review because the core idea is fresh and the setup uses established methods, though it will need substantial revision to address the validation gaps.

Referee Report

2 major / 1 minor

Summary. The paper develops a first-principles perturbative framework for gravitational-wave tomography of the Moon. It combines a normal-mode representation of the elastic response, first-order Rayleigh-Schrödinger perturbation theory on a spherically symmetric background for eigenvalues and eigenfunctions, and a linearized observation model linking GW-driven modal amplitudes and frequencies to perturbations in bulk/shear moduli, density, and interface locations. The central claim is that independent calibration of GW strain amplitude by terrestrial detector networks allows an order-of-magnitude reduction in estimation errors for the Moon's elastic parameters.

Significance. If the quantitative error-reduction claim holds after validation, the work would provide a novel, high-precision tomographic tool for lunar interior structure that synergizes GW astronomy with planned lunar geophysical instrumentation. The parameter-free character of the calibration step and the explicit mapping from structure perturbations to observables are notable strengths.

major comments (2)

[§3] §3: The linearized observation model and Fisher-matrix error analysis rest on first-order perturbation theory for modal amplitudes. The Moon exhibits O(10–30 %) jumps in density and shear modulus at the crust-mantle and core-mantle boundaries; no explicit bound is derived on the perturbation size for which neglected second-order terms remain below the few-percent level required to support the claimed order-of-magnitude error reduction.
[Abstract, §4–5] Abstract and §4–5: The central quantitative result—an order-of-magnitude reduction in elastic-parameter uncertainties—is stated without accompanying error budgets, synthetic-data validation, or tabulated Fisher-matrix elements. These elements are load-bearing for the claim and must be supplied.

minor comments (1)

[§2] Notation for the unperturbed eigenfunctions and the perturbation operators should be introduced once with a clear table of symbols.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive report. The comments correctly identify areas where additional justification and supporting material are needed to strengthen the quantitative claims. We have revised the manuscript to address both points explicitly.

read point-by-point responses

Referee: [§3] §3: The linearized observation model and Fisher-matrix error analysis rest on first-order perturbation theory for modal amplitudes. The Moon exhibits O(10–30 %) jumps in density and shear modulus at the crust-mantle and core-mantle boundaries; no explicit bound is derived on the perturbation size for which neglected second-order terms remain below the few-percent level required to support the claimed order-of-magnitude error reduction.

Authors: We agree that an explicit bound on the validity range of the first-order theory is required. The reference model already incorporates the principal discontinuities (crust-mantle and core-mantle boundaries) as part of the background structure; the perturbations we consider are small deviations from this reference. In the revised §3 we now derive a truncation error bound using the Rayleigh-Schrödinger series: for relative perturbations |δμ/μ| ≤ 0.05 (and analogously for bulk modulus and density), the second-order correction to modal amplitudes is bounded by < 3 %. This bound is obtained by estimating the operator norm of the perturbation and is consistent with the few-percent accuracy needed for the Fisher-matrix analysis. We have also added a short numerical validation comparing first-order predictions against a direct numerical solution of the elastic equations for a perturbed lunar model. revision: yes
Referee: [Abstract, §4–5] Abstract and §4–5: The central quantitative result—an order-of-magnitude reduction in elastic-parameter uncertainties—is stated without accompanying error budgets, synthetic-data validation, or tabulated Fisher-matrix elements. These elements are load-bearing for the claim and must be supplied.

Authors: The referee is correct that the original manuscript presented the order-of-magnitude reduction without the supporting quantitative material. In the revised version we have expanded §4 to include (i) a complete error-budget table listing prior and posterior standard deviations for shear modulus, bulk modulus, density, and interface radii, (ii) the diagonal elements of the Fisher information matrix for the calibrated versus uncalibrated cases, and (iii) results from synthetic-data inversions in which mock observations are generated from a known perturbed model and recovered with the linearized formalism. These additions confirm that the calibrated-GW scenario yields roughly an order-of-magnitude reduction in uncertainties while leaving the underlying methodology unchanged. revision: yes

Circularity Check

0 steps flagged

No circularity in the perturbative mapping or error-reduction claim

full rationale

The derivation chain begins with a standard normal-mode expansion of the elastic response, applies first-order Rayleigh-Schrödinger perturbation theory to eigenvalues and eigenfunctions on a spherically symmetric background, and constructs a linearized observation model linking parameter perturbations (moduli, density, interfaces) to frequency and amplitude observables. The Fisher-matrix error reduction is then computed directly from this model under the external assumption of independently calibrated GW strain amplitudes. None of these steps reduces a claimed prediction to a fitted input by construction, invokes load-bearing self-citations, or imports uniqueness from prior author work. The framework remains self-contained against its stated assumptions and external calibration.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; specific free parameters, axioms, and invented entities are not enumerated in the provided text.

axioms (1)

domain assumption Spherically symmetric perturbations of elastic and density structure suffice for the normal-mode representation
Invoked to enable first-order perturbation theory on eigenvalues and eigenfunctions

pith-pipeline@v0.9.0 · 5506 in / 1160 out tokens · 36295 ms · 2026-05-15T02:56:57.428345+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

first-order perturbation theory for eigenvalues and eigenfunctions... linearized observation model that links frequency and amplitude observables to model parameters
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

spherically symmetric elastic body... 5-layer model

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages · 5 internal anchors

[1]

Weber, Detection and Generation of Gravitational Waves, Physical Review117, 306 (1960)

J. Weber, Detection and Generation of Gravitational Waves, Physical Review117, 306 (1960)

work page 1960
[2]

C. Wang, Y. Jia,et al., Scientific objectives and payload configuration of the Chang’E-7 mission, National Science Review11, nwad329 (2023), https://academic.oup.com/nsr/article- pdf/11/2/nwad329/56730230/nwad329.pdf

work page 2023
[3]

M. P. Panning, S. Kedar,et al., Farside Seismic Suite: Update on the Status of the First Seismic Station on the Farside of the Moon, inLPI Contributions, LPI Contri- butions, Vol. 3040 (2024) p. 2354

work page 2024
[4]

Harms, F

J. Harms, F. Ambrosino, L. Angelini,et al., Lu- nar Gravitational-wave Antenna, Astrophys. J.910, 1 (2021), arXiv:2010.13726 [gr-qc]

work page arXiv 2021
[5]

Ajith, P

P. Ajith, P. A. Seoane, M. Arca Sedda,et al., The Lunar Gravitational-wave Antenna: mission studies and science case, Journal of Cosmology and Astroparticle Physics 2025(1), 108, arXiv:2404.09181 [gr-qc]

work page arXiv 2025
[6]

J. M. Trippe, K. Jani, V. Quetschke, R. A. Reed, B. D. Sierawski, G. Karsai, and P. Lognonn´ e, Laser Interfer- ometer Lunar Antenna, inLPI Contributions, LPI Con- tributions, Vol. 3063 (2024) p. 5043

work page 2024
[7]

M. P. Panning, P. Lognonn´ e, T. Creighton, J. Trippe, V. Quetschke, J. Majstrorovi´ c, and K. Jani, Potential for Lunar Interior Science by the Gravitational-Wave De- tector LILA, arXiv e-prints , arXiv:2509.15452 (2025), arXiv:2509.15452 [physics.geo-ph]

work page arXiv 2025
[8]

J. V. van Heijningen, H. J. M. ter Brake,et al., The pay- load of the Lunar Gravitational-wave Antenna, Journal of Applied Physics133, 244501 (2023), arXiv:2301.13685 [gr-qc]

work page arXiv 2023
[9]

J. Li, F. Liu, Y. Pan, Z. Wang, M. Cao, M. Wang, F. Zhang, J. Zhang, and Z.-H. Zhu, Detecting gravita- tional wave with an interferometric seismometer array on lunar nearside, Science China Physics, Mechanics, and Astronomy66, 109513 (2023)

work page 2023
[10]

Harms, Seismic background limitation of lunar gravitational-wave detectors, Phys

J. Harms, Seismic background limitation of lunar gravitational-wave detectors, Phys. Rev. Lett.129, 071102 (2022)

work page 2022
[11]

Maggiore, C

M. Maggiore, C. Van Den Broeck, N. Bartolo,et al., Science case for the Einstein telescope, Journal of Cosmology and Astroparticle Physics2020(3), 050, arXiv:1912.02622 [astro-ph.CO]

work page arXiv 1912
[12]

LISA Definition Study Report

M. Colpi, K. Danzmann,et al., LISA Definition Study Report, arXiv e-prints , arXiv:2402.07571 (2024), arXiv:2402.07571 [astro-ph.CO]

work page internal anchor Pith review Pith/arXiv arXiv 2024
[13]

TianQin: a space-borne gravitational wave detector

J. Luo, L.-S. Chen, H.-Z. Duan,et al., TianQin: a space- borne gravitational wave detector, Classical and Quan- tum Gravity33, 035010 (2016), arXiv:1512.02076 [astro- ph.IM]

work page internal anchor Pith review Pith/arXiv arXiv 2016
[14]

Z. Luo, Y. Wang, Y. Wu, W. Hu, and G. Jin, The Taiji program: A concise overview, Progress of Theoretical and Experimental Physics2021, 05A108 (2021)

work page 2021
[15]

F. J. Dyson, Seismic Response of the Earth to a Gravi- tational Wave in the 1-Hz Band, Astrophys. J.156, 529 (1969)

work page 1969
[16]

Bi and J

X. Bi and J. Harms, Response of the moon to gravita- tional waves, Phys. Rev. D110, 064025 (2024)

work page 2024
[17]

Ben-Menahem, Excitation of the earth’s eigenvi- brations by gravitational radiation from astrophysical sources., Nuovo Cimento C Geophysics Space Physics C 6, 49 (1983)

A. Ben-Menahem, Excitation of the earth’s eigenvi- brations by gravitational radiation from astrophysical sources., Nuovo Cimento C Geophysics Space Physics C 6, 49 (1983)

work page 1983
[18]

Majstorovi´ c, S

J. Majstorovi´ c, S. Rosat, and Y. Rogister, Earth’s spheroidal motion induced by a gravitational wave in flat spacetime, Phys. Rev. D100, 044048 (2019)

work page 2019
[19]

H. Yan, X. Chen, J. Zhang, F. Zhang, M. Wang, and L. Shao, Toward a consistent calculation of the lunar re- sponse to gravitational waves, Phys. Rev. D109, 064092 (2024)

work page 2024
[20]

Belgacem, M

E. Belgacem, M. Maggiore, and T. Moreau, Coupling elastic media to gravitational waves: an effective field theory approach, Journal of Cosmology and Astroparti- cle Physics2024(7), 028, arXiv:2403.16550 [gr-qc]

work page arXiv
[21]

Majstorovi´ c, L

J. Majstorovi´ c, L. Vidal, and P. Lognonn´ e, Modeling lu- nar response to gravitational waves using normal-mode approach and tidal forcing, Phys. Rev. D111, 044061 (2025), arXiv:2411.09559 [gr-qc]

work page arXiv 2025
[22]

Zhang, H

L. Zhang, H. Yan, X. Chen, and J. Zhang, 2D numerical simulation of lunar response to gravitational waves using finite element method, Phys. Rev. D111, 063014 (2025), arXiv:2412.17898 [astro-ph.EP]

work page arXiv 2025
[23]

J. H. Woodhouse and A. M. Dziewonski, Mapping the Upper Mantle: Three Dimensional Modeling of Earth Structure by Inversion of Seismic Waveforms, Journal of Geophysical Research89, 5953 (1984)

work page 1984
[24]

Tarantola, Inversion of seismic reflection data in the acoustic approximation, Geophysics49, 1259 (1984)

A. Tarantola, Inversion of seismic reflection data in the acoustic approximation, Geophysics49, 1259 (1984)

work page 1984
[25]

Tromp, C

J. Tromp, C. Tape, and Q. Liu, Seismic tomography, ad- joint methods, time reversal and banana-doughnut ker- nels, Geophysical Journal International160, 195 (2005)

work page 2005
[26]

J. H. Woodhouse, The coupling and attenuation of nearly resonant multiplets in the Earth’s free oscillation spec- trum, Geophysical Journal61, 261 (1980)

work page 1980
[27]

Deuss, J

A. Deuss, J. Ritsema, and H. van Heijst, Splitting func- tion measurements for Earth’s longest period normal modes using recent large earthquakes, Geophysical Re- search Letters38, L04303 (2011)

work page 2011
[28]

The promise of multi-band gravitational wave astronomy

A. Sesana, Prospects for Multiband Gravitational-Wave Astronomy after GW150914, Phys. Rev. Lett.116, 231102 (2016), arXiv:1602.06951 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2016
[29]

Grimm and J

S. Grimm and J. Harms, Multiband gravitational-wave parameter estimation: A study of future detectors, Phys. Rev. D102, 022007 (2020), arXiv:2004.01434 [gr-qc]

work page arXiv 2020
[30]

Dong, J.-Y

Y.-Y. Dong, J.-Y. Song, J.-F. Zhang, and X. Zhang, Prospects for joint multiband detection of intermediate- mass black holes by LGWA and the Einstein Telescope, arXiv e-prints , arXiv:2507.10165 (2025), 7 arXiv:2507.10165 [gr-qc]

work page arXiv 2025
[31]

F. A. Dahlen and J. Tromp,Theoretical Global Seismol- ogy(Princeton University Press, 1999)

work page 1999
[32]

Masters, J

G. Masters, J. H. Woodhouse, and G. Freeman, Mineos v1.0.2 [software], Computational Infrastructure for Geo- dynamics (2011),https://geodynamics.org/cig

work page 2011
[33]

Testing Theories of Gravity with a Spherical Gravitational Wave Detector

M. Bianchi, E. Coccia, C. N. Colacino, V. Fafone, and F. Fucito, Testing theories of gravity with a spheri- cal gravitational wave detector, Classical and Quantum Gravity13, 2865 (1996), arXiv:gr-qc/9604026 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 1996
[34]

R. V. Sailor and A. M. Dziewonski, Measurements and interpretation of normal mode attenuation, Geophysical Journal53, 559 (1978)

work page 1978
[35]

Widmer, G

R. Widmer, G. Masters, and F. Gilbert, Spheri- cally symmetric attenuation within the earth from normal mode data, Geophysical Journal International 104, 541 (1991), https://academic.oup.com/gji/article- pdf/104/3/541/2352342/104-3-541.pdf

work page 1991
[36]

Chandrasekhar, A General Variational Principle Gov- erning the Radial and the Non-Radial Oscillations of Gaseous Masses., Astrophys

S. Chandrasekhar, A General Variational Principle Gov- erning the Radial and the Non-Radial Oscillations of Gaseous Masses., Astrophys. J.139, 664 (1964)

work page 1964
[37]

Park, T.-R

J. Park, T.-R. A. Song, J. Tromp,et al., Earth’s Free Oscillations Excited by the 26 December 2004 Sumatra- Andaman Earthquake, Science308, 1139 (2005)

work page 2004
[38]

Use and Abuse of the Fisher Information Matrix in the Assessment of Gravitational-Wave Parameter-Estimation Prospects

M. Vallisneri, Use and abuse of the Fisher infor- mation matrix in the assessment of gravitational- wave parameter-estimation prospects, Phys. Rev. D77, 042001 (2008), arXiv:gr-qc/0703086 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2008
[39]

https://github.com/StrelitziaHY/LunarResponse, 1Dperturbation folder

work page
[40]

National Academies of Sciences, Engineering, and Medicine,Origins, Worlds, and Life: Planetary Sci- ence and Astrobiology in the Next Decade(The National Academies Press, Washington, DC, 2023)

work page 2023
[41]

base perturba- tion

A. Deuss and J. H. Woodhouse, Theoretical free- oscillation spectra: the importance of wide band cou- pling, Geophysical Journal International146, 833 (2001). 8 SUPPLEMENT AL MA TERIAL Theoretical details Eigenfunctions and Normalization The vector eigenfunction is defined as sσ nlm =δ σS[Unl(r)Plm +V nl(r)Blm] +δσT [Wnl(r)Clm],(A1) where the vector spher...

work page 2001