{"id":"5543bb5c-1caa-43fd-a126-9eae26d14d56","arxiv_id":"2508.16759","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Stars and dark matter in tidally stripped two-component galaxies lose mass at the same truncation energy, enabling dark matter constraints from stellar kinematics.","lead":"This paper uses idealized N-body simulations to show that stars and dark matter in a dwarf galaxy are tidally stripped at the same binding energy, regardless of their density profiles. If confirmed, the energy of stripped stars in tidal streams could be used to infer the dark matter content of disrupted galaxies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unbinding procedure's self-potential-only criterion could make the star/dark-matter truncation-energy equality an artifact.","rationale":"The reader's weakest_assumption identifies the iterative unbinding procedure as the load-bearing step, and I agree. The central claim is about measured truncation energies of the two components, and those measurements are outcomes of the bounded-remnant identification. Because the unbinding procedure neglects the host potential and assumes spherical symmetry, and because the two components have markedly different radial extents, the systematic error in the binding energy is likely to differ between stars and dark matter. If this error correlates with the truncation-energy measurement, the equality could be induced by the method rather than by the physics. The proposed test directly compares an alternative bound-remnant definition and would settle whether the equality persists when the host potential is included. The reader's verdict of CONDITIONAL already appropriately requires addressing this and other issues (overreach in 'regardless of initial profiles', unavailability of code/data, missing M200 energy-space plots), so I do not recommend changing the verdict. Secondary concerns about the limited parameter space and the definition of truncation at the 50% point are real but less load-bearing: even if the equality holds only for the tested cases, the paper's main phenomenology stands, whereas an unbinding artifact would invalidate the headline claim outright.","tokens_in":15193,"tokens_out":9901,"duration_ms":110006,"concrete_test":"Recompute the bound remnants for simulations M10_10A and M200_10A using an alternative binding criterion that includes the host potential—for example, the boosted-potential method (Stucker et al. 2021) or a Jacobi-energy criterion with the NFW host potential added to the satellite potential—and recompute the truncation energies Et,star and Et,dark as defined in §3.3. If the difference |Et,star − Et,dark| remains below 1% of Et for both simulations, the equality is robust; if the difference grows to the size of the reported 0.5% residuals or larger, the claim is an artifact of the self-potential-only unbinding procedure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that stars and dark matter are stripped to the same truncation energy rests entirely on the bound remnants identified by the iterative unbinding procedure in §2.4. That procedure uses only the satellite's self-potential (Eq. 10), explicitly neglecting the host potential, and assumes spherical symmetry. Because the stellar component is more concentrated than the dark matter in all models (rdark/rstar = 2 or 10), the error from neglecting the host tidal field is asymmetric: dark matter particles at large radii are more likely to be misclassified as self-bound, shifting the dark-matter truncation energy toward lower initial binding energies. The stellar component, being more tightly bound, is affected less. The observed equality to <0.5% could therefore be a property of the unbinding algorithm rather than a physical law. The paper provides no cross-check with an alternative bound-remnant definition (e.g., a Jacobi/tidal-radius criterion or the boosted-potential framework), and it does not quantify the non-sphericity of the remnants that would justify Eq. (10).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents 16 idealized N-body simulations of two-component (stellar + dark matter) satellites with gamma profiles embedded in a fixed NFW host potential. It reports that in initial binding-energy space, stars and dark matter are stripped to the same truncation energy to within typically 0.5%, across the explored inner slopes, mass/size ratios, and orbital configurations. It further finds that systems with a cuspy stellar component and a cored dark-matter halo evolve into dark-matter-deficient (DMD) galaxies and illustrates this with a model of NGC 1052-DF2. The authors advocate an energy-space framework for multi-component tidal stripping and suggest that observable stellar kinematics in tidal streams could constrain the dark-matter content and orbital history of disrupted dwarfs.","tokens_in":15486,"tokens_out":8842,"duration_ms":106610,"significance":"If the claimed energy-truncation equality holds beyond the simulated parameter space, it would be a genuinely useful simplifying principle: stripping of stars and dark matter would be organized by a single initial binding-energy threshold, independent of component density profile. The paper's strengths are that the simulations are clearly specified, equal-mass particles are used to avoid spurious heating, stability tests are shown, and the truncation-energy residuals are small and displayed. The central result is measured rather than fitted, so it is not circular in an obvious way. However, the strength of the claim currently exceeds the evidence: the bound-remnant definition is not cross-checked, the parameter grid is narrow, and the full survival functions are not compared between components. These issues are addressable but must be fixed before the universal statement can be accepted.","major_comments":[{"comment":"The iterative unbinding that defines the bound remnant uses only the satellite's spherical self-potential, explicitly neglecting the host potential. Since dark matter is more extended than the stars (rdark/rstar = 2 or 10), any bias from neglecting the tidal field and non-inertial frame is asymmetric between the components. The equality Et,star ≈ Et,dark in Fig. 9 could therefore be a property of the 'self-bound' definition rather than a physical stripping law. Please cross-check with a Jacobi/tidal-radius or boosted-potential binding criterion (e.g., Stücker et al. 2023), or quantify the host-tidal correction to Eq. (10). Without this, the central claim is not robustly established.","section":"§2.4, Eq. (10)"},{"comment":"The statement that stripping is identical 'regardless of initial profiles' is too broad for the explored grid: only γstar, γdark ∈ {0,1}, Mdark/Mstar ∈ {10,200}, rdark/rstar ∈ {2,10}, two orbits, and one host potential. All satellites are spherical and isotropic, and the M200 energy–circularity maps are not shown. The data support a statement such as 'for the range of two-component gamma-profile models explored here'. A universal law would require broader variation, including continuous slopes, anisotropy, orbital families, and host concentrations. This matters because the abstract's strongest claim is universality.","section":"Abstract and §3.3, Fig. 9"},{"comment":"The 'identical stripping' conclusion is based on a single scalar, the 50% truncation energy Et. However, Figs. 7 and 8 show a gradual transition in the survival fraction and a noticeable circularity dependence, and the residual analysis compares only Et. For the M200 runs, no energy–circularity maps are shown, so it is unknown whether the full survival function—not just the median contour—is identical for stars and dark matter. Please report a quantitative comparison of survival functions (e.g., the width and circularity dependence of the truncation transition) for each component, or explicitly restrict the claim to Et.","section":"§3.3, Figs. 7–9"}],"minor_comments":[{"comment":"The definition of tunit appears to have a typo: 'tunit = sqrt(r_s^3/(G M_sat^3))' should be sqrt(r_s^3/(G M_sat)); the printed expression has incorrect dimensions.","section":"§2.3"},{"comment":"The text says 'As with M10 01A' when discussing M10 10B, but the figure caption compares with M10 00A. Please clarify which simulation is the reference.","section":"§3.3, Fig. 8"},{"comment":"The x-axis is labeled 't' but the text says 'as a function of tunit'. Label the axis as 't/tunit' or note that t is in code units.","section":"Fig. 9"},{"comment":"The phrase 'The M200 simulations (not shown)' is misleading because M200 results appear in Fig. 9. Rephrase to indicate that only the energy–circularity maps are not shown.","section":"§3.3"},{"comment":"The DF2 case study finds the DMD phase lasts only about one orbital period (~200 Myr) before disruption. This transient nature deserves more discussion if the model is to support the existence of observed DMD galaxies.","section":"§4, Fig. 11"}],"recommendation":"major_revision","confidential_remarks":"This is a solid idealized simulation study with a clear and potentially important result, but the headline universality claim is currently stronger than the evidence. The requested cross-check of the bound-remnant definition and a more careful qualification of the parameter coverage should be achievable without redoing the entire simulation suite. No concerns about novelty; the two-component extension is a useful step beyond single-component energy-truncation models."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the main thing to know: the paper reports that in idealized N-body simulations of a star+DM satellite in a fixed host potential, both components are stripped to essentially the same truncation energy, defined by the energy at which half the initial particles remain. That equality is new for two-component systems, and the residual plots show it holding to better than ~0.5% across their runs. I buy that as a genuine extension of the single-component energy-truncation picture, and the simulations are clean: equal-mass particles, stability checks out to 5000 t_unit, and the gamma-profile family covers cuspy and cored stars and DM.\n\nThe soft spots are real but not fatal. The unbinding procedure in §2.4 uses only the satellite self-potential (Eq. 10), which is standard practice but not the only choice. The DM is more extended than the stars in every model, so if the self-potential-only criterion misclassifies some outer particles as bound, the two components could be affected asymmetrically. The paper itself notes the host potential is neglected but gives no cross-check with an alternative bound-remnant definition, say a boosted-potential or tidal-radius criterion. That's the load-bearing part of the measurement, so I'd want to see that check before calling the equality universal.\n\nSecond, the parameter grid is small: four gamma pairs, two mass/size ratios, two orbits. The phrase 'regardless of their initial profiles' outruns that. The M200 runs are mentioned but not shown in the main figures—they say the pattern is the same, but the reader can't verify. Also, the truncation-energy residuals have no formal error bars; 'typically less than 0.5 per cent' is a descriptive statement, not a quantified uncertainty.\n\nThe DF2 case study is a reasonable illustration, not a decisive test, and the discussion honestly lists missing baryonic physics.\n\nNet: this is a worthwhile paper for the tidal-stripping community, and the central two-component result deserves a careful referee. The authors need to either expand the parameter coverage or soften the 'regardless' claim, and they should test the binding criterion's sensitivity. I would send it to peer review.","headline":"Two-component tidal stripping in energy space looks like a real result, but 'regardless of initial profiles' is a stretch given the narrow grid and the untested binding criterion.","tokens_in":15910,"tokens_out":3298,"would_cite":true,"duration_ms":41534,"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":"Tidal stripping cuts stars and dark matter at the same energy.","keywords":["tidal stripping","dwarf galaxies","dark matter-deficient galaxies","energy truncation","N-body simulations","cusp-core structure","ultra-diffuse galaxies","stellar kinematics"],"falsifier":"Re-run the same orbits with a binding-energy criterion that includes the host's tidal potential (for example, the boosted-potential definition) and check whether stellar and dark-matter truncation energies still agree to within about one percent; if they separate, the matched cutoff is an artifact of the self-potential-only unbinding.","tokens_in":15125,"feed_emoji":"🌌","tokens_out":6825,"duration_ms":72855,"temperature":0.7,"pith_summary":"The paper uses idealized N-body simulations of two-component dwarf galaxies—stars embedded in dark matter—orbiting in a fixed host potential, and asks whether tidal stripping removes particles by anything other than how tightly bound they were at infall. It claims that, in energy space, stars and dark matter are stripped identically: both components lose particles down to the same truncation energy, even when their inner density profiles are very different (cored versus cuspy). If this holds, the energy distribution of stars stripped from a disrupted dwarf is a direct readout of the dark matter structure and orbital history of the parent system, without assuming a particular dark matter profile. The paper also argues that a galaxy with cored dark matter and cuspy stars naturally becomes dark-matter-deficient through tides, matching the properties of NGC 1052-DF2 and similar systems.","feed_headline":"Tidal stripping cuts stars and dark matter at the same energy","feed_subtitle":"Those equal energy cutoffs let shredded stars reveal a dwarf galaxy's dark matter mass.","key_machinery":"The load-bearing machinery is the initial energy-circularity plane. Each particle is tagged with its binding energy E normalized by the satellite's initial central potential, and its circularity epsilon_circ = L/L_max; at every snapshot the simulation keeps only particles with E < 0 computed iteratively from the satellite's own potential. The claim is carried by the truncation energy E_t—the energy at which 50% of the initial particles remain—which the paper shows is common to stars and dark matter in every run. Supporting this is the gamma profile, an (alpha,beta)=(1,4) double power law whose inner slope gamma encodes cuspiness, used to build all satellites.","core_discovery":"At the level of the paper's own claims, the discovery is a statement of conservation: tidal stripping removes particles by a threshold in initial binding energy, and that threshold is identical for the stellar and dark-matter components. Across 16 simulations spanning mass ratios of 10:1 and 200:1, two different orbits, and the four combinations of cored or cuspy stellar and dark-matter profiles, the truncation energy—defined as the energy at which half of the initially bound particles remain—agrees between components to better than about 0.5 percent. The same pattern appears in a tailored model of NGC 1052-DF2, where the surviving stars are the most tightly bound material at infall. Consequ","pith_inferences":["The equality of truncation energies suggests tidal stripping is a purely gravitational selection on initial binding energy that does not care what particle species carries it; a natural extension is that globular clusters or satellite galaxies embedded in the same halo would be stripped at the same E_t.","The paper's DMD timeline implies dark-matter-deficient dwarfs are a transient phase lasting about one orbit; if so, their observed abundance could constrain how often such orbits actually occur in groups and clusters, not just the formation mechanism.","A testable extension is to rerun the analysis with baryonic feedback or a live, evolving host potential; if the common E_t survives those changes, the inference from stellar streams to dark matter mass becomes robust in real galaxies.","One could also measure E_t directly in a stellar stream by fitting the energy cutoff of stream stars in the dwarf's rest frame and comparing it with the dark matter mass inferred from the stream's orbit—an observation that would test the paper's central claim."],"forward_implications":["Single-component energy-truncation models extend to two-component galaxies: one number, E_t, describes the stripping of both stars and dark matter.","From observed stellar kinematics alone one can estimate E_t and, combined with an independent estimate of the orbital tidal field, infer the bound dark matter mass of a disrupted dwarf.","Stellar velocity dispersions must remain below the local escape speed, giving a lower bound on E_t and hence on the bound dark matter mass of any self-bound satellite.","Cored dark matter with a cuspy stellar profile is a natural tidal pathway to dark-matter-deficient galaxies, with the DMD phase lasting roughly one orbital period before disruption.","If either component is cuspy, it partially shields the other from complete disruption, so cuspy systems survive longer and may be over-represented among observed remnants."],"supporting_citations":[{"why":"Introduces the energy-truncation model of subhalo stripping and the method for generating multi-component initial conditions used here.","marker":"Drakos et al. 2017"},{"why":"Supplies the iterative satellite-frame centering and self-potential approximation (their Equation 10) used to define the bound remnants.","marker":"Drakos et al. 2020"},{"why":"Shows the energy-truncation model reproduces cuspy and cored collisionless systems, the single-component result this paper extends to two components.","marker":"Drakos et al. 2022"},{"why":"Provides the energy-circularity binning used to display which particles survive stripping.","marker":"Choi et al. 2009"},{"why":"Establishes the universal outer-stripping behavior and justifies that the exact outer slope of the satellite profile is quickly erased.","marker":"Hayashi et al. 2003"},{"why":"Proposed that cored dark matter plus cuspy stars can produce dark-matter-deficient galaxies, the scenario this paper tests and supports.","marker":"Ogiya 2018"},{"why":"Provides the DF2 structural and orbital parameters and the roughly 9 Gyr accretion age used in the NGC 1052-DF2 case study.","marker":"Ogiya et al. 2022"},{"why":"Supplies the standard method for generating stable multi-component N-body initial conditions and flags numerical relaxation effects on central density.","marker":"Kazantzidis et al. 2004"},{"why":"Demonstrates that cored profiles lose central density and eventually disrupt, the contrast case for cuspy survival invoked in the discussion.","marker":"Peñarrubia et al. 2010"}],"fun_headline_variants":["Stars and dark matter get stripped at the same energy","Tidal stripping cuts stars and dark matter at identical energy","Shredded stars reveal dark matter's energy cutoff","Equal energy stripping links stars to dark matter halo","Stripped stars trace dark matter's energy threshold"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The iterative unbinding step keeps only particles with negative energy computed from the satellite's own potential, ignoring the host, and assumes this correctly identifies the physically bound remnant; if that step is biased, the matched truncation energies could be an artifact.","fun_headline_variants_meta":{"raw":{"variants":["Stars and dark matter get stripped at the same energy","Tidal stripping cuts stars and dark matter at identical energy","Shredded stars reveal dark matter's energy cutoff","Equal energy stripping links stars to dark matter halo","Stripped stars trace dark matter's energy threshold"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000144,"raw_usage":{"total_tokens":986,"prompt_tokens":694,"completion_tokens":292,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":438,"completion_tokens_details":{"reasoning_tokens":216}},"tokens_in":438,"tokens_out":292,"duration_ms":3425,"temperature":1.0,"reasoning_tokens":216,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T17:10:57.383655+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the same orbits with a binding-energy criterion that includes the host's tidal potential (for example, the boosted-potential definition) and check whether stellar and dark-matter truncation energies still agree to within about one percent; if they separate, the matched cutoff is an artifact of the self-potential-only unbinding.","supporting_citations":[{"cited_title":"E., Taylor, J","cited_arxiv_id":null,"evidence_quote":"Introduces the energy-truncation model of subhalo stripping and the method for generating multi-component initial conditions used here."}],"review_version":1}