{"id":"381aa9ee-13c7-4f9f-b1f1-362a6e15462b","arxiv_id":"2507.00410","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"For dynamical dark matter halos, the orbiting density profile at fixed mass is set by one scale, the halo radius, whose scatter shrinks from 16% to 11% when formation time is included.","lead":"A study of simulated dark matter halos shows that each halo's orbiting density profile at fixed mass can be described by a single scale, the halo radius, with roughly 16% scatter. Accounting for halo formation time reduces that scatter to about 11%, simplifying how halo structure relates to mass and assembly history.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The single-degree-of-freedom claim rests on an α∞–rh correlation fit to the same two-parameter fits that define both variables; without covariance or mock validation, the tight relation in Eq. (11) may be a fitting degeneracy rather than physical.","rationale":"The reader's weakest_assumption identifies exactly the same concern: the α∞–rh relation is measured on the same fits and imposed in the one-parameter refits, with no covariance or independent validation. I find this to be the most load-bearing issue because the paper's headline claim ('single dynamical variable') follows directly from interpreting a fitted correlation as physical. The paper does report that the one-parameter fits are 'very similar' in quality, but this is qualitative and cannot distinguish a true one-parameter family from a two-parameter model with a degeneracy. The concrete test above would settle it. I therefore agree with the reader's CONDITIONAL verdict; no verdict change is needed.","tokens_in":10357,"tokens_out":4546,"duration_ms":49468,"concrete_test":"Use the Hessian of the cost function (Eq. 5) at each halo's best-fit (ln rh, α∞) to obtain the 2D covariance matrix, and measure the direction of its principal eigenvector. If the median angle is close to dα∞/d ln R = 1 (Eq. 11), the correlation is consistent with a fitting degeneracy. Complement this with a mock test: generate profiles using Eq. (1) with independent Gaussian scatter in ln rh and α∞ (e.g., σ_ln rh = 0.16, σ_α∞ = 0.10) and the same particle noise, then run the identical fitting pipeline. If the recovered Δα∞–ln R relation has slope near 1 with small scatter, the pipeline alone produces the apparent one-parameter family, undercutting the paper's central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section III.B, the paper demonstrates a tight correlation between the two fitted parameters Δα∞ and ln R (Eqs. 10–11), and interprets this as evidence that individual halo profiles have one effective degree of freedom. The line is fit to the very same individual two-parameter fits that define both quantities (Eqs. 8–9), and the scatter in Fig. 2 is not decomposed into measurement error versus intrinsic scatter. The adopted profile model (Eq. 1) with ε fixed at 0.037 means the slope transition occurs at x ~ ε, which is outside the fitted radial range (r > 6×softening ≈ 0.09 h⁻¹ Mpc); only the regime where α(x) is near α∞ is constrained. If the likelihood for individual halos has a ridge along Δα∞ ≈ ln R + const (e.g., because Poisson noise in the outer profile shifts rh and α∞ together), the fitted points will populate that ridge even if the true profiles have two independent degrees of freedom. The subsequent one-parameter refits impose Eq. (11) and therefore cannot validate it. No independent simulation or mock-catalog test is provided. The central claim thus hinges on an unverified assumption that the fitted correlation is intrinsic.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the dark matter density profiles of dynamical halos in a single N-body simulation. It fits individual halo profiles with the two-parameter model of Salazar et al., identifies a tight correlation between the fitted slope offset Δα∞ and the scaled halo radius at fixed mass, adopts this correlation as a one-parameter model, measures the scatter in halo radius, and studies its dependence on formation time and accretion rate. The central claims are that individual profiles at fixed mass have one effective degree of freedom, the halo radius, and that the profile shape can be used to read off both mass and formation time.","tokens_in":10736,"tokens_out":8827,"duration_ms":101980,"significance":"If the single-degree-of-freedom claim is established, this would be a useful empirical simplification of halo structure: the orbiting density profile of dynamical halos would be a one-parameter family at fixed mass, with a modest additional dependence on formation time. The paper is transparent about its methods, makes the analysis code publicly available, and includes a useful comparison with the independent Diemer parameterization. The main limitation is that the central correlation is calibrated and then imposed on the same data, so the paper currently falls short of demonstrating that the correlation is physical rather than a fitting degeneracy.","major_comments":[{"comment":"The central claim that individual profiles have a single effective degree of freedom rests entirely on the tight Δα∞–ln R relation, but this relation is measured from the same two-parameter fits that define both variables, and the subsequent one-parameter refits in Section III.C impose Eq. (11) and therefore cannot validate it. The fixed value ε=0.037 places the slope transition at x∼0.037, below the fitted radial range r>0.09 h−1 Mpc, so the data constrain only the approach of α(x) to α∞ over the fitted interval; this is exactly the regime in which a degeneracy between rh and α∞ could masquerade as a physical correlation. I ask the authors to demonstrate, for example by fitting mock halos with independent rh and α∞, by cross-validating Eq. (11) on an independent subsample or simulation, or by reporting the individual-fit covariance and showing that its noise direction is not aligned with Eq. (11), that the correlation is not an artifact of the fitting procedure.","section":"III.B (Eq. 11), Fig. 2"},{"comment":"The intrinsic scatter σln rh|Morb = 0.16 is derived by fitting Var(ln rh/rh,st|N) = σ0^2 + k/N to the one-parameter refits, but the fit is not shown, no uncertainty on σ0 is quoted, and the variance model is not tested. Since these refits already assume Eq. (11), any systematic error in the assumed relation enters the quoted scatter. Please report the fitted k, the uncertainty on σ0, and validate the particle-noise decomposition with mocks or by checking stability across radial binning and the choice of δ in Eq. (5).","section":"III.C (Eq. 14)"},{"comment":"The formation-time relation is calibrated on the same high-mass sample used to define aRF via median(a60|Morb), and the result is summarized by the residual scatter σ = 0.11 without an uncertainty or residual diagnostics. The paper should show that the residuals are consistent with zero mean and constant variance as a function of mass and aRF, and quote an uncertainty on the 0.11 value. The abstract's wording that only a small fraction of the scatter is due to formation time should also be quantified relative to the variance: (0.16)^2 versus (0.11)^2 leaves roughly half the variance unexplained, which is not obviously a small fraction.","section":"IV (Eqs. 17-18)"},{"comment":"The internal consistency check between Eq. (20) and Eq. (21) yields a 0.029 difference between the right-hand side (−0.094) and the directly measured value (−0.064). The attribution of this offset to non-Gaussian tails is plausible, but it should be verified, for instance by showing that the tail correction has the required sign and size. As it stands, the model has a small unmodeled offset in a relation that is used to support the formation-time dependence of the halo radius.","section":"IV (Eq. 22)"}],"minor_comments":[{"comment":"The second bullet says that 'at fixed mass, the slope of the profile (α∞) is tightly correlated with halo mass'; this appears to be a typo, as the text and Fig. 2 show a correlation with halo radius (or scaled radius R), not with halo mass.","section":"VI, bullet list"},{"comment":"The phrase 'characterized the orbiting orbiting profiles of spherical halos' contains a duplicated word and should be corrected.","section":"VI, first comparison paragraph"},{"comment":"The manuscript mixes the spellings 'halos' and 'haloes'; please choose one and use it consistently.","section":"Throughout"},{"comment":"The best-fit value of k in Eq. (14) is not reported, which makes the variance model difficult to reproduce; please add the fitted value and its uncertainty.","section":"III.C (Eq. 14)"},{"comment":"The uncertainties on the fitted line come from jackknifing the box, but the number of jackknife realizations and the scatter around the line are not stated; reporting these would help readers assess the tightness of the correlation in Fig. 2.","section":"III.B (Eq. 10)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid empirical characterization, but the central single-parameter claim is not yet established because the defining correlation is calibrated and enforced on the same data. I would like to see a mock-based or split-sample validation before publication. The concern is testable and fixable within the scope of the manuscript, so major revision rather than rejection seems appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what you should know: this is a solid empirical paper with public code and transparent fitting methods, but the headline claim—that individual dynamical halo profiles at fixed mass live on a one-parameter family—is not established as strongly as the paper suggests. It rests on a correlation measured from the same fits that define both variables, with no independent simulation validation and no error bars on the headline scatter values.\n\nWhat's new: the observation that the two fitted profile parameters (inner slope and halo radius) are tightly correlated at fixed mass, and that the scatter in halo radius is ~16%, shrinking to ~11% when relative formation time is included. The comparison to Diemer's recent paper is honest and useful, and the code availability is a plus. The internal consistency check in Section IV is a nice touch, even if it's approximate.\n\nNow the soft spots, in proportion. The biggest is the single-degree-of-freedom claim. The authors fit each halo with two free parameters, plot the resulting alpha_inf against r_h, fit a line to that plot, then impose the line and refit everything with only r_h free. That procedure cannot validate the relation—it assumes it. They never quantify the scatter around the relation, and they don't propagate the covariance between the two fitted parameters. The stress-test concern about a fitting degeneracy (Poisson noise shifting r_h and alpha_inf together) is legitimate. It's not fatal, because Diemer's Figure 10 shows a similar correlation using a different profile model and different simulations, suggesting some degree of physical reality. But the paper doesn't lean on that as external validation; it just cites it as qualitative agreement. A referee should ask for a mock-catalog test or an independent simulation fit.\n\nSecond, the headline scatter numbers (16% and 11%) have no uncertainties. The variance model includes a particle-noise term, but the error on the intrinsic scatter is never quoted. For a paper whose main quantitative results are these two percentages, that's a real gap.\n\nThird, the formation-time analysis is restricted to halos with a60 > 0.4 and orbital masses above 10^14 h^-1 M_sun, so the 11% figure applies to a limited range. The paper acknowledges this, but it should be made clear that the low scatter claim is conditional.\n\nWho's this for: anyone modeling halo density profiles or using cluster profiles for mass calibration. It's a useful empirical contribution, but I'd treat the one-parameter compression as provisional until validated elsewhere.\n\nMy recommendation: send it to a serious referee. The question is worth asking, the execution is mostly careful, and the issues are addressable in revision. A competent referee will likely ask for cross-validation and error propagation, but that doesn't make the paper unworthy of review.","headline":"Plausible but unproven: the claim that individual dynamical halo profiles have one degree of freedom needs independent validation and covariance handling.","tokens_in":11202,"tokens_out":3599,"would_cite":false,"duration_ms":40616,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that individual dynamical halos of fixed mass have orbiting density profiles controlled by a single variable, the halo radius $r_{\\rm h}$, with halo radii scattered by about 16% at fixed mass and by about 11% after…","keywords":["dark matter halos","density profiles","dynamical halos","orbiting particles","halo radius","formation time","cosmological simulations","halo mass scatter"],"falsifier":"Fit every halo in an independent simulation with both $\\alpha_\\infty$ and $r_{\\rm h}$ free, then refit with $\\alpha_\\infty$ forced by Equation (12); if the forced fits degrade significantly, or if the free fits scatter about $\\Delta\\alpha_\\infty = 0.05 + \\ln R$ much more widely than the quoted errors, the single-degree-of-freedom claim is not supported.","tokens_in":10192,"feed_emoji":"🌌","tokens_out":16315,"duration_ms":142356,"temperature":0.7,"pith_summary":"This paper adopts the dynamical-halo definition in which a halo is the collection of orbiting particles around its self-generated potential, and asks how the orbiting density profile of such a halo depends on mass. It claims that at fixed mass the profile of an individual halo is controlled by a single dynamical variable, the halo radius $r_{\\rm h}$, because the fitted inner slope $\\alpha_\\infty$ is tied to $r_{\\rm h}$ through $\\Delta\\alpha_\\infty = 0.05 + \\ln R$. The halo radius is lognormally scattered by about 16% at fixed mass; including the relative formation time $a_{\\rm RF}$ through $r_{\\rm h,mod} = (1.91-a_{\\rm RF})\\,r_{\\rm h,st}(M_{\\rm orb})$ lowers the scatter to about 11%. If correct, a halo's mass and formation time can be read off from the shape of its orbiting profile alone.","feed_headline":"One radius sets a dark halo's density profile","feed_subtitle":"For dynamical halos at fixed mass, profile shape is set by one radius; formation time tightens scatter from 16% to 11%.","key_machinery":"The load-bearing object is the halo radius $r_{\\rm h}$, the length scale at which the exponential truncation $\\exp(-r^2/2r_{\\rm h}^2)$ turns on in the adopted profile model [9]. With the inner scaling $\\epsilon$ fixed at 0.037, the slope $\\alpha_\\infty$ is replaced by the deterministic relation $\\Delta\\alpha_\\infty = 0.05 + \\ln R$, where $R \\equiv r_{\\rm h}/r_{\\rm h,st}(M_{\\rm orb})$, turning the two-parameter fit into a one-parameter fit. Then $r_{\\rm h}$ is linked to the relative formation time $a_{\\rm RF} = a_{60}/{\\rm median}(a_{60}|M_{\\rm orb})$ by $r_{\\rm h,mod} = (1.91-a_{\\rm RF})\\,r_{\\rm h,st}(M_{\\rm orb})$, which tightens the predicted radius. The slope–radius relation is what makes the profile a single-degree-of-freedom family; the formation-time relation is what makes the radius partly predictable from accretion history.","core_discovery":"The central discovery is that the orbiting profiles of dynamical halos form an effectively one-parameter family at fixed mass. Fitting each halo with the truncated power-law model $\\rho_{\\rm orb}(x) = A(x/\\epsilon)^{-\\alpha(x)}\\exp(-x^2/2)$ yields two shape parameters, the asymptotic slope $\\alpha_\\infty$ and the halo radius $r_{\\rm h}$; the paper shows these are tightly correlated within a mass bin, $\\Delta\\alpha_\\infty = 0.05 + \\ln(r_{\\rm h}/r_{\\rm h,st}(M_{\\rm orb}))$. This correlation lets every halo's profile be described by $r_{\\rm h}$ alone, with the slope imposed by Equation (12). The halo radius at fixed orbiting mass is lognormally distributed with intrinsic scatter $\\sigma_{\\ln r_{\\rm h}|M_{\\rm orb}} = 0.16$; late-forming halos are more compact, and using the relative formation time through Equation (18) reduces the scatter to $\\sigma_{\\ln r_{\\rm h}|M_{\\rm orb},a_{\\rm RF}} = 0.11$. The paper takes these two scalings as evidence that the orbiting profile of a dynamical halo is a single-degree-of-freedom object whose extent encodes mass and formation history.","pith_inferences":["If the relation $\\Delta\\alpha_\\infty = 0.05 + \\ln R$ holds across redshifts and simulation resolutions, the orbiting profile could serve as a low-cost estimator of halo mass and formation time that bypasses full merger-tree construction.","Because the residual scatter after including formation time is still 11%, a natural next test is whether combining $a_{\\rm RF}$ with the recent accretion rate removes more of the scatter than either variable alone; the paper compares but does not combine the two.","If profile shape determines mass and formation time, then lensing or satellite-kinematics measurements of the inner slope may be able to infer halo radius, turning the single-degree-of-freedom claim into a route from observed profiles to halo assembly history.","The 16% scatter in halo radius at fixed mass implies a corresponding spread in the truncation scale of the orbiting profile, which should be propagated into stacked halo-model predictions of large-scale structure."],"forward_implications":["A halo's orbiting density profile at fixed mass is fully characterized by one number, $r_{\\rm h}$, so the profile shape carries no independent information beyond its spatial extent.","From a measured profile slope and radius, Equations (11) and (18) can be inverted to estimate both the orbiting mass and the relative formation time of the halo.","The intrinsic scatter in $r_{\\rm h}$ at fixed mass is about 16%; adding formation time reduces it to about 11%, setting the accuracy with which mass, and mass plus formation time, predict the profile extent.","Late-forming halos are more compact and have shallower inner slopes, extending the known concentration–formation-time connection to the orbiting profiles of dynamical halos.","The same trends appear with an alternative four-parameter profile parameterization, suggesting the single-degree-of-freedom behavior is not an artifact of the adopted fitting function alone."],"supporting_citations":[{"why":"It defines dynamical halos as collections of orbiting particles and supplies the algorithm used to tag particles as orbiting or infalling.","marker":"[1]"},{"why":"It supplies the fiducial profile model and the stacked-profile parameter values that the paper fits to individual halos.","marker":"[9]"},{"why":"It is the cosmological simulation from which all halo and particle data in this paper are taken.","marker":"[12]"},{"why":"It produces the initial halo catalog whose particles are then classified into orbiting and infalling sets.","marker":"[15]"},{"why":"It supplies the alternative four-parameter profile parameterization used in the comparison against the fiducial model.","marker":"[11]"},{"why":"It is the parallel analysis of orbiting profiles that the paper compares with, providing independent evidence for the radius–slope trend.","marker":"[17]"},{"why":"It establishes the link between halo formation history and internal structure that motivates the relative-formation-time analysis.","marker":"[16]"},{"why":"It provides the accretion-rate measurements used to compare accretion rate with formation time as predictors of halo radius.","marker":"[23]"}],"fun_headline_variants":["One radius controls dark halo density","Halo profiles collapse to a single radius","Dark halo shape from one dynamical variable","Formation time tightens halo profile scatter","Single radius sets dark halo orbiting profile"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the tight correlation between fitted inner slope and halo radius, $\\Delta\\alpha_\\infty = 0.05 + \\ln R$, reflects a physical single degree of freedom and not a degeneracy in the two-parameter fitting procedure, because the same fits that define the relation are then refit under it.","fun_headline_variants_meta":{"raw":{"variants":["One radius controls dark halo density","Halo profiles collapse to a single radius","Dark halo shape from one dynamical variable","Formation time tightens halo profile scatter","Single radius sets dark halo orbiting profile"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1407,"prompt_tokens":983,"completion_tokens":424,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":362}},"tokens_in":599,"tokens_out":424,"duration_ms":5829,"temperature":1.0,"reasoning_tokens":362,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:18:25.666484+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit every halo in an independent simulation with both $\\alpha_\\infty$ and $r_{\\rm h}$ free, then refit with $\\alpha_\\infty$ forced by Equation (12); if the forced fits degrade significantly, or if the free fits scatter about $\\Delta\\alpha_\\infty = 0.05 + \\ln R$ much more widely than the quoted errors, the single-degree-of-freedom claim is not supported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defines dynamical halos as collections of orbiting particles and supplies the algorithm used to tag particles as orbiting or infalling."},{"cited_title":"The Density Profile of Dynamical Halos","cited_arxiv_id":"2507.00410","evidence_quote":"It supplies the fiducial profile model and the stacked-profile parameter values that the paper fits to individual halos."},{"cited_title":"A dynamics-based density profile for dark haloes. I. Algorithm and basic results","cited_arxiv_id":"2112.03921","evidence_quote":"It produces the initial halo catalog whose particles are then classified into orbiting and infalling sets."},{"cited_title":"read off","cited_arxiv_id":null,"evidence_quote":"It is the parallel analysis of orbiting profiles that the paper compares with, providing independent evidence for the radius–slope trend."}],"review_version":1}