{"id":"45cb31e7-fa66-4de8-9882-cc94c2e20009","arxiv_id":"2608.13347","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Relaxing the simplifying assumptions in the excursion-set-theory derivation of the most probable outer density profile makes it match the theory's own numerics but disagree more strongly with N-body simulations.","lead":"A new test shows that the analytic profile used for the outskirts of galaxy clusters matches simulations only when derived with simplifying assumptions. A more careful derivation from excursion set theory fits the theory's own numbers but strays further from real simulations, so the profile should be treated as a fitting formula.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Divergence claim rests on uncontrolled comparison to N-body profiles from a prior publication; growth of divergence could be artifact of mismatched mass/radius/window conventions.","rationale":"The reader's verdict ACCEPT is reasonable given the paper's internal consistency: the analytic mode from Eq. (6) matches the numerically-resolved DD (Figs. 2, 5), and the paper explicitly frames the window-function explanation as a hypothesis for future work. However, the headline claim — that the US profile should only be used as an effective description because relaxing assumptions worsens agreement with simulations — rests on an external comparison to previously published N-body profiles. That comparison is not controlled for mass definition, radius mapping, or window function, all of which are known to differ between the excursion-set calculation and simulation-based profile extraction. If the β↔radius mapping or mass definition is mismatched, the apparent divergence and its growth could be an artifact rather than a physical consequence of relaxing assumptions. The full spherical-collapse stiffness (ρ̂ < 10) adds a secondary uncertainty, since the 'relaxed' profiles in the mode region still rely on the approximate conversion Eq. (7), but the primary load-bearing issue is the uncontrolled N-body comparison. A direct re-analysis of the simulations with consistent conventions would settle this concern. Therefore the appropriate verdict is CONDITIONAL: accept the theoretical and numerical contributions, but condition the central conclusion on the proposed controlled comparison.","tokens_in":13315,"tokens_out":3986,"duration_ms":40925,"concrete_test":"Re-extract density profiles from the Korkidis & Pavlidou (2024) N-body simulations (or equivalent) with controlled definitions: fix halo mass definition (e.g. M200c), convert β to radius using the same enclosed-mass relation used for the simulated profiles, and compare to modes from Eq. (6) and Eq. (17). If the divergence and its growth with relaxation persist, the conclusion stands; if the curves realign or the ordering changes, the reported divergence is an artifact of the comparison conventions.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central negative result — that relaxing assumptions increases divergence from N-body profiles — depends on the comparability of the analytic modes (Eqs. 6 and 17) with the simulated profiles imported from Korkidis & Pavlidou (2024, Fig. 5). The paper does not re-analyze the simulations; it reuses published profiles. The analytic curves in Fig. 6 are evaluated as functions of β, the clustering parameter, while the simulated profiles are functions of radius; the β↔radius mapping is not derived or tested in this work. The excursion-set calculation assumes a sharp-k-space window and a specific mass-variance S(m); N-body profiles are constructed with configuration-space top-hat filters and specific halo mass definitions. If the β-to-radius mapping or the mass definition does not match the analytic conventions, the horizontal offset between curves is uncontrolled, and the 'growing divergence' as assumptions are relaxed could reflect this offset rather than the physics. In addition, the full spherical-collapse conversion (Eq. 8) could not be evaluated in the mode region ρ̂ < 10 (Sec. 2.3 and Fig. 4), so the 'relaxed' profiles in Fig. 6 still use the approximate conversion Eq. (7); the claim that relaxing assumption (3) leaves the divergence unaffected is therefore based on extrapolation of the derivative matching at ρ̂ > 14.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper re-derives the most probable outer density profile from excursion set theory's double distribution (DD). It transforms the DD from linear overdensity to nonlinear matter density, derives a cubic polynomial for the mode of the transformed distribution (Eq. 6), and evaluates the mode numerically across mass and clustering parameter β. The paper then compares the resulting analytic profiles with N-body profiles taken from Korkidis & Pavlidou (2024, Fig. 5), finding that the analytic mode agrees with the numerically realized DD but diverges more strongly from the simulated profiles as simplifying assumptions are relaxed. The authors attribute this mismatch to the difference between the sharp k-space window function used in excursion set theory and the configuration-space top-hat filter used in simulations, and conclude that the universal-scaling (US) profile should be used only as a few-parameter effective description fitted to simulations.","tokens_in":13569,"tokens_out":5301,"duration_ms":52884,"significance":"If the negative result is established, it is a useful clarification for the field: the US profile is not a first-principles prediction of excursion set theory, and its success as a fit to simulations must be understood through the window-function/Markovian assumption. The paper's internal derivation is self-contained, the analytic mode is validated against a direct numerical evaluation of the DD (Figs. 2 and 5), and the analysis code is released, which are concrete strengths. However, the central comparison to N-body profiles is not controlled: the analytic curves are functions of β while the simulated profiles are functions of radius, and the β↔radius mapping is neither derived nor tested. The conclusion that the divergence grows as assumptions are relaxed is therefore not yet fully supported, and in addition the full spherical-collapse conversion could not be evaluated in the mode region, so assumption (3) was not actually relaxed in the main comparison.","major_comments":[{"comment":"The central claim that the corrected profile diverges from N-body profiles, and that this divergence grows as assumptions are relaxed, rests on an uncontrolled comparison: the analytic modes in Fig. 6 are computed as functions of the clustering parameter β (via Eqs. 6 and 17), while the simulated profiles are imported from Korkidis & Pavlidou (2024, Fig. 5) as functions of radius. The mapping from β to radius is not derived, stated, or tested anywhere in this paper. Without such a mapping, the horizontal offset between the analytic curves and the simulated profiles is uncontrolled, so the apparent growth of the divergence could be an artifact of mismatched mass, radius, or window conventions rather than a physical effect. Please provide the β↔radius mapping used to place the modes in Fig. 6, or redo the comparison with simulations analyzed under the same mass and window definitions as the analytic calculation.","section":"Sec. 3, Fig. 6"},{"comment":"The full spherical-collapse conversion, Eq. (8), could be evaluated only for ρ̂ ≳ 14 because of the stiffness of Eq. (11), while all modes discussed in Sec. 3 lie at ρ̂ < 10. Consequently, the 'relaxed' profiles shown in Fig. 6 still use the approximate conversion, Eq. (7), in exactly the region where the modes occur. The paper therefore does not actually relax assumption (3) in the comparison that supports the successive-relaxation claim. This limitation is acknowledged in the text, but it should be made explicit that the 'divergence grows as we relax successive assumptions' statement applies only to assumptions (1), (2), (4), and (5), and not to assumption (3), unless the full conversion is evaluated down to the mode region.","section":"Sec. 2.3 and Fig. 4"}],"minor_comments":[{"comment":"Equation (4) is typeset ambiguously as ρ̂_us = (1 − β − γ)^(−δc); it should be ρ̂_us = (1 − β^{−γ})^{−δc} (or an equivalent explicit form), since the power-law index γ enters as an exponent of β.","section":"Eq. (4)"},{"comment":"Equation (6) is a cubic polynomial in X, so it has up to three real roots; the manuscript does not state which root corresponds to the physical mode nor how that branch is selected in Figs. 5 and 6. This should be specified to make the mode calculation reproducible.","section":"Eq. (6) and Figs. 5–6"},{"comment":"The sentence in Sec. 3 that says the results 'rule out assumptions (2)–(4)' is inaccurate, because assumption (3), the spherical-collapse approximation, is not evaluated in the mode region; only assumptions (2), (4), and (5) are quantitatively addressed. This wording should be corrected.","section":"Sec. 3"},{"comment":"The Fig. 5 caption appears to have a truncated axis label; the horizontal axis should be explicitly labeled as β (and the vertical axis as the mode ρ̂) so that the figure is interpretable without reference to the main text.","section":"Fig. 5 caption"},{"comment":"The reference to Tanoglidis et al. (2016) is listed as an arXiv identifier only; if the paper has appeared in a journal, the full bibliographic information should be provided.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's internal consistency and code release are genuine strengths, and the derivation of Eq. (6) appears to be a useful, self-contained contribution. My recommendation is driven entirely by the two load-bearing issues in the simulation comparison: the unspecified β-to-radius mapping in Fig. 6, and the incomplete relaxation of the spherical-collapse assumption in the mode region. If those are addressed, the paper would be a solid contribution to A&A; as it stands, the central negative claim is not yet established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: a counterintuitive negative result, honestly handled. This paper removes simplifying assumptions from a widely-used analytic density profile and finds the corrected version fits N-body simulations worse, not better. That is a genuinely useful thing to know.\n\nWhat's new and good: the cubic mode equation for the properly transformed double distribution (Eq. 6 / Appendix A) is original, and it is checked against numerical evaluation of the DD in two independent ways (Figs. 2 and 5). The code is released. The paper also makes a clean methodological point: the first-principles derivation is internally consistent, yet its agreement with simulations is worse than the older, more approximate US profile. So the US profile should be treated as a few-parameter fitting formula, not a prediction from excursion set theory. That conclusion is clearly stated and actionable.\n\nSoft spots, in order of severity.\n\nFirst, the central comparison to N-body profiles is imported from Korkidis & Pavlidou (2024) rather than re-derived here. The analytic curves are functions of β, the clustering parameter; the simulated profiles are extracted under top-hat mass definitions. The mapping between β and radius, and the matching of mass conventions, is not derived or tested in this paper. If the alignment is off, the growing divergence as assumptions are relaxed could be partly a horizontal offset between curves rather than purely the physics. The authors argue that the window-function mismatch is the physics, but the comparison as presented cannot fully distinguish that from a technical artifact. A referee should ask for a direct re-analysis of the simulations with matched conventions, or at least an explicit statement of the mapping.\n\nSecond, the full spherical collapse conversion (Eq. 8) could only be evaluated for ρ̂ ≳ 14 because of stiffness; the modes sit at ρ̂ < 10. The paper discloses this clearly. The claim that relaxing assumption (3) leaves the divergence unchanged is therefore an extrapolation, albeit a reasonable one given the derivative matching at large ρ̂. Not fatal, but the conclusions read slightly stronger than the numerics support.\n\nThe Markovian-vs-correlated-window explanation is appropriately framed as a hypothesis, and the discussion of non-Markovian alternatives is well grounded.\n\nWho gains: anyone using the US profile as a fitting formula, and anyone working on excursion set theory. The paper deserves a serious referee. I would send it out, with the expectation that the comparison be made rigorous.","headline":"A solid and honest paper with a counterintuitive negative result, but the central comparison to N-body profiles borrows conventions from prior work and needs to be made rigorous before the quantitative divergence claim is taken at face value.","tokens_in":14080,"tokens_out":9729,"would_cite":true,"duration_ms":79527,"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":"Removing the simplifying assumptions from the excursion-set derivation makes the predicted profile agree with the underlying distribution but diverge more from N-body simulations, so the universal-scaling profile is an effective fit…","keywords":["outer density profile","excursion set theory","double distribution","turnaround radius","universal scaling profile","window function","Markovian random walk","spherical collapse"],"falsifier":"Re-extract the most probable outer density profile from the same N-body simulations using a filter whose Fourier increments are independent (or a smooth k-space filter matched to the analytic variance convention), holding mass definitions fixed, and check whether the corrected analytic modes now match; if the divergence persists under matched windows, the window-function explanation would be refuted.","tokens_in":13097,"feed_emoji":"🌌","tokens_out":14741,"duration_ms":126016,"temperature":0.7,"pith_summary":"The paper asks whether the 'universal-scaling' outer density profile of galaxy clusters—a formula derived from excursion set theory under a stack of simplifying assumptions—is still obtained when those assumptions are removed one by one. The authors re-derive the most probable profile from the double distribution, the joint probability of a collapsed object's mass and the density contrast of its surroundings, this time transforming properly to nonlinear density before taking the mode, keeping the mass-dependent term, and using a realistic matter variance. They find that each correction makes the analytic mode agree better with the numerically evaluated distribution but worse with the N-body simulation profiles that motivated the original formula. They attribute the growing mismatch to the window function: excursion set theory's sharp filter in Fourier space (k-space) produces independent, Markovian increments, while simulated profiles are built with a top-hat filter in real space whose increments are correlated. The paper concludes that the universal-scaling profile should be treated as a few-parameter effective description fitted to simulations, not as a first-principles prediction of excursion set theory.","feed_headline":"A more faithful excursion-set profile fits simulations worse","feed_subtitle":"The universal-scaling profile survives only as a fitted empirical description, not a first-principles prediction.","key_machinery":"The load-bearing object is the double distribution (DD), the joint probability for a collapsed object of mass $m$ to be embedded in a region of overdensity at the larger scale $\\beta m$, and in particular the mode of slices of that distribution, which is interpreted as the most probable outer density profile. The paper's technical core is a cubic polynomial, $2A'\\tilde\\delta_c X^3 + (2A'\\eta - 2B'\\tilde\\delta_c)X^2 - (2B'\\eta + 2\\tilde\\delta_c + \\tilde\\delta_c^2)X - \\eta(1+\\tilde\\delta_c)=0$ with $X=\\hat\\rho^{-1/\\tilde\\delta_c}$, that gives the mode after the distribution has been properly transformed from linear to nonlinear overdensity; it replaces the earlier shortcut of differentiating first and converting afterwards. Around this core, the machinery includes the choice of spherical-collapse conversion, the Bardeen transfer function for the mass variance $S(m)$, and the sharp k-space window whose independent increments make the density trajectory Markovian—the property that the paper identifies as the source of the mismatch with top-hat-filtered simulation profiles.","core_discovery":"On the paper's own terms, the central discovery is negative: the earlier simplifying assumptions were not innocent shortcuts. When the double distribution is transformed to nonlinear density before the mode is taken, when the dropped term that carries mass dependence is retained, and when the power-law variance is replaced by a Bardeen transfer-function variance, the resulting most probable profile continues to match the numerical mode of the double distribution—all internal consistency checks pass—but it departs from the N-body profiles, and the departure grows as each assumption is relaxed. The paper reads this as evidence that the discrepancy lies inside excursion set theory itself, in the Markovian random walk induced by the sharp k-space window function, which is not the window used to measure simulated profiles. It therefore concludes that the universal-scaling profile is not a first-principles prediction of the outer density profile; it is a functional form that works well as an effective description once its few parameters are fitted to simulated data.","pith_inferences":["Editorial extension: if the window-function diagnosis is correct, then re-extracting simulated profiles with a sharp-k-like filter that has independent Fourier increments should restore agreement with the corrected analytic modes; the paper does not run this test.","Editorial extension: the monotonic growth of the divergence suggests the simplifying assumptions and the Markovian approximation partially cancel each other's errors in the original profile, so fitting the universal-scaling profile to simulations may be absorbing both effects at once.","Editorial extension: because the full spherical-collapse conversion could not be evaluated for $\\hat\\rho \\lesssim 14$, the behavior of the correctly transformed mode in exactly the low-density regime where the main conclusions are drawn still rests on the approximate conversion.","Editorial extension: the same Markovian-versus-correlated mechanism, if real, should affect other excursion-set predictions such as conditional mass functions and clustering statistics, not only outer density profiles."],"forward_implications":["The turnaround-radius probe of dark energy remains usable, but its theoretical profile must be calibrated against N-body simulations for each mass range rather than taken as a parameter-free prediction.","A correct first-principles prediction of the outer density profile will require a non-Markovian formulation of excursion set theory with correlated density trajectories.","The universal-scaling profile's empirical success means its functional form and fitted parameters still carry information about the quasi-linear density field, even though its excursion-set derivation is not first-principles.","The correctly transformed mode of the double distribution can serve as an internal benchmark for numerical implementations of excursion set theory, independently of its agreement with simulations."],"supporting_citations":[{"why":"Derives the double distribution and supplies the approximate spherical-collapse conversion on which both the universal-scaling profile and the transformed mode rest.","marker":"Pavlidou & Fields (2005)"},{"why":"Proposes the universal-scaling profile and provides the N-body simulation profiles that the corrected modes are compared against.","marker":"Korkidis & Pavlidou (2024)"},{"why":"Identifies the turnaround radius with the first deviation between simulated and analytic outer profiles, motivating the physical use of the profile.","marker":"Korkidis & Pavlidou (2025)"},{"why":"Establishes the excursion-set random walk and the Markovian property of the sharp k-space window used in the derivation.","marker":"Bond et al. (1991)"},{"why":"Provides the transfer function used to compute the realistic mass variance.","marker":"Bardeen et al. (1986)"},{"why":"Offers the path-integral, non-Markovian formulation the paper proposes as the route to a rigorous re-derivation.","marker":"Maggiore & Riotto (2010)"},{"why":"Derives a Fokker-Planck solution with a mass-dependent collapse threshold, another candidate framework for the re-derivation.","marker":"Lapi & Danese (2020)"}],"fun_headline_variants":["Excursion-set profile loses predictive power when refined","Refined density profile diverges more from simulations","Universal profile degrades as assumptions are relaxed","Most probable profile only empirical, not derived","Better theory, worse fit to N-body simulations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The divergence verdict assumes that the published N-body profiles, extracted with a top-hat window and particular mass definitions, are directly comparable to the sharp-k analytic calculation; if the two window conventions are not commensurable, the measured divergence is distorted.","fun_headline_variants_meta":{"raw":{"variants":["Excursion-set profile loses predictive power when refined","Refined density profile diverges more from simulations","Universal profile degrades as assumptions are relaxed","Most probable profile only empirical, not derived","Better theory, worse fit to N-body simulations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000617,"raw_usage":{"total_tokens":2879,"prompt_tokens":972,"completion_tokens":1907,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":1838}},"tokens_in":588,"tokens_out":1907,"duration_ms":12945,"temperature":1.0,"reasoning_tokens":1838,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:12:40.913605+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-extract the most probable outer density profile from the same N-body simulations using a filter whose Fourier increments are independent (or a smooth k-space filter matched to the analytic variance convention), holding mass definitions fixed, and check whether the corrected analytic modes now match; if the divergence persists under matched windows, the window-function explanation would be refuted.","supporting_citations":[{"cited_title":"& Fields, B","cited_arxiv_id":null,"evidence_quote":"Derives the double distribution and supplies the approximate spherical-collapse conversion on which both the universal-scaling profile and the transformed mode rest."},{"cited_title":"& Pavlidou, V","cited_arxiv_id":null,"evidence_quote":"Proposes the universal-scaling profile and provides the N-body simulation profiles that the corrected modes are compared against."},{"cited_title":"& Pavlidou, V","cited_arxiv_id":null,"evidence_quote":"Identifies the turnaround radius with the first deviation between simulated and analytic outer profiles, motivating the physical use of the profile."},{"cited_title":"M., Bond, J","cited_arxiv_id":null,"evidence_quote":"Provides the transfer function used to compute the realistic mass variance."},{"cited_title":"& Riotto, A","cited_arxiv_id":null,"evidence_quote":"Offers the path-integral, non-Markovian formulation the paper proposes as the route to a rigorous re-derivation."},{"cited_title":"& Danese, L","cited_arxiv_id":null,"evidence_quote":"Derives a Fokker-Planck solution with a mass-dependent collapse threshold, another candidate framework for the re-derivation."}],"review_version":1}