{"id":"49ff0533-b04f-45cf-bc00-0f062c121ef3","arxiv_id":"2411.09733","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A new hydrodynamic code for fuzzy dark matter shows that its quantum pressure creates dark-matter cores in all simulated dwarf halos, whereas baryons create cores only at certain masses, and stellar properties stay close to cold dark matter.","lead":"This paper presents a new simulation code that combines fuzzy dark matter quantum pressure with full hydrodynamics, and applies it to zoom-in simulations of dwarf galaxies. The main result is that quantum-pressure cores and baryon-feedback cores appear at different halo masses and redshifts, while stellar observables remain almost indistinguishable from cold dark matter predictions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Untested SPH quantum-potential discretization (Eqs. 13–15) is load-bearing: the claimed FDM cores and CDM-like stellar observables could be numerical artifacts unless the quantum-pressure estimate converges at the simulated smoothing lengths.","rationale":"The reader's weakest_assumption is the load-bearing point: the SPH discretization of the quantum potential is the foundation for the core-formation and stellar-observable claims. I agree that the absence of convergence or wave-mechanics validation makes the central claim conditional. My emphasis is slightly narrower: the most concrete risk is not primarily the Madelung single-velocity approximation but the lack of evidence that the SPH Laplacian, evaluated at the actual smoothing lengths in these runs, converges to the true quantum pressure where it matters. The small sample size and the run-count inconsistency are real but secondary; they weaken generalization, not the mechanism. Since the reader already assigned CONDITIONAL with medium risk and my concern matches that assessment, I recommend no verdict change.","tokens_in":15240,"tokens_out":6693,"duration_ms":77201,"concrete_test":"Re-run one dwarf DMO system (e.g. M_2, m22=2) with (i) 8× higher mass resolution and (ii) the same resolution but a higher-order/MLS estimate of ∇²ρ in Eq. (14), keeping all other physics fixed. If the z=0 logarithmic slope α(0.4–1 kpc) and the core radius change by more than ~0.1 in slope or 20% in radius relative to the published run, the FDM core claim is resolution- and discretization-dependent. As an independent cross-check, evolve the same initial conditions with a Schrödinger–Poisson grid code at matched resolution and compare the density profile inside 2 kpc; agreement within ~10% would close the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Every new physical result follows from the SPH estimate of the FDM quantum pressure in §3.1. Equations (13)–(15) replace ∇²√ρ/√ρ with a smoothed-particle Laplacian of the density, with kernel width h_i set by neighbour counts. The paper offers no convergence test, no comparison with a direct Schrödinger–Poisson solver, and no quantitative demonstration that h_i resolves the de Broglie wavelength for m22=2 and m22=8 at the radii where cores form. In low-density regions h_i can exceed the de Broglie scale, while in the central region the discrete Laplacian of a noisy particle density is biased and oscillatory; either effect can suppress or exaggerate the quantum pressure responsible for the claimed cores. Because the same QP term also shapes DM velocity dispersion and the 'strongest effect wins' combination of baryonic and FDM core formation, the central conclusion that FDM cores appear in all dwarf haloes while stellar observables stay CDM-like is not yet decoupled from the discretization. The extra massive-halo run in Fig. 3, which shows no core, also indicates that 'all haloes' requires the dwarf-mass qualifier, but that is secondary to the numerical issue. As a pilot code paper, the missing validation limits the result to an unverified claim rather than established physics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents fuzzy-gasoline, a modified version of the gasoline2 smooth-particle-hydrodynamics code that adds the Fuzzy Dark Matter (FDM) quantum potential to cosmological hydrodynamical simulations. It applies the code to three NIHAO dwarf galaxies of different masses, runs each as CDM, FDM with m22=8, and FDM with m22=2, in both hydrodynamic and dark-matter-only versions, and adds one massive halo as a control. The main scientific claims are that baryons and FDM both create cored dark-matter profiles through distinct mechanisms, that FDM cores appear in all simulated haloes especially at low mass and high redshift, and that stellar observables such as star formation histories and stellar velocity dispersions remain close to CDM predictions. The paper is explicitly a pilot study and emphasizes the code-development aspect.","tokens_in":15522,"tokens_out":5515,"duration_ms":53969,"significance":"If the numerical method is valid, this would be an important step: it is the first cosmological hydrodynamical simulation suite that combines FDM quantum-pressure effects with a full baryonic model including cooling, star formation, and supernova feedback. The design of comparing hydrodynamical and dark-matter-only runs for the same initial conditions is clean, and the massive-halo control run provides a useful falsification check against a naive 'FDM always cores' interpretation. However, the central results currently rely on an unvalidated SPH discretization of the quantum potential and on only three dwarf haloes, with no convergence tests, no error bars, and no comparison to a Schrödinger-Poisson reference. The significance is therefore conditional on additional validation and on tempering the statistical claims.","major_comments":[{"comment":"The SPH estimate of the quantum potential is the only FDM microphysics in the code, and every new physical result follows from it. The paper provides no convergence test in mass resolution or in the neighbor number N_N(i), no comparison with a direct Schrödinger-Poisson solver, and no demonstration that the smoothing length h_i resolves the de Broglie wavelength at the radii where cores are claimed. In low-density regions h_i can exceed the de Broglie scale, while in the central region the discrete Laplacian of a noisy particle density may bias or oscillate, either suppressing or exaggerating quantum pressure. Because the same quantum-pressure term shapes the DM density profiles, the DM velocity dispersions, and the 'strongest effect wins' combination with baryonic core formation, the central claims are not yet decoupled from this discretization. Please add (i) convergence runs with different mass resolutions and N_N(i), (ii) a validation of at least one halo against a Schrödinger-Poisson solver or a semi-analytic soliton solution, and (iii) a quantitative comparison of h_i with the local de Broglie wavelength in the core region for m22=2 and m22=8.","section":"§3.1, Eqs. (13)–(15)"},{"comment":"The abstract and Section 3.2 state that the work uses 'more than 30 zoom-in simulations', but the paper actually describes three unique dwarf haloes (L, M, S) plus one massive control halo, each run in several model variants; the tables list 22 run entries in total. More importantly, all conclusions about 'all haloes' and about stellar observables being 'remarkably similar' to CDM rest on a single halo per mass bin, with no halo-to-halo scatter or error bars. This is a load-bearing limitation for the statistical wording of the central claims. Either the sample should be expanded or the conclusions and abstract should be reworded to state explicitly that these are pilot results based on one halo per mass.","section":"§3.2 and Tables 1, A1–A5"},{"comment":"The abstract claims that 'FDM-induced cores emerge in all haloes', but the massive control system in Fig. 3 shows no noticeable core induced by either FDM model or by baryons. The text later qualifies this by noting the massive system does not feature relevant deviations, but the abstract and conclusions do not carry that qualifier. The claim should be restricted to the simulated dwarf mass range, or the wording should otherwise be reconciled with Fig. 3.","section":"Abstract and Fig. 3"},{"comment":"In the lowest-mass system S, the m22=2 run forms zero stars and the text attributes this partly to the stellar component lying 'below resolution'. This makes the statement that FDM generally delays and suppresses star formation not fully decoupled from resolution effects at the low-mass end. Similarly, the 'remarkable similarity' of stellar observables to CDM is asserted for systems where star formation is near the resolution threshold. The paper should either show that the S_2 result is converged (for example, with a higher-resolution rerun) or explicitly state in the conclusions that this part of the result is tentative.","section":"§4.2, S system"}],"minor_comments":[{"comment":"There is a typo in the kernel argument: W(|r_j - r_j|, h_i) should presumably be W(|r_j - r_i|, h_i).","section":"Eq. (11)"},{"comment":"The phrase 'more than 30 zoom-in simulations' is inconsistent with the three galaxies listed in Table 1 and the run counting in Tables A1–A5; please clarify whether 'simulations' includes subhaloes or other variants, or correct the number.","section":"§3.2"},{"comment":"The statement that positive values of the slope alpha are 'dynamically unstable values due to noise' and are therefore omitted/shaded needs a reference or quantitative justification; as plotted, the apparent constant-core behavior of the FDM model is partly a consequence of this truncation.","section":"Fig. 4, right panel"},{"comment":"The introduction and methods mention black hole physics as part of the baryonic model, but the run tables list no black hole particles or black hole properties; please specify whether black holes form in these dwarf runs or are absent because of the halo masses.","section":"§3.2"},{"comment":"The phrase 'able to correctly simulate evolving astrophysical systems' is too strong for the approximate Madelung-SPH treatment used here; 'model' or 'approximately simulate' would be more accurate given the limitations discussed above.","section":"§5"},{"comment":"The data availability statement says data will be shared 'on reasonable request' but provides no repository link or code release; providing at least the analysis scripts and a documented version of fuzzy-gasoline would strengthen reproducibility.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is a promising pilot study with a clean hydro-vs-DMO design and a useful massive-halo control, and I see no circularity: no quantity is fit to data and the FDM physics is included as model input rather than imposed on the output. The main risk is that the SPH quantum-potential discretization has not been validated, and the sample size does not support the statistical wording used in the abstract. If the authors add a convergence/validation analysis and temper the claims to the actual sample, the paper could be suitable for publication; as it stands, the central claims are not yet established. The mismatch between 'more than 30 simulations' and the three unique haloes should also be corrected before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Fuzzy Gasoline is a genuinely new step for FDM simulations: it puts the SPH quantum-pressure routine from ax-gadget into GASOLINE2 and adds the full NIHAO baryonic physics, so you can follow galaxy formation with FDM through to z=0. On top of that, the paper makes a physical claim that I don't think is in the earlier FDM literature: baryon-driven and FDM-driven cores form at different redshifts and masses, and when both are present the final profile tracks the stronger effect rather than the sum. The massive-halo control run is a nice touch—it shows the effect turns off where it should.\n\nThe main thing to know is that the central result is not yet pinned down numerically. The quantum potential is computed by SPH estimates of the second derivative of the density (Eqs. 13–15), and there is no convergence test, no comparison to a known soliton solution or a Schrödinger–Poisson run, and no check that the particle smoothing length resolves the de Broglie wavelength in the cores. In low-density outskirts h_i can be large, and in the center the discrete Laplacian of a noisy density is biased. Either way, the FDM core that appears in 'all haloes' could be partly a smoothing artifact. That matters because the same QP term also shapes the velocity dispersion and the 'strongest effect wins' behavior.\n\nOther weaknesses are proportional. It's only one halo per mass, with no error bars—fine for a pilot but not enough for strong conclusions about the population. The abstract says 'thirty zoom-in simulations' but the text describes three halos; count the hydro and DMO runs and you get about twenty, not thirty. That looks like a typo, not a red flag. Code and data are only 'available on reasonable request,' which is weak by modern standards but not disqualifying for a pilot.\n\nThe authors are aware this is a pilot study and say so. If the SPH discretization is validated, the astrophysical conclusions are plausible and would be useful to both simulators and dwarf-galaxy observers. I'd send it to a referee, but the referee should insist on a convergence test and a check against a Schrödinger–Poisson soliton before the core claim is accepted.","headline":"New FDM+baryon simulation tool, but the core-forming result rests on an unvalidated SPH quantum-pressure term; deserves review with mandatory convergence tests.","tokens_in":16045,"tokens_out":3002,"would_cite":true,"duration_ms":29610,"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":"Fuzzy dark matter's quantum pressure creates cores in every dwarf dark matter halo, earlier and at lower masses than baryonic feedback can, while leaving stellar observables almost indistinguishable from cold dark matter predictions.","keywords":["fuzzy dark matter","dwarf galaxies","galaxy formation","hydrodynamical simulations","quantum potential","core formation","NIHAO","soliton"],"falsifier":"Run the same dwarf halo with the same initial conditions in FUZZY-GASOLINE and in a direct Schrödinger–Poisson grid solver and compare the dark matter density profile at $z=0$; if the core radius differs by more than the run-to-run scatter, the SPH quantum potential is not faithful. Alternatively, search for dark matter cores in dwarf galaxies at $z \\gtrsim 2$, where baryonic feedback is too slow to make cores; finding no cores at that epoch would contradict the paper's claim that FDM cores exist at all redshifts.","tokens_in":15042,"feed_emoji":"🌌","tokens_out":5843,"duration_ms":52736,"temperature":0.7,"pith_summary":"Using a new smoothed-particle hydrodynamics code that combines the fuzzy dark matter quantum potential with full baryonic physics, the paper simulates thirty dwarf galaxies with halo masses $10^9$ to $10^{11}$ solar masses under two FDM boson masses and under CDM. It claims that cores in dark matter density profiles can form in two distinct ways: FDM's quantum pressure makes a core in every halo, most prominently in low-mass systems at high redshift, while baryonic feedback makes cores later, only in haloes near $10^{10}$ to $10^{11}$ solar masses. When both mechanisms operate, the final profile follows the strongest of the two suppressions rather than adding them. Despite large differences in the dark sector, star formation histories and stellar velocity dispersion profiles remain very close to CDM, so FDM can alter dark matter structure without betraying itself in stellar observations.","feed_headline":"Fuzzy dark matter carves cores in every dwarf halo","feed_subtitle":"First baryonic FDM simulations show star formation matches CDM even though dark matter differs.","key_machinery":"The load-bearing piece is the quantum potential $Q = \\frac{\\hbar^2}{2 m_\\chi^2} \\frac{\\nabla^2 \\sqrt{\\rho}}{\\sqrt{\\rho}}$, which appears in the Madelung-transformed Euler equation as an extra repulsive force. The code FUZZY-GASOLINE evaluates this potential with SPH density and its first and second derivatives (Eqs. 13–15), and it symmetrizes the two-particle acceleration so energy is conserved by construction. This machinery turns wave mechanics into a force law for N-body particles, letting baryonic cooling, star formation, and feedback run alongside the FDM dynamics.","core_discovery":"The central discovery is that baryons and low-mass fuzzy dark matter are each independently capable of turning a cuspy halo into a cored one, and in different regimes: FDM-induced cores appear in all haloes, especially low-mass ones at high redshift, whereas baryon-driven cores appear only in a mass range around $10^{10}$–$10^{11}\\,M_\\odot$ and at low redshift. The combined simulation shows the dark matter profile is set by whichever effect suppresses the central density more strongly, not by a sum of the two. The authors also find the velocity dispersion profiles flatten and lower in similar ways under both effects, though for different reasons. Stellar observables—star formation histories and velocity dispersion—are only mildly affected even when the dark matter structure differs substantially, which makes FDM with boson masses $m_{22}=2$ and $8$ difficult to reject using stellar data alone.","pith_inferences":["If the SPH quantum-potential scheme holds up, the same Madelung machinery could be ported to other wave-like dark matter models (e.g., repulsive self-interactions), but the single-valued velocity field assumption will need direct comparison with full Schrödinger–Poisson solvers before those extensions are trusted.","The threshold effect—$m_{22}=8$ raising the virial mass of small haloes—implies a potentially observable signature in the abundance of dwarf satellites around Milky Way–like hosts: a systematic excess relative to CDM at the low-mass end, which would be a clean test independent of stellar profiles.","Because stellar observables are degenerate between FDM and CDM, constraints on the boson mass from stellar kinematics alone may be systematically biased; combining core sizes with the requirement that stars form at all could break the degeneracy."],"forward_implications":["High-redshift cores in dwarf galaxies, if observed, would be a direct FDM signature, because baryons cannot produce cores early enough.","FDM and baryonic core formation are not additive: in systems where both operate, the final dark matter profile matches the stronger of the two single effects, guiding how future simulations and analytic models combine them.","Stellar observables like star formation histories and velocity dispersions cannot cleanly tell FDM from CDM, so distinguishing the models requires measuring the dark matter distribution itself—e.g., through inner density slopes or kinematics of dark-matter-dominated systems.","In the smallest haloes, a heavier FDM boson ($m_{22}=8$) can actually increase the halo's virial mass and boost star formation compared to CDM because mass from unresolved substructures is redistributed into the main halo, whereas $m_{22}=2$ delays or quenches star formation entirely."],"supporting_citations":[{"why":"Supplies the SPH quantum-potential scheme that fuzzy-gasoline refines and symmetrizes.","marker":"Nori & Baldi 2018"},{"why":"Provides the NIHAO halo sample and the baryonic physics model for cooling, star formation, and feedback.","marker":"Wang et al. 2015"},{"why":"Defines fuzzy dark matter and the soliton core profile that motivates the core predictions.","marker":"Hu et al. 2000"},{"why":"Suggests the smoothed-particle approach for evaluating FDM dynamics that the implementation follows.","marker":"Mocz & Succi 2015"},{"why":"Independently proposes the same SPH route for the quantum potential in cosmological N-body codes.","marker":"Marsh 2015"},{"why":"The gasoline2 code framework into which the FDM routines are integrated.","marker":"Wadsley et al. 2017"},{"why":"Documents baryon-induced core formation in CDM dwarfs and the mass range where it operates.","marker":"Tollet et al. 2016"},{"why":"Explains the feedback-driven potential fluctuations that create baryonic cores.","marker":"Pontzen & Governato 2012"},{"why":"Provides the hybrid soliton-plus-NFW profile used to characterize FDM haloes.","marker":"Schive et al. 2014"}],"fun_headline_variants":["FDM and baryons each core dwarf halos, but differently","Core formation in dwarfs: FDM and baryons take separate paths","Fuzzy dark matter and baryons all erase cusps, but stars match CDM","Dwarfs: FDM and baryons both make cores, stars give no clue"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulation's cores and stellar similarities rest on the assumption that the SPH representation of the quantum potential faithfully reproduces true wave mechanics at the resolution used; if the smoothing length or the single-fluid Madelung approximation hides interference and vorticity, the cores and the CDM-like stellar properties would be numerical artifacts.","fun_headline_variants_meta":{"raw":{"variants":["FDM and baryons each core dwarf halos, but differently","Core formation in dwarfs: FDM and baryons take separate paths","Fuzzy dark matter and baryons all erase cusps, but stars match CDM","Dwarfs: FDM and baryons both make cores, stars give no clue"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000588,"raw_usage":{"total_tokens":2774,"prompt_tokens":975,"completion_tokens":1799,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":1714}},"tokens_in":591,"tokens_out":1799,"duration_ms":11272,"temperature":1.0,"reasoning_tokens":1714,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:21:54.218476+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same dwarf halo with the same initial conditions in FUZZY-GASOLINE and in a direct Schrödinger–Poisson grid solver and compare the dark matter density profile at $z=0$; if the core radius differs by more than the run-to-run scatter, the SPH quantum potential is not faithful. Alternatively, search for dark matter cores in dwarf galaxies at $z \\gtrsim 2$, where baryonic feedback is too slow to make cores; finding no cores at that epoch would contradict the paper's claim that FDM cores exist at all redshifts.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Independently proposes the same SPH route for the quantum potential in cosmological N-body codes."}],"review_version":1}