{"id":"22b4c95d-4c39-4f62-960e-b47101876c57","arxiv_id":"2608.05117","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A validated Cartesian numerical framework computes gravitoelastic correlation tensors for Rayleigh and body waves in half-space and full-space media, benchmarked against new analytical closed forms.","lead":"This paper derives analytical formulas for how seismic vibrations create gravitational noise on test masses, and validates a numerical program that computes these correlations for different wave types and detector geometries. It gives gravitational-wave observatory designers a tool to plan seismometer arrays for noise cancellation.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The derivative-free representation Eq. (14) is used for half-space body-wave benchmarks without proving that the S_infinity boundary integral in Eq. (12) vanishes; for propagating P/SV fields that integral is oscillatory and cannot be assumed zero.","rationale":"The reader's weakest assumption identifies exactly the load-bearing issue: the transition from Eq. (12) to Eq. (14) requires vanishing of the S_infinity integral, and for half-space body waves this is not justified. My reading of the manuscript confirms that this is not a minor technicality. The kernel-contraction formulation is one of the two integral formulations claimed to be validated, and the half-space body-wave benchmarks are central to the paper's stated scope. The analytical closed forms in Sec. 3.2 are themselves derived independently of Eq. (14) from the plane-wave acceleration amplitudes, so the concern does not necessarily invalidate the analytical reference; however, if the kernel-contraction numerical code is based on Eq. (14), the comparison could be misleading. The manuscript gives no argument that the non-evanescent P/SV correlations make the boundary term vanish, and the convergence sweeps over L_z and L_xy do not address the infinite-boundary term because they use a finite domain. I agree with the conditional verdict: the central numerical-validation claim is plausible but not yet established for half-space body waves in the kernel-contraction formulation. The lack of public code and the internal-only benchmarks are secondary concerns; the S_infinity condition is the single most load-bearing unproven step. No change to the reader's verdict is needed.","tokens_in":17170,"tokens_out":8071,"duration_ms":79875,"concrete_test":"Evaluate the boundary integral in Eq. (12) explicitly for the half-space body-wave correlation tensor built from Eqs. (58)-(61) with the reflection coefficients of Eq. (62). Numerically compute I_S_infinity(R, H) over the truncated domain boundary (bottom at z' = -H, lateral at rho' = R) for (R/lambda, H/lambda) = (5,5), (10,10), (20,20), (40,40) at a fixed observation point; if I_S_infinity extrapolates to a value comparable to the principal-value or local term, Eq. (14) is invalid and the kernel-contraction half-space body-wave benchmarks must be recomputed using the full Eq. (12) or the bulk-surface Eq. (7). If I_S_infinity tends to zero, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's derivative-free representation, Eq. (14), is obtained from Eq. (12) only by dropping the S_infinity boundary integral. The text states this is valid 'if the contribution at infinity vanishes,' but it never verifies that condition for half-space body waves. In those benchmarks the displacement correlation tensor is not evanescent: the P/SV depth profiles in Eqs. (58)-(61) contain propagating factors e^{±ik_z z'}, and the horizontal angular average leaves Bessel functions J_n(kappa*Delta_rho) that decay only as (kappa*Delta_rho)^{-1/2}. Consequently, the surface integral over the lower boundary z' -> -infinity is oscillatory rather than exponentially decaying, and the lateral boundary contribution is at most conditionally convergent. No proof, symmetry argument, or numerical check is given that this S_infinity term vanishes. Since the kernel-contraction numerical formulation used for half-space body waves appears to rely on Eq. (14), a nonzero S_infinity term would mean the numerical reference itself is missing a boundary contribution. The small validation errors reported in Sec. 5.2 would then not establish the central claim for the kernel-contraction half-space body-wave cases; they could instead reflect the same omission in both the numerical and closed-form expressions. The paper explicitly flags the condition, but the condition remains unexamined for the non-evanescent body-wave case.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Cartesian midpoint-grid numerical framework for computing gravitoelastic correlation tensors (cross-correlations between Newtonian acceleration at a test mass and seismic displacement at sensor positions) for Rayleigh waves and body waves in half-space and full-space homogeneous isotropic media. It derives closed-form analytical benchmarks, including a universal angular-averaged tensor, mode-specific displacement and acceleration profiles, incoherent P-SV mixtures, and leading-order spherical-cavity corrections. The numerical framework is then validated against these closed forms across all nine tensor components, with convergence studies over resolution, domain size, test-mass position, frequency, and cavity radius, and with spatial map comparisons. The central claim is that the numerical framework reproduces the analytical tensors with mean total errors between about 6e-4 and 1.3e-2, thereby providing a unified tool for Newtonian-noise modeling and sensor-array design.","tokens_in":17408,"tokens_out":15502,"duration_ms":134629,"significance":"If the validation is sound, this paper provides a useful, systematically benchmarked numerical tool for Newtonian-noise estimation in future gravitational-wave detectors such as the Einstein Telescope. The strengths are the thoroughness of the convergence studies (all nine tensor components, multiple geometries, frequency and cavity-radius sweeps), the well-defined normalized error metrics, and the explicit use of equipartition-based diffuse-field fractions that are not fitted to the numerical results. The paper also clearly separates exact full-space models from principal-value and approximate-cavity references, and it identifies regimes where finite-cavity and truncation errors matter. The main limitations are that no code or data are made publicly available, and two technical points in the derivation and in the scope of the cavity formula need clarification before the validation claims can be fully accepted.","major_comments":[{"comment":"Equation (14) is obtained from Eq. (12) by dropping the boundary integral at infinity, justified only 'if the contribution at infinity vanishes.' For the half-space body-wave benchmarks used in Section 5, the correlation tensor C_αj(r',r) is not evanescent: the depth profiles in Eqs. (58)-(61) contain propagating factors e^{±ik_z z'}, and after the angular average the Bessel functions J_n(κΔρ) decay only as (κΔρ)^{-1/2}. The manuscript never verifies that the S∞ term in Eq. (12) vanishes, yet Section 3.3.2 refers to Eq. (40) as the 'exact interior relation' without qualification. Please add a demonstration that the boundary integral vanishes (for example, using |J_n(κΔρ)| ≤ C(κΔρ)^{-1/2} to bound the hemisphere boundary integral by a constant times R^{-1/2} → 0), or state explicitly the additional assumption needed and restrict the exactness claim to the cases where the condition is established.","section":"Sec. 2.2.2 and Sec. 3.3.2"},{"comment":"Equation (41) is presented as the leading-order cavity-corrected closed form for both the kernel-contraction and bulk-surface models, but the text states that it 'is valid for Rayleigh waves and the full-space medium.' The numerical validation, however, applies Eq. (41) as the reference for half-space body-wave cavity configurations (e.g., the 'half-space body-wave cavity model' in Fig. 4 and the small-radius plateaus in Fig. 5). If the stated restriction is correct, those validation results are not meaningful; if the restriction is not intended, the sentence should be revised. Please clarify the domain of validity of Eq. (41) and ensure that the validation is consistent with that domain.","section":"Sec. 3.4 and Sec. 5.2.3"}],"minor_comments":[{"comment":"The statement that code and data are available 'upon request' rather than in a permanent repository limits reproducibility of the reported convergence numbers; please consider releasing the code or providing a detailed description of the data products.","section":"Data and code availability"},{"comment":"Please specify which integral formulation (bulk-surface or kernel-contraction) the reported mean-total-error ranges refer to, since both formulations are claimed to be verified and the reader cannot tell whether the ranges cover one or both.","section":"Sec. 5.2.1"},{"comment":"The grid discretization does not describe how the spherical cavity boundary is handled on the Cartesian grid (for example, whether cells whose centers lie inside the sphere are excluded with a staircase approximation); please specify this, since the discussion of non-monotonic convergence cites the 'Cartesian representation of the spherical cavity boundary' as a source of error.","section":"Appendix B"},{"comment":"The acceleration transfer functions B_pot(z0) and B_vert(z0) are listed without derivation; please provide a derivation or a precise reference to the standard half-space potential solutions from which they follow.","section":"Appendix D.2.2"},{"comment":"The projection formula contains c_xy + c_yx; if the gravitoelastic tensor is symmetric, this should be noted and the expansion written as 2c_xy, otherwise the asymmetry should be explained.","section":"Eq. (24)"}],"recommendation":"major_revision","confidential_remarks":"The S∞ concern raised in the external stress-test is real as a rigor gap but not as an error: the algebraic decay of the Bessel functions is sufficient to make the boundary integral vanish, so the condition in Eq. (14) actually holds for the half-space body-wave benchmarks. The more consequential issue is the apparent contradiction between the claimed validity of Eq. (41) and its use for half-space body-wave cavity cases; that needs to be resolved in the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Send this one out. The paper has a real contribution: Eq. (35) is a clean universal closed-form gravitoelastic tensor for horizontally isotropic fields, and Appendix D turns the mode-specific Rayleigh/body-wave cases into a small set of master-equation mappings. The Cartesian midpoint-grid framework is validated against those closed forms for Rayleigh waves and body waves, half-space and full-space, above- and below-ground, with and without a spherical cavity. That is directly useful for Newtonian-noise cancellation and sensor-array design for Einstein Telescope.\n\nThe validation work is solid. The error metrics are well defined, all nine tensor components are checked, and the convergence sweeps over ppw, domain size, frequency, and cavity radius are systematic. The maps at 0.5–2.4% mean discrepancy support the central claim. I also credit the paper for clearly labeling the cavity corrections as leading-order and for stating that scattering is not included.\n\nThe stress-test concern about the S∞ term in Eq. (14) does not survive a close read. It is true that the body-wave correlation tensor is not evanescent and the paper only says the surface term vanishes. But it does vanish: on the lower boundary the oscillatory Bessel factors and the 1/R^2 kernel give a contribution O(|z|−1/2), and on the lateral boundary the Fourier transform in z′ decays exponentially. The result is correct, but the derivation should say this. Right now the condition is asserted, and a referee should ask for a short proof or at least a numerical check.\n\nThe real weakness is reproducibility. The paper presents a numerical framework as its main deliverable, yet no code or data are released. “Available upon request” is not good enough for a validation paper. I would require a permanent repository before publication. The analytical benchmarks are internal—they come from the same physical model—so the validation is mostly about the quadrature and algebra, not about independent physics. The paper is honest about this, but it is still a limitation.\n\nFull-space formulas are mostly from Harms 2019; the new material is the half-space body-wave closed forms with incident P/SV mixtures and the cavity corrections. I would cite this for those formulas and for the master-equation mapping.\n\nRecommendation: accept for peer review, with minor-to-moderate requests: release code/data, and justify the S∞ boundary condition explicitly.","headline":"A solid, useful synthesis of gravitoelastic correlation tensors for Newtonian-noise modeling, held back mainly by missing code/data and one hand-waved boundary term.","tokens_in":17945,"tokens_out":7531,"would_cite":true,"duration_ms":66427,"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":"A Cartesian numerical framework reproduces the analytical gravitoelastic correlation tensors for Rayleigh waves and body waves in half-space and full-space media, with mean total errors from 6.17e-4 to 1.31e-2.","keywords":["Newtonian noise","gravity gradient noise","gravitoelastic correlation tensor","Rayleigh waves","body waves","noise cancellation","seismometer arrays","gravitational-wave detectors"],"falsifier":"Evaluate the surface term at infinity in the derivative-free representation for a propagating P/SV half-space field at the paper's nominal parameters; a nonzero result means the analytical body-wave closed forms are incomplete, and the reported numerical agreement would not certify the framework against the true gravitoelastic response.","tokens_in":16952,"feed_emoji":"🌊","tokens_out":8579,"duration_ms":77660,"temperature":0.7,"pith_summary":"The paper is trying to establish that the cross-correlations between seismic ground motion and the gravitational acceleration on a suspended test mass—the gravitoelastic correlation tensor—can be computed reliably with a Cartesian numerical integration scheme, given the statistics of the seismic field. To prove this, the authors derive exact and asymptotic analytical forms for Rayleigh waves and body waves in half-space and full-space media, including above-ground, underground, and spherical-cavity configurations, and show that the numerical integrals match them. If true, this gives future gravitational-wave observatories a practical way to design seismometer arrays and subtract seismic Newtonian noise without simulating the full elastic wavefield in the time domain.","feed_headline":"Newtonian-noise correlations computed to under 2 percent error","feed_subtitle":"Cartesian integration matches analytical tensors for Rayleigh and body waves, above and below ground, with and without cavities.","key_machinery":"The central object is the gravitoelastic correlation tensor, the cross-spectral tensor linking the Newtonian acceleration perturbation at the test mass to the seismic displacement at a seismometer. The argument runs on two integral representations: a bulk-surface integral over the two-point displacement correlation tensor, and a derivative-free kernel-contraction form obtained by moving the derivative onto the Newtonian kernel, which produces a principal-value volume integral plus a local term $(4\\pi/3)G\\rho_0 C_{ij}(\\mathbf{r}_0,\\mathbf{r})$. The numerical engine is a Cartesian midpoint grid with even cell counts and depth alignment, while the analytical benchmarks come from angular Bessel identities that convert azimuthal averages into $J_0$, $J_1$, and $J_2$ combinations, and a diffuse-field equipartition prescription fixes the P/SV power fractions in mixed-wave models.","core_discovery":"The central claim is that one Cartesian midpoint-grid integration framework reproduces the analytical gravitoelastic correlation tensors $c_{ij}(\\mathbf{r},\\omega)=\\langle\\delta a_i(\\mathbf{r}_0,\\omega)\\,\\xi_j^*(\\mathbf{r},\\omega)\\rangle$ for Rayleigh waves and for body waves in homogeneous half-space and full-space media, for test masses above and below ground and with or without a spherical source-exclusion cavity. The validation compares all nine tensor components over spatial maps and reports mean total Frobenius errors between $6.17\\times10^{-4}$ and $1.31\\times10^{-2}$, with peak-normalized map discrepancies below roughly nine percent in the worst case. On this evidence, the framework is a validated tool for Newtonian-noise cancellation studies and sensor-array design in future low-frequency gravitational-wave detectors.","pith_inferences":["Beyond the paper, the per-component error curves suggest a practical grid-design recipe: coarser horizontal grids may suffice once the xx and yy residuals dominate, while zz-sensitive configurations need finer vertical sampling.","The same Cartesian integrator could be fed with displacement correlation tensors computed for heterogeneous or topographic models, with the homogeneous benchmarks in this paper serving as a regression suite for such extensions.","If the surface-at-infinity term in the half-space derivation turns out not to vanish, the analytical body-wave benchmarks would need revision, but the numerical framework itself would remain testable against a corrected reference."],"forward_implications":["With a validated correlation tensor, Newtonian-noise cancellation studies can compute sensor-array correlations directly on a Cartesian grid, avoiding full time-domain wavefield simulations for benchmark design.","The same framework covers above-ground, underground, and spherical-cavity test-mass geometries, so one code can be applied across the configurations of a proposed detector site.","Because all nine tensor components are validated, off-diagonal correlations are available for multi-axis and arbitrary-arm-orientation array optimization, not just scalar noise power.","The frequency and cavity-radius sweeps delimit the regimes where the small-cavity principal-value approximation is reliable, telling users how small the cavity must be for the simple analytical model to apply."],"supporting_citations":[{"why":"Defines the gravitoelastic correlation tensor and supplies the bulk-surface integral, the full-space P/S kernels, and the Rayleigh secular equation used as analytical benchmarks.","marker":"[6]"},{"why":"Provides the free-surface reflection and conversion coefficients used to build the half-space body-wave displacement profiles.","marker":"[15]"},{"why":"Establishes the diffuse-field equipartition ratio used to fix the P/S power fractions in mixed-wave benchmarks.","marker":"[16]"},{"why":"Provides the diffuse-field spectral weighting used with [16] to set the half-space and full-space P fractions.","marker":"[17]"},{"why":"Formulates Newtonian-noise cancellation with gravitoelastic correlations, motivating the array-design target of the framework.","marker":"[8]"}],"fun_headline_variants":["Newtonian noise correlations: precise to <2% error","Gravitoelastic tensors validated for wave fields","Half-space and full-space Newtonian noise models unified","Analytical benchmarks for seismic Newtonian noise","Cartesian framework matches analytical noise tensors"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analytical half-space body-wave benchmarks are derived only after dropping a surface integral at infinity, and the paper does not prove that this term vanishes for propagating body waves.","fun_headline_variants_meta":{"raw":{"variants":["Newtonian noise correlations: precise to <2% error","Gravitoelastic tensors validated for wave fields","Half-space and full-space Newtonian noise models unified","Analytical benchmarks for seismic Newtonian noise","Cartesian framework matches analytical noise tensors"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000739,"raw_usage":{"total_tokens":3265,"prompt_tokens":877,"completion_tokens":2388,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":493,"completion_tokens_details":{"reasoning_tokens":2315}},"tokens_in":493,"tokens_out":2388,"duration_ms":16982,"temperature":1.0,"reasoning_tokens":2315,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:36:42.346518+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the surface term at infinity in the derivative-free representation for a propagating P/SV half-space field at the paper's nominal parameters; a nonzero result means the analytical body-wave closed forms are incomplete, and the reported numerical agreement would not certify the framework against the true gravitoelastic response.","supporting_citations":[{"cited_title":"Relativ.226 updated and expanded review","cited_arxiv_id":null,"evidence_quote":"Defines the gravitoelastic correlation tensor and supplies the bulk-surface integral, the full-space P/S kernels, and the Rayleigh secular equation used as analytical benchmarks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the free-surface reflection and conversion coefficients used to build the half-space body-wave displacement profiles."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the diffuse-field equipartition ratio used to fix the P/S power fractions in mixed-wave benchmarks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the diffuse-field spectral weighting used with [16] to set the half-space and full-space P fractions."},{"cited_title":"Quantum Grav.36145006","cited_arxiv_id":null,"evidence_quote":"Formulates Newtonian-noise cancellation with gravitoelastic correlations, motivating the array-design target of the framework."}],"review_version":1}