{"id":"ec3b008b-4e9c-41aa-a1cf-8cf65e81edbe","arxiv_id":"2411.14614","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A two-parameter model including soliton self-gravity shows that simulated fuzzy dark matter cores sit between the self-gravity and host-halo regimes, and that core-halo scatter depends on both halo concentration and core-intrinsic features.","lead":"Fuzzy dark matter halos are modeled as a wave-like soliton core inside a smooth outer halo, and the authors infer two numbers for each simulated halo that control the core's size. The paper argues that the scatter in the core-halo relation comes partly from the halo's concentration and partly from the core's own internal history, and demonstrates a way to estimate the dark matter particle mass from halo profiles.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Core-halo scatter claim may be circular: reconstructed B is a deterministic function of the same f=rt/rc used to bin the data","rationale":"The reader's verdict is CONDITIONAL and already flags the reconstruction assumptions and the tautological relation between f and B. My stress-test sharpens this: the load-bearing evidence for the headline claim is the correlation of core-halo scatter with B, but B is derived from the same f=rt/rc ratio used in the binning. This is not a fatal mathematical error in the model, and the FDM mass reconstruction (Section 5.3) and consistency checks in Appendix C are credible and independently checkable. However, without a demonstration that B is not merely a proxy for f, the paper's central claim about intrinsic soliton features remains unsupported. A concrete partial-correlation or mock test would settle this. The verdict should remain CONDITIONAL, with the condition being that the claimed B-driven scatter survives control for f.","tokens_in":20593,"tokens_out":3171,"duration_ms":35677,"concrete_test":"Using the same cosmological simulation data, compute the residuals of log10 rc after removing the median rc–Mh relation (in narrow Mh bins). Then compute the partial Spearman correlation between those residuals and B = log10(beta/beta_crit), controlling for log10(f). If the partial correlation becomes insignificant or changes sign, the apparent B-dependence is a proxy for the input f, undercutting the central claim. A complementary mock test: generate halos from the model with known alpha, beta and with f scattered independently of beta, run the reconstruction, and check whether the recovered B spuriously correlates with rc residuals.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that scatter in the core-halo relation is driven by intrinsic soliton features (beta) rests on the observed correlation between rc-Mh residuals and B = log10(beta/beta_crit) shown in Figures 10 and 13. But B is not an independent observable: for each halo, beta is obtained from Eq. (5.1), beta = alpha (1+c(rt/rc)^2)^8 / [(rt/rs)(1+rt/rs)^2], using exactly the measured ratio f = rt/rc (plus the reconstructed alpha and rs). Thus B is a monotonic function of f for fixed alpha and rs. Since f is itself an input that encodes the same transition-radius information used to define the core/halo decomposition, the claim that the rc–Mh scatter 'correlates with B' is largely equivalent to saying it correlates with the input f. The paper interprets this as evidence for local soliton dynamics, but it may simply reflect the fitting procedure or the assumed NFW outer profile. No quantitative test is provided to show that B carries information independent of f (e.g., a partial correlation or a null mock test). Without such a test, the central claim that intrinsic soliton features, rather than the concentration–mass scatter alone, explain the core-halo scatter is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models a fuzzy dark matter halo as a spherically symmetric Schrödinger-Poisson ground-state soliton embedded in an external NFW potential, characterized by two dimensionless parameters α (halo potential strength) and β (soliton central density). The authors solve the system numerically, derive the limiting behaviors in which either soliton self-gravity or the host halo dominates, define a critical βcrit, and obtain universal upper and lower bounds on core radius and core mass. They then reconstruct α, β, and the concentration parameter cvir for simulated halos using rc, rt, and Mh, concluding that simulated cores lie in the intermediate regime β/βcrit ~ O(1)-O(100). On this basis they argue that the scatter in the core-halo relation cannot be explained solely by the concentration-mass scatter but is also influenced by intrinsic soliton features, and they demonstrate a reconstruction of the FDM mass accurate to about 10%.","tokens_in":20864,"tokens_out":7406,"duration_ms":68709,"significance":"If the central scatter claim were established, the paper would provide a useful interpretive framework for FDM simulations and a practical model for observational core-halo analyses. The numerical implementation appears internally consistent, the two limiting cases are checked against the numerical solutions, the universal bounds in Sec. 3.3 are a useful byproduct, and the FDM mass reconstruction is an explicit, falsifiable application. However, the main scatter conclusion currently rests on a definitional correlation between B and the input ratio f=rt/rc, and the paper does not provide a quantitative test showing that B carries information independent of f. The core-halo scatter claim therefore needs substantial additional support before it can be accepted.","major_comments":[{"comment":"The central claim that the core-halo scatter is significantly influenced by intrinsic soliton features rests on the correlation between rc-Mh residuals and B = log10(beta/beta_crit). But B is not an independent observable: for each halo, beta is obtained from Eq. (5.1) using exactly the measured f = rt/rc, with alpha and rs fixed by the reconstruction, and beta_crit is a function of alpha. In the limit-B regime where xc ~ alpha^{-1/3}, the ratio beta/beta_crit reduces to a function of f alone, so the strong f-B correlation shown in Fig. 10 and the color trends in Fig. 13 are substantially built in by construction. The authors should provide a partial-correlation analysis of rc-Mh residuals with B after controlling for f, or a null mock test with synthetic scatter in f and cvir, to demonstrate that beta carries information beyond the input ratio f. Without such a test, the abstract's conclusion is not established.","section":"Sec. 5.1-5.2, Eq. (5.1), Figs. 10 and 13"},{"comment":"The reconstruction solves a single equation r_th_c(cvir) = rc for the concentration parameter, but the paper does not establish that this solution is unique or that the inferred parameters are stable to the assumed NFW form and to the fitting formula (A.1). If r_th_c(cvir) is non-monotonic or has multiple crossings in the sampled halo range, the reconstructed alpha, beta, and B depend on which branch is selected. The authors should report the number of solutions found and the sensitivity of the reconstructed parameters to the measurement uncertainties in rc and rt, for example through a bootstrap or Monte Carlo propagation of the input data.","section":"Sec. 5.1, Eq. (5.2)"},{"comment":"The FDM mass reconstruction treats Mc as an additional independent input, but the text states that in the simulation data 'Mc is related to rc' and that rt was determined using Eq. (3.8), which is strictly valid only in limit A. If Mc was derived from rc in this way, Eq. (5.6) is not an independent consistency condition and the ~10% agreement in Fig. 14 is partly by construction. The authors should clarify whether Mc is an independent simulation measurement or is derived from rc; in the latter case, they should demonstrate that the mass reconstruction works using only (Mh, rt, rc).","section":"Sec. 5.3 and Appendix C"}],"minor_comments":[{"comment":"The text reads 'pawer-law behavior'; this should be 'power-law behavior'.","section":"Sec. 1"},{"comment":"After Eq. (3.8), 'parsection' appears to be a typo for 'parsec'.","section":"Sec. 3.1.1"},{"comment":"The heading 'Maximam core radius and minimam core mass' contains two typos; it should read 'Maximum core radius and minimum core mass'.","section":"Sec. 3.3"},{"comment":"The text says 'Newton-Rapthon method'; this should be 'Newton-Raphson method'.","section":"Sec. 4.1"},{"comment":"After Eq. (4.10), 'reosonable approximation' should be 'reasonable approximation'.","section":"Sec. 4.2"},{"comment":"The labels for limit A and limit B appear to be swapped: Eq. (3.12) and Eq. (3.14) belong to limit B, while Eq. (3.5) and Eq. (3.7) belong to limit A.","section":"Figure 3 caption"},{"comment":"The dashed line in panel (c) should presumably correspond to rlim_c,max rather than rlim_c,min, to be consistent with Eq. (3.35).","section":"Figure 5 caption"},{"comment":"The sentence 'independent of neither β nor cvir' should be 'independent of both β and cvir' or equivalently 'dependent on neither β nor cvir'.","section":"Sec. 6"},{"comment":"The text says 'as a bybroduct'; this should be 'as a byproduct'.","section":"Sec. 5.1"},{"comment":"The text says 'used for the reconstructuction'; this should be 'used for the reconstruction'.","section":"Sec. 5.3"}],"recommendation":"major_revision","confidential_remarks":"The paper's numerical model and mass reconstruction are solid and publishable in principle, but the headline claim about the origin of core-halo scatter needs a nontrivial statistical test to separate B from f. The current evidence is largely a consistency relation. I would recommend major revision rather than rejection, because the missing test is well defined and could be supplied within the scope of the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the genuinely new content is the numerical solution of the Schrödinger-Poisson system with soliton self-gravity plus an external NFW potential, organized around two dimensionless parameters (alpha, beta). The limiting cases reduce to known results—limit A is the Schive soliton, limit B is Taruya-Saga—but the interpolation between them is new, and the parameter reconstruction from simulation data is a useful tool. The punchline that survives scrutiny is the mass check: using Mh, rc, rt, and Mc, they recover the fiducial m_psi ≈ 8e-23 eV to about 10%, and Appendix C shows reconstructed core masses track simulation values to roughly 30%. That gives the model real predictive content.\n\nSecond, the broader claim—that scatter in the core-halo relation is significantly influenced by intrinsic features of the soliton core rather than only concentration scatter—does not survive a careful reading. The issue is Eq. (5.1). beta is reconstructed from exactly the measured ratio f = rt/rc, together with alpha and rs. The correlation between f and B shown in Figure 10 is therefore largely built in by construction; it cannot independently support the claim that B, rather than f itself, drives the rc-Mh scatter. Likewise, A is a deterministic function of cvir and Mh, so the A-cvir correlation in Figure 9 is also partly bookkeeping. A partial correlation at fixed f, or a mock test with a null model where the outer profile is pure NFW, would be needed to show B carries information independent of the input profile. No such test is given. So the central abstract claim is not established as presented.\n\nThe equilibrium assumption is a further soft spot: real simulated cores are not stationary and spherical, and the paper's own model assumes they are. If cores are breathing or perturbed, reconstructed beta is a fit parameter, not a direct measurement of local dynamics. This is acknowledged implicitly by the setup, but the impact on the conclusion is not quantified.\n\nWhat the paper does well: clean asymptotic analysis, honest appendix checks, and a credible numerical pipeline. The qualitative result that soliton self-gravity matters in the intermediate regime between the two limits holds up.\n\nWho should read it: people working on FDM core-halo relations and observers estimating m_psi from dwarf density profiles. It deserves a serious referee. I would recommend accepting with major revision: add the mock or partial-correlation analysis, reframe the scatter claim, and state clearly what remains model-dependent.","headline":"Solid two-parameter core-halo model and a useful FDM mass reconstruction, but the central scatter claim is not established because the key parameter is reconstructed from the same f=rt/rc it is correlated against.","tokens_in":21352,"tokens_out":3788,"would_cite":true,"duration_ms":40365,"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 scatter in fuzzy dark matter's core–halo relation comes from the soliton core's own self-gravity and local dynamics, not just from halo concentration scatter.","keywords":["fuzzy dark matter","soliton core","core-halo relation","Schrödinger-Poisson equation","NFW halo","ultralight dark matter","soliton self-gravity","concentration-mass relation"],"falsifier":"Track a single simulated soliton core over a dynamical time and reconstruct $\\beta$ at many snapshots: if $\\beta$ fluctuates by orders of magnitude or jumps at merger events while the profile still fits the ground-state template, then the reconstructed scatter reflects non-equilibrium dynamics rather than an intrinsic, quasi-steady core property. Alternatively, measure the density profile's deviation from the fitting form in eq. (3.21); a systematic mismatch at fixed $(\\alpha,\\beta)$ would show the reconstruction procedure is biased.","tokens_in":20398,"feed_emoji":"🌌","tokens_out":8907,"duration_ms":68635,"temperature":0.7,"pith_summary":"This paper sets out a minimal model of a fuzzy dark matter halo in which the central soliton core is a ground state of the Schrödinger–Poisson system living inside an external Navarro–Frenk–White halo potential. The model is controlled by two dimensionless numbers: $\\alpha$, which measures the host halo's gravitational pull on the core, and $\\beta$, which measures the soliton's self-gravity through its central density. Applying the model to the soliton cores seen in cosmological and merger simulations, the authors reconstruct these parameters for each halo and find that simulated halos sit between the two extreme regimes, with $\\beta/\\beta_{\\rm crit}\\sim O(1)$ to $O(10^2)$, so both self-gravity and the host halo matter. Their central claim is that the scatter in the core–halo relation cannot be blamed only on the scatter in the halo's concentration–mass relation; intrinsic features of the soliton core, potentially from local dynamics at the halo center, contribute substantially. They also show that the fuzzy dark matter particle mass can be recovered from halo profile data to about 10% accuracy, which gives observational applications.","feed_headline":"Soliton self-gravity scatters fuzzy dark matter cores","feed_subtitle":"A two-parameter model shows simulated FDM halos need both soliton and halo gravity, and recovers the particle mass to ~10%.","key_machinery":"The load-bearing object is a stationary, spherically symmetric ground-state solution of the Schrödinger–Poisson equations for a soliton sitting in an external NFW potential (eqs. 2.17–2.18). The solution is labeled by two dimensionless parameters: $\\alpha$ (host halo strength, eq. 2.16) and $\\beta$ (central wavefunction amplitude squared, eq. 2.19). The transition between the two limiting behaviors occurs at $\\beta_{\\rm crit}(\\alpha)\\simeq 1.92\\,\\alpha^{4/3}$; in one limit the soliton's self-gravity dominates, in the other the NFW potential dominates. The reconstruction uses the measured transition radius to write $\\beta$ as a function of the concentration parameter, then solves self-consistently for the value of $c_{\\rm vir}$ whose predicted core radius matches the simulated one, using a fitting formula for $x_c(\\alpha,\\beta)$ from appendix A. The same machinery yields universal upper and lower bounds on core radius and core mass for a given halo mass and FDM particle mass.","core_discovery":"On the paper's own terms, the discovery is a corrected picture of the core–halo relation in fuzzy dark matter. Previous work approximated the core as a bound state governed mostly by the host halo or mostly by its own gravity, and attributed the scatter of the relation to the scatter in halo concentration. This paper solves the full stationary Schrödinger–Poisson system with both potentials and reconstructs $(\\alpha, \\beta, c_{\\rm vir})$ from the measured core radius $r_c$, transition radius $r_t$, and halo mass $M_h$. The reconstructed halos occupy the intermediate regime $\\beta/\\beta_{\\rm crit}\\sim 1$–$10^2$, and the reconstructed $\\beta$ scatters widely even at fixed halo mass, with no strong correlation with concentration. The authors therefore claim that the diversity of the core–halo relation reflects not only the halo's assembly history, encoded in $c_{\\rm vir}$, but also the soliton's own local state, and they show that the scatter is correlated with both the concentration parameter and the ratio $f=r_t/r_c$.","pith_inferences":["If the reconstructed $\\beta$ scatter is really tracing local soliton dynamics, then tracking individual simulated cores across time should show $\\beta$ wandering on a dynamical timescale; that is a direct simulation test of the paper's interpretation.","Applying the same reconstruction pipeline to snapshots before and after a major merger would quantify how much of the $\\beta$ scatter is merger-driven rather than quasi-steady intrinsic variation.","Replacing the NFW external potential with a baryonic density profile—an extension the paper mentions—would turn the model into a tool for predicting how baryonic feedback changes the core–halo relation and its scatter.","The universal bounds could be turned into a null test: an observed core larger than $r_{c,\\rm max}$ for an assumed particle mass would rule out that mass."],"forward_implications":["Most simulated fuzzy dark matter halos lie in the regime $\\beta/\\beta_{\\rm crit}\\sim 1$–$10^2$, so any realistic core–halo model must include both soliton self-gravity and the host halo potential.","The scatter in the core–halo relation splits into two channels: the concentration parameter $c_{\\rm vir}$ correlates with $\\alpha$, while the ratio $f=r_t/r_c$ correlates with $\\beta/\\beta_{\\rm crit}$.","Reconstructed concentration parameters fall below the CDM concentration–mass relation, roughly between two existing FDM predictions, consistent with suppressed small-scale structure.","Using the core mass as an extra input, the model recovers the simulation's FDM particle mass ($8\\times10^{-23}$ eV) to within about 10%.","For a given FDM mass and halo mass, the core radius has an upper bound and the core mass a lower bound that do not depend on $\\beta$ or $c_{\\rm vir}$, giving observational criteria."],"supporting_citations":[{"why":"Establishes the original core–halo scaling and the soliton density profile that this model reproduces in the self-gravity-dominated limit.","marker":"[14]"},{"why":"Supplies the simulation dataset (core radius, transition radius, halo mass) and the definition of the transition radius used for parameter reconstruction.","marker":"[20]"},{"why":"Supplies the cosmological fuzzy dark matter simulation whose halos are reconstructed in section 5.","marker":"[25]"},{"why":"Provides the analytical limit-B solution (NFW-dominated core) and the concentration-dependent core–halo relation that this paper extends by adding soliton self-gravity.","marker":"[31]"},{"why":"Previous quantitative claim that concentration–mass scatter explains core–halo scatter; the paper's central result is a challenge to it.","marker":"[30]"},{"why":"The CDM concentration–mass relation used as a comparison baseline for the reconstructed concentration parameters.","marker":"[32]"},{"why":"Supplies the numerical shooting method for the Schrödinger–Poisson ground state used in section 4.","marker":"[38]"},{"why":"One of the fuzzy dark matter concentration–mass models used to compare the reconstructed concentration values.","marker":"[39]"}],"fun_headline_variants":["Soliton self-gravity dictates fuzzy dark matter core-halo scatter","Two potentials explain fuzzy dark matter core-halo scatter","Soliton self-gravity is key to fuzzy dark matter core scatter","Fuzzy dark matter core-halo scatter traced to soliton self-gravity","Fuzzy dark matter mass reconstruction via soliton-halo model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reconstruction assumes that each simulated soliton is a stationary, spherically symmetric ground state of the Schrödinger–Poisson system with an external NFW potential, and that the transition radius marks the crossing where the fitted soliton and NFW profiles have equal density; if real cores are oscillating or out of equilibrium, the inferred $\\beta$ and its scatter could be artifacts of the model.","fun_headline_variants_meta":{"raw":{"variants":["Soliton self-gravity dictates fuzzy dark matter core-halo scatter","Two potentials explain fuzzy dark matter core-halo scatter","Soliton self-gravity is key to fuzzy dark matter core scatter","Fuzzy dark matter core-halo scatter traced to soliton self-gravity","Fuzzy dark matter mass reconstruction via soliton-halo model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00107,"raw_usage":{"total_tokens":4537,"prompt_tokens":1054,"completion_tokens":3483,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":3392}},"tokens_in":670,"tokens_out":3483,"duration_ms":20985,"temperature":1.0,"reasoning_tokens":3392,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:05:56.080154+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track a single simulated soliton core over a dynamical time and reconstruct $\\beta$ at many snapshots: if $\\beta$ fluctuates by orders of magnitude or jumps at merger events while the profile still fits the ground-state template, then the reconstructed scatter reflects non-equilibrium dynamics rather than an intrinsic, quasi-steady core property. Alternatively, measure the density profile's deviation from the fitting form in eq. (3.21); a systematic mismatch at fixed $(\\alpha,\\beta)$ would show the reconstruction procedure is biased.","supporting_citations":[{"cited_title":"Analytical approach to core-halo structure of fuzzy dark matter","cited_arxiv_id":"2208.06562","evidence_quote":"Provides the analytical limit-B solution (NFW-dominated core) and the concentration-dependent core–halo relation that this paper extends by adding soliton self-gravity."}],"review_version":1}