{"id":"be64a3d4-e2d6-46d0-bc05-66f67bf3773f","arxiv_id":"2501.16317","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Stellar 3D shapes closely track dark matter halo shapes in simulated dwarf galaxies that lack a stellar disk, making stellar shape a plausible proxy for halo shape in low-mass dwarfs.","lead":"This paper measures the 3D shapes of stars and dark matter in 80 simulated dwarf galaxies and finds that in galaxies without a stellar disk, the stellar shape closely tracks the dark matter halo shape. It matters because measuring only stars in real dwarf galaxies could become a way to infer dark matter halo shapes and test dark matter models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Low-mass nondisky sample may measure shapes inside constant-density DM cores, making stellar-DM agreement an artifact; need density-slope check at 2Reff.","rationale":"The reader's weakest assumption is that the shape-recovery pipeline at 2Reff is unbiased. I agree this is the load-bearing premise, but I identify a more specific and more directly testable failure mode: near-constant-density DM cores that extend beyond 2Reff in the lowest-mass nondisky galaxies. The paper asserts, without measuring core radii, that 2Reff is outside a typical DM core. If that assertion is false for a substantial part of the low-mass sample, then both stellar and DM axis ratios are recovered in a regime where the shape is ill-defined, and the apparent stellar-DM agreement could be an artifact rather than evidence for the physical proxy claim. This does not overturn the paper's other results (e.g., disky galaxies differing from their halos, observed stellar-shape validation), so the appropriate disposition remains conditional: the central claim should be accepted only after the core-radius/density-slope check is passed. I therefore keep the reader's CONDITIONAL verdict unchanged while sharpening the condition that must be met.","tokens_in":26185,"tokens_out":7029,"duration_ms":77227,"concrete_test":"For every galaxy, compute spherically averaged DM and stellar density profiles and measure the local log-slope alpha = dln(rho)/dln(r) at 2Reff, or fit a cored profile to obtain r_core. Define a 'well-defined shape' subsample with |alpha| at 2Reff above a threshold (e.g., |alpha| > 0.3, or 2Reff > 1.5 r_core), following Fischer & Valenzuela (2023). Recompute Figures 3, 4, and the KS tests on this subsample. If the nondisky stellar-DM Q/S agreement and near-unity SDM/S* ratios weaken or disappear, the central claim is an artifact of ill-defined fits in constant-density cores; if the agreement persists, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that stellar shapes trace DM halo shapes in nondisky dwarfs rests on shape measurements at 2Reff (Section 2.3). The paper itself cites Fischer & Valenzuela (2023) that shape measurements are undefined in regions of near-constant density, but it only asserts that 2Reff lies outside a typical DM core ('a higher radius than a typical DM core at approximately the size of the galaxy half-light radius') without measuring core radii. For the lowest-mass nondisky galaxies (M* ~ 10^6 M_sun), Reff can be well below 0.3 kpc while feedback-generated cores in dwarfs extend to ~0.5-1 kpc, so 2Reff may fall inside a flat-density core. Inside such a core the iterative ellipsoid fit is ill-conditioned: the local density gradient is too weak to define isodensity surfaces, and the recovered axis ratios for both stars and DM are dominated by Poisson noise and smoothing rather than by physical shape. This would produce the illusion of strong stellar-DM shape agreement exactly in the low-mass regime where the paper claims the proxy is most reliable, and could explain the KS p=0.49 for Q and near-unity SDM/S* at low masses. The paper's own Figure 16 shows shapes becoming rounder toward the center, consistent with this artifact. The concern is testable.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper measures intrinsic 3D shapes of the stellar and dark-matter components in 80 simulated dwarf galaxies from the Marvel, DC Justice League, and Massive Dwarfs zoom-in suites, spanning stellar masses of about 1e6 to 1e10 Msun. Shapes are derived from iterative inertia-tensor fits in radial shells, evaluated at twice the effective radius (2 Reff), and summarized by axis ratios Q=B/A and S=C/A and by triaxiality. The central claim is that, for galaxies without a stellar disk, the stellar distribution closely tracks the dark-matter halo shape, so stellar shapes can serve as a proxy for DM shapes in low-mass dwarfs; the paper also reports that disks produce measurable stellar-DM shape differences, that triaxiality tracks between the two components, that axes are generally aligned, that two supernova feedback implementations give similar shapes, and that recent mergers with ratios greater than about 4 do not strongly perturb the measured shapes. The results are validated against observed dwarf shapes from Kado-Fong et al. (2020).","tokens_in":26406,"tokens_out":7022,"duration_ms":72693,"significance":"If the central claim is correct, the paper provides a practical route to infer DM halo shapes from stellar photometry in dwarf galaxies, which would be valuable for discriminating CDM from SIDM and for understanding baryonic effects on halo structure. The study has notable strengths: it uses multiple simulation suites with different resolutions and feedback models, it reports per-galaxy comparisons as well as ensemble KS tests, it includes an observational validation, and it states that the analysis code and figure data will be publicly available. The main caveat, discussed below, is whether the shape measurement at 2 Reff for the lowest-mass nondisky galaxies is well defined in regions that may have near-constant DM density; if that concern is not resolved, the inferred stellar-DM shape agreement could be partly an artifact of the measurement pipeline.","major_comments":[{"comment":"The load-bearing claim that stellar shapes trace DM shapes in nondisky dwarfs depends on shape measurements at 2 Reff, but the paper does not measure the local DM density slope or core radius at 2 Reff for the low-mass sample. Section 2.3 asserts that 2 Reff is 'a higher radius than a typical DM core' and cites Read et al. (2016) and Fitts et al. (2017, 2019), yet no density profile or core radius measurement is presented. For the lowest-mass galaxies (M* ~ 1e6 Msun, where Section 3.3 reports 2 Reff values as small as 0.6 kpc), feedback-generated cores on the order of 0.5-1 kpc would place 2 Reff inside a near-constant-density region. In such a region the iterative ellipsoid fit is ill-conditioned, and the recovered axis ratios for both stars and DM may be dominated by Poisson noise and smoothing rather than by physical shape, biasing the two components toward spurious agreement. Figure 16, which shows shapes becoming rounder toward the center, is consistent with this concern. The authors should report d log rho_DM/d log r at 2 Reff (or equivalently, core radii) for the nondisky sample and demonstrate that the KS result (p=0.49 for Q) and the near-unity SDM/S* ratios at low masses are not driven by galaxies with flat central density slopes, for example by repeating the analysis after excluding systems with slope near zero at 2 Reff.","section":"Section 2.3 / Section 3.2"}],"minor_comments":[{"comment":"The conclusion that mergers with ratios greater than about 4 do not perturb galaxy shapes is based on two merger events with ratios 4.4 and 5.6, as the text acknowledges. The abstract and conclusions state this as a general result ('a dwarf galaxy's shape is largely unperturbed by recent mergers (with merger ratios >4)'); this should be rephrased to indicate the small number of events and the limited range of merger ratios probed.","section":"Section 3.5 / Abstract"},{"comment":"The claim that shape measurements are 'robust to different implementations of baryonic feedback' rests on a comparison of eight galaxies in one simulation volume. The authors note the small sample in Section 4.2, but the abstract and summary conclusion do not carry this caveat; the wording should be softened to reflect the limited statistical power.","section":"Section 4.2 / Abstract"},{"comment":"In the final paragraph of Appendix D, the text says 'Q = B/A and S = B/A values closer to one'; the second expression should be S = C/A.","section":"Appendix D"},{"comment":"The slopes quoted for QDM/Q* and SDM/S* versus stellar mass do not include a discussion of whether the individual galaxy ratios are consistent with unity within the estimated measurement uncertainties. Adding representative error bars or a statement about typical per-galaxy uncertainties would make the low-mass 'near-unity' claim easier to evaluate.","section":"Section 3.2 / Figure 4"},{"comment":"The statement that '82% of our sample is aligned according to our expectations from shapes' uses shape categories that are derived from the same axis-ratio measurements used to define the expected alignment patterns, so this is a consistency check rather than an independent confirmation; the wording should make that explicit.","section":"Section 3.4"}],"recommendation":"major_revision","confidential_remarks":"The core-density concern in Section 2.3/3.2 is the main reason for the major revision recommendation. If the authors can show, from their own density profiles, that 2 Reff for the nondisky low-mass sample lies outside a flat-density core, or that the stellar-DM agreement persists when such systems are excluded, the central claim would be on much firmer ground. The merger and feedback conclusions are secondary and should be presented with explicit sample-size caveats in the abstract and conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The useful new result is the systematic, multi-suite comparison of intrinsic 3D shapes of stars and DM in dwarfs from 1e6 to 1e10 Msun, with the mass-dependent stellar/DM axis ratio ratios, triaxiality slopes, alignment fractions, and robustness checks. The central claim—nondisky dwarf stellar shapes trace DM halo shapes—is well supported: KS p=0.49 for Q, near-unity SDM/S* at low mass, triaxiality slope 0.99±0.08 with R2=0.77, and alignment consistent with shape categories. They also validate against Kado-Fong et al., which is a genuine check rather than a fit.\n\nWhat the paper does well: it measures shapes with an iterative inertia tensor, enforces particle thresholds and convergence, reports slopes and scatter, and includes radial-profile appendices. The disk classification is checked visually and kinematically; only 14% sit near the boundary. Code and figure data are public. Simulations are proprietary, but the paper says as much and offers collaboration.\n\nSoft spots, in proportion. The merger claim is the weakest: it rests on two galaxies with ratios 4.4 and 5.6, and the abstract overstates it more than the text does. The feedback null is based on eight matched galaxies in one volume; useful, but underpowered for a strong claim of robustness. The three suites differ in resolution and cosmology; they do check environment and see no difference, which mitigates this.\n\nThe core-radius concern is worth taking seriously. They cite Fischer & Valenzuela that shape measurements are undefined in near-constant-density regions, and they assert 2Reff lies outside typical DM cores, but they do not measure density slopes or core radii for their sample. For the lowest-mass galaxies this is checkable. It is not obviously fatal—their nondisky low-mass shapes are prolate/triaxial, not spherical, so they are not simply reporting round objects—but a density-slope measurement at 2Reff would settle it. I would want that added before citing this as the definitive proxy result.\n\nCircularity is not an issue: shapes are measured and compared, not fitted to the conclusion. The citation pattern is appropriate; they engage EDGE, FIRE, Chua et al., and Valenzuela, and the novelty claim holds.\n\nWho is this for: observers who want to use dwarf stellar shapes to infer DM shapes, and simulators comparing CDM vs SIDM predictions. It deserves a serious referee, and I would cite it with the core-radius caveat noted.","headline":"A solid, careful simulation study showing nondisky dwarf stellar shapes track DM halo shapes across 1e6-1e10 Msun; side claims on mergers and feedback are softer than the abstract implies.","tokens_in":27038,"tokens_out":3432,"would_cite":true,"duration_ms":34754,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d","98.62.Gq"],"model":"deepseek-v4-flash","headline":"In diskless dwarf galaxies, the stellar distribution mirrors the dark matter halo's 3D shape, making stellar shapes a practical proxy for halo shapes and a route to testing dark matter models.","keywords":["dwarf galaxies","dark matter halo shapes","galaxy intrinsic shapes","triaxiality","cold dark matter simulations","cosmological zoom-in simulations","stellar disks","supernova feedback"],"falsifier":"Recompute the shapes of the same simulated galaxies with a different measurement pipeline — for instance the standard (non-reduced) inertia tensor, a different radial binning scheme, or much stricter per-bin particle thresholds — and repeat the Kolmogorov–Smirnov tests: if the $Q$ distributions of nondisky stars and dark matter are no longer statistically indistinguishable (p=0.49) and the per-galaxy $S_{\\rm DM}/S_*$ ratios drift away from unity, the claimed correspondence is a product of the fitting procedure rather than the physics. An observational check would be deep imaging of a few dozen field dwarfs below $10^{7.5}\\,M_\\odot$ whose deprojected stellar shapes come out systematically rounder or more oblate than the prolate, triaxial halo shapes the simulations predict.","tokens_in":25943,"feed_emoji":"🌌","tokens_out":16922,"duration_ms":129185,"temperature":0.7,"pith_summary":"The paper sets out to establish that the 3D shape of a dwarf galaxy's stars can stand in for the 3D shape of its dark matter halo — a claim that matters because halo shape is one of the few predicted differences between cold dark matter and self-interacting dark matter. Using 80 dwarf galaxies from three cosmological zoom-in simulation suites spanning stellar masses $10^6$–$10^{10}\\,M_\\odot$, the authors compare axis ratios and triaxiality of stars and dark matter measured at twice the effective radius. They find that in galaxies without a stellar disk the stellar and dark-matter axis-ratio distributions are statistically indistinguishable, whereas in disky galaxies the disk flattens the stellar distribution and breaks the correspondence. If correct, the result lets observers use plain starlight in low-mass, diskless dwarfs to infer the shape of the underlying dark matter halo and thereby discriminate between dark matter models.","feed_headline":"For diskless dwarfs, starlight reveals the dark halo's shape","feed_subtitle":"In 80 simulated dwarfs, stellar and dark matter shapes match when no disk forms, so stars can probe dark matter.","key_machinery":"The load-bearing object is the shape tensor $S_{ij} = (1/M)\\sum_k m_k r_{k,i} r_{k,j}$ — the moment-of-inertia tensor of particle positions — evaluated in iteratively fitted ellipsoidal shells whose eigenvalues give the principal axes $A \\geq B \\geq C$ and hence the axis ratios $Q = B/A$ and $S = C/A$. After smoothing the axis-ratio profiles $Q(r)$ and $S(r)$ with 3rd- to 5th-order polynomials, the shapes are read off at twice the effective radius ($2 R_{\\rm eff}$), the radius where self-interacting-dark-matter sphericalization would be detectable. Adaptive radial binning with a floor of 5000 star particles sets the resolution limits, and the triaxiality parameter $T = (1-Q^2)/(1-S^2)$ places each galaxy between oblate ($T < 1/3$), triaxial ($1/3 < T < 2/3$), and prolate ($T > 2/3$). The disk classification — thin ($S_* < 0.4$) and circular ($Q_* > 0.65$) — separates the cleanly dark-matter-tracing population from the disk-dominated one.","core_discovery":"The paper claims that stellar shape follows dark matter halo shape across its 80 simulated dwarf galaxies. The correspondence is strongest where baryons are dynamically subdominant: in nondisky galaxies, which dominate below stellar mass about $10^{7.5}\\,M_\\odot$, a two-sample Kolmogorov–Smirnov test finds the stellar and dark-matter distributions of the intermediate-to-major axis ratio $Q$ statistically indistinguishable (p=0.49) and the minor-to-major ratio $S$ only modestly different (p=0.015), with per-galaxy $S_{\\rm DM}/S_*$ ratios near unity. Stellar triaxiality tracks halo triaxiality with slope $0.99\\pm0.08$ ($R^2=0.77$) in nondisky galaxies, while disky galaxies deviate from the one-to-one relation because the disk makes stars flatter ($S_{\\rm DM}/S_*$ averages $2.92\\pm0.60$). The authors further argue that disk formation above about $10^{7.5}\\,M_\\odot$ starts to round and flatten the dark matter itself, that stellar and dark-matter axes are strongly aligned (most clearly for the minor axis in disky systems), that swapping between two supernova feedback implementations leaves both shapes unchanged, and that recent mergers with ratios above about 4 move the axis ratios by no more than about 0.05. The paper concludes that stellar shape measurements are a reliable tool for inferring dark matter halo shapes in dwarf galaxies.","pith_inferences":["If stellar shape is a faithful halo-shape tracer, then the sphericalization that self-interacting dark matter produces in inner halos should appear as a measurable roundness trend in the stellar axis-ratio profiles $Q(r)$ and $S(r)$ of low-mass diskless dwarfs toward the center — a test the paper gestures at but does not run.","The $10^{7.5}\\,M_\\odot$ disk-formation threshold suggests a practical survey design rule: shape-based dark-matter probes should be restricted to diskless dwarfs below this mass, while shape distributions of higher-mass dwarfs should be treated as baryon-contaminated.","The proxy claim could be stress-tested by forward-modeling the observational pipeline — projecting the simulated galaxies to 2D, measuring ellipticities as a survey would, and applying standard deprojection techniques — to verify that the stellar–DM shape correlation survives the inference process before trusting it on real data.","If the proxy holds, the population distribution of dwarf galaxy shapes becomes a dark-matter-model diagnostic: the fractions of prolate, triaxial, and oblate field dwarfs are set by structure formation in CDM and would be shifted by dark-matter self-interactions, so shape surveys could constrain the self-interaction cross section."],"forward_implications":["Observed stellar shapes in nondisky dwarf galaxies below about $10^{7.5}\\,M_\\odot$ can be read directly as dark matter halo shapes, giving observers a handle on halo shape without kinematics or lensing.","Measuring stellar triaxiality in a sample of diskless dwarfs can discriminate between the prolate halos that cold dark matter produces and the rounder halos that self-interacting dark matter predicts.","The stellar disk is the main thing that breaks the stellar–DM correspondence: above $10^{7.5}\\,M_\\odot$ the disk makes stars much flatter than the halo, so dark-matter inference should be restricted to diskless galaxies.","Shape measurements appear insensitive to the supernova feedback implementation, so dwarf-galaxy shapes do not encode subgrid baryonic physics and results from different simulation suites can be combined.","Recent mergers with ratios above 4 change both stellar and dark-matter axis ratios by less than about 0.05, so merger activity does not need to be accounted for when interpreting dwarf galaxy shapes."],"supporting_citations":[{"why":"Supplies the iterative ellipsoid-fitting algorithm whose converged eigenvalues define the principal axes and axis ratios.","marker":"M. Tomassetti et al. (2016)"},{"why":"Provides the adaptive radial binning strategy and particle-count floor that the shape measurements adopt.","marker":"M. Zemp et al. (2011)"},{"why":"Origin of the shape tensor definition used to measure the ellipsoidal shells.","marker":"B. Allgood et al. (2006)"},{"why":"The observational dwarf-shape distribution the simulations are validated against and the main observational comparator.","marker":"E. Kado-Fong et al. (2020)"},{"why":"Shows shape fits converge to good accuracy even with few hundred particles and supplies the low-mass dark-matter-only halo shapes the nondisky dwarfs resemble.","marker":"C. A. Vera-Ciro et al. (2014)"},{"why":"Establishes that baryonic feedback leaves dwarf halo shapes largely unchanged while rounding higher-mass halos, anchoring the mass-threshold interpretation.","marker":"K. T. E. Chua et al. (2022)"},{"why":"Shows gas-poor ultrafaint dwarfs keep prolate dark matter halos, the low-mass behavior this paper extends to the full dwarf range.","marker":"M. D. A. Orkney et al. (2023)"},{"why":"Describes the superbubble supernova feedback model used in the Massive Dwarfs suite, central to the feedback-robustness comparison.","marker":"B. W. Keller et al. (2014)"},{"why":"Directly compares blastwave and superbubble feedback in the Marvel suite and finds similar results, justifying the combination of feedback models.","marker":"B. Azartash-Namin et al. (2024)"}],"fun_headline_variants":["Starlight traces dark halos in dwarf galaxies","In diskless dwarfs, stars reveal dark matter shapes","Dwarf galaxy stars mirror dark matter halo shapes","Diskless dwarfs: stellar shapes match dark halos","Stars in dwarf galaxies trace dark matter shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole stellar–DM comparison rests on the assumption that the iterative shape-tensor fitting, with its chosen radial bins, particle-count floors, and polynomial smoothing, recovers the true 3D shapes of stars and dark matter equally well at twice the effective radius; if the measurement pipeline biases the two components differently, the apparent shape correspondence could be partly artificial.","fun_headline_variants_meta":{"raw":{"variants":["Starlight traces dark halos in dwarf galaxies","In diskless dwarfs, stars reveal dark matter shapes","Dwarf galaxy stars mirror dark matter halo shapes","Diskless dwarfs: stellar shapes match dark halos","Stars in dwarf galaxies trace dark matter shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000351,"raw_usage":{"total_tokens":2030,"prompt_tokens":1180,"completion_tokens":850,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":796,"completion_tokens_details":{"reasoning_tokens":775}},"tokens_in":796,"tokens_out":850,"duration_ms":7686,"temperature":1.0,"reasoning_tokens":775,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:33:02.406099+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the shapes of the same simulated galaxies with a different measurement pipeline — for instance the standard (non-reduced) inertia tensor, a different radial binning scheme, or much stricter per-bin particle thresholds — and repeat the Kolmogorov–Smirnov tests: if the $Q$ distributions of nondisky stars and dark matter are no longer statistically indistinguishable (p=0.49) and the per-galaxy $S_{\\rm DM}/S_*$ ratios drift away from unity, the claimed correspondence is a product of the fitting procedure rather than the physics. An observational check would be deep imaging of a few dozen field dwarfs below $10^{7.5}\\,M_\\odot$ whose deprojected stellar shapes come out systematically rounder or more oblate than the prolate, triaxial halo shapes the simulations predict.","supporting_citations":[{"cited_title":"2016,MNRAS ,458, 4477","cited_arxiv_id":null,"evidence_quote":"Supplies the iterative ellipsoid-fitting algorithm whose converged eigenvalues define the principal axes and axis ratios."},{"cited_title":"Y., Gnedin, N","cited_arxiv_id":null,"evidence_quote":"Provides the adaptive radial binning strategy and particle-count floor that the shape measurements adopt."},{"cited_title":"E., Huang, S., et al","cited_arxiv_id":null,"evidence_quote":"The observational dwarf-shape distribution the simulations are validated against and the main observational comparator."},{"cited_title":"A., Sales, L","cited_arxiv_id":null,"evidence_quote":"Shows shape fits converge to good accuracy even with few hundred particles and supplies the low-mass dark-matter-only halo shapes the nondisky dwarfs resemble."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that baryonic feedback leaves dwarf halo shapes largely unchanged while rounding higher-mass halos, anchoring the mass-threshold interpretation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows gas-poor ultrafaint dwarfs keep prolate dark matter halos, the low-mass behavior this paper extends to the full dwarf range."},{"cited_title":"W., Wadsley, J., Benincasa, S","cited_arxiv_id":null,"evidence_quote":"Describes the superbubble supernova feedback model used in the Massive Dwarfs suite, central to the feedback-robustness comparison."}],"review_version":1}