{"id":"6bf12344-0dd2-4da6-848b-4c94280cb0e0","arxiv_id":"2505.21713","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Satellite mean stellar mass densities anti-correlate with cluster-centric distance, with a transition radius near 0.5 R_vir that could be used as a photometric probe of the host tidal field.","lead":"Satellite galaxies in clusters are stripped by the host's tidal force, and this study finds that their average stellar mass density rises toward the cluster centre, with a characteristic transition radius around half the cluster virial radius. If correct, this offers a purely photometric way to map the gravitational field of galaxy clusters without needing spectroscopy.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Transition-radius calibration is degenerate with the assumed pre-infall stellar density distribution: the toy model's IC scale (2×10^5 M⊙ kpc^-3) sets where Eq. (4) crosses |τ|, and the paper never quantifies how ℜ⋆ shifts under the 'different initial conditions' it says were tested.","rationale":"The paper's central observational claim is that the stellar-density transition radius ℜ⋆ is a photometric diagnostic of the host's tidal field. The physical mechanism (Eqs. 1-4) is standard tidal truncation, and the TNG/TNG50/EAGLE simulations plus the Fornax/Virgo data provide mutually supporting evidence that mean satellite densities rise inward. I credit the main trend as credible; the null tests in Fornax/Virgo are suggestive, though the combined 5σ statement is weakened by the non-significant MW/M31 p-values and by the arbitrary metric choice. The most load-bearing step is not the existence of the density-radius trend, but the calibration of ℜ⋆ to the host gravitational field. That calibration runs through the pre-infall density scale of the satellite population: in the toy model the transition appears where |τ(r)| crosses the population's initial mean density; if that density scale differs between environments, the same |τ(r)| would yield a different ℜ⋆. The paper asserts that different ICs were tested but does not report how ℜ⋆ responded; without that, a reader cannot tell whether the observed near-universality of ℜ⋆/Rvir is a property of the host field or an artifact of assuming similar satellite initial densities. This is precisely the concern raised by the reader, and I agree with it. The conditional verdict is appropriate: the main phenomenon is well supported, but the 'photometric diagnostic probe of host gravitational field' claim needs a quantitative demonstration of insensitivity to the initial density distribution.","tokens_in":34686,"tokens_out":6010,"duration_ms":66187,"concrete_test":"Take the TNG50 Fornax-type and TNG100 Virgo-type clusters; for each satellite, identify the subhalo progenitor at infall and record its mean stellar density within 2h⋆ (or the density at 2Rvir before stripping). Then rescale each satellite's present-day density by the ratio of the population-mean infall density to a reference value (e.g. ±0.5 dex) and re-run the moving-average/ℜ⋆ measurement of Sect. 5.1 in 3D and projection. If ℜ⋆ shifts by more than 0.1 Rvir, the transition radius is not a clean host-field probe; if it is stable, the initial-density degeneracy is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that ℜ⋆ 'probes the gravitational field strength of the host halo' requires that ℜ⋆ is set by the host tidal field rather than by the density scale of the satellite population before infall. In the toy model, Eq. (4) (ρ_tr ≥ |τ|) makes stripping begin where the host tidal field crosses the satellite's mean density; the population's pre-infall log-normal is centred at ρ_tr,init = 2×10^5 M⊙ kpc^-3 (Sect. 3.1). Changing that centre changes the crossing radius, so ℜ⋆ shifts unless the host profile also changes. The paper only says 'tested different initial conditions (IC) as well' with no quantitative result, leaving the degeneracy unquantified. The TNG simulations mitigate this because their infall densities are self-consistent, and the observed ℜ⋆/Rvir values roughly match TNG; however, the observed samples (Virgo, Fornax, Local Group) have no direct measure of their own pre-infall stellar density distribution. Since the proposed photometric diagnostic requires predicting ℜ⋆ from |τ| alone, a factor-of-two variation in the initial density scale (plausible given the log-normal width and different host masses) would move the inferred tidal field strength by a comparable factor, breaking the calibration. This is the load-bearing assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that satellite galaxies' mean stellar mass densities anti-correlate with cluster-centric distance because tidal stripping removes low-density outer material, leaving denser remnants near cluster centres. Using a semi-analytical toy model, IllustrisTNG and EAGLE simulations, and observed satellite populations in Fornax, Virgo, the Milky Way, and Andromeda, the authors propose a 'transition radius' Rstar ≈ 0.5 R_vir, inside which the moving-average stellar density profile rises toward the centre. They present this as a photometric diagnostic for the host halo's gravitational field strength, supported by null-hypothesis tests claiming a combined ~5sigma detection.","tokens_in":34941,"tokens_out":7024,"duration_ms":79922,"significance":"If the central claim is robust, the proposed relation would be a valuable, purely photometric tool for probing the tidal fields and dark-matter distributions of galaxy clusters, where spectroscopic kinematics for many dwarf satellites are expensive or impossible. The paper has real strengths: it combines three complementary approaches, uses public cosmological simulations at multiple resolutions and with two different subgrid models, explicitly checks projection effects, sample completeness, mass segregation, and gas content, and performs a null-hypothesis shuffle test for the observed samples. The TNG/EAGLE simulation results are self-consistent and the Fornax/Virgo observations are suggestive. However, the toy model's transition radius is essentially set by the assumed pre-infall satellite density scale, and the observed values of Rstar/Rvir in Table 1 do not cluster tightly at 0.5, so the headline universality and the diagnostic calibration are not yet established at the level the abstract claims.","major_comments":[{"comment":"The toy-model transition radius is set by the assumed initial density scale, not determined independently by the host potential. Equation (4) starts stripping where the satellite's mean density equals |tau|; with the initial total mean densities drawn from a log-normal centered at rho_tr,init = 2e5 M_sun kpc^-3, the location where the moving-average density profile rises is largely determined by that adopted centre and by the host profile. The text says 'tested different initial conditions (IC) as well' but reports no quantitative result for how Rstar shifts. This is load-bearing because the proposed photometric diagnostic requires Rstar to map to |tau| without knowing the pre-infall density distribution; the paper should quantify the shift in Rstar for a factor-of-two variation in the log-normal centre, or explicitly restrict the diagnostic to differential trends within a single population.","section":"Sect. 3.1, Eq. (4)"},{"comment":"The measurements reported do not support a universal value Rstar ≈ 0.5 R_vir as stated in the abstract. The observed values are Fornax [ii] = 0.36 R_vir, Virgo [ii] = 0.39 R_vir, MW [i] = 0.36 R_vir versus [ii] = 0.54 R_vir, and M31 [i]/[ii] = 0.49/0.46 R_vir, while the simulated values span roughly 0.40-0.58 R_vir depending on simulation and projection. The Fornax and Virgo offsets from 0.5 R_vir are much larger than the quoted bootstrap uncertainties, and the two metrics for the Milky Way differ by 40 kpc, suggesting systematics not captured by those uncertainties. The abstract and conclusions should state the measured range and the dominant systematic uncertainties rather than presenting a single 0.5 R_vir value.","section":"Table 1 and Sect. 5.1"},{"comment":"The observational comparison cannot currently distinguish a host-tidal-field calibration from environmental differences in pre-infall satellite structure. The observed samples provide no measurement of the pre-infall mean stellar mass density distribution for Virgo, Fornax, the MW, or M31, so the assumption that the toy-model initial log-normal is universal across these environments is untested. The null-hypothesis test demonstrates that the Fornax and Virgo moving-average minima are unlikely to be random (p = 4% and 0.1%), but it does not test whether the radius of the minimum is determined by the host tidal field rather than by the satellites' initial density structure. This caveat should be stated explicitly, or the authors should use the simulations to split satellite populations by infall density and show that Rstar/R_vir remains unchanged.","section":"Sects. 4.3 and 5.2"},{"comment":"The composite significance statement should be qualified. The individual p_NH values are 4% (Fornax), 0.1% (Virgo), 17% (MW), and 26% (M31), so the MW and M31 samples are individually consistent with the null hypothesis at the usual significance levels. The combined probability 1.8e-6 is dominated by Virgo, and the paper does not justify treating the four environments as independent, equally powered tests of the same null hypothesis. I recommend reporting the four individual p-values and the combined statistic with its assumptions, and avoiding the unqualified phrase 'highly significant' for the full observed sample.","section":"Sect. 5.2, Conclusions item IV"}],"minor_comments":[{"comment":"There is a typo: 'ellitpical' should be 'elliptical'.","section":"Introduction, first paragraph"},{"comment":"The moving-average definition sums over j = -q/2 to +q/2 but does not specify how bins near the radial boundaries are treated; please clarify the boundary convention.","section":"Eq. (17)"},{"comment":"For reproducibility, please list the adopted virial radius for each system directly in the table (or in a machine-readable footnote) and give Rstar/R_vir with the same asymmetric uncertainties as the kpc values.","section":"Table 1"},{"comment":"The polynomial order 7 is described as 'optimized,' but no optimization criterion or sensitivity test is given; please state the criterion and show that the measured Rstar is stable for polynomial orders 5-9.","section":"Sect. 5.1, Eq. (20)"},{"comment":"The phrase 'All are moving averages are calculated' should read 'All moving averages are calculated.'","section":"Captions of Figs. 4, A.1, A.2"},{"comment":"A data/code availability statement is missing; the paper relies on public simulations and catalogs, but the delorean orbital code and the analysis scripts used to produce Rstar are not linked to a repository or version.","section":"Global"}],"recommendation":"major_revision","confidential_remarks":"For the editor: this is a promising paper with a strong multi-method approach, and the simulation and observational evidence for a density-distance trend is real. The main problems are overclaims: the toy-model Rstar is calibrated by the assumed initial density scale, the observed Rstar/R_vir values do not converge tightly on 0.5, and the composite null-hypothesis significance is dominated by one system. I believe these issues are addressable with additional analysis or with careful reframing, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The central empirical result is real: in Virgo and Fornax the mean stellar mass density of satellites rises inward within a transition radius, and the trend is reproduced in TNG and EAGLE. That is a useful extension of the Local Group density-pericentre work to cluster scales, using good photometric samples. The null-hypothesis shuffling test is a nice touch, and the projection analysis is careful. But the paper sells the transition radius too hard. Observed R*/Rvir are 0.35-0.39 in Fornax/Virgo, 0.46-0.49 in M31, and 0.36-0.54 in the MW — not the abstract's ~0.5 Rvir. And the 'photometric diagnostic that probes the host gravitational field' is the load-bearing claim that isn't yet supported.\n\nWhy: Eq. (4) makes stripping begin where the satellite's initial mean density equals the local tidal field. The toy model sets that initial density from a log-normal centred at 2e5 M_sun/kpc3, and the paper says only that 'different initial conditions were tested' without reporting how R* shifts. If the pre-infall stellar density scale varies by environment or host mass, the crossing radius moves, and R* no longer isolates the halo field. The TNG/EAGLE clusters partly rescue this because their infall densities are self-consistent, and the observed R* do land near the TNG values. But the observed samples have no direct measure of their own pre-infall density distributions, so the calibration claim is a factor-of-two away from being demonstrated. The authors should either quantify that shift or reframe R* as a trend tracer rather than a calibrated field-strength probe.\n\nTwo smaller things. The combined 5sigma from p=4%, 0.1%, 17%, 26% is just the product of four p-values; that needs a stated combination method and an independence argument, and it hides that the Local Group systems are individually insignificant. And no toy-model code is released, which matters because the interpretation leans on it.\n\nNet: the density-distance relation itself is solid and worth citing. The universal 0.5 Rvir and the diagnostic calibration are overreaching. Send it to a referee, but the referee should push for a corrected abstract, a quantitative IC-degeneracy test, and a proper meta-analysis of the p-values.","headline":"A genuinely useful density-distance result for cluster satellites, but the 'universal 0.5 Rvir transition' and the 'photometric probe of halo gravity' claim outrun what the analysis actually shows.","tokens_in":35618,"tokens_out":3453,"would_cite":true,"duration_ms":38371,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a satellite galaxy's mean stellar mass density records the tidal field of its host cluster, with a transition near half the virial radius.","keywords":["satellite galaxies","galaxy clusters","tidal stripping","stellar mass density","dwarf galaxies","density-distance relation","tidal field diagnostics"],"falsifier":"Measure $\\Re_{\\star}$ in a cluster whose tidal field is independently known from satellite kinematics or weak lensing, and compare it with the radius where that tidal field equals the initial stellar density scale assumed in the model; the claim fails if the photometric transition sits elsewhere or shifts when the satellite sample is changed.","tokens_in":34412,"feed_emoji":"🌊","tokens_out":7722,"duration_ms":74621,"temperature":0.7,"pith_summary":"This paper argues that the mean stellar mass density of a satellite galaxy population is a readout of the gravitational tidal field of its host cluster. In a toy stripping model, in cosmological galaxy simulations, and in observed satellites around Fornax, Virgo, the Milky Way, and M31, the density profile is flat or slightly rising in the outskirts and turns upward inside a transition radius $\\Re_{\\star}\\approx0.5\\,R_{\\mathrm{vir}}$ where the host tidal field begins to dominate. The paper claims this stellar-density–distance relation ($\\bar{\\rho}^{\\star}-r$) works as a photometric tracer of the host tidal field, because it uses only stellar masses and half-light radii and needs no satellite kinematics. If correct, the transition radius gives a practical way to probe the dark matter distribution of distant clusters from imaging alone.","feed_headline":"Satellite galaxy cores expose cluster tides to 0.5 Rvir","feed_subtitle":"Photometry alone could map the dark matter pull of clusters too distant for orbital measurements.","key_machinery":"The load-bearing object is the moving average of the mean stellar mass density within twice the stellar half-mass radius, $\\langle\\bar{\\rho}^{\\star}_{2h_{\\star}}\\rangle(r)$, computed with a logarithmic filter over the satellite population. It is tied to the host through the Jacobi density identity $\\bar{\\rho}_J=|\\tau|$, the condition that material remains bound only where the satellite's mean density exceeds the local tidal field magnitude, and through the relation $|\\tau|\\propto 1/r$ for a spherically symmetric host. Tidal stripping shrinks $h_{\\star}$, raising the central stellar mass density, and the transition radius $\\Re_{\\star}$ is defined empirically either as the minimum of the moving average or as the zero of its radial derivative in the range $0.2<R/R_{\\mathrm{vir}}<1$.","core_discovery":"On its own terms, the paper's discovery is that the tidal truncation condition $\\bar{\\rho}_{\\mathrm{tr}}\\ge|\\tau|$, which ties a satellite's bound mean density to the host tidal field magnitude, imprints a feature in the population-average stellar density profile. The moving average of $\\bar{\\rho}^{\\star}_{2h_{\\star}}$ declines gently from the virial radius inward, reaches a minimum at $\\Re_{\\star}\\approx0.5\\,R_{\\mathrm{vir}}$, and then rises steeply toward the cluster centre, tracking the tidal field profile of an NFW-like host. The same feature appears in the idealized orbital model, in hydrodynamical cosmological simulations over a host mass range $10^{12}$–$10^{15.5}\\,M_\\odot$, and in all four observed systems, with projection lowering the measured radius by roughly 15%. The paper concludes that $\\Re_{\\star}$ marks the radius inside which the host tidal field has significantly processed stellar components and can serve as a photometric diagnostic of host gravitational field strength.","pith_inferences":["A testable extension not in the paper: stack many clusters and compare the photometric $\\Re_{\\star}$ with a tidal field mapped independently by weak lensing or satellite kinematics; agreement would validate the method while disagreement would localize the assumption.","Applying the same tracer to globular cluster systems, where internal kinematics are measurable in nearby galaxies, could calibrate the photometric diagnostic against dynamical tidal fields.","If pre-infall mean stellar densities differ systematically between environments, $\\Re_{\\star}$ would no longer isolate host properties; comparing clusters with similar virial mass but different merger histories would expose this dependence."],"forward_implications":["The $\\bar{\\rho}^{\\star}-r$ transition can be measured from photometry alone, extending tidal-field studies to clusters beyond the Local Group where satellite orbits cannot be measured.","Inside $\\Re_{\\star}$ the stellar components of surviving satellites have been tidally processed, so cluster-centric trends in galaxy size–mass relations must be interpreted with this trimming in mind.","The consistent $0.4$–$0.6\\,R_{\\mathrm{vir}}$ location across host masses implies the feature can benchmark the virial radii of clusters and groups.","The null-hypothesis test puts the joint chance of the observed transition appearing randomly in four systems near $2\\times10^{-6}$, supporting the reality of the feature."],"supporting_citations":[{"why":"Supplies the tidal truncation equations (Jacobi radius and truncation density) on which the toy model is built.","marker":"Binney & Merrifield (1998)"},{"why":"Provides the NFW host density profile used to compute the tidal field and virial radii.","marker":"Navarro et al. (1996)"},{"why":"Gives the abundance-matching relation used to set satellite dark matter halo masses before infall.","marker":"Moster et al. (2013)"},{"why":"Establishes the earlier Milky Way density–pericentre relation that this work extends to stellar densities and cluster-centric distance.","marker":"Kaplinghat et al. (2019)"},{"why":"Supplies the Milky Way and M31 satellite catalogues used for the Local Group comparison.","marker":"Putman et al. (2021)"},{"why":"Supplies the Virgo cluster satellite catalogue used for the observed density profile.","marker":"Kim et al. (2015)"},{"why":"Supplies the deep Fornax cluster satellite catalogue with stellar masses used in the observational analysis.","marker":"Eigenthaler et al. (2018)"},{"why":"Provides the Fornax cluster mass model used to define its tidal field profile.","marker":"Schuberth et al. (2010)"},{"why":"Provides tidal-track simulations showing that half-mass radii shrink during stripping, supporting the use of $h_{\\star}$ as a tracer.","marker":"Errani et al. (2022)"}],"fun_headline_variants":["Cluster tides leave a 0.5 Rvir mark in satellite stars","Satellite star densities reveal cluster tide frontier","Photometric radius maps dark matter halo strength","Tidal cutoff at 0.5 Rvir reshapes satellite galaxies","Cluster tide edge found in satellite stellar densities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument isolates the host tidal field only if satellites that have not yet been stripped have the same distribution of mean stellar mass densities in every environment; otherwise the radius where the host tidal field crosses that density scale is set partly by the satellites' initial structure, not by the host.","fun_headline_variants_meta":{"raw":{"variants":["Cluster tides leave a 0.5 Rvir mark in satellite stars","Satellite star densities reveal cluster tide frontier","Photometric radius maps dark matter halo strength","Tidal cutoff at 0.5 Rvir reshapes satellite galaxies","Cluster tide edge found in satellite stellar densities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00029,"raw_usage":{"total_tokens":1738,"prompt_tokens":1027,"completion_tokens":711,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":643,"completion_tokens_details":{"reasoning_tokens":633}},"tokens_in":643,"tokens_out":711,"duration_ms":7383,"temperature":1.0,"reasoning_tokens":633,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:25:07.054306+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $\\Re_{\\star}$ in a cluster whose tidal field is independently known from satellite kinematics or weak lensing, and compare it with the radius where that tidal field equals the initial stellar density scale assumed in the model; the claim fails if the photometric transition sits elsewhere or shifts when the satellite sample is changed.","supporting_citations":[{"cited_title":"P., Naab, T., White, S","cited_arxiv_id":null,"evidence_quote":"Gives the abundance-matching relation used to set satellite dark matter halo masses before infall."},{"cited_title":"E., Zheng, Y., Price-Whelan, A","cited_arxiv_id":null,"evidence_quote":"Supplies the Milky Way and M31 satellite catalogues used for the Local Group comparison."},{"cited_title":"2010, A & A, 513, A52","cited_arxiv_id":null,"evidence_quote":"Provides the Fornax cluster mass model used to define its tidal field profile."}],"review_version":1}