{"id":"12b89a7a-eabd-4e51-91a7-6aa1e344bf74","arxiv_id":"2605.30367","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Gaia quasar proper motions give a Solar System acceleration amplitude of 5.72 μas yr^{-1} with vector components (0.40, -5.09, -2.40), consistent with prior work but with 1.5-2.5 times larger credible intervals after marginalizing higher multipoles.","lead":"The paper analyzes the dipole pattern in proper motions of distant quasars observed by Gaia to infer the acceleration of the Solar System relative to the quasar rest frame. This yields an updated measurement with wider uncertainties after accounting for higher-order angular patterns in the data.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.3","headline":"Cross-correlations may miss unmodeled systematics that bias the marginalised dipole","rationale":"The reader's weakest_assumption is exactly the load-bearing step; the abstract-only limitation noted in the rationale reinforces that the claim cannot be accepted without verification of the cross-correlation completeness. No other internal inconsistency is visible from the given material.","tokens_in":1847,"tokens_out":360,"duration_ms":19328,"concrete_test":"Inject a mock systematic field into the Quaia catalogue that is correlated with an unmodelled map (e.g., a redshift-binned density variation orthogonal to the supplied stellar-density and scanning-law templates) at the 10 % level of the observed higher-multipole power; re-run the SBI pipeline and check whether any dipole component shifts by >1σ relative to the reported (0.40, -5.09, -2.40) values.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (widened but still consistent dipole, no redshift dependence, kinematic origin) rests on the assertion that cross-correlations with the Gaia scanning law, stellar density, and stellar proper-motion maps are sufficient to diagnose and remove all non-kinematic contributions to the higher multipoles before simulation-based marginalisation. If any residual power remains (e.g., from redshift-dependent quasar selection, colour-dependent astrometric systematics, or scanning-law harmonics not captured by the stellar maps), it can still couple into the dipole even after joint inference, because the pseudo-C_ℓ estimator and SBI prior do not automatically null those modes. The abstract provides no quantitative test of completeness (e.g., null tests on injected systematics or power-spectrum residuals after subtraction), leaving this the weakest link in the argument.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper measures the Solar System's acceleration vector via the dipole in quasar proper motions using Gaia EDR3 and Quaia catalogues. It characterizes the full angular power spectrum with the pseudo-C_ℓ estimator and uses simulation-based inference (SBI) to jointly constrain the dipole (g_x, g_y, g_z) and higher multipoles, applying cross-correlations with Gaia scanning strategy, stellar density, and stellar proper-motion maps to diagnose systematics. The inferred acceleration is consistent with prior results but with credible intervals widened by factors of 1.5–2.5; the best estimate from Quaia is (0.40^{+0.70}_{-0.70}, -5.09^{+0.54}_{-0.54}, -2.40^{+0.55}_{-0.58}) μas yr^{-1} (amplitude 5.72 μas yr^{-1}), with no significant redshift dependence, supporting a kinematic origin.","tokens_in":2027,"tokens_out":731,"duration_ms":15665,"significance":"If the systematic control holds, the result provides a more conservative and robust determination of the acceleration by properly marginalizing higher-multipole degeneracies that previous analyses appear to have neglected, yielding wider but more reliable uncertainties. The absence of redshift dependence adds supporting evidence for the kinematic interpretation. The use of SBI for joint inference and explicit cross-correlation diagnostics are methodological strengths that could be adopted more broadly in astrometric dipole studies.","major_comments":[{"comment":"Abstract and framework description: the assertion that cross-correlations with the Gaia scanning law, stellar density, and stellar proper-motion maps are sufficient to diagnose and remove all non-kinematic contributions to higher multipoles is load-bearing for the central claim of a clean kinematic dipole. No quantitative completeness test (e.g., power-spectrum residuals after subtraction, null tests on injected systematics, or recovery fractions for simulated non-kinematic modes) is described, leaving open the possibility that residual power couples into the dipole even after SBI marginalization.","section":"Abstract, framework description"},{"comment":"Results section (Quaia best estimate): the reported credible intervals already incorporate higher-multipole marginalization, but without an explicit comparison table or figure showing the dipole posterior with versus without the higher-multipole model (or with versus without the cross-correlation cleaning step), it is difficult to quantify how much of the factor 1.5–2.5 widening is due to each component.","section":"Results"},{"comment":"Methods (pseudo-C_ℓ + SBI): the pseudo-C_ℓ estimator and SBI prior must be shown to be insensitive to any unmodeled redshift-dependent selection or colour-dependent astrometric systematics that could survive the stellar-map cross-correlations; a dedicated injection-recovery test on the full pipeline would directly address this.","section":"Methods"}],"minor_comments":[{"comment":"Notation: the amplitude is quoted as 5.72_{-0.52}^{+0.53} while the vector components use asymmetric errors; a brief statement on how the amplitude posterior is derived from the vector components would improve clarity.","section":"Abstract"},{"comment":"Figure clarity: the cross-correlation maps and power spectra should include explicit labels for which multipole ranges are used in the SBI fit versus those used only for diagnostics.","section":"Figures"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading and constructive comments. We address each major comment below and describe the revisions we will implement.","responses":[{"response":"We agree that quantitative validation strengthens the central claim. We will add a new subsection to the Methods describing power-spectrum residuals after modeled subtraction and injection-recovery tests using simulated non-kinematic modes. These tests will quantify residual coupling into the dipole after SBI marginalization and support the diagnostic power of the cross-correlations.","revision_made":"yes","referee_comment":"[Abstract, framework description] Abstract and framework description: the assertion that cross-correlations with the Gaia scanning law, stellar density, and stellar proper-motion maps are sufficient to diagnose and remove all non-kinematic contributions to higher multipoles is load-bearing for the central claim of a clean kinematic dipole. No quantitative completeness test (e.g., power-spectrum residuals after subtraction, null tests on injected systematics, or recovery fractions for simulated non-kinematic modes) is described, leaving open the possibility that residual power couples into the dipole even after SBI marginalization."},{"response":"We will add a new figure to the Results section comparing the dipole posterior under four cases: full model, without higher-multipole marginalization, without cross-correlation cleaning, and without both. This will explicitly quantify the contribution of each element to the reported widening of the credible intervals.","revision_made":"yes","referee_comment":"[Results] Results section (Quaia best estimate): the reported credible intervals already incorporate higher-multipole marginalization, but without an explicit comparison table or figure showing the dipole posterior with versus without the higher-multipole model (or with versus without the cross-correlation cleaning step), it is difficult to quantify how much of the factor 1.5–2.5 widening is due to each component."},{"response":"We will add a dedicated injection-recovery test to the Methods section. Mock catalogues will be generated with injected redshift-dependent selection and colour-dependent astrometric systematics, processed through the full pseudo-C_ℓ + SBI pipeline, and the recovered dipole parameters examined for bias after cross-correlation cleaning.","revision_made":"yes","referee_comment":"[Methods] Methods (pseudo-C_ℓ + SBI): the pseudo-C_ℓ estimator and SBI prior must be shown to be insensitive to any unmodeled redshift-dependent selection or colour-dependent astrometric systematics that could survive the stellar-map cross-correlations; a dedicated injection-recovery test on the full pipeline would directly address this."}],"tokens_in":1668,"tokens_out":552,"duration_ms":26692,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The core result here is that fitting the full proper motion power spectrum and marginalizing over higher multipoles with simulation-based inference increases the credible intervals on the acceleration vector by 1.5–2.5 times compared with earlier dipole-only analyses, while the central value stays consistent. Their Quaia-based numbers are (0.40, -5.09, -2.40) μas yr⁻¹ with amplitude 5.72 μas yr⁻¹ and no detected redshift trend.\n\nThey do two things that are worth noting. First, they actually run the joint inference instead of assuming the dipole is isolated. Second, they apply the same pipeline to both the Gaia EDR3 quasar list and the Quaia catalogue and show the results line up. That is straightforward and useful for anyone who needs a more conservative uncertainty on the local acceleration.\n\nThe weakest part is the claim that cross-correlations with the scanning law, stellar density, and stellar proper motions are enough to remove non-kinematic power from the higher multipoles. The abstract gives no numbers on residual power after subtraction or on null tests with injected systematics, so it is not clear how much leakage into the dipole remains possible. Without those checks the marginalised intervals could still be optimistic.\n\nThe work is incremental but cleanly executed on public data. It is mainly for people who already follow Gaia astrometry and galactic dynamics papers and want the updated error budget. The methods are standard enough that a referee could evaluate them quickly. I would send it to peer review.","headline":"The paper widens the error bars on the Solar System acceleration dipole by marginalizing higher multipoles with SBI, but the cross-correlation checks for systematics look incomplete.","tokens_in":2490,"tokens_out":390,"would_cite":false,"duration_ms":14872,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Quasar proper motions reveal Solar System acceleration at 5.72 μas per year, with uncertainties 1.5-2.5 times larger after accounting for higher multipoles.","keywords":["quasar proper motions","solar system acceleration","Gaia catalogue","proper motion dipole","angular power spectrum","simulation-based inference","astrometric systematics"],"falsifier":"Detection of statistically significant redshift dependence in the measured acceleration amplitude or a result inconsistent with independent determinations of Solar System motion relative to the cosmic microwave background.","tokens_in":2741,"feed_emoji":"🌌","tokens_out":725,"duration_ms":16826,"temperature":0.7,"pith_summary":"The paper measures the dipole in quasar proper motions observed by Gaia to infer the acceleration of the Solar System relative to the distant quasar frame. It models the full angular power spectrum using the pseudo-C_ℓ approach and simulation-based inference to fit the dipole jointly with higher multipoles, then uses cross-correlations with scanning strategy and stellar maps to check for systematics. Applied to Gaia EDR3 and Quaia catalogues, the acceleration matches earlier values but with substantially wider credible intervals, and shows no redshift dependence.","feed_headline":"Quasar motions give Solar acceleration 5.72 μas/yr with wider errors","feed_subtitle":"Higher multipoles enlarge credible intervals 1.5-2.5 times yet preserve the dipole and show no redshift trend","key_machinery":"The pseudo-C_ℓ formalism paired with simulation-based inference that jointly constrains the dipole and higher multipole power in the proper motion field while using cross-correlations to diagnose systematics.","core_discovery":"The acceleration of the Solar System is (0.40^{+0.70}_{-0.70}, -5.09^{+0.54}_{-0.54}, -2.40^{+0.55}_{-0.58}) μas yr^{-1} with amplitude 5.72 μas yr^{-1} from the Quaia catalogue. This value is consistent with prior determinations, yet the credible intervals widen by factors of 1.5 to 2.5 once higher-multipole degeneracies are marginalised, indicating earlier uncertainty estimates were optimistic. The signal exhibits no significant redshift dependence, supporting its kinematic origin.","pith_inferences":["Similar higher-multipole marginalisation may be needed in other large-scale astrometric dipole studies to avoid underestimating errors.","The widened intervals could shift how this acceleration measurement is combined with galactic dynamics or local group motion models.","If non-kinematic residuals remain after the cross-correlation tests, they would most likely appear as excess power at specific multipoles tied to the scanning law."],"forward_implications":["Earlier published uncertainties on the Solar System acceleration were optimistic by factors of 1.5-2.5.","The absence of redshift dependence strengthens the case that the dipole is kinematic rather than systematic.","Joint modelling of dipole and higher multipoles is required for unbiased dipole inference in future astrometric catalogues.","The framework can be reapplied to larger or deeper quasar samples to tighten the acceleration constraints."],"fun_headline_variants":["Quasar proper motions measure solar acceleration of 5.72 μas/yr","Higher multipoles increase solar acceleration error bars by 1.5-2.5x","Solar accel 5.72 μas/yr from Quaia quasars shows no redshift trend","Gaia quasar data widens solar system acceleration uncertainties"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"Cross-correlations with Gaia scanning strategy, stellar density, and stellar proper motion maps suffice to identify and remove non-kinematic contributions to higher multipoles that could bias the dipole.","fun_headline_variants_meta":{"raw":{"variants":["Quasar proper motions measure solar acceleration of 5.72 μas/yr","Higher multipoles increase solar acceleration error bars by 1.5-2.5x","Solar accel 5.72 μas/yr from Quaia quasars shows no redshift trend","Gaia quasar data widens solar system acceleration uncertainties"]},"model":"grok-4.3","cost_usd":0.006645,"raw_usage":{"total_tokens":3147,"prompt_tokens":763,"num_sources_used":0,"completion_tokens":83,"cost_in_usd_ticks":66449500,"prompt_tokens_details":{"text_tokens":763,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2301,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":763,"tokens_out":83,"duration_ms":17009,"temperature":1.0,"reasoning_tokens":2301,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-30T18:34:10.953986+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Detection of statistically significant redshift dependence in the measured acceleration amplitude or a result inconsistent with independent determinations of Solar System motion relative to the cosmic microwave background.","supporting_citations":[],"review_version":1}