{"id":"7811eb09-538f-4679-9d84-62d0fa849070","arxiv_id":"2509.07470","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Satellites that form most of their stars early keep dark matter cusps, while late or continuous star formers develop oscillating cores, so the dwarf diversity problem may be explained by feedback history and tides within cold dark matter.","lead":"This paper re-simulates the same spiral galaxy with five different star formation recipes to study its orbiting dwarf satellites. It finds that how early satellites form their stars, plus tidal pulls from the host, decides whether they keep dark matter cusps or flatten into cores, potentially solving the dwarf galaxy 'diversity problem' within standard dark matter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No resolution-convergence test for γ in tidally stripped subhalos; the oscillating core/cusp regime could be numerical, undermining the t90%–γ claim.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the inner slope γ of tidally stripped subhalos may not be numerically converged or physically meaningful. This is the critical link in the argument because the central claim—that a stable cusp corresponds to early star formation while late star formation produces oscillating core/cusp behavior—is supported by γ values and their time evolution. If those γ measurements are compromised, the diversity claim loses its quantitative support. I considered other potential weaknesses: the small sample (8–9 subhalos per run, with matched halos across runs), the post-hoc choice of t90% thresholds, and possible confounding between star formation timing and orbital history. However, these would weaken the statistical strength or generalizability of the claim, whereas an unverified resolution dependence could invalidate the measurement itself. The paper even flags the risk in Section 4, so this is not an invented objection. The proposed concrete test—a higher-resolution rerun or a resolution-based quality cut on the γ samples—would directly settle whether the oscillations and the t90%–γ trend are physical. Since the reader's verdict CONDITIONAL already reflects this uncertainty, my stress-test does not change the verdict; it sharpens the condition under which the claim would be accepted.","tokens_in":27065,"tokens_out":6732,"duration_ms":77087,"concrete_test":"Re-run at least one hydro simulation (e.g., Mffϵ100-MecFB, which exhibits cores and oscillations) with 8× better mass resolution (m_p ≈ 2.4e4 M_sun) in the zoom-in region, and compare the γ(t) tracks and the t90%–γ relation against the fiducial run. If the tracks and relation are unchanged, the concern is resolved. A cheaper complementary test: for every subhalo snapshot, compute the number of bound particles within r<1 kpc and the tidal radius; restrict the analysis to snapshots with N_bound(r<1 kpc) > 100 and r_tidal > 2 kpc, then check whether the oscillations and the t90% split persist. If they disappear, the claim is resolution- or boundary-driven rather than physical.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on measuring the inner dark matter slope γ within 1 kpc and tracking its time evolution for subhalos that are often heavily tidally stripped. Section 4 cites Borukhovetskaya et al. (2022) showing that insufficient resolution can artificially inflate remnant sizes and distort tidal evolution, but the paper presents no convergence test for its own γ measurements. Subhalos are selected with M_sub > 80 m_p (m_p ≈ 1.9e5 M_sun), giving only ~500–5000 particles even before stripping; after tidal mass loss, the bound population within 1 kpc can be much smaller. If resolution biases γ or its oscillations, the t90%–γ correlation and the claimed 'fluctuating regime' would be artifacts. This is the most load-bearing weakness because the unique part of the paper—the oscillating inner slopes driven by feedback and tides—rests entirely on these γ(t) tracks, whereas the population-level survival trends are less affected by this specific numerical issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes subhalo properties in the Mochima suite, a set of zoom-in simulations of the same Milky-Way-like host run with a dark-matter-only baseline and five baryonic prescriptions that vary the star formation and supernova feedback models. A refined phase-space subhalo identification method is introduced, and the authors use it to study subhalo survival, stellar content, mass spectra, and the inner dark-matter slope gamma of the resolved satellites. The central claim is an emergent correlation between star formation timing and inner structure: galaxies that formed 90% of their stellar mass more than 7 Gyr ago tend to retain stable cusps at z=0, while galaxies with late or extended star formation show oscillating inner slopes that alternate between core-like and cusp-like values. This is compared with Local Group dwarf data, and the authors argue that the observed diversity of inner profiles can arise within LCDM from the interplay of feedback history and tidal environment.","tokens_in":27327,"tokens_out":6950,"duration_ms":84648,"significance":"If correct, the t90%-gamma relation would be a valuable physical explanation for the dwarf galaxy diversity problem, connecting star formation timing and tidal interactions to core/cusp structure. The controlled same-host suite is a genuine strength: it isolates the effect of subgrid baryonic physics while holding initial conditions fixed, and the population-level trends (host concentration, potential depth, stellar binding, survival) are plausible and well aligned with external results. The paper also makes a good-faith comparison with independent Local Group observations and cites its own limitations explicitly. However, the central claim currently rests on a small, non-independent sample and on gamma measurements in heavily stripped subhalos for which no resolution convergence test is provided. The result is therefore promising and potentially important, but not yet established at the level claimed.","major_comments":[{"comment":"The central claim rests on gamma(t) tracks in tidally stripped subhalos, but no resolution-convergence test for gamma is presented. The paper itself cites Borukhovetskaya et al. (2022) [110] and warns that insufficient resolution can inflate remnant sizes, yet the only safeguard is the selection threshold M_sub > 80 m_p (Section 2.1), with m_p ~= 1.9e5 Msun. After tidal stripping, the number of particles within 1 kpc can be far smaller, and the oscillating core/cusp regime could be numerical rather than physical. A dedicated convergence test, or at least a conservative particle-number cut with a repeat of the t90%-gamma analysis, is required before the oscillating regime can be interpreted as a real dynamical phenomenon.","section":"Section 4 / Figure 6-9"},{"comment":"The split at t90% = 7 Gyr appears to be chosen after inspecting the figure and is not statistically validated. The text also uses a different criterion ('more than 10% of the stellar mass assembled within the last 5 Gyr'), which is not equivalent to t90% < 7 Gyr. Please specify the operational definition, report how many objects fall on each side of the threshold, and test the significance of the correlation (e.g., a rank correlation with a permutation test). Without this, the qualitative visual impression in Figure 7 is not sufficient support for a central claim.","section":"Section 4 / Figure 7"},{"comment":"The resolved sample is only ~8-9 subhalos per baryonic run, and many of these are the same Lagrangian halos re-simulated with different feedback models, so the effective number of independent systems is small. The shaded evolutionary tracks in Figure 7 compress 8 Gyr of gamma(t) into a 1 Gyr interval in t90% space, visually multiplying the number of points and potentially making the correlation look stronger than it is for z=0 objects. The z=0 relation should be shown separately with proper error bars, and the compressed tracks should not be treated as independent data points.","section":"Section 3.3 / Figure 7"}],"minor_comments":[{"comment":"The text says halo F 'continues forming stars until approximately 5 Myr ago'; given the lookback-time axes in Figure 6, this should presumably be '5 Gyr ago'.","section":"Section 3.3"},{"comment":"Draco and WLM appear twice in the inner panel because they are drawn from different data sources. Please distinguish the duplicates with different markers or list them explicitly in the caption.","section":"Figure 7"},{"comment":"The phase-space selection uses several empirically chosen thresholds (d6D < 10, 90 and 135 degree angles). No sensitivity test is shown for these choices. Since the main population trends are robust, this is not blocking, but a brief robustness statement would help.","section":"Appendix A"},{"comment":"The subhalo identification method is described in detail, but there is no quantitative validation against other finders or against known recovery/incompleteness curves for the adopted set of thresholds. A comparison with a standard finder in one run would strengthen the methodological claims.","section":"Section 2.1"},{"comment":"The SHMR comparison notes that simulated stellar masses are 'systematically higher' than observational constraints, but no offset or scatter is quantified. Reporting the median offset and its run-to-run range would make the statement more precise.","section":"Section 3.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its limitations and the controlled suite is well suited to the question. The main issue is that the headline t90%-gamma claim needs stronger numerical and statistical support: a convergence test for gamma in stripped subhalos and a pre-specified, non-repeated-measure analysis of the correlation. If those can be supplied in a revision, this could become a solid JCAP paper. If not, the oscillating-regime claim should be softened substantially."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this paper uses the Mochima suite — the same host halo re-simulated five ways plus a DMO baseline — to argue that the inner DM slope of dwarf satellites tracks the timing of star formation. Galaxies that formed 90% of their stellar mass early keep cusps; galaxies with extended or recent star formation show cores that oscillate back toward cusps. The t90%-gamma trend is plausible, consistent with EDGE and with Local Group data, and the extension to satellites — where tides add time-dependent fluctuations — is genuinely new. I'd take the population-level result seriously.\n\nThe paper does several things well. The subhalo identification method (potential subtraction plus phase-space selection) is a real improvement over single-snapshot spherical cuts, and the authors test it against tidal-radius estimates. The population-level findings on host concentration, baryonic potential depth, and subhalo survival are clearly presented and honestly connected to earlier work. They also cite EDGE and flag their own resolution caveat, which is the right instinct.\n\nNow the soft spots, in rough order of importance. First, the unique claim — oscillating slopes in tidally stripped subhalos — rests on gamma(t) tracks measured within 1 kpc for systems with few hundred to few thousand particles before stripping, and fewer after. The paper cites Borukhovetskaya et al. on resolution inflating remnant sizes, but does not run a convergence test for gamma in these heavily stripped systems. That is the load-bearing weakness: if the oscillations are numerical, the fluctuating regime collapses, even though the early-SFH-to-cusp trend survives. Second, the sample is small — 8-9 subhalos per run — so the scatter in Figure 7 is driven by a handful of objects per feedback model. Third, the 7 Gyr threshold looks post-hoc; it is not derived or tested against other splits. These are not fatal, but they prevent the central claim from being more than conditional.\n\nWho benefits: anyone working on dwarf galaxy structure, satellite survival, or subgrid feedback. The paper deserves a serious referee; I would send it out and ask for convergence tests on gamma, an explicit sensitivity analysis of the threshold, and a candid discussion of how many particles are actually inside 1 kpc at late times. With those, the oscillating-slope result could become solid. As it stands, the paper is a useful, honest contribution to a live debate, but the headline claim needs another round of scrutiny.","headline":"Controlled re-simulation suite shows star-formation timing separates cusps from cores in satellites, with a plausible but not yet proven oscillatory regime driven by tides and late feedback.","tokens_in":27843,"tokens_out":1467,"would_cite":true,"duration_ms":21788,"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":"Satellite dwarfs that finish forming stars early keep steep dark matter cusps, while late or bursty star formation makes their inner dark matter slope oscillate between core and cusp — a diversity produced within cold dark matter, not again","keywords":["dark matter cusps and cores","dwarf satellite galaxies","star formation history","baryonic feedback","subhalo survival","tidal stripping","cosmological zoom-in simulations","core-cusp diversity"],"falsifier":"Run the same simulations at higher mass and force resolution and check whether gamma oscillations in late-star-forming satellites persist; if they vanish, the fluctuating regime is numerical. Observationally, measure resolved star formation histories and dark matter inner slopes for a larger sample of Milky Way and Andromeda satellites: if late-forming dwarfs consistently show steep cusps, or early-forming dwarfs show cores, the t90% relation fails.","tokens_in":26985,"feed_emoji":"🌌","tokens_out":11308,"duration_ms":120594,"temperature":0.7,"pith_summary":"This paper uses a suite of cosmological zoom-in simulations of one Milky-Way-like spiral galaxy, rerun six times with different treatments of gas, star formation, and supernova feedback, to ask what controls the dark matter content of its dwarf satellite galaxies. It argues that the observed diversity of inner dark-matter profiles is not a failure of cold dark matter: it is what the standard model predicts once star formation history and tidal interactions are included. The central result is a timing relation: satellites that assembled 90% of their stellar mass more than 7 Gyr ago end up at redshift zero with a stable cusp, while satellites that formed more than 10% of their stars in the last 5 Gyr have inner slopes that oscillate between core-like and cusp-like values on gigayear timescales. Along the way the paper shows that baryon-rich subhalos resist tidal disruption better, and that host halo concentration, itself modulated by feedback, sets the survival of low-mass dark subhalos. If the timing relation holds, the 'diversity problem' becomes a problem of when and how stars formed, not a sign that dark matter needs new physics.","feed_headline":"Star-formation timing decides dwarf dark matter cusps or cores","feed_subtitle":"Satellites forming 90% of their stars early keep steep cusps; late, bursty star formation makes inner slopes oscillate.","key_machinery":"The two key quantities are gamma, the logarithmic slope of the dark matter density profile measured within 1 kpc, and t90%, the lookback time by which a galaxy assembled 90% of its stellar mass. The argument follows each resolved subhalo over the last 8 Gyr, tracking gamma alongside enclosed dark matter, stellar and gas mass and the internal and external parts of the gravitational potential, so cusp formation and erosion are tied to specific events in the star formation history and orbit. This tracking is made possible by a subhalo identification method that reconstructs a local gravitational potential for each clump and selects bound particles with an energy criterion, a directional constra","core_discovery":"The paper claims that the inner dark matter slope of a dwarf satellite, gamma = d ln rho / d ln r inside 1 kpc, is set jointly by star formation history and tidal environment rather than by halo mass or cosmological initial conditions alone. Galaxies whose star formation stalls early keep a stable gravitational potential and therefore retain a steep cusp, sometimes deepening it by adiabatic contraction. Galaxies with recent or extended star formation experience repeated supernova-driven gas outflows that temporarily flatten the cusp into a core, followed by re-contraction as stars accumulate; tidal shocks at pericentric passage add another source of fluctuation. The result is a fluctuating r","pith_inferences":["A sharper observational test would use SFH measurements for a larger sample of dwarfs beyond the Local Group, comparing t90% to inner slopes; the paper's handful of Local Group dwarfs is suggestive but small.","If the oscillating regime is physical, theoretical predictions of dark matter annihilation or gravitational lensing in dwarf satellites should use time-averaged central densities rather than the instantaneous gamma.","A decisive numerical check is to rerun the same suite at higher resolution: the early-truncation cusps should remain steep and stable, while the late-star-forming oscillations should persist; if they disappear, the fluctuating regime is at least partly a resolution artifact."],"forward_implications":["Cusps and cores are not a bimodal halo property: the same satellite can pass through both states, so single-epoch inner-slope measurements for late-forming dwarfs should be interpreted with caution.","Surveys of Local Group satellites can test the relation directly: dwarfs with early truncated star formation should consistently show steep inner slopes, while late-forming dwarfs should show much larger scatter.","Because baryon-rich subhalos are harder to destroy, the satellite stellar mass function is coupled to feedback physics; simulations with the strongest feedback produce fewer, more massive surviving satellites.","Host halo concentration, changed by baryonic feedback, becomes a population-level predictor: more concentrated hosts destroy more low-mass dark subhalos, which links the host's inner dark matter profile to its satellite demographics."],"supporting_citations":[{"why":"Defines the Mochima suite: the same galaxy re-simulated in one dark-matter-only and five baryonic runs that form the entire data set.","marker":"[46]"},{"why":"Provides the host dark matter structure analysis and the profile-fitting method used to measure inner slopes gamma.","marker":"[52]"},{"why":"Supplies the observed cumulative star formation histories and t90% values for Local Group dwarf galaxies used in the comparison.","marker":"[99]"},{"why":"Establishes the observed diversity of dwarf rotation curves that the paper aims to reproduce within cold dark matter.","marker":"[14]"},{"why":"Supplies the mechanism by which repeated supernova-driven gas outflows transform dark matter cusps into cores.","marker":"[90]"},{"why":"Shows that continued central star formation can drive re-contraction and reformation of cusps, supporting the oscillating regime.","marker":"[42]"},{"why":"A high-resolution isolated-dwarf study showing that cusp versus core outcome depends on the timing of star formation, the closest prior result to the t90 relation.","marker":"[109]"},{"why":"Warns that insufficient resolution inflates remnant sizes and distorts tidal evolution, the caveat limiting interpretation of heavily stripped low-mass subhalos.","marker":"[110]"}],"fun_headline_variants":["Early star formation locks in dark matter cusps","Bursty star formation wiggles dwarf dark matter cores","Star formation timing rules dwarf dark matter slopes","Dark matter cusps vs cores: star formation history decides","Dwarf dark matter diversity traced to star formation timing"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The whole timing result rests on the premise that the measured inner dark matter slope within 1 kpc is a numerically converged, physically meaningful quantity even in small, heavily tidally stripped satellites; the paper flags that insufficient resolution can inflate remnant sizes and distort tidal evolution, so part of the oscillating regime could be numerical.","fun_headline_variants_meta":{"raw":{"variants":["Early star formation locks in dark matter cusps","Bursty star formation wiggles dwarf dark matter cores","Star formation timing rules dwarf dark matter slopes","Dark matter cusps vs cores: star formation history decides","Dwarf dark matter diversity traced to star formation timing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000185,"raw_usage":{"total_tokens":1175,"prompt_tokens":780,"completion_tokens":395,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":317}},"tokens_in":524,"tokens_out":395,"duration_ms":4368,"temperature":1.0,"reasoning_tokens":317,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T22:08:27.530749+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same simulations at higher mass and force resolution and check whether gamma oscillations in late-star-forming satellites persist; if they vanish, the fluctuating regime is numerical. Observationally, measure resolved star formation histories and dark matter inner slopes for a larger sample of Milky Way and Andromeda satellites: if late-forming dwarfs consistently show steep cusps, or early-forming dwarfs show cores, the t90% relation fails.","supporting_citations":[{"cited_title":"Cosmological simulations of the same spiral galaxy: the impact of baryonic physics","cited_arxiv_id":"2004.06008","evidence_quote":"Defines the Mochima suite: the same galaxy re-simulated in one dark-matter-only and five baryonic runs that form the entire data set."},{"cited_title":"Cosmological simulations of the same spiral galaxy: connecting the dark matter distribution of the host halo with the subgrid baryonic physics","cited_arxiv_id":"2301.06189","evidence_quote":"Provides the host dark matter structure analysis and the profile-fitting method used to measure inner slopes gamma."},{"cited_title":"The formation of cores in galaxies across cosmic time -- the existence of cores is not in tension with the LCDM paradigm","cited_arxiv_id":"2310.13055","evidence_quote":"Shows that continued central star formation can drive re-contraction and reformation of cusps, supporting the oscillating regime."}],"review_version":1}