{"id":"7bd16ac6-ffba-4ba5-9ed5-875a8b2ede22","arxiv_id":"2505.03085","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The paper proposes that Kepler's supernova progenitor may come from a merged satellite galaxy and estimates that alien Type Ia supernovae occur at roughly 2 to 5 per 100,000 years.","lead":"This paper asks whether some of the Milky Way's Type Ia supernovae come from stars stolen from smaller galaxies, and argues that Kepler's supernova of 1604 is a candidate. It estimates that such 'alien' SNe Ia happen every 20,000 to 60,000 years, a small but real fraction of all explosions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Alien-candidate claim rests on unresolved knot proper-motion systematics: the two published catalogs differ by ~3.3σ and flip the substructure association, so 'extragalactic origin' is not yet secure.","rationale":"The reader's weakest-assumption pick is the correct one. Section 2.1 assumes the optical knots are dense nitrogen-rich circumstellar material moving with the progenitor; the two proper-motion catalogs disagree at about 3.3σ, and the paper itself says the difference likely comes from using different knot samples. This is a systematic effect, not a random error, so displaying both sets side by side in Figures 1–3 overstates the precision of Kepler's (E, Lz). The comparison to Malhan et al. accretion boxes changes qualitatively: vK77 places Kepler inside several substructure boxes, while Bv91 removes it from all of them. The broader statement that Kepler is not in the disk/bulge in-situ boxes would also apply to an ordinary halo star; it is not a specific accretion signature. A clean expanding-shell fit to the knots would settle whether the assumed systemic velocity is correct. I do not see a comparable internal problem in the rate calculation. The DFA approximation is explicitly derived, the mass-loss constant 0.4 is justified in Appendix A, and the paper acknowledges the estimate is a lower bound. The two methods share inputs and are thus not fully independent, but the agreement mostly checks the DFA algebra rather than the physical rates; the dominant uncertainty is the adopted Macc, which the paper states. Therefore the rate component supports a modest 'small but possibly detectable' conclusion. The Kepler candidate, however, is not yet secure because of the knot proper-motion systematics. The paper's own Conclusion (Section 5, item 1) concedes the substructure association is inconclusive, so the authors are appropriately careful. Since the reader's CONDITIONAL verdict already reflects this, I recommend no change to the verdict.","tokens_in":21041,"tokens_out":12417,"duration_ms":143983,"concrete_test":"Fit an expanding-shell model to the 50 knots of Bandiera & van den Bergh (1991) using the radial velocities from Blair et al. (1991) plus both proper-motion sets, and solve for the systemic proper motion as a free parameter; if the fitted systemic PM differs from both vK77 and Bv91 by more than the quoted errors, the assumption that the knot motions trace the progenitor is false and the kinematic-based alien classification is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is Section 2.1's assumption that the optical knots' proper motions trace the progenitor's space motion. The two published catalogs, vK77 (pmDE = 10.9 ± 1.8 mas/yr) and Bv91 (pmDE = 4.84 ± 0.49 mas/yr), disagree at about 3.3σ, and the paper itself attributes the discrepancy to different knot samples. This is a systematic, not merely statistical, uncertainty. In Section 3.2 the inferred (E, Lz) position shifts enough that vK77 overlaps Sagittarius, Cetus, G-E, L/W, and Helmi, while Bv91 overlaps none of the accretion boxes; Section 5 explicitly concedes that the substructure association is inconclusive. The weaker claim 'not in the Malhan in-situ boxes' would also be true for an ordinary high-energy halo star, so it does not by itself establish an accreted origin. If the knots are ejecta or shocked CSM rather than progenitor material, the Kepler-based part of the paper loses its observational anchor. The alien-SNe-Ia rate estimate is a separate, transparent lower-bound calculation and is not affected by this concern.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces the notion of 'alien SNe Ia' produced by stars that were accreted into the Milky Way through dry mergers, and it presents two main results. First, using Gaia DR3 data for stars around Kepler's supernova remnant, the authors compare the progenitor's kinematics and dynamics in (E, Lz) and action space with nearby stars and with known substructures, concluding that the progenitor is kinematically and dynamically anomalous and may not belong to the in-situ Milky Way population, thus proposing it as a possible alien SN Ia candidate. Second, the paper estimates the recent rate of alien SNe Ia from seven known merger remnants (Kraken, Gaia-Enceladus-Sausage, Helmi streams, Sequoia, Sagittarius, Wukong/LMS-1, and Cetus) using two methods: a galactic chemical evolution calculation with GalCEM, and a newly introduced delta-function approximation (DFA). The GalCEM method gives 1.5e-5 to 5.0e-5 yr^-1 (0.9 to 3.0 events in the past 60 kyr), while the DFA gives 3.1e-5 yr^-1 with 1.84 events; both are honestly labeled lower bounds under the assumption of no post-merger star formation.","tokens_in":21278,"tokens_out":4992,"duration_ms":58245,"significance":"If the Kepler association were secure, this would be the first identified candidate for a supernova progenitor of extragalactic origin in the Milky Way, and the paper would open a new observational window on accreted stellar populations. The rate calculation is a transparent forward estimate that uses literature values for Macc and zacc, a plausible DTD family, and public, reproducible code; the lower-bound character is clearly stated. The new DFA method is clearly derived and its approximations are enumerated, which is useful even if the method is approximate. However, the Kepler candidate claim rests on the assumption that optical knots trace the progenitor's space motion, and the two available proper-motion catalogs are mutually inconsistent at about 3.3 sigma, flipping the substructure association. The rate methods are also not fully independent because they share the same satellite masses, accretion times, and DTD inputs, so their consistency should not be overinterpreted as independent validation.","major_comments":[{"comment":"The kinematic case for Kepler's progenitor being an alien SN Ia depends entirely on the assumption that the optical knots' proper motions trace the progenitor's space motion. The two published measurements disagree strongly: van den Bergh & Kamper (1977) give pmDE = 10.9 +/- 1.8 mas/yr while Bandiera & van den Bergh (1991) give pmDE = 4.84 +/- 0.49 mas/yr, a difference of roughly 3.3 sigma. As the paper itself shows in Section 3.2, adopting the vK77 values places Kepler's progenitor in or near several accretion-event boxes (Sagittarius, Cetus, G-E, L/W, Helmi), whereas adopting the Bv91 values places it in none of the accretion boxes. Since the knots could in principle be shocked ejecta or circumstellar material rather than unshocked progenitor material, the conclusion that the progenitor is not in-situ is not secured. The paper should either provide independent evidence for the knot-progenitor connection or explicitly frame the entire Kepler analysis as conditional on a still-unresolved systematic choice, with the abstract and title adjusted accordingly.","section":"Section 2.1 and Section 3.2"},{"comment":"The claim that Kepler's progenitor is 'not in situ' is weakened by the lack of a statistical control sample of halo stars. Showing that the progenitor lies outside the Malhan et al. (2022) in-situ boxes is not by itself evidence for an accreted origin, because an ordinary high-energy halo star with a large vertical excursion would also fall outside those boxes. The nearby-star comparison in Figure 1 is dominated by disk stars, so a clean separation from them is expected for any halo-like orbit. To make the 'alien' classification load-bearing, the authors should quantify how often field halo stars occupy the same (E, Lz) region, for example by drawing a matched control sample from Gaia and computing the probability that an in-situ star would appear as anomalous as Kepler's progenitor under the same selection and orbit-integration pipeline.","section":"Section 3.2 and Figure 1"},{"comment":"The statement that 'Methods I and II are relatively independent, so we consider our results to be reliable' is not justified. Method II uses the same Macc and zacc values for the same seven satellites and adopts the same power-law DTD family from Maoz & Graur (2017) that is used in one branch of Method I; both methods also assume dry mergers with no post-merger star formation. Given these shared inputs, the consistency of the two rate estimates is largely a check on the algebraic approximations of the DFA, not an independent confirmation of the astrophysical result. The paper should soften this claim and present the two methods as complementary estimates sharing common systematic uncertainties, with the dominant uncertainty being the satellite masses and DTD normalization.","section":"Section 4.2"}],"minor_comments":[{"comment":"There are several typographical and formatting issues, including 'Kranken' in the discussion of Malhan et al. (2022) (should be 'Kraken'), the inline notation 'd4.5' in the Introduction, and inconsistent use of 'top/bottom panels' in the Figure 3 caption when the figure appears to be arranged in rows; these should be cleaned up before publication.","section":"Throughout"},{"comment":"In the DFA derivation, the choice of the harmonic mean for DTD(Delta t_har) is reasonable, but the associated error budget in Appendix B would benefit from being stated more concretely; for instance, the claim that the '1 sigma dynamic range of Macc is about 5' should be reconciled with the quoted 0.3 dex systematic uncertainty, which corresponds to a factor of about 2, so that the reader can verify the relative sizes of the systematic terms.","section":"Section 4.1.2 and Appendix B"},{"comment":"The discussion correctly acknowledges that the substructure association is inconclusive, but this caveat appears only after the abstract and Section 3.2 have already stated the 'alien' hypothesis in stronger terms; the authors should ensure the abstract matches the level of certainty expressed in the conclusions.","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The paper has a solid, reproducible rate calculation that is clearly labeled as a lower bound, and the DFA derivation is a useful methodological contribution. The main weakness is that the Kepler candidate claim, which is central to the paper's framing and title, rests on unresolved proper-motion systematics and on an assumption about the nature of the optical knots. This is fixable by reframing the Kepler result as explicitly conditional and by adding a control-sample analysis, but it requires substantive revision rather than copy-editing. The rate section could stand as a separate contribution if the Kepler claim were removed, but as written the two parts should be brought to a consistent level of uncertainty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a decent paper that introduces a genuinely new idea—using the known Milky Way merger inventory to estimate a present-day rate of SNe Ia from accreted stars—and it does that part carefully. The Kepler candidate is the weaker half, and the paper knows it. The two published proper-motion catalogs disagree at about 3.3 sigma, and the substructure association flips depending on which you use. The authors are upfront that the association is inconclusive. So the Kepler part should be read as a well-hedged candidate, not a detection.\n\nWhat's actually new: combining the Kruijssen et al. / Naidu et al. accretion masses and epochs with SN Ia DTDs to get an 'alien' rate. The delta-function approximation (Eq. 8) is a compact closed-form tool, and the error budget is thoughtful—they show the main uncertainty is Macc, not the DTD or mass-loss details. The rate estimates, 1.5e-5 to 5.0e-5 per year (0.9–3.0 events per 60 kyr), are transparent forward calculations and honestly labeled lower bounds. No fitting to the signal; no circularity.\n\nSoft spots. First, the Kepler proper-motion systematics. vK77 and Bv91 disagree in pmDE by ~3.3 sigma, and the overlap with accreted substructures changes qualitatively. The paper says new astrometry is needed, which is correct, but that means the 'alien' label on Kepler is not secure. Also, being outside the Malhan in-situ boxes does not by itself prove an accreted origin; a high-energy halo star would look similar. Second, the two rate methods share the same Macc, zacc, and DTD family, so they are not fully independent. The paper calls them 'relatively independent,' which oversells a little—they differ in how SFH is handled, not in the underlying inputs. That doesn't break the rate estimate, but the agreement between methods is less impressive than it looks. Third, no code or data shipped; the paper relies on public code (GalCEM, Galpy) and literature tables. Not a dealbreaker, but releasing the orbit-integration and rate-calculation scripts would help reproducibility.\n\nFor whom: people working on SN Ia rates, Galactic archaeology, or SNR demographics. It defines a new population at a quantifiable rate, and the DFA formula may be useful elsewhere. It deserves a serious referee; I would send it to review, expecting the Kepler claims to be softened further and the method-independence caveat addressed. I'd probably cite the rate estimate and the DFA approximation but treat the Kepler association as tentative.","headline":"Solid, honest rate estimate for 'alien' SNe Ia from known mergers; the Kepler candidate is suggestive but the proper-motion systematics keep it from being more than a candidate.","tokens_in":21867,"tokens_out":2731,"would_cite":true,"duration_ms":27454,"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 Kepler's supernova had an extragalactic progenitor, and that 'alien' Type Ia supernovae from Milky Way mergers occur at a detectable rate of 1.5e-5 to 5e-5 per year.","keywords":["Type Ia supernovae","supernova remnants","Milky Way mergers","galactic chemical evolution","stellar kinematics","Kepler's supernova","Gaia DR3","delay-time distribution"],"falsifier":"Measure the proper motion of Kepler's optical knots with sub-milliarcsecond precision using astrometric monitoring over a decade; if the resulting space motion places the progenitor inside the (E, Lz) region that Malhan et al. (2022) assign to in-situ disk or bulge populations, the claim that Kepler is an alien SN Ia would be falsified. Detecting a surviving companion star with velocity consistent with in-situ populations would likewise contradict the accreted-origin scenario.","tokens_in":20789,"feed_emoji":"💥","tokens_out":11591,"duration_ms":96581,"temperature":0.7,"pith_summary":"This paper argues that some Type Ia supernovae in the Milky Way originate from stars that were born in dwarf galaxies later swallowed by our own, and it names these events 'alien SNe Ia.' Using Gaia DR3 data, the authors show that the progenitor of Kepler's supernova (SN 1604) has higher total energy and different angular momentum than its neighboring stars and known in-situ populations, suggesting it may be an immigrant from a disrupted satellite. They then estimate how often such alien supernovae occur today by combining the accretion times and stellar masses of seven known mergers with two independent methods, obtaining a recent rate of roughly $1.5\\times10^{-5}$ to $5.0\\times10^{-5}$ per year. If correct, this means roughly one to three of the supernova remnants currently visible in the Milky Way should be 'alien,' offering a new way to read the Galaxy's merger history from its explosive relics.","feed_headline":"Kepler's supernova may be from a swallowed galaxy","feed_subtitle":"The paper counts one to three 'alien' Type Ia supernovae exploding in the Milky Way in the last 60,000 years.","key_machinery":"The analysis runs on two pieces of machinery. First, orbit integration in the Galpy package, using two Milky Way potentials (MWPotential2014 and McMillan17), converts the position and velocity of Kepler's progenitor and of 3,507 nearby Gaia DR3 stars into integrals of motion: total energy $E$, angular momentum $L_z$, and actions $(j_r,j_p,j_z)$. Comparing these against the $(E,L_z)$ boxes that Malhan et al. (2022) assign to in-situ and accreted substructures is what produces the 'alien' kinematic signature. Second, the rate estimate combines the accretion time $z_{\\rm acc}$ and stellar mass $M_{\\rm acc}$ of each of the seven satellites with the delay-time distribution (DTD) of SNe Ia: method I uses the GalCEM one-zone chemical evolution code with two DTDs (Greggio 2005; Maoz & Graur 2017), while method II, introduced in this paper, is the delta-function approximation, which replaces the unknown star-formation history with a single 'harmonic' delay time $\\Delta t_{\\rm har}$ and a constant stellar mass loss of 40%, yielding $N \\approx (M_{\\rm acc}/0.6)\\,{\\rm DTD}(\\Delta t_{\\rm har})\\,t_{\\rm SNR}$. This machinery lets the authors estimate the alien SN Ia rate using only the mass and accretion time of each merger.","core_discovery":"The paper's central claim is that Kepler's supernova progenitor is kinematically and dynamically distinct from the in-situ Milky Way stellar population: in $(E,L_z)$ space it lies far from the disk and bulge regions identified by Malhan et al. (2022), and in action space it has larger radial and vertical actions than its surrounding stars. The authors interpret this as evidence that the progenitor was accreted into the Milky Way from a disrupted satellite galaxy, making SN 1604 a candidate 'alien' Type Ia supernova. They further claim that alien SNe Ia from seven known accreted galaxies (Kraken, Gaia-Enceladus-Sausage, the Helmi streams, Sequoia, Sagittarius, Wukong/LMS-1, and Cetus) should occur at a present-day rate of $1.5\\times 10^{-5}$ to $5.0\\times 10^{-5}$ yr$^{-1}$ by galactic chemical evolution modeling, or $3.1^{+1.8}_{-1.1}\\times 10^{-5}$ yr$^{-1}$ by a new delta-function approximation, corresponding to 0.9 to 3.0 events over the last 60 kyr. These are lower bounds because the models assume dry mergers with no post-merger star formation.","pith_inferences":["A direct test of the paper's logic would be to apply the same $(E,L_z)$ and action-space comparison to other historical Type Ia remnants, such as Tycho's SNR, whose in-situ association is usually assumed; an unexpected outlier there would suggest that alien progenitors are common enough to be found among the small set of historical supernovae.","If Kepler were ever confirmed as alien (for example, through the kinematics of a surviving companion star), its known age and the host merger's accretion time would pin down a single delay time for one SN Ia, providing an individual anchor for the delay-time distribution of accreted stellar populations.","The delta-function approximation could be turned into a forecasting tool for other galaxies: given a measured merger history from stellar halos, one could predict the present-day fraction of 'alien' SNe Ia in external galaxies and compare with resolved stellar-population studies.","A chemical follow-up prediction follows from the paper's own logic: an alien SNR should show nucleosynthetic or circumstellar abundance patterns reflecting the lower metallicity of its dwarf-galaxy origin, distinguishing it from in-situ remnants of similar age."],"forward_implications":["Roughly one to three of the supernova remnants visible today should be 'alien' within the canonical 60 kyr SNR lifetime, so remnant surveys can expect to find a few kinematic outliers like Kepler.","Kepler's unusual high-velocity, plane-escaping motion and its asymmetric shell are naturally explained if its progenitor was accreted, without invoking special binary evolution.","The delta-function approximation gives a quick estimate of recent SN Ia rates for any accreted galaxy from just its stellar mass and accretion time, usable beyond the Milky Way.","Because the estimates assume dry mergers with no post-merger star formation, the true alien SN Ia rate is likely higher, strengthening the case that a detectable fraction of SNe Ia are extragalactic in origin.","The Milky Way's merger history should leave observable traces not only in stars but also in the remnants of thermonuclear supernovae."],"supporting_citations":[{"why":"Supplies one of the two proper-motion measurements of Kepler's optical knots; with this value Kepler's progenitor has about a 1% escape probability.","marker":"van den Bergh & Kamper 1977"},{"why":"Supplies the other, more recent proper-motion measurement; with this value escape is nearly impossible and associations with substructures change.","marker":"Bandiera & van den Bergh 1991"},{"why":"Argues the dense, nitrogen-rich knots in Kepler trace the progenitor's motion and derives its high velocity above the Galactic plane.","marker":"Bandiera 1987"},{"why":"Provides the (E, Lz) ranges for in-situ populations and accreted substructures used to classify Kepler as non-in-situ.","marker":"Malhan et al. 2022"},{"why":"Supplies the accretion redshifts and stellar masses for five of the seven mergers considered.","marker":"Kruijssen et al. 2020"},{"why":"Supplies the truncation redshifts and stellar masses for Wukong/LMS-1 and Cetus.","marker":"Naidu et al. 2022"},{"why":"Provides the power-law delay-time distribution with slope -1.07 and normalization used in both rate methods.","marker":"Maoz & Graur 2017"},{"why":"Provides the fiducial single-degenerate delay-time distribution implemented in GalCEM for method I.","marker":"Greggio 2005"},{"why":"Provides the GalCEM code used to evolve each satellite galaxy and compute SN Ia rates in method I.","marker":"Gjergo et al. 2023"},{"why":"Supplies the 60 kyr typical radio lifetime of SNRs used to convert rates into expected event numbers.","marker":"Frail et al. 1994"}],"fun_headline_variants":["Kepler's supernova may be an alien from a merged galaxy","Milky Way mergers spawn alien supernovae","Kepler's supernova: an immigrant star","Alien supernovae from Milky Way's past mergers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the dense, nitrogen-rich optical knots in Kepler's remnant trace the motion of the progenitor star itself—not shocked ejecta or circumstellar material—and the two published proper-motion measurements of those knots disagree, so the 'alien' kinematic signature is fragile.","fun_headline_variants_meta":{"raw":{"variants":["Kepler's supernova may be an alien from a merged galaxy","Milky Way mergers spawn alien supernovae","Kepler's supernova: an immigrant star","Alien supernovae from Milky Way's past mergers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000938,"raw_usage":{"total_tokens":4103,"prompt_tokens":1132,"completion_tokens":2971,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":748,"completion_tokens_details":{"reasoning_tokens":2906}},"tokens_in":748,"tokens_out":2971,"duration_ms":18265,"temperature":1.0,"reasoning_tokens":2906,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:59:53.498399+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the proper motion of Kepler's optical knots with sub-milliarcsecond precision using astrometric monitoring over a decade; if the resulting space motion places the progenitor inside the (E, Lz) region that Malhan et al. (2022) assign to in-situ disk or bulge populations, the claim that Kepler is an alien SN Ia would be falsified. Detecting a surviving companion star with velocity consistent with in-situ populations would likewise contradict the accreted-origin scenario.","supporting_citations":[{"cited_title":"1991, ApJ, 374, 186, doi: 10.1086/170108","cited_arxiv_id":null,"evidence_quote":"Supplies the other, more recent proper-motion measurement; with this value escape is nearly impossible and associations with substructures change."},{"cited_title":"1987, ApJ, 319, 885, doi: 10.1086/165505","cited_arxiv_id":null,"evidence_quote":"Argues the dense, nitrogen-rich knots in Kepler trace the progenitor's motion and derives its high velocity above the Galactic plane."}],"review_version":1}