{"id":"6a9f2ea9-bcfd-43e4-abe5-5d7e59c0ed46","arxiv_id":"2412.13701","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Adding the octupole moment of galaxy cross-correlations to the dipole improves forecasted constraints on the local position invariance violation parameter alpha by about 11% at smin = 15 Mpc/h.","lead":"This paper forecasts how well future galaxy surveys can test a fundamental physics principle, local position invariance, by using an extra pattern in galaxy clustering. It finds that combining the octupole with the dipole improves the expected constraints by about 11% on conservative scales.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 11% improvement is computed from an octupole model (Eq. 5) that has not been validated against N-body simulations; because the claimed gain is modest, even moderate model error could erase it. This matches the paper's own stated limitation, so conditional acceptance is the right call.","rationale":"I read the paper as making a Fisher-forecast claim: adding the octupole moment to the dipole improves the projected constraint on the LPI parameter alpha by 6% (smin=5) and 11% (smin=15). The forecast machinery is internally consistent: Eq. (6) uses the full covariance including off-diagonal dipole-octupole terms, Eq. (8) is a standard Gaussian covariance, and the numerical values in Fig. 4 follow from the stated model. The dipole part of the model is supported by earlier simulation work (Refs. [30,33]), which is real independent evidence. The weakest load-bearing element is the octupole prediction itself: Eq. (5) has not been validated in N-body simulations, and the authors say so in the Conclusions. This is not a manufactured objection; it is the exact condition needed for the 11% number to be reliable. My independent assessment agrees with the reader's weakest_assumption. I considered other possible concerns, such as line-of-sight definition, Gaussian covariance, and the Gaussian prior on Roff, but these are either standard, subdominant, or common to the dipole-only forecast, so they do not change the relative improvement. The honest finding is that the central claim is conditional on the octupole model being accurate. Since the reader already assigned CONDITIONAL with MODERATE confidence, no verdict change is needed. If the N-body check later shows the octupole amplitude differs substantially, the verdict should move to REJECT; if it confirms the model, the paper can be accepted.","tokens_in":1215,"tokens_out":958,"duration_ms":80824,"concrete_test":"Apply a sensitivity test to the Fisher forecast: rescale the alpha-dependent (Delta_epsilon_NL) part of the octupole in Eq. (5) by factors 0.8 and 1.2, and recompute sigma_alpha for all survey pairs and for smin=15 Mpc/h. If the average improvement relative to dipole-only drops below ~5% or reverses, the unvalidated octupole amplitude is load-bearing. Additionally (or instead), measure the octupole in an N-body simulation with the same halo occupation and off-centering model and compare its amplitude and scale dependence to Eq. (5) over s=5-30 Mpc/h; the simulation measurement is the decisive check.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim is the ~11% (smin=15 Mpc/h) and ~6% (smin=5 Mpc/h) improvement in sigma_alpha from adding the octupole to the dipole. The forecast ratio sigma_dipole+octupole/sigma_dipole is determined by the alpha-sensitivity vector of the octupole relative to the dipole in data space. The octupole model, Eq. (5), is taken from Ref. [33] and combines the gravitational-redshift term proportional to Delta_epsilon_NL with a wide-angle Doppler term proportional to Delta_b (s/d). While the dipole version of this model has been checked against N-body measurements, the octupole has not; the authors explicitly state in the Conclusions that 'a verification of octupole prediction requires a more proper simulation setup.' Because the octupole is a subdominant, higher-order signal, its amplitude and scale dependence are more sensitive than the dipole to the modeling of off-centering (Roff), small-scale velocity dispersion, and the linear-theory power spectrum over 5-30 Mpc/h. If the true octupole alpha-response differs by ~20-30% from Eq. (5), the 11% average improvement could shrink below significance or change sign for some survey pairs. Thus the central numerical claim rests on an unvalidated ingredient.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter extends a dipole-only Fisher forecast of the LPI-violating parameter alpha in galaxy cross-correlations to include the octupole moment. Using the analytic asymmetric-correlation model of Eq. (5), a Gaussian, plane-parallel covariance (Eq. 8), and survey specifications for DESI, Euclid, PFS, and SKA1/2, the authors find that adding the octupole to the dipole improves the marginalized 1-sigma error on alpha by 6% for smin=5 Mpc/h and 11% for smin=15 Mpc/h on average. The combined constraints are sigma_alpha ~ 0.028 and 0.14 for the two scale cuts. The authors present the smin=15 Mpc/h result as the conservative case and explicitly state in the Conclusions that the octupole model has not yet been verified with N-body simulations.","tokens_in":12413,"tokens_out":7769,"duration_ms":67865,"significance":"If the forecast is robust, this is a useful, low-cost forecast: it shows that a higher-order odd multipole can tighten an equivalence-principle test without new data. The Fisher formalism is transparent, the survey setup is clearly described, and the authors are commendably explicit about the main limitation. The central idea is worth publishing, but the numerical size of the claimed improvement rests on an unvalidated octupole model, and the combined constraints assume independence between overlapping survey pairs. These issues are addressable but are load-bearing rather than cosmetic.","major_comments":[{"comment":"The central improvement claim rests on the octupole model in Eq. (5), which the authors state in the Conclusions has not been verified with N-body simulations ('a verification of octupole prediction requires a more proper simulation setup'). Because the octupole is a subdominant signal and its alpha-response is not tested against simulations, a moderate error in its amplitude or scale dependence could erase or even reverse the 11% improvement. I ask the authors to either provide a simulation-level validation of the octupole model on the scales used (5-30 Mpc/h), or to demonstrate quantitatively that the forecast is robust to plausible mismodeling (e.g., by perturbing the octupole alpha-response within a range allowed by off-centering, velocity dispersion, and nonlinear power-spectrum uncertainties). Without one of these, the quoted 6-11% improvements are conditional on an unvalidated ingredient.","section":"Eq. (5) and Conclusions"},{"comment":"The covariance matrix in Eq. (8) is derived under the Gaussian and plane-parallel approximations, while the signal model in Eq. (5) includes wide-angle corrections. For the smin=5 Mpc/h case, the included separations lie in a regime where non-Gaussian and nonlinear contributions to the covariance are expected to be important; the authors caution that these 'could affect the constraints,' but they still quote sigma_alpha ~ 0.028 as a headline result. For smin=15 Mpc/h the scales are more benign, but the model still depends on nonlinear halo-potential and velocity-dispersion terms. Please quantify the impact of non-Gaussian covariance, or clearly demote the smin=5 Mpc/h result and present smin=15 Mpc/h as the headline.","section":"Eq. (8) and Results"},{"comment":"The combined constraints sigma_alpha ~ 0.028 and 0.14 are obtained by summing inverse variances from different cross-correlation pairs. These pairs are not independent: they share the same survey volume and, in several cases, one of the two galaxy samples (e.g., SKA2 appears in DESI-BGS x SKA2, DESI-LRG x SKA2, DESI-ELG x SKA2, and SKA2 x Euclid). The cross-covariance between the different pair measurements is nonzero, so simple inverse-variance addition can overstate the combined precision. Please either construct a joint Fisher matrix that includes cross-pair covariances, or present the pair-by-pair results as separate forecasts without an 'All combined' number.","section":"Results (Fig. 4)"}],"minor_comments":[{"comment":"The parameter alpha does not appear explicitly in Eq. (2) or Eq. (5); please state how the LPI-violating parameter enters the model (presumably through the gravitational-redshift term in epsilon_NL) so that the Fisher derivatives d xi / d alpha are unambiguous.","section":"Model and Forecast formalism"},{"comment":"The 'average improvement of 11%' is not defined; please specify whether it is the unweighted mean, inverse-variance weighted mean, or median of the ratios in the bottom panels of Fig. 4, and state over which survey pairs the average is taken.","section":"Results"},{"comment":"There are minor formatting and notation issues: 'smax = 30, Mpc/h' and 'smin = 5Mpc /h' should be typeset consistently, and the notation Roff,X/Y in the parameter vector theta is confusing (presumably Roff,X and Roff,Y).","section":"Setup"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the paper's central idea is worth publishing, but the headline improvement and the combined constraints are currently conditional on an unvalidated octupole model and on an independence assumption between overlapping survey pairs. These are fixable in a revision; I would not reject. The authors' explicit admission of the missing N-body verification is commendable but does not itself supply the missing validation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead this one if you care about equivalence-principle tests from galaxy clustering. It is a straightforward Fisher forecast that adds the octupole moment of the cross-correlation function to the dipole-only analysis of Saga et al. 2023. The new number is an 11% average improvement in the forecasted error on the LPI-violating parameter alpha when smin = 15 Mpc/h, and 6% at smin = 5 Mpc/h. That is the honest takeaway: a modest, not transformative, gain.\n\nWhat the paper does well: the Fisher formalism is laid out cleanly, the covariance matrix is specified, and the nuisance parameters (bias, off-centering) are handled with priors taken from prior work. The dipole-only results reproduce Ref. [42], which is a good sanity check. The authors also present the scale-cut dependence honestly, distinguishing the aggressive smin = 5 Mpc/h from the conservative smin = 15 Mpc/h. They explicitly note that non-Gaussian and nonlinear contributions could affect the smaller-scale results. No overfitting, no magical thinking.\n\nThe soft spot is the one the authors themselves flag in the Conclusions: the octupole model, Eq. (5), has not been verified with N-body simulations. The dipole model has been checked; the octupole has not. Since the claimed improvement is only 11%, even a moderate 20-30% error in the octupole's alpha-response could shrink that gain below significance. This is a real limitation, and the stress-test note is right that the central numerical claim rests on an unvalidated ingredient. It is not a fatal flaw for a forecast, but it means the forecast should be read as conditional on the model.\n\nThe self-citation pattern is not a problem here: the model and survey setup come from the authors' prior work, but the dipole has been simulation-tested, and the forecast is a legitimate extrapolation. The paper is honest about what is new and what is not.\n\nWho is this for? Specialists in relativistic effects in galaxy clustering who are planning LPI or gravitational-redshift measurements with DESI, PFS, Euclid, or SKA. They will find the numbers useful. General cosmologists can skip it.\n\nRecommendation: send it to peer review. The octupole validation issue should be explicitly raised as a required caveat or a point to address in revision, but it is not grounds for rejection. The paper is a competent, clearly written forecast that knows its own limits.\n\nBest,\n[You]","headline":"Clean, honest Fisher forecast extending the LPI dipole analysis to the octupole; the 11% gain is real but rests on an octupole model the authors themselves have not yet validated with simulations.","tokens_in":12950,"tokens_out":1636,"would_cite":false,"duration_ms":15925,"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":"Forecast: adding the octupole moment of the galaxy cross-correlation function to the dipole improves expected constraints on the LPI-violating parameter $\\alpha$ by 11% for conservative scale cuts ($s_{\\rm min}=15\\,{\\rm Mpc}/h$) and by 6%…","keywords":["local position invariance","gravitational redshift","galaxy cross-correlation function","octupole moment","dipole moment","Fisher forecast","Einstein equivalence principle","redshift-space distortions"],"falsifier":"Use two galaxy populations in an N-body simulation with realistic halo potentials, measure the octupole moment of their cross-correlation between 5 and 30 Mpc/h, and compare it with Eq. (5); if the measured octupole differs from the model by more than the forecast's 1-sigma band, the predicted 11% improvement would not be realized on those scales.","tokens_in":1649,"feed_emoji":"🔭","tokens_out":1999,"duration_ms":69064,"temperature":0.7,"pith_summary":"The paper argues that a higher-order statistical pattern in how two different galaxy populations are distributed, the octupole moment, carries extra information about a possible violation of local position invariance, the idea that non-gravitational physics is the same everywhere. It forecasts how well upcoming surveys would measure the violation parameter $\\alpha$, which scales deviations of the gravitational redshift from general relativity. The central result is that adding the octupole to the dipole tightens the expected measurement by about 11 percent when only larger separations are used, and by 6 percent when small separations are included. Because this improvement requires no additional data, it makes higher-order multipoles worth including in future equivalence-principle tests.","feed_headline":"Galaxy octupole moment sharpens equivalence-principle test by 11%","feed_subtitle":"Adding the octupole to dipole clustering tightens the gravitational-redshift bound by 11 percent, at no extra cost.","key_machinery":"The central object is the asymmetric part of the redshift-space cross-correlation function $\\xi^{XY}_{\\rm asym}$, expanded in odd Legendre multipoles; the dipole ($\\ell=1$) and octupole ($\\ell=3$) are the two terms used. Its model, Eq. (5), combines a plane-parallel piece from gravitational redshift, generated by a nonlinear halo potential, with standard Doppler terms and a leading $(s/d)$ wide-angle correction. The argument is carried by a Fisher matrix over five parameters, namely $\\alpha$, the off-centering radius $R_{\\rm off}$, and the two galaxy biases, with an analytic Gaussian covariance and survey-overlap volumes. The octupole's key property is that it lacks the $b_X b_Y$ term that dominates the dipole's real-space contribution, giving it a complementary bias dependence and therefore additional constraining power.","core_discovery":"On the paper's own terms, the discovery is that the octupole moment of the cross-correlation function between two galaxy samples is a useful and independent probe of the LPI-violating parameter $\\alpha$. This moment, which arises from the Doppler effect mixing with gravitational redshift from nonlinear halo potentials plus a leading wide-angle correction, has a different dependence on galaxy bias than the dipole does, so it breaks parameter degeneracies rather than merely repeating the dipole. In Fisher forecasts for surveys like DESI, Euclid, PFS, and SKA, adding the octupole to the dipole improves the expected $1\\sigma$ constraint on $\\alpha$ by an average of 6% for $s_{\\rm min}=5\\,{\\rm Mpc}/h$ and 11% for $s_{\\rm min}=15\\,{\\rm Mpc}/h$; combining all survey pairs yields $\\sigma_\\alpha \\approx 0.028$ and $\\approx 0.14$, respectively. The authors present this as the first use of the small-scale octupole in tests of local position invariance.","pith_inferences":["If N-body simulations confirm the octupole model, a natural next step is to include the $\\ell=5$ triakontadipole, which the paper neglects, in the same Fisher forecasts.","The octupole's sign behavior, which unlike the dipole shows no sign flip, could be used as a consistency check to separate a genuine $\\alpha$ signal from systematics.","Because the improvement is largest when small scales are excluded, the real gain depends on nonlinear modeling; better simulations could make aggressive scale cuts usable and yield improvements larger than the reported 6%.","The same Fisher setup can be converted into a likelihood for actual overlapping survey data, making the forecast directly testable with real measurements."],"forward_implications":["Expected LPI constraints improve by about 11% for $s_{\\rm min}=15\\,{\\rm Mpc}/h$ and 6% for $s_{\\rm min}=5\\,{\\rm Mpc}/h$ when the octupole is added to the dipole, without requiring new observations.","Cross-correlations between low-bias galaxy samples, such as DESI-BGS with SKA1/2, benefit the most from the octupole.","Combining all considered survey pairs gives $\\sigma_\\alpha \\approx 0.14$ for conservative scales and $\\sigma_\\alpha \\approx 0.028$ for aggressive scales.","The octupole alone can constrain $\\alpha$ below roughly 0.5 for some survey pairs when small scales are used.","Future equivalence-principle tests based on galaxy clustering should include odd multipoles beyond the dipole."],"supporting_citations":[{"why":"Supplies the analytic model of the dipole and octupole moments, including Eq. (5), on which the forecasts are built.","marker":"[33]"},{"why":"Establishes the dipole-only forecast pipeline, parameter priors, and survey-overlap setup that this paper extends to the octupole.","marker":"[42]"},{"why":"Provides the N-body measurements of the dipole that the analytic model is known to reproduce.","marker":"[30]"},{"why":"Develops the nonlinear halo-potential model for gravitational redshift and velocity dispersion used in the correlation function.","marker":"[32]"},{"why":"A recent independent study of the octupole via the wide-angle effect, which this paper contrasts with its small-scale octupole.","marker":"[50]"},{"why":"Defines the DESI survey samples (BGS, LRG, ELG) whose number densities and biases enter the forecasts.","marker":"[37]"},{"why":"Defines the Subaru PFS [O II] ELG sample used in the cross-correlation forecasts.","marker":"[36]"},{"why":"Defines the Euclid H-alpha emitter sample used in the cross-correlation forecasts.","marker":"[38]"},{"why":"Defines the SKA1 and SKA2 HI galaxy samples used in the cross-correlation forecasts.","marker":"[39]"}],"fun_headline_variants":["Octupole galaxy clustering tightens equivalence-principle test by 11%","Octupole boosts gravitational-redshift test of equivalence principle by 11%","Adding octupole to dipole sharpens cosmology test of Einstein equivalence","Octupole galaxy moment improves local position invariance test by 11%","Higher-order odd multipole of galaxy clustering sharpens LPI test 11%"],"cache_read_input_tokens":15104,"weakest_assumption_plain":"The octupole's predicted signal, taken from an analytic model that the paper itself says has not yet been verified with a proper simulation setup, is accurate on the 5-30 Mpc/h scales used; if it is not, the forecasted improvement changes.","fun_headline_variants_meta":{"raw":{"variants":["Octupole galaxy clustering tightens equivalence-principle test by 11%","Octupole boosts gravitational-redshift test of equivalence principle by 11%","Adding octupole to dipole sharpens cosmology test of Einstein equivalence","Octupole galaxy moment improves local position invariance test by 11%","Higher-order odd multipole of galaxy clustering sharpens LPI test 11%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000652,"raw_usage":{"total_tokens":3032,"prompt_tokens":1031,"completion_tokens":2001,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":647,"completion_tokens_details":{"reasoning_tokens":1901}},"tokens_in":647,"tokens_out":2001,"duration_ms":12657,"temperature":1.0,"reasoning_tokens":1901,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:52:19.421977+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use two galaxy populations in an N-body simulation with realistic halo potentials, measure the octupole moment of their cross-correlation between 5 and 30 Mpc/h, and compare it with Eq. (5); if the measured octupole differs from the model by more than the forecast's 1-sigma band, the predicted 11% improvement would not be realized on those scales.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the analytic model of the dipole and octupole moments, including Eq. (5), on which the forecasts are built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the dipole-only forecast pipeline, parameter priors, and survey-overlap setup that this paper extends to the octupole."},{"cited_title":"Breton, Y","cited_arxiv_id":null,"evidence_quote":"Provides the N-body measurements of the dipole that the analytic model is known to reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Develops the nonlinear halo-potential model for gravitational redshift and velocity dispersion used in the correlation function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"A recent independent study of the octupole via the wide-angle effect, which this paper contrasts with its small-scale octupole."},{"cited_title":"Takada, R","cited_arxiv_id":null,"evidence_quote":"Defines the Subaru PFS [O II] ELG sample used in the cross-correlation forecasts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the SKA1 and SKA2 HI galaxy samples used in the cross-correlation forecasts."}],"review_version":1}