{"id":"dcb31e67-8f56-4c3e-b8d6-ff8e1574d5ff","arxiv_id":"1908.08056","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Calcium-strong transients may originate from rare binary stars formed dynamically in globular clusters and ejected before exploding, matching their observed wide offsets from host galaxies.","lead":"This paper argues that calcium-strong transients, a faint and puzzling class of stellar explosions, may be born in globular clusters and then kicked out before they explode. The proposal matches the explosions' very large distances from their host galaxies and predicts a hidden population of such blasts near galaxy centers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on an unquantified GC dynamical-formation rate, and the paper's own §3.4 mass budget already demands an implausibly high efficiency (~1 CaST per 8.4 M_sun of GC mass), so the rate assumption is the load-bearing weakness.","rationale":"The reader's weakest assumption is exactly the one I consider load-bearing. The paper is honest about the gap, but the dynamical birth rate is what separates a causal GC scenario from a radial-distribution coincidence. Section 3.4's efficiency estimate is the strongest internal evidence that the scenario is strained: requiring nearly half of all GC WDs to participate in CaSTs over a Hubble time, or 3% within a 1 Gyr window, is a severe demand for channels that are rare in the field. The proposed mass-segregation boost of 2-7x is real but partial, and dense core-collapsed clusters are a minority of the GC population, so the effective denominator is smaller than the total GC mass. The ejection/hardening argument is physically plausible, and the radial-distribution match is suggestive but not unique, since old metal-poor stars match as well. I would not move to reject because the authors explicitly flag the rate calculation as future work and keep the alternative old-metal-poor-star channel open; a conditional acceptance that requires the rate calculation is the right level. If a future N-body/population-synthesis calculation does not reproduce even ~1/10 of the required efficiency, the GC explanation should be abandoned in favor of the alternative the paper itself preserves.","tokens_in":20266,"tokens_out":14787,"duration_ms":146626,"concrete_test":"Run a cluster Monte Carlo / N-body + binary population synthesis calculation for representative GCs spanning core densities n~10^4-10^7 pc^-3 and masses 10^5-10^6 M_sun (e.g., with CMC or MOCCA), evolving each cluster for a Hubble time and counting He+O/Ne WD mergers, He WD+NS disruptions, and He-burning-star+WD interactions that occur after binary-hardening ejection. Convert the yield to a mass-normalized rate and compare to the required η ≈ 1 event per 8.4 M_sun (or 1 per 120 M_sun for the 1 Gyr-delay variant). If the simulated yield is an order of magnitude or more below this benchmark, the GC production scenario cannot supply the observed CaST rate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (three rare binary channels formed dynamically in GCs, then hardened and ejected) requires that GC-specific dynamical formation rates supply the observed volumetric CaST rate. This is the one assumption the paper explicitly refuses to quantify: §3.3 states that the authors 'proceed under the assumption that CaST progenitors are formed dynamically in GCs as hard binaries at a rate larger than, or at least competitive with, that in the field' and defers the calculation to future N-body/binary population synthesis work. The internal rate check in §3.4 makes the burden concrete: an event rate of 1.21e-5 yr^-1 Mpc^-3 over a Hubble time from a GC mass density of 1.4e6 M_sun Mpc^-3 requires η ≈ 1 CaST per 8.4 M_sun of GC mass, i.e., ~40% of all GC WDs must end up in CaSTs (0.3 WDs per M_sun). Even the favorable 'only last 1 Gyr' variant requires ~1/120 M_sun and ~3% of GC WDs, before multiplying by the 2-7x mass-segregation correction and by the fact that only a subset of GCs have core densities near the fiducial n=1e7 pc^-3 needed for rapid hardening. Since the field rate for the He WD+NS channel is already several orders of magnitude too low before GC enhancement, the needed enhancement is not a small factor. A measured dynamical birth rate that falls short would leave only a radial-distribution coincidence, not a progenitor solution.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies calcium-strong transients (CaSTs), a class of faint, spectroscopically peculiar explosions with high nebular Ca/O ratios, large galactocentric offsets, and a relatively high volumetric rate. The authors first show that the CaST radial distribution is qualitatively similar to the radial distribution of globular clusters (GCs) and of old, metal-poor stars. They then propose that CaSTs originate from one of three binary channels—He+O/Ne white dwarf mergers, He white dwarf disruptions by neutron stars, or accretion from low-mass He-burning stars onto white dwarfs—that are rare in the field but can be formed dynamically in GCs. Binary hardening both raises the interaction rate and ejects the binaries before explosion, explaining the absence of GCs at transient sites. The paper also quantifies the required CaST production efficiency from GC mass and derives a falsifiable prediction of a comparable CaST rate in nuclear star clusters. The central claim is conditional on an assumed, uncalculated dynamical formation rate, and the paper is transparent about this limitation.","tokens_in":20653,"tokens_out":4402,"duration_ms":47764,"significance":"If the GC-production scenario is correct, the paper would connect CaSTs to dynamical binary formation in dense stellar systems and provide a new, testable prediction for nuclear star clusters. The paper makes several useful contributions: it assembles the observational constraints on CaSTs, derives a theoretical GC radial distribution from independent empirical scaling relations, identifies three specific progenitor channels that warrant further study, and explicitly quantifies the rate-efficiency tension that any GC scenario must overcome. The falsifiable NSC prediction and the honesty about the missing dynamical calculation are strengths. The main significance is as a hypothesis-shaping paper rather than a demonstrated progenitor model; the central claim rests on an assumption that the authors themselves defer to future work.","major_comments":[{"comment":"The central assumption of the GC scenario is not quantified. The paper states that \"we proceed under the assumption that CaST progenitors are formed dynamically in GCs as hard binaries at a rate larger than, or at least competitive with, that in the field\" and defers the calculation to future N-body and population-synthesis work. This is load-bearing because §2.4 quotes field WD+NS merger rates of only 3e-16 to 3e-15 yr^-1 Msun^-1, so the required GC enhancement is not a small factor. As written, the paper demonstrates that the GC scenario is consistent with the radial distribution only if an uncalculated formation rate happens to be favorable; it does not demonstrate that the proposed channels supply the observed CaST rate. A quantitative requirement on the dynamical birth rate for each channel, or a reframing of the claim as a strictly conditional hypothesis, is needed.","section":"Section 3.3"},{"comment":"The rate-efficiency budget derived in §3.4 is severe and, as the authors acknowledge, not resolved by the paper. The constant-Hubble-time case requires η ≈ 1 CaST per 8.4 Msun of GC mass, which with 0.3 WDs per Msun implies roughly half of all GC white dwarfs participating; the 1-Gyr-delay case gives about 3% before applying the 2–7× mass-segregation boost and the restriction to the subset of GCs with densities near 10^7 pc^-3 needed for rapid hardening. Since the hardening timescale in Eq. (3) is evaluated at n = 10^7 pc^-3, typical GCs will be much less efficient, so the required per-cluster efficiency is even higher. This budget should be converted into a required dynamical formation rate for each candidate channel and compared with existing constraints from cataclysmic variables, X-ray binaries, and double-WD populations, rather than left as a caveat.","section":"Section 3.4, Eq. (6)"},{"comment":"The claimed consistency between the GC radial distribution and the CaST distribution is assessed only by visual comparison of cumulative distributions. No significance test (e.g., Kolmogorov–Smirnov) is given, the gold sample contains only eight objects, and the paper itself notes in §3.1 that the Palomar Transient Factory is biased against recovering CaSTs close to their hosts, making the gold distribution an upper limit. With these issues, the visual agreement is suggestive but not a quantitative empirical constraint. The paper should either provide a test statistic or explicitly label the match as qualitative for both Figure 1 and Figure 2.","section":"Section 3.1, Figure 1"}],"minor_comments":[{"comment":"The uncertainty in the GC system effective radius–halo mass relation, including the alternative slope from Hudson & Robison (2018) noted in the footnote, is not propagated into the cumulative distributions in Figures 1 and 2; showing a range of Sérsic indices and scaling-relation parameters would make the robustness of the conclusion clearer.","section":"Section 3.1"},{"comment":"The quoted Frohmaier et al. rate is asymmetric (+1.13/-0.39), but the derived efficiency η is presented as a single value; reporting the corresponding range in η would help the reader see how sensitive the mass-budget argument is to the uncertain rate.","section":"Section 3.4"},{"comment":"The statement that overall TNG100 stellar density profiles match observed stacked galaxies is used to support the reliability of age and metallicity binning, but this is a weaker validation than a direct comparison of halo stellar populations; the sentence should be phrased as a modeling assumption.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"This is an interesting and readable hypothesis paper, but the main claim is not yet demonstrated: the dynamical formation rate is assumed rather than computed, and the paper's own rate budget shows a strong tension. I would support publication after a revision that either quantifies the required formation rates and compares them to existing dynamical constraints, or explicitly repositions the paper as a constraints-and-predictions study rather than a progenitor solution. The NSC prediction is a genuinely useful falsifiable corollary and should be retained."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper gives the CaST community the right target. It shows quantitatively, using the halo mass function, the M_GCs/M_halo relation, and the Forbes (2017) size relation, that the projected radial distribution of globular clusters in the local universe matches the observed CaST distribution. It then proposes a concrete mechanism: rare WD binaries formed dynamically in GCs, hardened by three-body encounters, ejected before interaction, and exploding as CaSTs. That combination is genuinely new. The old-metal-poor-star matching was already in Yuan et al. (2013) and Perets (2014), and the three binary channels come from earlier work. The new pieces are the GC-distribution construction and the hardening/ejection framing, plus the nuclear star cluster prediction.\n\nThe paper is honest about its limits. Section 3.3 explicitly states that the dynamical birth rate of the proposed binaries is assumed, not computed. Section 3.4 shows the burden is stiff: constant production over a Hubble time requires about one CaST per 8.4 solar masses of GC mass, which translates to roughly half of all cluster white dwarfs. Even the more favorable 1 Gyr variant needs about 3% of GC WDs before applying the 2-7x mass segregation correction, and it only works for clusters dense enough to harden binaries. The authors themselves rule out the He WD + NS channel because too few neutron stars are retained in GCs. So the central claim is not yet demonstrated. It is a plausibility argument with a clear falsifiable prediction (CaSTs in nuclear star clusters) and a clear required calculation (N-body plus binary population synthesis formation rates). That is the right way to frame it.\n\nWhere are the soft spots? The radial distribution comparison is visual; there is no significance test for matching 8 or 14 objects to the model curves, and the gold sample is acknowledged to be biased. That is a real but minor issue, because the silver sample and the host galaxy properties line up reasonably. The NSC prediction is interesting but not yet testable, given how hard it is to recover CaSTs in galaxy cores, and no quantitative expected rate is given beyond \"similar rate.\"\n\nWho is this for? Observers working on Ca-rich gap transients and theorists doing binary population synthesis in dense stellar systems. It deserves a serious referee. I would send it out, and I would push the authors to quantify the dynamical formation rate in future work, but I would not require that calculation as a condition for publishing this scenario.","headline":"This paper makes a real contribution by constructing the local-universe globular cluster radial distribution and arguing that dynamical formation plus hardening/ejection in GCs can explain the CaST population, but the central formation-rate assumption is still unquantified and the scenario is a plausible framework, not a demonstrated solution.","tokens_in":21224,"tokens_out":2333,"would_cite":true,"duration_ms":24009,"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":"This paper argues that calcium-strong transients come from rare binaries forged and then ejected by globular clusters.","keywords":["calcium-strong transients","calcium-rich supernovae","globular clusters","white dwarf mergers","tidal disruption","binary hardening","nuclear star clusters","transient astronomy"],"falsifier":"A search capable of finding faint transients in the bright cores of galaxies that finds no nuclear-star-cluster CaST population at a rate comparable to the observed outskirt rate would falsify the dynamical-production scenario. A direct dynamical calculation showing that the three channels form in globular clusters far too slowly to meet the required mass budget would do the same.","tokens_in":20003,"feed_emoji":"💥","tokens_out":9188,"duration_ms":83164,"temperature":0.7,"pith_summary":"Calcium-strong transients (CaSTs) are faint, fast explosions with unusually high calcium-to-oxygen nebular line ratios, and they appear tens of kiloparsecs from their host galaxies far more often than ordinary supernovae. This paper argues that their extreme locations and relatively high rate point to a production site inside globular clusters: rare binaries -- helium white dwarfs merging with oxygen/neon white dwarfs, helium white dwarfs tidally disrupted by neutron stars, or low-mass helium-burning stars dumping material onto white dwarfs -- are formed dynamically in dense clusters and then ejected by binary-hardening recoil before they explode. The ejected binaries inherit the globular cluster radial distribution, explaining both the large offsets and the absence of any cluster at the explosion site. The paper shows that old, metal-poor stars also match the distribution, but no concrete progenitor scenario exists for that alternative. If the cluster scenario is correct, CaSTs should also be produced near the centers of galaxies in nuclear star clusters, a prediction that future surveys can test.","feed_headline":"Calcium-strong transients may be born in globular clusters","feed_subtitle":"Rare binaries ejected from dense star clusters would explain their huge offsets and predict flares near galaxy centers.","key_machinery":"The load-bearing object is the globular cluster itself, treated as a dynamical factory: its high stellar density creates rare hard binaries through captures and exchanges, and the same environment then hardens them. The quantitative machinery is the balance between the binary-hardening recoil and the cluster's gravitational pull. Repeated encounters shrink a hard binary's orbit and give the binary a growing recoil velocity; once that velocity exceeds the cluster escape speed, the binary is ejected while its orbit is still wide enough that gravitational-wave inspiral has not yet brought the components into contact. The paper estimates that the ejection separation and the gravitational-wave-dominated separation are comparable for its three candidate channels, so ejection before interaction is plausible for at least some systems. The other load-bearing component is the mass budget: combining the observed volumetric CaST rate with the local mass density of globular clusters requires roughly one CaST per 8.4 solar masses of cluster mass if production has been steady for a Hubble time, or one per 120 solar masses if it has been active for only the last gigayear.","core_discovery":"The central claim is that the observed population of calcium-strong transients is best understood as the product of dynamically formed binaries that are ejected from globular clusters before interacting. The paper matches the projected galactocentric radial distribution of the observed CaSTs to the theoretical radial distribution of globular clusters in the local universe, built from the halo mass function and a power-law relation between cluster system size and halo mass. It then identifies three binary channels whose field rates are too low but whose formation and interaction rates could be enhanced in dense cluster environments: He plus O/Ne white dwarf mergers, tidal disruption of He white dwarfs by neutron stars, and stable accretion from low-mass He-burning stars onto white dwarfs. Binary hardening simultaneously shrinks the orbits, raising the interaction rate, and eventually gives recoil kicks that exceed the cluster escape speed, ejecting the binary before mass transfer or merger. That is why no globular cluster is seen at a CaST site despite the cluster origin. The paper explicitly flags that the required production efficiency is high: if CaSTs have been produced at a constant rate for a Hubble time, nearly half of all white dwarfs in globular clusters must participate, and this budget appears to rule out the neutron-star channel because few neutron stars are retained in clusters.","pith_inferences":["If the globular cluster origin is right, the CaST rate should scale with the total globular cluster mass of a galaxy or halo, so more massive clustered environments should contribute disproportionately; this could be tested by comparing host halo masses of future CaST samples.","The tight mass budget hints that globular clusters may have been more massive in the past or that only a subpopulation of dense, core-collapsed clusters produces most CaSTs; identifying which clusters contribute could sharpen predictions for the nuclear-star-cluster rate.","The same ejection mechanism may apply to other compact-object transients with anomalously extended radial distributions, so a systematic search for hostless explosions tracing globular-cluster halos could reveal whether CaSTs are one example of a broader dynamical channel.","A direct N-body calculation of the three binary channels in realistic cluster models would turn the currently assumed formation-rate enhancement into a measured rate; if it comes out below field rates, the old-metal-poor-star alternative would need a concrete progenitor to remain viable."],"forward_implications":["CaSTs should continue to appear far from their hosts, tracing the globular cluster radial distribution; the observed sample should converge toward the broader silver-sample distribution rather than the gold-sample one as surveys improve.","Deep imaging at CaST sites should usually reveal no host globular cluster, because the progenitors are ejected before exploding, though events born in clusters with unusually high escape velocities may still be found inside clusters.","Nuclear star clusters should produce CaSTs at a rate comparable to the observed outskirt rate, and the current lack of such detections is a selection effect that core-sensitive searches can test.","The helium-white-dwarf plus neutron-star channel is unlikely to be the main route, because the number of neutron stars retained in globular clusters is too small to provide the required event rate.","Explosion models of the remaining two channels should produce low radioactive yields and ejecta dominated by intermediate-mass elements, matching the faint, calcium-dominated spectra of CaSTs."],"supporting_citations":[{"why":"Defines the prototype SN 2005E and proposes helium-shell detonation as the explosion mechanism the paper builds on.","marker":"Perets et al. 2010"},{"why":"Establishes the calcium-rich gap transient class and its photometric, spectroscopic, and host-galaxy properties.","marker":"Kasliwal et al. 2012"},{"why":"Provides the volumetric CaST rate that sets the required production efficiency and quantifies detection incompleteness in galaxy cores.","marker":"Frohmaier et al. 2018"},{"why":"Reports deep searches finding no globular cluster at CaST sites, motivating the ejection mechanism.","marker":"Lunnan et al. 2017"},{"why":"First suggested that globular cluster or old-metal-poor stellar distributions match the CaST radial distribution.","marker":"Yuan et al. 2013"},{"why":"Supplies the relation between globular cluster system size and halo mass used to build the predicted GC radial distribution.","marker":"Forbes 2017"},{"why":"Provides the halo mass function used to weight globular cluster systems by halo mass.","marker":"Tinker et al. 2008"},{"why":"Gives the recoil-ejection criterion for hardened binaries that underlies the ejection argument.","marker":"Sigurdsson & Phinney 1993"},{"why":"Quantifies neutron star retention in globular clusters, used to reject the He WD plus NS channel.","marker":"Ivanova et al. 2008"},{"why":"Shows the field WD-plus-neutron-star merger rate is far too low, motivating dynamical formation in clusters.","marker":"Toonen et al. 2018"}],"fun_headline_variants":["Calcium-strong transients may hail from globular clusters","CaSTs likely ejected from dense clusters before exploding","Globular clusters may be the birthplaces of CaSTs","Cluster-spawned binaries may explain calcium transients","CaSTs produced in globular clusters, then flung out"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The assumption that the three candidate binaries are actually formed dynamically inside globular clusters at rates at least competitive with field formation is not quantified in the paper, and without it the radial-distribution match is only a coincidence.","fun_headline_variants_meta":{"raw":{"variants":["Calcium-strong transients may hail from globular clusters","CaSTs likely ejected from dense clusters before exploding","Globular clusters may be the birthplaces of CaSTs","Cluster-spawned binaries may explain calcium transients","CaSTs produced in globular clusters, then flung out"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000975,"raw_usage":{"total_tokens":4222,"prompt_tokens":1101,"completion_tokens":3121,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":717,"completion_tokens_details":{"reasoning_tokens":3038}},"tokens_in":717,"tokens_out":3121,"duration_ms":24422,"temperature":1.0,"reasoning_tokens":3038,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:51:10.395257+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A search capable of finding faint transients in the bright cores of galaxies that finds no nuclear-star-cluster CaST population at a rate comparable to the observed outskirt rate would falsify the dynamical-production scenario. A direct dynamical calculation showing that the three channels form in globular clusters far too slowly to meet the required mass budget would do the same.","supporting_citations":[{"cited_title":"P., et al","cited_arxiv_id":null,"evidence_quote":"First suggested that globular cluster or old-metal-poor stellar distributions match the CaST radial distribution."}],"review_version":1}