{"id":"ed5c280c-515b-4b5f-b9e2-97c69dfd0cfa","arxiv_id":"2607.13971","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Mean and variance of the primordial spin factor are fit to a universal function of the initial density-potential cross-correlation coefficient, with fitted parameters claimed constant across cosmologies and smoothing scales.","lead":"A heuristic fitting formula connects the distribution of a 'primordial spin factor' to the cross-correlation coefficient between initial density and potential fields. The authors claim it is universal across cosmologies, but their wCDM test is degenerate and the relation is a fit rather than a derivation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed universality across cosmologies is untested: the wCDM runs differ from ΛCDM only in dark energy, which is dynamically negligible at z=99, so all four simulations share the same early-universe Gaussian initial conditions.","rationale":"The reader's verdict is CONDITIONAL with high confidence, and my analysis supports that conditionality. The most load-bearing weakness is not precisely the reader's stated 'only q' assumption: for isotropic Gaussian fields, the normalized one-point covariance of the density and potential Hessians is indeed fixed up to the single cross-correlation q, so τ̄ and Sτ may plausibly depend only on q. The stronger, unambiguous problem is that the four cosmologies in Tables 1–2 are not cosmologically distinct at z=99. Dark energy is dynamically negligible at that redshift for all four models, so the simulations differ only in late-time background evolution and share the same primordial power spectrum. Consequently, the paper's central empirical claim of universality across cosmologies is not supported by the presented data. A concrete test with genuinely different primordial spectra would settle whether the fitted parameters are truly universal or merely characteristic of the ΛCDM power spectrum shape. The reader's rationale does mention the wCDM degeneracy, so there is partial agreement, but the formal weakest_assumption is somewhat different and arguably less likely to be the true point of failure. The correct verdict remains CONDITIONAL: the empirical relation is promising but the universality claim requires reframing or additional tests.","tokens_in":11271,"tokens_out":15225,"duration_ms":156124,"concrete_test":"Run the same Hessian pipeline at z=99 on Gaussian initial-condition simulations with the same σ8, Ωm, Ωb, and h but substantially different primordial power spectra — e.g., ns=0.8 and ns=1.2, or a running spectral index — using the same nine (Rf1,Rf2) smoothing pairs. Fit Eqs. (10)–(11) and compare the six parameters to Tables 1–2. If any parameter shifts by more than the quoted statistical errors, the universal relation does not hold for arbitrary early-universe power spectra, and the claim must be restricted to the tested ΛCDM-like spectrum. Alternatively, for a fast analytic check, generate Gaussian realizations of δ and Φ with different Pδ(k) shapes, compute τ̄(q) and Sτ(q) directly from the Hessian fields, and test whether all curves collapse onto a single q-only relation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eqs. (10)–(11) describe τ̄(q) and Sτ(q) 'universally valid ... regardless of the smoothing scales for both of the cosmologies.' The empirical support is the constancy of the fitted parameters in Tables 1–2 across four simulation sets. But at z=99, the dark-energy density fraction is ≤ ~3×10^-3 for ΛCDM and even smaller for the wCDM models (for w=-1.5 it scales as a^{4.5} and is ~10^-9 of matter). The expansion history at that epoch is therefore matter-dominated for all four models, and the linear transfer functions are effectively identical given the shared Ωm, Ωb, h, ns, and σ8. The wCDM runs are not four different early universes; they are one Gaussian initial power spectrum repeated with different late-time background parameters. Thus the parameter constancy demonstrates scale-independence of the fit for a single Pδ(k) shape, not cosmology-independence. The heuristic argument in Section 2's first bullet that only q matters is not a proof, and even if the Gaussian one-point Hessian distribution is q-only (which is plausibly true by isotropy), the specific six-parameter functional form is fitted, not derived, and no test with a genuinely different primordial spectrum (different ns, running, cutoff, or non-Gaussianity) is presented. The abstract's universality claim therefore rests on an untested dependence on the shape of Pδ(k).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a heuristic analytic relation between the initial density-potential cross-correlation coefficient q and the distribution of the primordial spin factor τ. The authors assume, in Section 2, that the mean and variance of τ depend only on q, and then fit two three-parameter functional forms, Eqs. (10)-(11), to measurements from four Multiverse N-body simulations at z=99 (two ΛCDM and two wCDM backgrounds) over nine combinations of smoothing scales. They report that the fitted parameters are robustly constant across these cases and conclude that the formula is universal, suggesting that q could be reconstructed from observed galaxy size distributions via the τ distribution.","tokens_in":11679,"tokens_out":9871,"duration_ms":102941,"significance":"The idea of connecting the early-universe density-potential cross-correlation to observable galaxy spin statistics is interesting and potentially useful. The numerical analysis is based on large simulations and a systematic set of smoothing scales, and the paper is clearly written. However, the central advertised result is not established: the 'analytic expression' is a six-parameter fitting function, not a derivation, and the same measurements are used both to determine and to validate it. The wCDM rows in Tables 1-2 are essentially degenerate with the ΩΛ=0.74 ΛCDM case at z=99, so they do not test cosmology-independence. If the claims are appropriately narrowed and an independent validation is added, the empirical relation could be a useful calibration, but as it stands the 'universality' claim is overreaching.","major_comments":[{"comment":"The 'analytic expression' is a phenomenological fit with six free parameters (three for τ̄ and three for Sτ). The same simulation data are used to set these parameters and then to demonstrate agreement in Figures 3 and 6; there is no holdout, out-of-sample prediction, or independent test. Consequently the abstract's statement that 'we prove that this analytic expression is universally valid' is not supported by the evidence. The authors should either add a genuine predictive test (e.g., a cosmology, smoothing scale, or power-spectrum shape not used in the fit) or explicitly reframe the result as an empirical calibration for the tested simulations.","section":"Section 3, Eqs. (10)-(11), Tables 1-2"},{"comment":"The assumption that τ̄ and Sτ depend on q alone is load-bearing but is only argued heuristically; no test is shown that fixes q while varying the individual auto-correlations σδ² and σΦ². Moreover, q is a correlation coefficient bounded by |q|≤1 via the Cauchy-Schwarz inequality applied to Eqs. (5)-(7), yet the best fit gives qc=1.049 in Table 1 and qs,c=0.979 in Table 2. The qc value lies outside the allowed domain, and the fitted Sτ formula becomes negative near q=1, directly contradicting the stated expectation that the function should vanish as q→1. The fitting domain and positivity constraints need to be addressed, or an explicit explanation is required for why the extrapolated qc>1 is physically meaningful.","section":"Section 2, first bullet, and Eqs. (4)-(7), Tables 1-2"},{"comment":"The wCDM runs do not provide an independent test of cosmology-independence. At z=99, dark energy is dynamically negligible for all four backgrounds; with Ωde=0.74 and w=-0.5 or -1.5, the early-universe expansion history and linear transfer function are essentially the same as for the ΩΛ=0.74 ΛCDM case. The exact agreement of the fitted parameters for these three rows to the quoted precision suggests that the underlying initial density fields are effectively identical or that the wCDM variation has no dynamical effect at the analyzed epoch. The only genuine shape variation is between the ΩΛ=0.74 and 0.64 runs, which differ in Ωm and thus in the matter-radiation equality scale. To support the claimed universality, the authors need to test initial power spectra with genuinely different shapes (e.g., different ns, running, cutoff, or non-Gaussianity) and to clarify whether independent random","section":"Section 3, Tables 1-2"},{"comment":"The proposed application—constraining early-universe physics by reconstructing q from galaxy size distributions—rests entirely on the q-only ansatz and on the specific functional forms of Eqs. (10)-(11). Since neither is derived, the discussion should clearly distinguish a testable hypothesis from an established result. At minimum, the authors should either provide a direct analytical derivation of the q-only dependence from the structure of the 12×12 covariance matrix or explicitly state that this is an empirical ansatz that has so far been checked only for a limited family of power spectra and smoothing scales.","section":"Section 4, Discussion"}],"minor_comments":[{"comment":"The words 'prove' and 'universally valid' are too strong for a six-parameter fit validated on the same data it was fitted to. Consider 'propose', 'show empirically', or 'calibrate'.","section":"Abstract and throughout"},{"comment":"Please clarify the matter density parameter for each cosmology. For flat ΛCDM with ΩΛ=0.74, Ωm=0.26, while ΩΛ=0.64 gives Ωm=0.36; the text does not state this explicitly and the later parenthetical 'equivalently, Ωm=0.26' may confuse readers.","section":"Section 3, paragraph 1"},{"comment":"The y-axis label '10k' appears garbled; the figure should clearly state whether the plotted quantities are kθ and kθ², and the tick labels should be readable.","section":"Figure 3 and Figure 6"},{"comment":"Fitting qc in the 'extrapolated range of q, i.e., beyond unity' is problematic because q is a correlation coefficient; this issue should be flagged in the main text and the domain of validity of Eqs. (10)-(11) stated explicitly.","section":"Footnote 1"},{"comment":"The integral notation is garbled ('Z s', missing differentials in places). Please rewrite using conventional differential notation with explicit integration domains.","section":"Eq. (2)"},{"comment":"Two Moon & Lee 2025 entries are cited in the text as 2025a and 2025b, but the reference list does not carry the 'a'/'b' suffixes. Please add them for unambiguous citation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has a useful empirical core, but the advertised universality claim is substantially stronger than what the data support. The identical fitted parameters for the ΩΛ=0.74, w=-0.5, and w=-1.5 runs should be investigated: it strongly suggests those runs do not provide independent early-universe realizations. A revised version that removes the 'prove' language, adds an out-of-sample test with a genuinely different power-spectrum shape, and addresses the q>1 fitting issue could be publishable; in the present form the central claim is not established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — this is a classic case of an empirical fit dressed up as a derived universal law. The new element is a six-parameter fitting formula, Eqs. (10)–(11), that links the mean and variance of the primordial spin factor tau to the density-potential cross-correlation coefficient q. That relation is genuinely new and, within the tested set-up, it appears to work: the Gamma distribution fits the simulated p(tau) nicely across nine smoothing-scale combinations, and the best-fit parameters are statistically constant across those scales. The proposed route from observable galaxy sizes to q, via the reconstruction of Moon & Lee (2025b), is a clever, inexpensive diagnostic if it holds up.\n\nThe problem is the word 'universal' and the word 'prove' in the abstract. The wCDM simulations are not tests of different early universes. At z=99, dark energy is dynamically irrelevant: for w=-1.5 it scales as a^{4.5} and is ~10^{-9} of the matter density. With the same Omega_m, Omega_b, h, n_s, and sigma8, the linear transfer functions are effectively identical. So all four simulation sets are one Gaussian P(k) shape repeated with different late-time parameters. The constancy of Tables 1 and 2 therefore demonstrates scale-independence for a single power spectrum, not cosmology-independence. The claim that q alone determines tau is an assumption stated in Section 2, not a proof; the 12-dimensional integration might marginalize other correlations, but the paper does not show it.\n\nThere are also two smaller but real issues. The 'analytic expression' has six free parameters and is fitted to the same measurements it is then validated against; that is legitimate as a parameterised empirical description, but it is not a predictive test. And the fitted threshold q_c is above unity (1.049 for tau-bar; the footnote says it is extrapolated beyond the physical range). Since q is a correlation coefficient, q>1 is unphysical, so the functional form has a regime that cannot be checked.\n\nWhat would fix this: reframe the claim as an empirical relation for Gaussian initial conditions with a standard CDM-like power spectrum; test against simulations with different n_s, a running index, a cutoff scale, or local-type non-Gaussianity; and ideally release the code and fitting data. The paper deserves a serious referee—the idea is worth taking seriously and the numerical analysis is careful—but the universality claim should not survive in its current form.","headline":"Useful empirical fit, but 'universal' is not earned—the wCDM runs do not actually change the early-universe power spectrum.","tokens_in":12102,"tokens_out":2937,"would_cite":false,"duration_ms":30980,"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":"The mean and variance of the primordial spin factor are determined solely by the density-potential cross-correlation coefficient, a universality that holds across smoothing scales in both ΛCDM and wCDM cosmologies.","keywords":["primordial spin factor","density-potential cross-correlation","tidal torque theory","Gamma distribution","halo angular momentum","early universe probes","large-scale structure","cosmological degeneracies"],"falsifier":"Run a cosmological simulation with a markedly different matter power spectrum shape (for instance, a changed spectral index or a bump in the spectrum) but with q matched to a ΛCDM case at the same smoothing scales; if the measured τ̄ and Sτ deviate from the universal curve, the q-only dependence is false. Alternatively, measure τ̄ and Sτ at fixed q but varying σ8 in the same cosmology — any significant shift would contradict the claim.","tokens_in":11168,"feed_emoji":"🌀","tokens_out":4150,"duration_ms":39974,"temperature":0.7,"pith_summary":"The paper tries to establish that the full distribution of the primordial spin factor — a measure of misalignment between the density and potential Hessian principal axes that sets the angular momentum of dark matter halos — is controlled by a single number, the cross-correlation coefficient q between the initial density and potential fields. Using N-body simulations across four cosmologies and nine smoothing-scale combinations, the authors show that the mean and variance of τ follow a universal analytic function of q alone, with fitted parameters independent of dark energy density, equation of state, and smoothing scales. A sympathetic reader would care because earlier work showed that the τ distribution can be reconstructed from observable galaxy size distributions, so this relation would let observers measure q from galaxies and, through q, probe early-universe physics without the usual cosmological parameter degeneracies.","feed_headline":"Galaxy spins hinge on one cosmic cross-correlation number","feed_subtitle":"A universal formula ties the spin factor's mean and variance to the density-potential correlation, turning galaxy sizes into a probe of the","key_machinery":"The central object is the correlation coefficient q, defined as the density-potential cross-correlation normalized by the square roots of the density and potential auto-correlations at two smoothing scales, and the heuristic functional form that expresses τ̄ and Sτ through q via a threshold parameter and power-law index. The work it does is to compress a formally 12-dimensional Gaussian integration over density and potential Hessian components into a q-only dependence, making the distribution of τ fully reconstructable from a single number.","core_discovery":"The authors propose heuristic formulas, Eqs. (10)–(11), in which both the mean τ̄ and variance Sτ of the primordial spin factor decrease from a plateau value as an exponential-with-power-law function of q, with a threshold value of q below which the decline is mild and above which it is steep. They find that three best-fit parameters in each formula — a normalization, a threshold, and a power-law index — remain constant, within errors, across ΛCDM cosmologies with different dark-energy densities and wCDM cosmologies with different equations of state, and across nine combinations of the two smoothing scales. Together with the Gamma-distribution form of p(τ), this means the entire spin-factor","pith_inferences":["The near-identical fitted parameters across four cosmologies hint that the functional form may be derivable from the Gaussian statistics of the smoothed fields alone; if so, the constants would be exact, not heuristic, and the relation could be tested with fully analytic calculations.","A direct testable extension is to vary the smoothing filter shape (e.g., top-hat instead of Gaussian) or the spectral index: if the τ–q relation shifts while q is held fixed, the universality claim is bounded to Gaussian-filtered Gaussian fields.","One can also test the relation at higher redshift or in simulations with non-Gaussian initial conditions; a deviation would signal either non-Gaussianity or a breakdown of the q-only assumption, either of which is cosmologically interesting."],"forward_implications":["Given the earlier result that the τ distribution can be reconstructed from the observed galaxy size distribution, the universal τ–q relation allows the density–potential cross-correlation coefficient q to be inferred from the same observables.","Because q is independent of the amplitude of the initial density fluctuations, this diagnostic is expected to suffer less from the standard cosmological parameter degeneracies than other large-scale-structure probes.","The universality over dark-energy density and equation of state means the relation is robust for a range of viable cosmologies, not just the standard one.","If q can be constrained on galactic scales, it opens a new probe of early-universe physics such as scale-dependent primordial non-Gaussianity, a running spectral index, or early dark energy."],"fun_headline_variants":["Universal law ties galaxy spin to one cross-correlation","Galaxy spin distribution boils down to a single number","Three fixed constants link spin to cosmic correlation","Spin factor's mean and variance follow one universal curve","Galaxy sizes may reveal early universe from spin correlation"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing assumption is that the mean and variance of the spin factor depend only on the cross-correlation coefficient q, meaning the 12-dimensional Gaussian integration marginalizes away every other property of the density and potential fields; if the shape or amplitude of their power spectra matters independently of q, the universal formula fails.","fun_headline_variants_meta":{"raw":{"variants":["Universal law ties galaxy spin to one cross-correlation","Galaxy spin distribution boils down to a single number","Three fixed constants link spin to cosmic correlation","Spin factor's mean and variance follow one universal curve","Galaxy sizes may reveal early universe from spin correlation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000451,"raw_usage":{"total_tokens":2134,"prompt_tokens":799,"completion_tokens":1335,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":1261}},"tokens_in":543,"tokens_out":1335,"duration_ms":10447,"temperature":1.0,"reasoning_tokens":1261,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T03:08:43.524366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a cosmological simulation with a markedly different matter power spectrum shape (for instance, a changed spectral index or a bump in the spectrum) but with q matched to a ΛCDM case at the same smoothing scales; if the measured τ̄ and Sτ deviate from the universal curve, the q-only dependence is false. Alternatively, measure τ̄ and Sτ at fixed q but varying σ8 in the same cosmology — any significant shift would contradict the claim.","supporting_citations":[],"review_version":1}