REVIEW 3 major objections 5 minor 53 references
Template-free symbolic regression finds that CMB data prefer an inverse-k primordial oscillation over standard linear and log templates.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-11 11:22 UTC pith:4YR7BDNU
load-bearing objection Template-free SR recovers the same inverse-k form on Planck and SPA; the result is real but weak (~1σ) and conditioned on a fixed baseline cosmology. the 3 major comments →
Inverse-k Primordial Oscillations from a Symbolic Regression Search
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Both Planck and the combined Planck+ACT+SPT-3G datasets independently select an inverse-k oscillation, cos(B/k) or sin(B/k) with B ≃ 4 Mpc^{-1}, as the leading low-complexity feature on the symbolic-regression Pareto front. On the combined data this inverse-k template yields the largest likelihood improvement over a featureless power-law spectrum (best-fit Δχ² ≃ −12.6) and a weak one-sigma preference for non-zero amplitude, whereas the conventional linear and logarithmic templates remain consistent with zero.
What carries the argument
Symbolic regression (PySR) that searches the space of analytic expressions for the feature f(k) multiplying the standard power-law primordial spectrum, using the CMB χ² as the sole objective while background and nuisance parameters are held fixed; the Pareto front of complexity versus fit improvement then isolates the inverse-k oscillation as the simplest competitive form.
Load-bearing premise
All cosmological background and nuisance parameters are locked to their best-fit power-law values during the search, so the recovered feature form is conditioned on that fixed baseline rather than being re-optimized jointly with the rest of the cosmology.
What would settle it
A full MCMC re-analysis that frees the background and nuisance parameters together with the inverse-k amplitude, frequency and phase, and checks whether the preference for non-zero amplitude and the superiority over linear and log templates survive on the same SPA data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper performs a template-free search for oscillatory features in the primordial power spectrum using symbolic regression (PySR) with CMB χ² as the objective. Both Planck and the combined SPA (Planck+ACT+SPT-3G) datasets independently recover an inverse-k oscillation of the form cos(B/k) or sin(B/k) with B ≃ 4 Mpc^{-1} as the leading low-complexity expression (complexity ~8) on the Pareto front. The authors then promote this form to a three-parameter template and compare it via MCMC to the standard linear and logarithmic oscillation templates on SPA, finding that inverse-k yields the largest improvement (Δχ² ≃ -12.6) and a weak 1σ preference for nonzero amplitude, while the other templates remain consistent with zero. A heuristic link to standard-clocks signals with p = -1 is noted.
Significance. If the inverse-k preference survives joint cosmological and nuisance variation, the work supplies a concrete, data-driven candidate feature beyond the usual linear/log templates and demonstrates that symbolic regression can surface interpretable analytic forms in CMB cosmology without hand-crafted templates. The dual-dataset convergence of the C8 expressions, the residual plots matching ACT/SPT polarization structure, and the explicit Pareto-front tables are genuine strengths. The result is of interest for early-universe model building (NEC-violating or phantom-like phases) and for future small-scale CMB experiments, even though the present preference remains only weak.
major comments (3)
- Methods and Supplemental Material: the entire SR search freezes all background cosmological and nuisance parameters at their power-law best-fit values, so the χ² objective varies only f(k). The claim that this is harmless is imported from the authors’ prior template analysis [27] and is not re-validated inside the SR pipeline. Because the headline result is that both datasets independently select cos(B/k) (or sin(B/k)) as the leading low-complexity feature, a joint MCMC or profile-likelihood check that frees at least the main ΛCDM parameters (and key ACT/SPT nuisances) while retaining the recovered inverse-k form is load-bearing; without it the “template-free discovery” remains conditioned on a fixed baseline.
- Table III and footnote [44]: the three templates share the same flat prior ω_X ∈ [1,100], yet g_lin, g_log and g_inv map this interval onto physically inequivalent frequency ranges. The paper itself notes that the comparison is therefore heuristic. Because the ranking of inverse-k over linear/log (Table IV, Fig. 3) is a central claim, either re-run the MCMC with priors that equalize the effective frequency content or quantify how much of the Δχ² advantage is prior-driven.
- SR results and Supplemental Material (resolution estimate): higher-complexity Planck expressions containing 1/k^{2} terms are discarded as numerically unreliable because they oscillate faster than the CLASS sampling (k_per_decade = 400). The same sampling criterion should be applied quantitatively to the C8 inverse-k solutions that carry the main claim; an explicit statement of the maximum reliable B (or a higher-resolution re-evaluation of Eqs. 2–3) is needed to confirm that the recovered B ≃ 4 Mpc^{-1} is not itself resolution-limited at the lowest k that affect the multipoles shown in Fig. 2.
minor comments (5)
- Abstract and Introduction: the phrase “model-independent” is too strong; the search is template-free within a fixed operator set and a fixed baseline cosmology. Soften to “template-free within a defined hypothesis class”.
- Fig. 1 caption and axis label: “2 relative to the power-law PPS” appears truncated; restore “Δχ^{2}” for clarity.
- Tables I–II versus V–VI: the main-text tables list only selected Pareto points; a one-sentence pointer that the full fronts appear in the Supplemental Material would help readers.
- Discussion: the standard-clocks interpretation (p = -1) is presented as heuristic, which is appropriate, but a brief remark on whether the recovered amplitude and phase are consistent with existing NEC-violating model predictions would strengthen the paragraph.
- Eq. (1) and subsequent expressions: units of k (Mpc^{-1}) are stated once; repeating them in the captions of Tables I–II would avoid ambiguity for readers who start from the tables.
Circularity Check
Mild self-citation for the fixed-parameter approximation; the inverse-k form itself emerges from free SR search rather than by construction.
specific steps
-
self citation load bearing
[Methods / Comparison with standard templates paragraph; citation [27]]
"It has been shown that the variation of the background and nuisance parameters has a negligible effect on the feature constraints [27]."
The SR search and subsequent MCMC ranking that produce the headline inverse-k selection are performed with all cosmological and nuisance parameters fixed to the power-law best-fit. The assertion that this freeze does not alter the preferred analytic form or Δχ² ranking rests only on a citation to the same authors' prior paper; without that self-citation the fixed-baseline χ² objective could in principle reshape the residual that SR attributes to f(k)≃A cos(B/k). The citation is therefore load-bearing for treating the recovered form as fully data-driven, though it does not insert the inverse-k expression by construction.
full rationale
The paper's central result is that PySR, run independently on Planck and SPA with a shared operator set and CMB χ² loss, selects inverse-k oscillations (C8 expressions in Tables I–II / Eqs. 2–3) as the leading low-complexity feature; that form is then re-parameterized as a three-parameter template and ranked against linear/log templates via MCMC. This is not definitional circularity, fitted-input-as-prediction, uniqueness smuggling, or renaming of a known result: the analytic expression is discovered by search, not inserted. The only load-bearing self-citation is the claim that freezing background and nuisance parameters at the power-law best-fit (Methods; Supplemental Material) has negligible effect on feature constraints, justified solely by the authors' prior work [27]. That is a standard methodological approximation in the literature and does not force the inverse-k functional form itself; the SR Pareto fronts and subsequent Δχ² ranking remain independent content. Score 2 reflects one minor self-citation that is not load-bearing for the claimed discovery. No other circular steps are present.
Axiom & Free-Parameter Ledger
free parameters (5)
- Inverse-k frequency B (SR C8) =
≃4.1–4.2 Mpc^{-1}
- Feature amplitude A_X (template MCMC) =
A_inv best-fit 0.0203; posterior 0.0163^{+0.0099}_{-0.0071}
- Frequency ω_X and phase φ_X =
e.g. inverse best-fit ω=83.32, φ=2.24
- Fixed ΛCDM + nuisance baseline =
best-fit power-law values (not re-listed)
- SR complexity cutoff / Pareto selection (C8) =
complexity 8 highlighted
axioms (6)
- domain assumption Primordial spectrum is P_R(k)=P_{R,0}(k)[1+f(k)] with power-law P_{R,0} and k_*=0.05 Mpc^{-1}.
- ad hoc to paper CMB χ² from CLASS with fixed cosmology is a sufficient objective for discovering the preferred analytic f(k).
- ad hoc to paper Operator set {+,-,×,/,^,exp,log,sin,cos,tanh}, no nested trig, max size 25, defines the hypothesis class.
- domain assumption CLASS sampling with 400 k-points per decade adequately resolves inverse-k phases for B≃4 over the main CMB-sensitive range.
- domain assumption A k^{1/p} oscillation can be generated by a massive field in a∝t^p (standard-clocks), so inverse-k suggests p=−1.
- standard math Standard arithmetic/trigonometric identities and likelihood comparison under flat priors.
invented entities (1)
-
Inverse-k primordial feature template f_inv∝cos(ω k_*/k + φ)
no independent evidence
read the original abstract
Oscillatory features in the primordial power spectrum, potential signatures of new physics in the early universe, are usually searched for using fixed templates. In this work, we perform a template-free search for primordial features using symbolic regression. We find that both Planck and the combined Planck+ACT+SPT-3G datasets independently select an inverse-$k$ oscillation, $\cos(B/k)$ with $B\simeq4\,\mathrm{Mpc}^{-1}$, as the leading low-complexity feature. Comparing this inverse-$k$ template with standard linear and logarithmic oscillating templates, we find that it fits the data best, showing a weak preference for a non-zero amplitude. Our results show that symbolic regression as a powerful machine learning technique can provide an interpretable, model-independent approach to cosmological discovery.
Figures
Reference graph
Works this paper leans on
-
[1]
A. H. Guth, Phys. Rev. D23, 347 (1981)
1981
-
[2]
A. D. Linde, Phys. Lett. B108, 389 (1982)
1982
-
[3]
Albrecht and P
A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett.48, 4 0.00 0.01 0.02 0.03 0.04 0.05 0.06 AX Inverse Linear Log FIG. 3. Marginalized posteriors of the oscillation amplitude AX for the linear, logarithmic and inverse-ktemplates on the SPA dataset. 1220 (1982)
1982
-
[4]
A. A. Starobinsky, Phys. Lett. B91, 99 (1980)
1980
-
[5]
N. Aghanimet al.(Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]
Pith/arXiv arXiv 2020
-
[6]
A. A. Starobinsky, JETP Lett.55, 489 (1992)
1992
-
[7]
J. A. Adams, B. Cresswell, and R. Easther, Phys. Rev. D64, 123514 (2001), arXiv:astro-ph/0102236
Pith/arXiv arXiv 2001
-
[8]
X. Wang, B. Feng, M. Li, X.-L. Chen, and X. Zhang, Int. J. Mod. Phys. D14, 1347 (2005), arXiv:astro- ph/0209242
arXiv 2005
-
[9]
Gong, JCAP07, 015 (2005), arXiv:astro- ph/0504383
J.-O. Gong, JCAP07, 015 (2005), arXiv:astro- ph/0504383
arXiv 2005
-
[10]
A. Slosaret al., Bull. Am. Astron. Soc.51, 98 (2019), arXiv:1903.09883 [astro-ph.CO]
Pith/arXiv arXiv 2019
-
[11]
Ach´ ucarroet al., (2022), arXiv:2203.08128 [astro- ph.CO]
A. Ach´ ucarroet al., (2022), arXiv:2203.08128 [astro- ph.CO]
Pith/arXiv arXiv 2022
-
[12]
X. Chen and C. Ringeval, JCAP08, 014 (2012), arXiv:1205.6085 [astro-ph.CO]
Pith/arXiv arXiv 2012
-
[13]
X. Chen, A. Loeb, and Z.-Z. Xianyu, Phys. Rev. Lett. 122, 121301 (2019), arXiv:1809.02603 [astro-ph.CO]
Pith/arXiv arXiv 2019
-
[14]
R. Easther, W. H. Kinney, and H. Peiris, JCAP05, 009 (2005), arXiv:astro-ph/0412613
Pith/arXiv arXiv 2005
-
[15]
P. D. Meerburg, D. N. Spergel, and B. D. Wandelt, Phys. Rev. D89, 063536 (2014), arXiv:1308.3704 [astro- ph.CO]
Pith/arXiv arXiv 2014
-
[16]
P. D. Meerburg and D. N. Spergel, Phys. Rev. D89, 063537 (2014), arXiv:1308.3705 [astro-ph.CO]
Pith/arXiv arXiv 2014
-
[17]
H. Peiris, R. Easther, and R. Flauger, JCAP09, 018 (2013), arXiv:1303.2616 [astro-ph.CO]
Pith/arXiv arXiv 2013
-
[18]
R. Easther and R. Flauger, JCAP02, 037 (2014), arXiv:1308.3736 [astro-ph.CO]
Pith/arXiv arXiv 2014
-
[19]
M. Aich, D. K. Hazra, L. Sriramkumar, and T. Souradeep, Phys. Rev. D87, 083526 (2013), arXiv:1106.2798 [astro-ph.CO]
Pith/arXiv arXiv 2013
-
[20]
Y. Akramiet al.(Planck), Astron. Astrophys.641, A10 (2020), arXiv:1807.06211 [astro-ph.CO]
Pith/arXiv arXiv 2020
-
[21]
F. Beutler, M. Biagetti, D. Green, A. Slosar, and B. Wal- lisch, Phys. Rev. Res.1, 033209 (2019), arXiv:1906.08758 [astro-ph.CO]
Pith/arXiv arXiv 2019
-
[22]
M. Ballardini, F. Finelli, F. Marulli, L. Moscardini, and A. Veropalumbo, Phys. Rev. D107, 043532 (2023), arXiv:2202.08819 [astro-ph.CO]
Pith/arXiv arXiv 2023
-
[23]
T. Mergulh˜ ao, F. Beutler, and J. A. Peacock, JCAP08, 012 (2023), arXiv:2303.13946 [astro-ph.CO]
Pith/arXiv arXiv 2023
-
[24]
Louiset al.(Atacama Cosmology Telescope), JCAP 11, 062 (2025), arXiv:2503.14452 [astro-ph.CO]
T. Louiset al.(Atacama Cosmology Telescope), JCAP 11, 062 (2025), arXiv:2503.14452 [astro-ph.CO]
Pith/arXiv arXiv 2025
-
[25]
Calabreseet al.(Atacama Cosmology Telescope), JCAP11, 063 (2025), arXiv:2503.14454 [astro-ph.CO]
E. Calabreseet al.(Atacama Cosmology Telescope), JCAP11, 063 (2025), arXiv:2503.14454 [astro-ph.CO]
Pith/arXiv arXiv 2025
-
[26]
E. Camphuiset al.(SPT-3G), Phys. Rev. D113, 083504 (2026), arXiv:2506.20707 [astro-ph.CO]
Pith/arXiv arXiv 2026
-
[27]
Z.-Y. Peng and Y.-S. Piao, Phys. Rev. D113, 023544 (2026), arXiv:2507.17276 [astro-ph.CO]
Pith/arXiv arXiv 2026
-
[28]
S. K. Nerval, R. Hlozek, H. T. Jense, and J. R. Bond, (2026), arXiv:2606.28310 [astro-ph.CO]
Pith/arXiv arXiv 2026
-
[29]
M. Braglia, X. Chen, and D. K. Hazra, JCAP06, 005 (2021), arXiv:2103.03025 [astro-ph.CO]
Pith/arXiv arXiv 2021
-
[30]
M. Braglia, X. Chen, and D. K. Hazra, Eur. Phys. J. C 82, 498 (2022), arXiv:2106.07546 [astro-ph.CO]
Pith/arXiv arXiv 2022
-
[31]
M. Braglia, X. Chen, and D. K. Hazra, Phys. Rev. D 105, 103523 (2022), arXiv:2108.10110 [astro-ph.CO]
Pith/arXiv arXiv 2022
-
[32]
M. Braglia, X. Chen, D. K. Hazra, and L. Pinol, JCAP 03, 014 (2023), arXiv:2210.07028 [astro-ph.CO]
Pith/arXiv arXiv 2023
-
[33]
C. Petretti, M. Braglia, X. Chen, D. Kumar Hazra, and S. Paban, JCAP06, 035 (2025), arXiv:2411.03459 [astro- ph.CO]
Pith/arXiv arXiv 2025
-
[34]
Schmidt and H
M. Schmidt and H. Lipson, Science324, 81 (2009)
2009
-
[35]
Interpretable machine learning for sci- ence with pysr and symbolicregression.jl,
M. Cranmer, “Interpretable machine learning for sci- ence with pysr and symbolicregression.jl,” (2023), arXiv:2305.01582 [astro-ph.IM]
Pith/arXiv arXiv 2023
-
[36]
S.-M. Udrescu and M. Tegmark, Sci. Adv.6, eaay2631 (2020), arXiv:1905.11481 [physics.comp-ph]
Pith/arXiv arXiv 2020
-
[37]
A. Sousa-Neto, C. Bengaly, J. E. Gonzalez, and J. Alcaniz, Phys. Dark Univ.50, 102108 (2025), arXiv:2502.10506 [astro-ph.CO]
Pith/arXiv arXiv 2025
-
[38]
S. M. Koksbang and A. Heinesen, (2026), arXiv:2604.05822 [astro-ph.CO]
Pith/arXiv arXiv 2026
-
[39]
D. J. Bartlett, H. Desmond, and P. G. Ferreira, IEEE Trans. Evol. Comput.28, 964 (2024), arXiv:2211.11461 [astro-ph.CO]
Pith/arXiv arXiv 2024
-
[40]
D. Blas, J. Lesgourgues, and T. Tram, JCAP07, 034 (2011), arXiv:1104.2933 [astro-ph.CO]
Pith/arXiv arXiv 2011
-
[41]
See Supplemental Material at [URL will be inserted by publisher] for details of the symbolic-regression setup, the resolution estimate, and the full Pareto-front results
-
[42]
C. Stahl, D. Werth, and V. Poulin, Phys. Rev. D111, 123514 (2025), arXiv:2502.02571 [astro-ph.CO]
Pith/arXiv arXiv 2025
-
[43]
R. Calderon, T. Simon, A. Shafieloo, and D. K. Hazra, JCAP01, 057 (2026), arXiv:2504.06183 [astro-ph.CO]
Pith/arXiv arXiv 2026
-
[44]
Note that the MCMC results depend on the choice of priors, in particular onω X, which is not equivalent across the three templates; the comparison should therefore be regarded as heuristic
-
[45]
J. Torrado and A. Lewis, JCAP05, 057 (2021), arXiv:2005.05290 [astro-ph.IM]
Pith/arXiv arXiv 2021
-
[46]
Y.-S. Piao and Y.-Z. Zhang, Phys. Rev. D70, 063513 (2004), arXiv:astro-ph/0401231
Pith/arXiv arXiv 2004
-
[47]
M. Baldi, F. Finelli, and S. Matarrese, Phys. Rev. D72, 083504 (2005), arXiv:astro-ph/0505552
Pith/arXiv arXiv 2005
-
[48]
Y. Cai and Y.-S. Piao, Phys. Rev. D103, 083521 (2021), 5 40 20 0 20 40 DTT [ K2] Inverse Linear Log Planck ACT SPT-3G 2 0 2 DTE [ K2] 500 1000 1500 2000 2500 0.5 0.0 0.5 DEE [ K2] FIG. 4. CMB residuals of the best-fit linear, logarithmic and inverse-ktemplates relative to the no-feature SPA baseline, compared with the Planck (cropped), ACT and SPT-3G data...
Pith/arXiv arXiv 2021
-
[49]
G. Ye, M. Zhu, and Y. Cai, JHEP02, 008 (2024), arXiv:2312.10685 [gr-qc]
Pith/arXiv arXiv 2024
-
[50]
Y. Cai, M. Zhu, and Y.-S. Piao, Phys. Rev. Lett.133, 021001 (2024), arXiv:2305.10933 [gr-qc]
Pith/arXiv arXiv 2024
-
[51]
S. Pan, Y. Cai, and Y.-S. Piao, Eur. Phys. J. C84, 976 (2024), arXiv:2404.12655 [astro-ph.CO]
Pith/arXiv arXiv 2024
-
[52]
Adeet al.(Simons Observatory), JCAP02, 056 (2019), arXiv:1808.07445 [astro-ph.CO]
P. Adeet al.(Simons Observatory), JCAP02, 056 (2019), arXiv:1808.07445 [astro-ph.CO]
Pith/arXiv arXiv 2019
-
[53]
Inverse-k Primordial Oscillations from a Symbolic Regression Search
K. Abazajianet al., (2019), arXiv:1907.04473 [astro- ph.IM]. 6 Supplemental Material for “Inverse-k Primordial Oscillations from a Symbolic Regression Search” Details of the SR search We write the primordial power spectrum (PPS) as PR(k) =P R,0(k) [1 +f(k)] (6) andf(k) is the primordial feature we search for. We usekin Mpc −1 as the input variable ofPySR....
Pith/arXiv arXiv 2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.