{"id":"1e79d644-4fed-457d-93b2-c3222be4a513","arxiv_id":"2607.06692","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Second and third moments of joint weak-lensing, tSZ, and X-ray maps improve primordial non-Gaussianity constraints by a factor of ~2 relative to lensing alone while breaking baryon and foreground degeneracies.","lead":"Multi-wavelength maps of weak lensing, thermal SZ, and X-rays can tighten constraints on early-Universe particle physics by a factor of two over lensing alone. The work shows how joint moments of these fields self-calibrate astrophysical nuisances while probing primordial non-Gaussianity.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"The factor-of-two gain is measured only under the paper’s own baryon SAM + halo-painting + transfer-function pipeline; residual model misspecification could re-open the very degeneracies that produce the gain.","rationale":"The reader correctly isolates the semi-analytic multi-wavelength pipeline as the weakest link. The forecast is internally consistent, the code and simulations are public, and the paper is transparent about resolution, IA, and foreground limitations. No internal contradiction or numerical instability undermines the reported Fisher ratios under the stated model. The load-bearing concern is therefore precisely the one the reader named: whether that model is complete enough for the degeneracy-breaking (and thus the factor-of-two) to survive when the true sky differs from the SAM. A single controlled misspecification test would settle whether the concern lands. Verdict remains CONDITIONAL; no stronger or weaker adjustment is warranted.","tokens_in":45211,"tokens_out":608,"duration_ms":5859,"concrete_test":"Re-run the full Fisher pipeline of Section 4 (same moments, same scale cuts, same extended SAM) after injecting one controlled misspecification: e.g., a Rankine–Hugoniot shock feature in the gas density at ~R_200c (as discussed in A.2) or a 50/50 cool-core/NCC split that alters only the core X-ray emissivity. If the marginalized σ(f_NL) for the “All” combination degrades by ≳10% relative to the baseline in Fig. 3, the claimed factor-of-two improvement is not robust to the residual systematics the paper itself flags.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (factor-of-two improvement on f_NL from adding tSZ/X-ray moments while marginalizing the extended baryon SAM) rests on the assumption that BaryonForge profiles, CIB transfer functions (Eq. A.9–A.10), radio HOD/LF (Eq. A.11–A.13), AGN LFs, and the two-parameter metallicity profile (Eq. A.19) fully span the correlated astrophysical systematics that couple the three fields. Section 4.5 and Appendices A.2–A.3 openly list missing physics (accretion shocks, cool-core/NCC bimodality, extended IA beyond NLA, CO, filtering). Because the Fisher matrix (Eq. 4.1) and the degeneracy-breaking shown in Fig. 4 are computed entirely inside this model family, any residual freedom outside it can re-correlate f_NL with the nuisance directions that the multi-probe combination is claimed to self-calibrate. The factor-of-two number is therefore conditional on model completeness at the precision of the forecast, not merely on statistical noise.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper forecasts constraints on several primordial non-Gaussianity (PNG) templates using second and third moments of weak-lensing convergence (LSST Y10), thermal SZ (enhanced SO + Planck), and X-ray count-rate (eROSITA DR3-like) maps. It builds multi-wavelength synthetic skies from the Ulagam N-body lightcones with Aarambam PNG initial conditions, applies a consistent baryon semi-analytic model (BaryonForge) across density, pressure, and X-ray emissivity, and correlates CIB, radio, and AGN foregrounds with the simulated density field. Using Fisher matrices with Hartlap correction and sequential marginalization over cosmology, IA, and an extended baryon SAM, the central claim is that the joint (WL+tSZ+X-ray) analysis recovers a factor-of-two improvement on f_NL relative to WL alone while self-calibrating nuisance parameters. Public code (Vaanam, BaryonForge, Aarambam) and simulation products are released.","tokens_in":45582,"tokens_out":1469,"duration_ms":18679,"significance":"If the forecast holds under realistic residual systematics, the work would elevate PNG constraints toward the multi-probe status already enjoyed by LCDM analyses, and would demonstrate that tSZ and X-ray moments can both add PNG information and break baryon/IA degeneracies that otherwise degrade WL-only forecasts by factors of 2–3. The public multi-wavelength map-maker and the consistent baryon SAM across three observables are concrete community assets. The numerical-convergence tests (Appendix C), Hartlap-corrected covariances, and explicit sequential marginalization over a large nuisance set are strengths relative to many Fisher forecasts in the literature. The result is therefore of clear interest for Stage-IV multi-wavelength cosmology even if the precise factor-of-two number is model-conditional.","major_comments":[{"comment":"Abstract and §4.2 / Fig. 3: The headline claim of a factor-of-two improvement on f_NL (joint vs WL-only, after full nuisance marginalization) is computed entirely inside the BaryonForge + transfer-function CIB (Eqs. A.9–A.10) + radio HOD/LF (Eqs. A.11–A.13) + AGN LF + two-parameter metallicity (Eq. A.19) model family. Section 4.5 and Appendices A.2–A.3 list missing physics (accretion shocks, cool-core/NCC bimodality, extended IA beyond NLA, CO, map filtering) that can re-correlate f_NL with the same nuisance directions the multi-probe combination is said to self-calibrate. The abstract and conclusions should state explicitly that the factor-of-two gain is conditional on completeness of this model family at the forecast precision, and should quantify (or bound) how residual freedom outside the varied parameters would degrade the gain, rather than only arguing O(10%) effects from amplitude","section":null},{"comment":"§3.3 and Appendix A.3: The X-ray forward model is described as “more simplistic” than the WL and tSZ pipelines and as requiring further investigation (masking, background subtraction, cool-core bimodality) before data application. Yet the “All” and “All++” configurations in Fig. 3 and the degeneracy-breaking narrative in Fig. 4 assign non-negligible weight to X-ray auto- and cross-moments. Either (i) demonstrate that the reported joint improvement remains stable when the X-ray cluster component is down-weighted or removed, or (ii) clearly separate the tSZ-driven gain from the X-ray-driven gain so that the central claim does not rest on the least mature map model.","section":null},{"comment":"§2.1 and §4.5: Maps are limited to NSIDE=1024 and halo catalogs to M_200c ≳ 10^14 M_⊙/h. The text correctly notes that this underestimates tSZ/X-ray signal and small-scale WL power, but the Fisher information (Eq. 4.1) and the scale-cut tests (θ > 10′, 25′) are still presented as representative of LSST Y10 / SO / eROSITA. Because the PNG signal is argued to live in peaks and one-halo regimes, a quantitative estimate of how much the factor-of-two ratio would change under higher resolution (or a statement that the ratio is conservative) is needed for the claim to be load-bearing rather than illustrative.","section":null}],"minor_comments":[{"comment":"Table 1: Several baryon SAM parameters (e.g. γ, δ, μ_θej) appear only in the extended set; a short cross-reference to A24 for their physical meaning would help readers who do not have that paper open.","section":null},{"comment":"Fig. 2: Planck comparison lines are useful; please state in the caption whether they are 68% marginalized constraints on the same templates (and same mass parameters ν, μ) as used here.","section":null},{"comment":"§2.3: The sparse fourth-moment subset used in “All++” is defined in prose; a compact equation or bullet list would make the datavector reproducible without re-reading the paragraph.","section":null},{"comment":"Appendix D: The correlated Poisson sampling is a nice technical contribution; a one-sentence pointer in the main text (§3.3 or §4) would help readers who skip the appendices understand why X-ray noise does not artificially inflate Fisher information.","section":null},{"comment":"Typos / notation: “propogating” → “propagating” (footnote 7); “vaccuum” → “vacuum” (NBD2 subsection); “derivate” → “derivative” (§4.1); “metallicty” → “metallicity” (§4.2).","section":null},{"comment":"§4.4 / Fig. 5: Clarify whether the Planck prior is applied as a diagonal Gaussian on (σ_8, Ω_m) only, or includes the full Planck covariance; this affects the quoted 30% degradation.","section":null}],"recommendation":"minor_revision","confidential_remarks":"Single-author thesis-style paper with heavy self-citation to the author’s prior simulation and baryon infrastructure papers; the novelty is in the multi-wavelength PNG forecast rather than in new simulation methods per se. That is fine for the journal if the forecast is clearly scoped as conditional. I do not see circularity or fabrication issues. Fit to astro-ph.CO / multi-probe cosmology is good."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The headline result is a clean factor-of-two improvement on f_NL (across Local, Equilateral, QSF, SI/SII, NBD2) when second and third moments of tSZ and X-ray are added to LSST Y10 lensing, while marginalizing cosmology, NLA IA, and an extended baryon SAM plus foreground amplitudes. That is new relative to the single-probe or WL+tSZ forecasts in the literature, and it is measured with a large Ulagam suite, Hartlap-corrected covariances, and explicit convergence tests in Appendix C.\n\nWhat the paper does well is the infrastructure. The same BaryonForge profiles feed density, pressure, and X-ray emissivity; CIB, radio, and AGN are painted from the density shells so the correlations are not free parameters; Vaanam, BaryonForge, Aarambam, and the light-cones are public. Degeneracy breaking on Mc, theta_ej, and alpha_nt is shown directly in the Fisher contours (Fig. 4). The sequential marginalization (Fig. 2–3) and the LCDM-prior tests (Fig. 5) are careful. Self-citations supply the simulation and baryon machinery, not the target claim.\n\nThe soft spot is exactly the one the stress-test flags, and the paper itself flags it in 4.5 and A.2–A.3: the factor-of-two lives inside this model family. Missing pieces (accretion shocks, cool-core bimodality, TATT-level IA, CO, filtering, NSIDE=1024, M>10^14 cut) could re-open some of the degeneracies the multi-probe combination is claimed to self-calibrate. That does not make the number fake; it makes it conditional on model completeness at the precision of the forecast. The author is transparent about this and does not oversell.\n\nThis is for people who already care about PNG forecasts or multi-wavelength baryon self-calibration. The math and citation pattern look solid; free parameters are listed and varied. I would send it to referees without hesitation. Engage if you work on PNG or on joint WL/SZ/X-ray pipelines; the public code alone is useful.","headline":"Solid multi-probe Fisher forecast showing a factor-of-two PNG gain from WL+tSZ+X-ray moments under a consistent baryon SAM; the number is real inside the model family, residual misspecification is the main caveat.","tokens_in":46160,"tokens_out":571,"would_cite":true,"duration_ms":6920,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Joint moments of weak lensing, tSZ, and X-ray maps tighten primordial non-Gaussianity constraints by a factor of two over lensing alone.","keywords":["primordial non-Gaussianity","weak lensing","thermal Sunyaev-Zeldovich","X-ray","baryon semi-analytic models","higher-order moments","multi-wavelength cosmology","Fisher forecasts"],"falsifier":"Apply the same multi-wavelength moment pipeline to real overlapping maps (for example DES/DECADE lensing + SPT/SO Compton-y + eROSITA counts) and test whether the recovered PNG posteriors tighten by the predicted factor of two relative to the lensing-only baseline after the same nuisance marginalization.","tokens_in":46090,"feed_emoji":"🌌","tokens_out":888,"duration_ms":7761,"temperature":0.7,"pith_summary":"Primordial non-Gaussianities leave imprints on the initial density field that later shape the statistics of peaks and massive halos. This paper shows that those imprints remain measurable in the second and third moments of three late-time fields: weak-lensing convergence, the thermal Sunyaev–Zeldovich effect, and X-ray photon counts. By forward-modeling all three observables from the same N-body lightcones with a shared semi-analytic baryon prescription and correlated foregrounds, the author finds that a joint analysis recovers twice the constraining power on the amplitude of several PNG templates relative to lensing alone, even after marginalizing over cosmology, intrinsic alignments, and an extended set of baryonic nuisance parameters. The multi-wavelength combination also self-calibrates many of those nuisance parameters by breaking degeneracies that pure lensing cannot resolve. The result matters because it converts existing and near-term multi-wavelength surveys into a practical route for probing particle physics of the early Universe on scales that the CMB and galaxy clustering have not yet reached.","feed_headline":"Lensing + tSZ + X-ray doubles PNG constraints","feed_subtitle":"Joint moments of three fields self-calibrate baryons and beat pure lensing by a factor of two.","key_machinery":"Simulation-based Fisher forecasts that paint consistent semi-analytic baryon profiles onto shared N-body lightcones, generate correlated multi-wavelength maps (weak lensing, Compton-y, X-ray counts), and extract second- and third-order aperture moments.","core_discovery":"Using second and third moments of the joint weak-lensing, thermal SZ, and X-ray fields yields a factor-of-two improvement in constraints on the amplitude of several primordial non-Gaussianity templates relative to a pure lensing analysis, while still marginalizing over cosmology, intrinsic alignments, and an extended baryon semi-analytic model.","pith_inferences":["Because the improvement is driven by degeneracy breaking around massive-halo peaks, other peak-sensitive probes such as cluster counts or kSZ may yield comparable multiplicative gains when added under the same consistent baryon model.","The public multi-wavelength map maker lowers the barrier for other groups to test whether higher-order statistics or machine-learned summaries can extract still more PNG information from the same simulated skies.","If residual model error in the baryon or foreground prescriptions re-introduces degeneracies, the claimed factor-of-two gain would shrink first for the most equilateral-like templates, which rely most heavily on intermediate and small scales."],"forward_implications":["PNG constraints competitive with or better than Planck become accessible from existing multi-wavelength survey overlaps once moments of tSZ and X-ray are co-analyzed with lensing.","Baryonic and foreground nuisance parameters that degrade pure-lensing forecasts can be self-calibrated by the joint data vector rather than left as free systematics.","The same map-making pipeline can be extended to higher resolution, larger sky area, and additional summary statistics to push further gains.","Models whose PNG signatures peak on small scales (for example certain non-Bunch–Davies templates) gain relatively more from the multi-probe combination."],"fun_headline_variants":["Lensing+tSZ+X-ray moments double PNG constraints","Joint 2nd-3rd moments of three fields double PNG reach","tSZ and X-ray double PNG limits over pure lensing","Multi-field moments self-calibrate baryons, double PNGs","Weak lensing plus SZ/X-ray doubles primordial NG power"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The semi-analytic baryon model and the adopted prescriptions for CIB, radio, AGN, and metallicity fully capture the correlated astrophysical systematics at the precision needed for unbiased forecasts.","fun_headline_variants_meta":{"raw":{"variants":["Lensing+tSZ+X-ray moments double PNG constraints","Joint 2nd-3rd moments of three fields double PNG reach","tSZ and X-ray double PNG limits over pure lensing","Multi-field moments self-calibrate baryons, double PNGs","Weak lensing plus SZ/X-ray doubles primordial NG power"]},"model":"grok-4.5","effort":"low","cost_usd":0.005586,"raw_usage":{"total_tokens":1507,"prompt_tokens":816,"num_sources_used":0,"completion_tokens":95,"cost_in_usd_ticks":55860000,"prompt_tokens_details":{"text_tokens":816,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":596,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":816,"tokens_out":95,"duration_ms":4932,"temperature":1.0,"reasoning_tokens":596,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T06:50:29.274082+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Apply the same multi-wavelength moment pipeline to real overlapping maps (for example DES/DECADE lensing + SPT/SO Compton-y + eROSITA counts) and test whether the recovered PNG posteriors tighten by the predicted factor of two relative to the lensing-only baseline after the same nuisance marginalization.","supporting_citations":[],"review_version":2}