{"id":"4b96a45b-2cda-43d6-b989-db9ab3f6fe80","arxiv_id":"2411.14598","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"TXS 0506+056's IceCube neutrinos are consistent with production in an X-ray-obscured core within 10-100 Schwarzschild radii, like NGC 1068, but only if the proton luminosity is near Eddington during flares.","lead":"Astrophysicists argue that the neutrinos IceCube saw from the active galaxy TXS 0506+056 may be made near its central black hole, in a hot X-ray bright region called a corona, rather than in the bright jet that gives the galaxy its blazar label. If right, the same core mechanism that powers neutrinos from NGC 1068 could explain TXS 0506+056, but only if its cosmic ray power approaches the Eddington limit, an extreme requirement.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The consistency argument for the dominant 2014-2015 burst assigns the 2017 flaring-state Swift-XRT flux as the coronal X-ray luminosity, despite this being an epoch-mismatched lower limit on an unknown coronal fraction; since the required proton luminosity scales inversely with this input, the…","rationale":"The reader's verdict and weakest-assumption identification are on target. I confirm that the single most load-bearing issue is the treatment of the X-ray luminosity. I sharpen it: the adopted LX is not only an unconstrained coronal fraction but also from a different epoch (2017 flare) than the dominant neutrino burst (2014-2015). Because Lp scales inversely with LX, this input directly controls whether the Eddington limit is violated by a factor of about 8 or by a much larger factor. The paper's own caveat that the observed flux 'serves as a lower limit' and that the coronal region alone is 'rather unknown' is acknowledged, but the numerical results do not propagate this uncertainty, and the central numbers (e.g., Lp = 3.0e47 erg/s for the 158-day burst) are presented without this caveat in the figures. The requested check—obtaining contemporaneous X-ray coverage or, failing that, bounding the coronal fraction—is the minimal experiment that would settle whether the 2014 burst can be accommodated by the corona model. The verdict remains CONDITIONAL: the model is a legitimate 'may' hypothesis, but it should not be treated as established until the X-ray input is anchored to the correct epoch and component.","tokens_in":13556,"tokens_out":10033,"duration_ms":95610,"concrete_test":"Construct contemporaneous X-ray light curves of TXS 0506+056 over 2014-2015 using Swift/XRT, MAXI, and Swift/BAT archival data, focusing on the 110-day burst window. If a 0.3-10 keV or 15-50 keV flux (or a bound on it) can be measured for that epoch, recompute Lp for the 158-day burst using that flux in Eqs. (2)-(6) and compare with Ledd. If the 2014 X-ray flux was less than one-third of the 2017 flaring-state value, the required Lp exceeds about 2e47 erg/s (more than about 5 Ledd at the adopted black-hole mass) and the core-corona explanation for the dominant burst loses its claimed consistency. If the flux was comparable or higher, the Eddington concern is alleviated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the 2014-2015 neutrino burst originates in the TXS 0506+056 corona rests on the value assigned to the coronal X-ray luminosity LX, which sets the photopion opacity in Eqs. (2)-(3) and, through Eq. (6), the proton luminosity Lp required to produce the observed Lnu. The adopted value, LX = 8.5e44 erg/s, is derived from the Swift-XRT 0.3-10 keV mean flaring-state flux cited to Keivani et al. (2018b), an epoch centered on the 2017 IC-170922A follow-up. The paper explicitly cautions that this observed flux is only a lower limit to the intrinsic flux and that the coronal fraction is unknown, then proceeds to treat it as the coronal luminosity. For the dominant 158-day 2014 burst there is no contemporaneous X-ray measurement; the corona could have been brighter or fainter in 2014 by an unknown factor. Since Lp scales as Lnu R EX / LX, a factor-of-3 decrease in the 2014 coronal LX would raise the required Lp from about 8 Ledd to about 24 Ledd, moving the model from 'challenging' to implausible, while a factor-of-3 increase would bring it below Eddington. The paper's Eddington-level requirement is therefore a reflection of an unconstrained input, not a robust prediction of the model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that the high-energy neutrinos observed by IceCube from the blazar TXS 0506+056, including the dominant 2014-2015 burst and the 2017 IC-170922A event, may originate in the obscured core of the active galactic nucleus rather than in the blazar jet. The authors apply a simple dimensional analysis for photopion production in an X-ray corona: the photopion opacity depends on the coronal X-ray luminosity, the emission radius (10-100 Schwarzschild radii), and the X-ray photon energy, and the required proton luminosity is inferred from the observed neutrino luminosity. They find that the proton luminosity needed for TXS 0506+056 approaches or exceeds the Eddington luminosity during the flares, which they identify as a challenging condition, while for NGC 1068 the requirement is sub-Eddingding. The paper also connects the inferred source density to the diffuse neutrino flux and argues that such core-corona sources can account for the diffuse astrophysical neutrino background. The analysis is qualitative and explicitly acknowledges the lower-limit nature of the observed X-ray flux and the lack of propagation of some uncertainties.","tokens_in":13826,"tokens_out":6552,"duration_ms":58112,"significance":"If the result holds, the paper would strengthen the emerging picture that neutrino emission from active galactic nuclei is powered by their cores rather than by relativistic jets, unifying the interpretation of NGC 1068 and TXS 0506+056. The dimensional analysis is transparent, and the paper is honest about its limitations, including the lower-limit status of the X-ray flux and the unpropagated mass and distance uncertainties. However, the quantitative claim that the 2014-2015 burst requires Eddington-level proton luminosity is not robust, because it depends on an assumed coronal X-ray luminosity that is unconstrained at the relevant epoch and because the L_nu-L_X relation in Fig. 3 is partially circular. These issues would need to be addressed before the central conclusion can be considered reliable.","major_comments":[{"comment":"The central quantitative conclusion that the 2014-2015 neutrino burst requires Eddington-level proton luminosity depends on assigning the observed 2017 flaring-state Swift-XRT flux to the coronal X-ray luminosity L_X of the 2014 epoch. The paper itself states that this observed flux is only a lower limit and that the coronal fraction is unknown, and there is no contemporaneous X-ray observation of the 2014 burst. Since tau_p_gamma is proportional to L_X and the required L_p is proportional to L_nu R E_X / L_X, a factor-of-3 uncertainty in the 2014 coronal L_X changes L_p from about 8 L_edd to about 24 L_edd (or down to about 3 L_edd), which moves the conclusion from \"challenging\" to implausible or to sub-Eddington. Please present the required proton luminosity as a function of L_X and carry out a sensitivity study over the plausible range of coronal fraction and epoch variability.","section":"Section 2, X-ray luminosity paragraph; Eqs. (3) and (6)"},{"comment":"The L_nu-L_X relation in Eq. (6) and Figure 3 is not an independent test of the model, because the proton luminosity L_p that fixes the normalization of the line is itself derived from the measured neutrino luminosity via Eq. (4), using the same opacity and radius assumptions. A source will therefore lie on the plotted band by construction if the same parameters are used in both equations. To make the comparison meaningful, the figure should show how the relation varies with an independently constrained L_p (for example, from accretion power or from an Eddington-ratio assumption) rather than with the inverse-inferred L_p.","section":"Section 2, Eq. (6) and Figure 3"},{"comment":"The stated decision not to propagate the uncertainties in the TXS black hole mass (3.1+29.9-2.7 x 10^8 M_sun) and luminosity distance is problematic for the central claim, because the Eddington ratio L_p/L_edd for the 158-day burst is the quantity that determines whether the scenario is viable. Since L_edd is proportional to the black hole mass, the factor-of-10 mass uncertainty alone spans a range from sub-Eddington to strongly super-Eddington proton power. Please propagate these uncertainties, or provide a separate argument showing that the conclusion is robust to them.","section":"Section 2, Table 1 and the paragraph following it"}],"minor_comments":[{"comment":"The citation \"Fiorillo et al. 2024a\" appears twice in the same sentence of the introduction, and the reference list contains duplicate entries for Fiorillo et al. with different author orders; please consolidate.","section":"Introduction, references"},{"comment":"The text refers to the \"hour-long 2017 burst\", but earlier in the paper the 2017 emission is described as a subdominant flare associated with the single neutrino IC-170922A, and its duration is not otherwise specified; please clarify whether the burst duration is really an hour or whether this is a typo.","section":"Section 3, paragraph on the 2017 burst"},{"comment":"The y-axis of Figure 1 appears to include both tau_p_gamma and tau_gamma_gamma, but the caption only mentions photopion opacities; please clarify what is plotted and define the gamma-gamma opacity line in the caption.","section":"Figure 1"},{"comment":"The numerical coefficient 0.7 in Eq. (6) is quoted without a derivation; please verify that it correctly follows from Eq. (4) in the small-opacity limit with the stated inelasticity and parameter choices.","section":"Equation (6)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a qualitative plausibility argument that fits the scope of the journal. The main concern is that the quantitative consistency check relies heavily on an assumed X-ray luminosity and on parameters inferred from the very neutrino data it seeks to explain; the authors should be encouraged to reframe the central claim as a conditional statement contingent on L_X and the black hole mass, and to provide a sensitivity analysis. The duplicate Fiorillo references and the unclear Figure 1 caption need editorial fixes. I would not reject the paper because the core idea is worth publishing once the load-bearing uncertainties are properly presented."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, honestly hedged consistency argument, not a demonstration. The new bit is the first explicit dimensional-consistency estimate of the corona model for TXS 0506+056, and the resulting constraint that the 2014 flare needs near-Eddington proton luminosity. The paper is worth reading and worth sending to a referee.\n\nWhat it does well: the equations are transparent, the analogies to NGC 1068 and NGC 4151 are apt, and the authors clearly state when they are assuming rather than measuring. They flag that the observed X-ray flux is a lower limit, that the coronal fraction is unknown, and that distance and Schwarzschild-radius errors are not propagated. That is good scientific hygiene.\n\nThe soft spots are real and the stress-test note points at the biggest one. The LX used for TXS is derived from the 2017 flaring-state Swift-XRT mean flux, and that epoch is attached to IC-170922A, not to the 2014 burst that dominates the neutrino flux. There is no contemporaneous X-ray measurement for 2014. Since the required proton luminosity scales as Lnu R EX / LX, a factor of three downward shift in the 2014 coronal X-ray luminosity moves the required Lp from roughly 8 Ledd to 24 Ledd. That turns 'challenging but possible' into 'probably not happening.' The authors' claim that the Eddington requirement is challenging is fair, but it is worth saying that the requirement is a reflection of an unconstrained input, not a prediction of the model.\n\nThere is also some circularity in the Lnu-LX relation in Eq. 6 and Figure 3: Lp is inferred from the measured Lnu through Eq. 4 and then used to draw the correlation. That makes Figure 3 an illustration of the model's normalization rather than an independent check. The proportionality itself was already argued by Kun et al. and Neronov et al., so the novelty is in applying the dimensional estimate to TXS, not in establishing the relation.\n\nNone of this kills the paper. The core-origin idea has independent support from gamma-ray obscuration during the 2014 burst and from the optical flash after IC-170922A, and the authors are appropriately careful to say 'may' rather than 'does.' But the central consistency check depends on parameters that are currently unmeasured, and the LX epoch mismatch is a genuine weakness.\n\nWho should read it: anyone working on AGN neutrino production or IceCube source interpretation. It deserves peer review. My recommendation is to send it out, with a request that the authors address the X-ray epoch mismatch, propagate the black-hole mass and distance uncertainties, and reframe Eq. 6/Figure 3 as a consistency relation rather than a prediction.","headline":"Plausible but not proven: the corona argument for TXS 0506+056 is transparent and honest, but the Eddington requirement hinges on an unmeasured, epoch-mismatched X-ray luminosity.","tokens_in":14447,"tokens_out":2481,"would_cite":true,"duration_ms":23346,"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":"TXS 0506+056's neutrinos may come from its hidden core, not its jet.","keywords":["TXS 0506+056","NGC 1068","neutrino astronomy","active galactic nuclei","coronal neutrino production","photopion opacity","gamma-ray obscuration"],"falsifier":"A NuSTAR or equivalent hard-X-ray monitoring campaign across a future TXS 0506+056 neutrino flare could settle the claim: if the intrinsic coronal X-ray luminosity during the flare is measured to be several times lower than the Swift-XRT flaring value of $2.3\\times10^{-12}\\,\\mathrm{erg\\,cm^{-2}\\,s^{-1}}$, the Eddington-luminosity requirement fails; conversely, a simultaneous GeV–TeV gamma-ray excess during a neutrino burst would falsify the gamma-ray-obscured core picture.","tokens_in":13272,"feed_emoji":"🔭","tokens_out":8808,"duration_ms":78358,"temperature":0.7,"pith_summary":"This paper argues that the high-energy neutrinos observed from the active galaxy TXS 0506+056, including the dominant 2014–2015 burst and the 2017 flare that produced IC-170922A, are produced not in its relativistic jet but in an obscured X-ray corona within 10–100 Schwarzschild radii of its central black hole. The same corona mechanism that explains the neutrino flux from NGC 1068 is applied to TXS 0506+056; both sources are X-ray bright and gamma-ray-obscured, so the neutral-pion gamma rays accompanying neutrino production are absorbed before escaping. Because TXS 0506+056 is roughly a hundred times farther away, the proton luminosity in its corona must approach the Eddington luminosity during the neutrino flares, a demanding but not impossible condition. If correct, this would unify the origin of neutrinos from the two best-established extragalactic neutrino sources and would make X-ray brightness and gamma-ray obscuration, rather than blazar classification, the primary guides for finding neutrino-producing active galaxies.","feed_headline":"TXS 0506+056's neutrinos may come from its hidden core","feed_subtitle":"If true, the same corona mechanism would unify the two clearest extragalactic neutrino sources yet seen.","key_machinery":"The carrying mechanism is the corona-disk model: a dense X-ray corona of radius $R$ around the black hole, with X-ray photon energy density $u_X = L_X/(4\\pi c R^2)$. The load-bearing identity is the photopion opacity formula $\\tau_{p\\gamma} \\simeq 70\\,(v_{\\mathrm{esc}}/c)^{-1}(R/R_S)^{-1}(E_X/1\\,\\mathrm{keV})^{-1}(L_X/L_{\\mathrm{edd}})$, which converts an observed X-ray luminosity and black-hole mass into a prediction for neutrino production efficiency. The companion relation $E_\\nu^2 Q_\\nu \\simeq \\tfrac{3}{8} f_{p\\gamma} E_p^2 Q_p$, with $f_{p\\gamma}=1-e^{-\\tau_{p\\gamma}}$, translates proton luminosity into neutrino luminosity and, in the large-opacity limit, reduces to $L_\\nu \\sim \\tfrac{3}{8}L_p$; these two steps carry the argument from observed fluxes to the required Eddington-level proton power.","core_discovery":"The central claim is a quantitative consistency argument: using the observed soft X-ray luminosity of TXS 0506+056 (about $8.5\\times10^{44}\\,\\mathrm{erg\\,s^{-1}}$ in the 0.3–10 keV band) and a black hole mass near $3\\times10^8\\,M_\\odot$, the photopion opacity $\\tau_{p\\gamma}$ reaches unity at emission radii $R/R_S \\sim 10$, just as it does for NGC 1068, and remains around 0.1 at $R/R_S \\sim 100$. In that corona, protons interacting with X-ray photons produce pions whose neutrinos match the hard $\\sim E^{-2}$ spectra and neutrino luminosities of both the 9.5-year time-integrated emission and the 158-day 2014 burst, while the same corona absorbs the accompanying gamma rays, explaining why the 2014 burst showed no elevated gamma-ray flux. The price is that the required proton luminosity during the flares approaches the Eddington luminosity of about $3.9\\times10^{46}\\,\\mathrm{erg\\,s^{-1}}$, roughly twenty times the adopted X-ray luminosity; the paper treats this as a challenge that may require super-Eddington activity or a proton-dominated corona.","pith_inferences":["Left implicit in the paper is a direct prediction: during future neutrino flares from radio-loud active galaxies, simultaneous hard-X-ray observations should show a coronal flare in the tens-of-keV band, while GeV–TeV gamma rays should dip rather than rise.","The Eddington-level proton luminosity requirement could be tested by searching for coronal proton dominance: if the corona's energy density is proton-dominated, X-ray reflection and reprocessing signatures should differ from those of electron-dominated coronae.","The suggested UV-photon origin for the 2017 single neutrino implies prompt optical brightening should track the neutrino alert on hour timescales, so systematic optical monitoring of neutrino alerts could discriminate between X-ray and UV target photons.","The population argument, applied to TXS-like episodic sources, implies that the duty cycle of Eddington-level flares, not the instantaneous source density alone, controls the diffuse flux; measuring flare rates across a larger sample of X-ray-bright blazars would sharpen this picture."],"forward_implications":["The dominant 2014–2015 burst and the 9.5-year time-integrated neutrino luminosity of TXS 0506+056 are consistent with production in an X-ray corona at $R \\sim 10\\text{--}100\\,R_S$, with proton luminosity at or near Eddington during the flaring interval.","The corona that produces the neutrinos also absorbs the neutral-pion gamma rays that accompany them, explaining why the 2014 burst showed no elevated gamma-ray flux and why the source is gamma-ray-obscured at the relevant times.","X-ray luminosity and gamma-ray obscuration, rather than blazar classification, become practical markers for selecting neutrino-producing active galaxies; the paper shows NGC 1068, NGC 4151, and TXS 0506+056 all sit close to the same $L_\\nu$–$L_X$ scaling.","If the corona power is limited by $L_p \\lesssim L_X$, the model caps the corona-powered neutrino luminosity of TXS 0506+056 at roughly $3.4\\times10^{44}\\,\\mathrm{erg\\,s^{-1}}$, so the full observed flux requires exceeding that bound during flares.","A population of NGC 1068-like sources with local density near $10^{-6}\\,\\mathrm{Mpc^{-3}}$ can account for the diffuse astrophysical neutrino flux, while TXS-like episodic sources contribute if their Eddington-level flares are on for about five percent of the time."],"supporting_citations":[{"why":"Reports the 2014–2015 neutrino burst of 19 events on a background of fewer than six; it supplies the dominant neutrino luminosity the model must reproduce.","marker":"Aartsen et al. 2018"},{"why":"Reports IC-170922A and the flaring-blazar association; it supplies the 9.5-year spectrum and the 158-day flare spectra used as constraints.","marker":"IceCube Collaboration et al. 2018"},{"why":"Identifies NGC 1068 as a neutrino-emitting active galaxy and anchors the corona mechanism that the paper applies to TXS 0506+056.","marker":"Abbasi et al. 2022"},{"why":"Provides the Swift-XRT 0.3–10 keV flaring-state flux and the NuSTAR hard-X-ray detection, the principal X-ray input for the opacity calculation.","marker":"Keivani et al. 2018b"},{"why":"Supplies the black hole mass estimate for TXS 0506+056 from the bulge luminosity, fixing the Eddington luminosity scale.","marker":"Padovani et al. 2019"},{"why":"Gives the redshift z ≈ 0.34 and luminosity distance 1774 Mpc that set the distance penalty in converting flux to luminosity.","marker":"Paiano et al. 2018"},{"why":"Provides the detailed core-corona modeling of NGC 1068 whose τpγ ≈ 1 result the paper's dimensional analysis reproduces.","marker":"Murase 2022"},{"why":"Supplies the intrinsic 2–10 keV X-ray flux and torus-model parameters used to normalize the comparison sources NGC 1068 and NGC 4151.","marker":"Ricci et al. 2017"}],"fun_headline_variants":["TXS 0506+056's neutrinos may share NGC 1068's hidden-core origin","Core corona may explain TXS 0506+056 and NGC 1068 neutrinos","TXS 0506+056 neutrinos trace to core, not jet"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the Swift-XRT 0.3–10 keV flux observed from TXS 0506+056 is comparable to the intrinsic X-ray luminosity of its corona, even though the paper concedes that the observed flux is only a lower limit and the coronal fraction is unknown.","fun_headline_variants_meta":{"raw":{"variants":["TXS 0506+056's neutrinos may share NGC 1068's hidden-core origin","Core corona may explain TXS 0506+056 and NGC 1068 neutrinos","TXS 0506+056 neutrinos trace to core, not jet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001437,"raw_usage":{"total_tokens":5810,"prompt_tokens":981,"completion_tokens":4829,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":4756}},"tokens_in":597,"tokens_out":4829,"duration_ms":30264,"temperature":1.0,"reasoning_tokens":4756,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:07:26.196580+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A NuSTAR or equivalent hard-X-ray monitoring campaign across a future TXS 0506+056 neutrino flare could settle the claim: if the intrinsic coronal X-ray luminosity during the flare is measured to be several times lower than the Swift-XRT flaring value of $2.3\\times10^{-12}\\,\\mathrm{erg\\,cm^{-2}\\,s^{-1}}$, the Eddington-luminosity requirement fails; conversely, a simultaneous GeV–TeV gamma-ray excess during a neutrino burst would falsify the gamma-ray-obscured core picture.","supporting_citations":[{"cited_title":"2018, Science, 361, eaat1378","cited_arxiv_id":null,"evidence_quote":"Reports IC-170922A and the flaring-blazar association; it supplies the 9.5-year spectrum and the 158-day flare spectra used as constraints."},{"cited_title":"2018, The Astrophysical Journal Letters, 854, L32","cited_arxiv_id":null,"evidence_quote":"Gives the redshift z ≈ 0.34 and luminosity distance 1774 Mpc that set the distance penalty in converting flux to luminosity."},{"cited_title":"J., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the intrinsic 2–10 keV X-ray flux and torus-model parameters used to normalize the comparison sources NGC 1068 and NGC 4151."}],"review_version":1}