{"id":"5a5b4fe2-c3a5-4325-826b-542f92957490","arxiv_id":"2501.03330","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Neutrinos from hyperaccreting Thorne-Zytkow objects should be detectable by Super-Kamiokande, IceCube, JUNO, and Hyper-Kamiokande out to the Small Magellanic Cloud.","lead":"This paper calculates the neutrino signals that Thorne-Zytkow objects, stars with a neutron star buried inside a giant envelope, should emit if they are accreting matter at extreme rates. It finds that current and upcoming neutrino detectors could spot these signals from sources as far away as the Small Magellanic Cloud, giving astronomers a new way to confirm these objects.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Detectability claims hinge on a steady-state Ye=0.5 disk model whose inner neutrino opacity and electron fraction are not checked; a direct re-derivation is needed.","rationale":"The reader's weakest_assumption already identifies steady-state modeling and fixed Ye=0.5 as the weak spot; I partially agree. My concern sharpens the issue: the most dangerous single point is not disk versus spherical geometry per se, but the neutrino-opacity treatment at the inner boundary combined with constant Ye=0.5. The paper is transparent about its simplifications and labels the flux an upper limit, which is honest; however, it does not quantify how much the SMC reach would shrink if the inner disk is optically thick or if Ye drops to beta-equilibrium values. The proposed test is feasible with existing steady-disk codes and would settle whether the flux normalization is robust. I do not see a fatal internal inconsistency or circular fitting: the local TZO rate enters only in the diffuse-background section and is not fitted to the detection claim, and the detector response follows standard published treatments. The central claim is plausible but conditional on an unquantified modeling assumption, so CONDITIONAL is the appropriate verdict. A REJECT would be unwarranted because the paper's own caveats and the absence of obvious algebraic errors give it genuine value as a feasibility study motivating targeted searches; an ACCEPT would be premature because the headline distance contours are not bracketed by systematic uncertainties.","tokens_in":20341,"tokens_out":3389,"duration_ms":30931,"concrete_test":"Recompute the neutrino luminosity and spectrum for Mdot=10^6 Msun/yr with a two-zone opacity treatment: (1) evaluate the neutrino optical depth tau_nu(r) from r_NS outward using the Fig. 5 density/temperature profiles and standard charged-current opacities; (2) in Eq. (A9) replace the lower integration limit with the radius where tau_nu ~ 2/3; (3) repeat the emission calculation with Ye=0.5 replaced by self-consistent beta-equilibrium Ye(r) from a steady disk code (e.g., Popham et al. 1999, Di Matteo et al. 2002). If the 5 kpc flux changes by less than 30%, the SMC reach stands; if it drops by more than a factor of a few, the Fig. 3 burst contours and the VX Sgr/HV 2112/HV 11417 statements need revision.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central detectability claim (Super-K/IceCube reaching SMC for Mdot=10^6 Msun/yr in Fig. 3) rests on the accretion-disk model of Appendix A.1. Two load-bearing assumptions are unchecked. First, the self-similar inner-disk solution (Zhang & Dai 2008) is applied down to ~20 km and the neutron star surface at 10 km is treated as opaque, but the neutrino optical depth of the inner disk is not computed; if the region inside ~20 km is neutrino-opaque, the effective decoupling radius is larger than r_NS and Eq. (A9) overestimates the escaping flux. Second, beta-process emission is computed with fixed Ye=0.5 (App. A.3), although beta equilibration in neutrino-dominated accretion flows drives Ye to ~0.1-0.3, reducing the free nucleon abundance and charged-current rates; the authors explicitly flag that they do not account for the change in the electron fraction induced by beta processes. Because S/sqrt(S+B) scales linearly with flux, an order-of-magnitude overestimate would move the SMC burst contours down to Galactic distances. The paper correctly labels its flux an upper limit, but that label does not quantify the sensitivity of the headline distance claims.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes neutrinos as a new messenger for identifying Thorne–Żytkow objects (TŻOs). It models hypercritical accretion onto a neutron star in two steady-state regimes: a spherical settling solution for accretion rates 10^-4 to <10^4 M_sun/yr and an accretion-disk solution for 10^4 to 10^6 M_sun/yr. Thermal neutrino emission from beta processes, pair annihilation, and other channels is computed, with MSW flavor conversion applied, and event rates are estimated for IceCube, Super-Kamiokande, Hyper-Kamiokande, and JUNO. The central detection claims are that TŻO bursts with Mdot ~ 10^6 M_sun/yr could be seen up to the Small Magellanic Cloud by Super-Kamiokande and IceCube, that steady TŻOs with Mdot ≳ 0.1 M_sun/yr could be detected beyond the SMC after one year, and that existing DSNB limits already exclude the most extreme diffuse TŻO scenarios. The paper is explicitly framed as an upper-limit feasibility study.","tokens_in":20609,"tokens_out":9360,"duration_ms":89924,"significance":"If the estimates are robust, the paper would open a genuinely new observational window: neutrinos could reveal the neutron-star core inside TŻO candidates that are electromagnetically ambiguous, and the proposed searches are concrete and implementable with operating or near-term detectors. The study's strengths are that it uses independent literature source models rather than fitting a target signal, it provides detector-specific event rates and signal-to-noise contours, it makes falsifiable predictions for named candidates (VX Sgr, HV 2112, HV 11417), and it extends the analysis to the diffuse background with a direct comparison to DSNB constraints. However, the headline detectability claims are only as strong as the underlying steady-state emission models, and several load-bearing assumptions are not quantitatively checked; the paper acknowledges some of these but does not calibrate their effect on the distance reach.","major_comments":[{"comment":"The disk-scenario flux integral in Eq. (A9) is truncated at a 'decoupling radius' R_dec, but the manuscript never defines R_dec or computes the neutrino optical depth of the inner accretion disk. The model adopts the Zhang & Dai (2008) self-similar solution down to about 20 km and treats only the neutron-star surface as opaque, yet neutrino-dominated accretion flows are commonly optically thick inside tens of km. If the actual decoupling radius is larger than r_NS, Eq. (A9) overestimates the escaping flux, and because S/sqrt(S+B) scales linearly with flux, the SMC reach shown in Fig. 3 (left) could shrink to Galactic distances. Please provide R_dec from the adopted disk solution, or compute the neutrino optical depth explicitly and demonstrate that the radii dominating the emission integral are transparent.","section":"§A.1, §A.3, Eq. (A9)"},{"comment":"The electron fraction is fixed at Y_e = 0.5 and the beta-process Q-value is set to zero, with no accounting for the change in Y_e induced by beta equilibration. This is explicitly stated in Appendix A.3, but its quantitative impact is not assessed. In neutrino-dominated accretion flows, beta equilibrium typically drives Y_e to values around 0.1–0.3, which changes the nucleon abundances entering the charged-current emissivities and alters the IBD-relevant anti-nu_e flux. An order-of-magnitude overestimate would move the burst contours in Fig. 3 from the SMC down to Galactic distances. Please quantify the response of the emitted anti-nu_e spectrum and luminosity to a representative beta-equilibrated Y_e (or to a range of Y_e values), or otherwise show that the SMC detection claim is robust to this uncertainty.","section":"§A.1, §A.3"},{"comment":"For the disk/burst scenario, the manuscript states that the density profile beyond the disk is not modeled and that 'we cannot test if flavour evolution is adiabatic,' yet the same MSW description is then applied as for the spherical envelope. Because the IBD event rate depends directly on the anti-nu_e flux, this flavor-conversion assumption is load-bearing for Table 1 and Fig. 3 (left). Please estimate the range of possible anti-nu_e fluxes for the disk geometry, or at least state explicitly how the assumed conversion scheme changes the event rates relative to no conversion or to a maximal-mixing alternative.","section":"§A.4"}],"minor_comments":[{"comment":"There is an inconsistency in the accretion-rate threshold for beta-process dominance: §3 says beta processes dominate for Mdot > 10^4 M_sun/yr, while §A.3 says they dominate 'for disks accreting at rates Mdot ≤ 10^4 M_sun/yr'; the inequality in A.3 appears to be a typo and should be corrected.","section":"§3 vs §A.3"},{"comment":"The threshold accretion rates quoted for the named candidates are not internally consistent: VX Sgr is said to require Mdot ≳ 10^-2 M_sun/yr in §4.3 but Mdot ≳ 10^-3 M_sun/yr in §6, and for HV 2112/HV 11417 §4.3 quotes one-year Super-K and JUNO constraints of 2.4 and 3.2 M_sun/yr, respectively, while §6 states detection above 0.5 M_sun/yr; these numbers should be reconciled.","section":"§4.3 vs §6"},{"comment":"The horizontal axis label reads 'log10 Ye [g cm^-3]', but Y_e is a dimensionless electron fraction; the label should presumably read 'log10 rho [g cm^-3]' or 'log10 rho_e' to match the plotted quantity.","section":"Fig. 5, bottom panel"},{"comment":"The sentence explaining the two assumptions of the Chevalier model is self-contradictory: saying the postshock flow is adiabatic 'relies on the fact that neutrino energy losses are not negligible' is confusing, since adiabaticity normally requires the opposite; please clarify the intended meaning.","section":"§A.2"},{"comment":"The quantity r is first described as 'the sum of the hit count per DOM' and then used in Eq. (B12) as if it were a rate in the deadtime factor 1/(1+r tau); please define r unambiguously as a rate and state the units.","section":"§B.1, Eq. (B12)"},{"comment":"The caption contains a typo: 'Hyper-Kamionade' should be 'Hyper-Kamiokande'.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-written feasibility study, and the topic is well matched to the journal's astroparticle/heavy-ion scope. My main concern is that the two highest-impact distance claims rest on quantities (the inner-disk neutrino decoupling radius and the electron fraction) that are not computed but are explicitly acknowledged as simplifications; a careful revision that quantifies these effects would make the claims much more credible. I also note that the rate normalization for the local TŻO abundance relies in part on a paper coauthored by one of the present authors, and the DSNB background model uses the same group's earlier work; this is not improper, but the authors should ensure the cited local rate is not treated as more secure than the literature warrants."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this before reading: the paper is a genuine first — a systematic TZO neutrino emission and detection study with a clean split into burst and steady accretion scenarios, concrete sensitivity estimates for IceCube, Super-Kamiokande/Hyper-Kamiokande, and JUNO, and direct application to named candidates (VX Sgr, HV 2112, HV 11417). The detection claim is conditional, not a detection, but it deserves referee time.\n\nWhat is actually good: the authors did the legwork that earlier generic accretion-disk neutrino studies did not do for TZOs. Fluxes, durations, event rates, and signal-to-noise contours follow from stated steady-state models and standard cross sections; the paper is internally consistent. It also draws a sharp and useful distinction between transient and steady signals, gives tailored search strategies for each, and uses current DSNB limits to exclude extreme TZO accretion scenarios. The upper-limit framing is honest. The citation pattern is unremarkable; the TZO local rate is an uncertain external input, not a fitted target.\n\nSoft spots, in order. The stress-test note mostly lands. The fixed Ye = 0.5 assumption is explicit in Appendix A.3, and the authors admit they do not account for beta-driven changes in electron fraction; in neutrino-dominated accretion flows Ye is expected to drift toward 0.1–0.3, which can suppress charged-current emission. The decoupling radius Rdec appears in Eq. (A9) but no value or opacity check for the inner disk is provided; if the inner region near ~20 km is neutrino-opaque, the escaping flux is overestimated. Since the headline figures scale linearly with flux, an order-of-magnitude systematics could shrink the SMC contours to Galactic distances. The paper labels fluxes as upper limits, but that label does not quantify this sensitivity. The flavor-conversion treatment for transient sources is an extrapolation and is flagged as such. None of this kills the qualitative conclusion — VX Sgr at 1.5 kpc remains testable with large statistics — but the quantitative horizons in Fig. 3 are softer than the prose suggests.\n\nBottom line: worth a serious referee. The referee should ask for a check of the inner-disk neutrino optical depth and a short study of the Ye dependence, or at least explicit systematic bands on Fig. 3. I would read it, cite it, and mention it to people planning TZO searches.","headline":"First systematic TZO neutrino study, useful and worth refereeing, but the SMC reach claims rest on unquantified disk-model simplifications (fixed Ye = 0.5, uncomputed decoupling radius).","tokens_in":21129,"tokens_out":3258,"would_cite":true,"duration_ms":33923,"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":"Neutrinos can identify the hidden neutron-star cores of Thorne-Zytkow objects.","keywords":["Thorne-Zytkow objects","hypercritical accretion","neutrino emission","neutron stars","Super-Kamiokande","IceCube","diffuse neutrino background"],"falsifier":"A year-long targeted search with gadolinium-loaded Super-Kamiokande or JUNO for electron antineutrinos in the 9.3–31.3 MeV window from VX Sgr (at 1.5–1.7 kpc) that finds no excess above the modeled background would falsify the claim that VX Sgr is a steady TZO accreting at ≳ $10^{-2}$ solar masses per year.","tokens_in":20092,"feed_emoji":"🔭","tokens_out":6177,"duration_ms":49319,"temperature":0.7,"pith_summary":"Thorne-Zytkow objects—stars with a neutron star buried inside a giant envelope—have been predicted for decades but are hard to identify because their light resembles ordinary red supergiants. This paper argues that hypercritical accretion onto the neutron star heats the surrounding gas to 1–10 MeV, producing a copious flux of neutrinos in the 1–100 MeV range. It predicts that current and next-generation neutrino observatories (Super-Kamiokande, IceCube, JUNO, Hyper-Kamiokande) can detect these neutrinos from sources as far as the Small Magellanic Cloud, and that the existing TZO candidates VX Sgr, HV 2112, and HV 11417 can be tested with data already being taken. If detected, neutrinos would positively identify the neutron-star core and measure the accretion rate, something electromagnetic observations alone cannot do.","feed_headline":"Neutrinos could expose hidden neutron-star cores","feed_subtitle":"Current detectors could test Thorne-Zytkow candidates out to the Small Magellanic Cloud.","key_machinery":"The central machinery is a pair of steady-state accretion models stitched together by accretion rate: for Ṁ between $10^{-4}$ and $10^{4}$ solar masses per year, a spherically symmetric shocked envelope with power-law density ρ(r) ∼ $r^{-3}$ and temperature T(r) ∼ $r^{-1}$ (following Chevalier 1989); for Ṁ from $10^{4}$ to $10^{6}$ solar masses per year, a steady advection-dominated/neutrino-dominated accretion disk (following Di Matteo et al. 2002 and Zhang & Dai 2008) with a self-similar inner region attached to the neutron star. Neutrino emission is dominated by electron-positron pair annihilation in the spherical case and by $\\beta$ processes (p + e^- ⇌ νe + n and n + e^+ ⇌ ν̄e + p) in the disk case, with MSW resonant flavor conversion in the envelope. A simple energy-budget argument, τ ≈ GM_NS Δm/(r_NS Lν), sets the maximum signal duration and splits the phenomenology into bursts versus steady sources.","core_discovery":"The paper establishes that for accretion rates between $10^{-4}$ and $10^{6}$ solar masses per year, the region around the neutron star reaches temperatures of roughly 1 to 10 MeV, where $\\beta$ processes and electron-positron pair annihilation copiously produce neutrinos of all flavors. In the transient 'TZO burst' regime (accretion rates ≳ $10^{4}$ solar masses per year), modeled with a steady accretion disk, the neutrino emission lasts up to about $10^{4}$ seconds; in the steady regime ($10^{-4}$ to $10^{4}$ solar masses per year), modeled with spherical accretion, the signal lasts from months to thousands of years. The paper computes the resulting fluxes, applies MSW flavor conversion, and derives detection horizons: Super-Kamiokande and IceCube are already sensitive to bursts beyond the Small Magellanic Cloud, and Super-Kamiokande, JUNO, and Hyper-Kamiokande can detect steady emission for accretion rates ≳ 0.1 solar masses per year at that distance. It also shows that the cumulative diffuse TZO neutrino flux can rival the diffuse supernova neutrino background, and that current DSNB limits already exclude the most extreme accretion and formation scenarios.","pith_inferences":["If hypercritical accretion is episodic rather than steady, the predicted neutrino luminosity and spectra could shift by orders of magnitude, so time-dependent simulations would sharpen or weaken the detection horizons.","A coincident detection of a luminous merger-driven transient with a neutrino burst would provide a direct probe of TZO formation channels, an idea the paper raises but does not develop.","Future MeV diffuse neutrino measurements may need to treat TZOs as a background component in DSNB searches, not just as a standalone signal."],"forward_implications":["A positive neutrino detection from VX Sgr, HV 2112, or HV 11417 would confirm a neutron-star core and constrain its accretion rate.","Super-Kamiokande and IceCube can already test TZO bursts lasting ≳ 0.1 s at accretion rates of 10^6 solar masses per year, out to and beyond the Small Magellanic Cloud.","Steady TZOs accreting above roughly 0.1 solar masses per year would be detectable by Super-Kamiokande, JUNO, and Hyper-Kamiokande within a year of data taking, even beyond the Small Magellanic Cloud.","The diffuse neutrino flux from all TZOs may overlap with the diffuse supernova neutrino background, and current DSNB limits already rule out the most extreme TZO accretion and formation rates."],"supporting_citations":[{"why":"It predicts the existence of Thorne-Zytkow objects, the class of sources this paper proposes to identify through neutrinos.","marker":"Thorne & Zytkow 1975, 1977"},{"why":"It supplies the steady spherical accretion model used for accretion rates below 10^4 solar masses per year.","marker":"Chevalier 1989"},{"why":"It supplies the steady accretion-disk model used for the high-accretion-rate TZO burst scenario.","marker":"Di Matteo et al. 2002"},{"why":"It provides the self-similar solution for the innermost disk region around the neutron star.","marker":"Zhang & Dai 2008"},{"why":"It sets the local TZO rate used for the diffuse flux calculation.","marker":"Nathaniel et al. 2024"},{"why":"It provides the IceCube event-rate and background model used for burst sensitivity.","marker":"Abbasi et al. 2011"},{"why":"It provides Super-Kamiokande's gadolinium-loaded detection efficiency and background rates for steady signals.","marker":"Abe et al. 2024"},{"why":"It defines JUNO's target mass and detection capabilities used for steady TZO searches.","marker":"An et al. 2016"}],"fun_headline_variants":["Neutrinos could expose Thorne-Zytkow objects","Neutrino detection may unmask Thorne-Zytkow stars","Hunting Thorne-Zytkow objects with neutrinos","Neutrinos to unveil hidden neutron-star cores"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions assume steady accretion at a fixed electron fraction Ye = 0.5 onto a 1.4 solar-mass, 10 km neutron star; if hypercritical accretion is intermittent or the electron fraction changes, the neutrino luminosity and spectral shape could differ by orders of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["Neutrinos could expose Thorne-Zytkow objects","Neutrino detection may unmask Thorne-Zytkow stars","Hunting Thorne-Zytkow objects with neutrinos","Neutrinos to unveil hidden neutron-star cores"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000542,"raw_usage":{"total_tokens":2662,"prompt_tokens":1078,"completion_tokens":1584,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":694,"completion_tokens_details":{"reasoning_tokens":1516}},"tokens_in":694,"tokens_out":1584,"duration_ms":14224,"temperature":1.0,"reasoning_tokens":1516,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:52:52.319358+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A year-long targeted search with gadolinium-loaded Super-Kamiokande or JUNO for electron antineutrinos in the 9.3–31.3 MeV window from VX Sgr (at 1.5–1.7 kpc) that finds no excess above the modeled background would falsify the claim that VX Sgr is a steady TZO accreting at ≳ $10^{-2}$ solar masses per year.","supporting_citations":[{"cited_title":"2002, Astrophys","cited_arxiv_id":null,"evidence_quote":"It supplies the steady accretion-disk model used for the high-accretion-rate TZO burst scenario."},{"cited_title":"Population synthesis of Thorne-\\.Zytkow objects: Rejuvenated donors and unexplored progenitors in the common envelope formation channel","cited_arxiv_id":"2407.11680","evidence_quote":"It sets the local TZO rate used for the diffuse flux calculation."},{"cited_title":"2011, Astron","cited_arxiv_id":null,"evidence_quote":"It provides the IceCube event-rate and background model used for burst sensitivity."}],"review_version":1}