{"id":"de35a744-a011-4051-90cc-b20286b2f227","arxiv_id":"2412.17898","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A 2D spectral element simulation of the Moon's response to gravitational waves qualitatively reproduces semi-analytical normal-mode results, demonstrating feasibility of numerical lunar GW response modeling.","lead":"This paper simulates, for the first time with a 2D finite-element method, how the Moon vibrates when a gravitational wave passes through it, and compares the results with older analytical calculations. The goal is to enable realistic simulations that could help design lunar seismometer arrays for gravitational-wave detection.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fig. 8 agreement is obtained only after a post-hoc 0.06 dex frequency shift that is neither derived from the cylindrical/spherical eigenmode analysis nor mode-by-mode; without an independent 2D benchmark the central feasibility claim remains unproven.","rationale":"The reader's weakest-assumption identification is accurate and is the same as mine: the 2D-to-3D comparison with a free frequency shift is the weakest link. I add specificity: Appendix A's own numbers imply a mode-dependent offset of up to ~0.16 dex for the fundamental, not the uniform 0.06 dex used in Fig. 8, so the shift is essentially a fitting parameter. This makes the validation weaker than it appears, but it does not make the paper's conclusion false—merely not yet tightly supported. The paper is a useful first step, clearly describes the simulation, and compares against an independent normal-mode solver, so rejection would be too harsh. The appropriate verdict is conditional acceptance pending an independent 2D analytic benchmark and quantification of the 2D/3D mapping. This matches the reader's verdict, so no adjustment is needed.","tokens_in":12074,"tokens_out":6325,"duration_ms":62416,"concrete_test":"Run the same FEM setup on a homogeneous cylindrical Moon with the Appendix A parameters, drive it with the Dyson force (m=2), and compare the extracted Tr(f), Th(f) with the exact 2D analytic solution (eigenfrequencies from Eq. A7 and eigenfunction amplitudes from Eq. A1) without any frequency shift. If the first several peaks agree in frequency to <1% and in amplitude to <10%, the FEM implementation is validated in 2D and the Fig. 8 shift can be interpreted as a genuine 2D/3D difference. If they do not agree, the 0.06 dex translation is masking a numerical error and the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the FEM simulation is feasible for computing lunar GW response rests primarily on Fig. 8, where the FEM-derived Tr(f) and Th(f) are compared with semi-analytical MINEOS results after 'horizontally translated by about 0.06 dex.' This translation is the load-bearing element: it is the only step that brings the first few resonant peaks into coincidence. The paper does not derive this shift from the 2D/3D comparison in Appendix A. There, a homogeneous spherical Moon (l=2) has fundamental eigenfrequency 0.971 mHz while a homogeneous cylindrical Moon (m=2) has 0.666 mHz, a ratio of 1.46 (~0.16 dex), with different ratios for higher modes (1.868 vs 1.970; 3.161 vs 3.087; etc.). Thus the expected cylindrical/spherical offset is mode-dependent and generally larger than 0.06 dex. Using a single fitted constant shift over the whole 1–20 mHz band can absorb a systematic frequency error in the FEM implementation, so the qualitative agreement in Fig. 8 does not by itself validate the Dyson-force implementation. The angular fits in Fig. 5 only check that the response has the expected quadrupolar pattern, not that frequencies or absolute amplitudes are correct. The feasibility claim therefore hinges on an asserted, not derived, equivalence between the 2D cylindrical model and the 3D spherical solution.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a two-dimensional spectral-element simulation of the lunar seismic response to gravitational waves, extending SPECFEM2D with a Dyson-type force density. The authors build a global 2D cylindrical Moon model with 140,576 fourth-order elements, drive it with a broad-band Gaussian source-time function, and extract radial and horizontal response functions Tr(f) and Th(f) in the 1–20 mHz band. They compare these with semi-analytical normal-mode response functions computed from the same lunar model, finding qualitative agreement only after shifting the FEM spectra horizontally by about 0.06 dex. They also analyze angular-resolution convergence and give a semi-analytical comparison of spherical and cylindrical eigenfrequencies in Appendix A. The paper concludes that FEM is feasible for lunar gravitational-wave response calculations and lays groundwork for 3D simulations.","tokens_in":12247,"tokens_out":5091,"duration_ms":52942,"significance":"If the validation is accepted, this is a useful proof-of-concept: it is the first global FEM/SEM calculation of the lunar GW response, it uses a long 100,000 s simulation, the model and code are publicly available on GitHub, and the angular-resolution test in Fig. 7 is a genuine convergence check. The paper is also honest about the limitations of the 2D plane-strain model. However, the central validation is weaker than the conclusion claims: the only quantitative comparison with the semi-analytical solution is made after applying a constant, unexplained horizontal frequency shift, while Appendix A itself shows that the cylindrical/spherical frequency offset is mode-dependent. Because the feasibility claim rests on that comparison, the current manuscript does not yet firmly establish that the FEM implementation reproduces the correct frequency content and amplitudes of the lunar response.","major_comments":[{"comment":"The central validation claim depends on the comparison in Fig. 8, but the FEM curves are shown only after being 'horizontally translated by about 0.06 dex'. This shift is not derived from Appendix A. For the homogeneous models in Table I, the spherical/cylindrical frequency ratios are mode-dependent: log10(0.971/0.666) = 0.164 for the first mode, log10(1.868/1.970) = -0.023 for the second, and 0.059 for the fourth. A single constant shift therefore cannot compensate for the geometrical difference in a principled way, and the agreement in Fig. 8 could absorb a systematic error in the FEM frequency scale. To support the feasibility claim, the authors should either benchmark the 2D FEM against an independent 2D cylindrical analytical solution, or derive the frequency mapping between their heterogeneous 2D model and the 3D spherical model rather than fitting it post hoc.","section":"§V, Fig. 8; Appendix A, Table I"},{"comment":"The angular fits in Fig. 5 test only that the response has the expected quadrupolar angular pattern, not that the eigenfrequencies or absolute amplitudes are correct. The comparison of Tr(f) and Th(f) in Fig. 8 is presented on logarithmic axes after a horizontal translation, with no quantitative error metric. The authors acknowledge visible deviations at high-frequency peaks and spectral troughs, but the conclusion 'agree generally well' is not backed by a criterion such as per-mode relative frequency error or a normalization-free amplitude comparison. A quantitative measure of agreement is needed because the paper explicitly states that the FEM results are the basis for concluding that the simulation is feasible for GW response calculations.","section":"§IV, Figs. 5–6; §V, Fig. 8"},{"comment":"The Dyson force is implemented as 4320 discrete point forces at 1° azimuthal spacing, but the paper does not describe how the forces are registered to the finite-element mesh at the sharp interfaces where the shear modulus changes. The authors themselves note that 'the way of constructing the grids and adding the force at these places could affect the numerical results' and defer this to future work. Since the Dyson force is concentrated at interfaces, the radial placement of the force relative to element boundaries could affect the excitation amplitudes of the modes. A radial-resolution or interface-registration test would strengthen the validity of the extracted Tr and Th amplitudes.","section":"§III, Fig. 1c; §V second paragraph"}],"minor_comments":[{"comment":"There are several typographical issues: 'aross' should be 'across', 'course grid' should be 'coarse grid', 'per unite strain' should be 'per unit strain', and 'Beside' at the start of the paragraph before Fig. 8 should be 'Besides'.","section":"Throughout"},{"comment":"The polarization tensor epsilon_ij is written as [[1,1,0],[1,-1,0],[0,0,0]], which contains an off-diagonal 1; for a linear plus polarization one would expect diag(1,-1,0). Please correct the matrix or clarify the intended convention.","section":"Eq. (1)"},{"comment":"The text says 'we calculate both m = 2 and m = 2.5'. For a cylindrical coordinate solution to be single-valued in the azimuthal angle, m should be an integer; m = 2.5 arises only from the Bessel-order relation with spherical Bessel functions and should be described as a formal correspondence rather than as a physical angular quantum number.","section":"Appendix A"},{"comment":"Several reference entries have inconsistent formatting, such as stray spaces in 'arXiv e-prints ,' and incomplete journal information for Ref. [30] and Ref. [41]. The reference list should be cleaned up.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of the journal and addresses a timely topic for lunar gravitational-wave detectors. The core problem is that the validation relies on a fitted 0.06-dex frequency shift, while the paper's own Appendix A shows that the 2D-to-3D frequency relation is mode-dependent. I do not think this is grounds for rejection, because a 2D analytical benchmark or a derived mapping is within reach, but the central claim needs to be re-supported before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is the first spectral-element simulation of the Moon's seismic response to gravitational waves, and it mostly works, but the central feasibility claim is carried by a post-hoc 0.06 dex frequency shift that the paper never derives. The authors are transparent about the shift, so this is an honest weakness, not a hidden one.\n\nWhat is genuinely new: they implement the Dyson force density in SPECFEM2D, run a 100,000 s simulation, test angular resolution convergence (three resolutions agree within 0.1%), and compare with MINEOS normal-mode solutions. The appendix comparing spherical and cylindrical eigenfrequencies is a useful contribution, and the model files appear to be on GitHub. The angular pattern of the response at the first three resonances matches the expected quadrupolar form, which is a meaningful sanity check.\n\nThe soft spot is Fig. 8. The FEM curves are shifted horizontally by 0.06 dex to align the first few peaks with MINEOS. That shift is not derived from the 2D/3D eigenmode comparison in Appendix A—in fact, that appendix shows the offset is mode-dependent and generally larger (0.16 dex for the fundamental). So the constant shift is absorbing a real model mismatch, and the agreement after shifting does not independently validate the force implementation. The horizontal response also deviates visibly at spectral troughs, and there is no uncertainty quantification. The abstract's talk of 'realistic' simulations overstates what is demonstrated; no topography or heterogeneity is actually included.\n\nThe circularity worry is minor. MINEOS is an independent solver, and the FEM results are not forced to match it, so the comparison is meaningful even though it is not conclusive.\n\nThis paper is for researchers building numerical tools for lunar GW detection and lunar seismology. It deserves peer review, but it needs a major revision: either derive the 2D-to-3D mapping or validate against a genuine 2D analytic solution, add uncertainty estimates, and temper the abstract. I would accept it conditionally.","headline":"First SEM simulation of lunar GW response, but validation rests on an underived frequency shift.","tokens_in":12875,"tokens_out":3490,"would_cite":true,"duration_ms":33428,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes that a two-dimensional finite-element simulation of the Moon reproduces semi-analytical gravitational-wave response functions, making numerical lunar GW seismology with realistic structure a feasible next step.","keywords":["gravitational waves","lunar seismology","spectral element method","finite element method","lunar response functions","normal modes","gravitational wave detector","force density"],"falsifier":"Run a full three-dimensional spectral-element simulation of the same homogeneous or layered lunar model, feed it the same gravitational-wave force density and source time function, and compare its first radial resonant peak with the spherical normal-mode prediction (about 0.971 mHz for a homogeneous Moon) and with the shifted 2D peak; if the shifted 2D result and the 3D result do not coincide, the 2D-to-3D mapping and the paper's feasibility claim are refuted.","tokens_in":11771,"feed_emoji":"🌕","tokens_out":10415,"duration_ms":88814,"temperature":0.7,"pith_summary":"Previous calculations of the Moon's response to gravitational waves used analytical normal-mode solutions that assume an ideal spherical Moon, so they cannot capture rough topography or lateral changes in the lunar crust and interior. This paper reports the first two-dimensional high-order finite-element (spectral element) simulation of the whole Moon driven by a gravitational wave, with the wave force applied where elastic rigidity changes sharply across layer boundaries. Comparing the simulated radial and horizontal response functions with semi-analytical normal-mode results over 1–20 mHz, the authors find general agreement in the positions and order of magnitude of the first resonant peaks once the numerical curves are shifted horizontally by about 0.06 dex to account for the difference between a 2D cylindrical model and a 3D spherical one. They conclude that the finite-element method is feasible for lunar gravitational-wave response calculations and that it lays the groundwork for future three-dimensional simulations with realistic lunar structure. The comparison also shows that the 2D model tracks radial response better than horizontal response, especially in the spectral troughs.","feed_headline":"2D lunar model reproduces gravitational-wave response","feed_subtitle":"Finite-element simulation matches semi-analytical curves, opening the path to 3D lunar GW seismology.","key_machinery":"The simulation is carried out with a two-dimensional spectral element method, a high-order finite-element technique, on a global lunar model built from 140,576 fourth-order elements with a mesh that grows coarser with seismic wave speed; the grid resolves frequencies up to about 0.2 Hz. The gravitational wave enters through the force density $\\vec{f}=\\nabla\\mu\\cdot h$, where $\\mu$ is the shear modulus and $h$ is the GW strain tensor, so the force acts almost entirely at interfaces where rigidity changes; in practice this is implemented as 4,320 point forces spaced 1° in azimuth. A Gaussian-wavelet source time function with a flat spectrum up to 20 mHz excites many frequencies in one 100,000-second run, and the response functions are obtained by Fourier transforming the surface displacement at each azimuth and fitting the angular dependence with the equatorial-plane form of the $\\ell=2$ response.","core_discovery":"On its own terms, the paper establishes that in an idealized symmetric Moon model, the response functions obtained from the two-dimensional spectral-element simulation agree generally well with previous semi-analytical normal-mode solutions. The radial response function $T_r(f)$ and the horizontal response function $T_h(f)$ both reproduce the order of magnitude and the locations of the first few resonant peaks in the 1–20 mHz band after a horizontal translation of about 0.06 dex, which compensates for the difference between the cylindrical Moon represented by the 2D grid and the spherical Moon assumed by the analytical solution. The simulated angular dependence also matches the pattern expected when the $\\ell=2$ spherical response is restricted to the equatorial plane. These points of agreement are what the authors read as evidence that their finite-element implementation is a feasible way to compute lunar gravitational-wave response and can serve as the foundation for 3D simulations.","pith_inferences":["If the 0.06-dex frequency shift reflects a systematic 2D-to-3D mapping rather than free calibration, then cheap 2D runs could be used to survey resonant bands across candidate lunar interior models before expensive 3D runs are undertaken.","A decisive test the paper leaves open is a full 3D spectral-element simulation of the same symmetric lunar model; direct comparison of its peaks with the shifted 2D curve would either validate the shift as physical or expose it as an artifact of the cylindrical approximation.","The cylindrical eigenfrequency calculation in the appendix suggests the offset depends on mode order and on the chosen azimuthal index, so a constant 0.06-dex shift may be a poor approximation for strongly heterogeneous or aspherical lunar models.","Because the force density lives at sharp shear-modulus interfaces, the simulation's sensitivity to grid refinement at those boundaries is testable by rerunning with finer crustal meshes and checking whether the peak positions and the horizontal-response mismatch move."],"forward_implications":["Sensitivity estimates for proposed lunar gravitational-wave detectors can be cross-checked against finite-element response curves rather than relying only on normal-mode models.","The same solver can be pointed at a Moon model containing craters, fractured regolith, and laterally heterogeneous layers, since FEM does not require spherical symmetry.","Quantitative predictions from the 2D model should be treated as more trustworthy for radial motion than for horizontal motion, especially at frequencies between the resonant peaks.","With more memory and longer runs, the current grid supports pushing the simulation band upward from 20 mHz toward 200 mHz, where the lunar GW detection band remains largely unexplored."],"supporting_citations":[{"why":"Provides the normal-mode Green's-function formalism for a sphere's response to a GW, including the force-density formulation the paper builds on.","marker":"[16]"},{"why":"Gives the previous semi-analytical lunar response functions and the lunar model that the FEM simulation is designed to match.","marker":"[17]"},{"why":"Supplies the tidal-force variant of the normal-mode response calculation that informs the semi-analytical reference curves.","marker":"[19]"},{"why":"Introduces the spectral-element method in seismology that the 2D simulation is based on.","marker":"[33]"},{"why":"Computes the eigenfrequencies and quality factors that define the semi-analytical reference response functions.","marker":"[35]"},{"why":"Provides the 2D spectral-element seismic wave code that the authors extend with the GW force terms.","marker":"[37]"},{"why":"Documents the numerical formulation of the spectral-element wave propagation used in the simulation.","marker":"[38]"},{"why":"Defines the gradient-of-shear-modulus force density $\\nabla\\mu\\cdot h$ used as the GW driving term.","marker":"[40]"},{"why":"Derives the eigenfrequency determinants for spherical and cylindrical elastic bodies used to explain the 2D-to-3D frequency offset.","marker":"[47]"}],"fun_headline_variants":["2D Moon simulation reproduces GW response","Finite-element Moon echoes gravitational-wave theory","2D lunar model validates gravitational-wave response","2D simulation paves way to 3D lunar GW seismology","Moon's GW response simulated in 2D finite-element model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The validation rests on treating a two-dimensional cylindrical Moon as a proxy for the three-dimensional spherical Moon and on applying a constant 0.06-dex frequency shift after the simulation; if that geometric mapping does not hold, the agreement seen in the comparison does not establish feasibility.","fun_headline_variants_meta":{"raw":{"variants":["2D Moon simulation reproduces GW response","Finite-element Moon echoes gravitational-wave theory","2D lunar model validates gravitational-wave response","2D simulation paves way to 3D lunar GW seismology","Moon's GW response simulated in 2D finite-element model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001718,"raw_usage":{"total_tokens":6763,"prompt_tokens":877,"completion_tokens":5886,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":493,"completion_tokens_details":{"reasoning_tokens":5810}},"tokens_in":493,"tokens_out":5886,"duration_ms":37777,"temperature":1.0,"reasoning_tokens":5810,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:09:01.022638+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a full three-dimensional spectral-element simulation of the same homogeneous or layered lunar model, feed it the same gravitational-wave force density and source time function, and compare its first radial resonant peak with the spherical normal-mode prediction (about 0.971 mHz for a homogeneous Moon) and with the shifted 2D peak; if the shifted 2D result and the 3D result do not coincide, the 2D-to-3D mapping and the paper's feasibility claim are refuted.","supporting_citations":[{"cited_title":"Majstorovi´ c, S","cited_arxiv_id":null,"evidence_quote":"Provides the normal-mode Green's-function formalism for a sphere's response to a GW, including the force-density formulation the paper builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the previous semi-analytical lunar response functions and the lunar model that the FEM simulation is designed to match."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Computes the eigenfrequencies and quality factors that define the semi-analytical reference response functions."},{"cited_title":"Komatitsch and J.-P","cited_arxiv_id":null,"evidence_quote":"Provides the 2D spectral-element seismic wave code that the authors extend with the GW force terms."},{"cited_title":"Komatitsch and J","cited_arxiv_id":null,"evidence_quote":"Documents the numerical formulation of the spectral-element wave propagation used in the simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the gradient-of-shear-modulus force density $\\nabla\\mu\\cdot h$ used as the GW driving term."},{"cited_title":"What can we learn about GW Physics with an elastic spherical antenna?","cited_arxiv_id":"gr-qc/0006102","evidence_quote":"Derives the eigenfrequency determinants for spherical and cylindrical elastic bodies used to explain the 2D-to-3D frequency offset."}],"review_version":1}