{"id":"46230515-db2c-4638-9ef7-2dd7f23b72b1","arxiv_id":"2502.01738","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":8,"one_line_summary":"Even with extreme assumptions, the corona of TXS 0506+056 cannot produce enough neutrinos to match IceCube, leaving the jet as the preferred neutrino origin.","lead":"Using a model of magnetic reconnection in the corona of the blazar TXS 0506+056, the authors compute how many neutrinos the corona should emit and find the number is too small to explain IceCube's detections. The result strengthens the traditional view that the blazar's jet, not its core, produces the observed neutrinos.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central conclusion hinges on the coronal X-ray luminosity estimate; if LX is instead at the observed Swift-XRT/NuSTAR level, the HL scenario could reach the IceCube flux, so the claim is contingent on the jet-dominated interpretation of the observed X-rays.","rationale":"The reader's weakest_assumption identifies the coronal X-ray luminosity inference as the load-bearing step, and my review agrees. The paper's internal logic is sound: given a coronal LX, the AM3 calculations and the energy-budget argument (ηp+ηX<1) show that even optimistic proton loading cannot match the observed neutrino flux. The EL scenario's failure against the 2014/15 flare is also robust because it already saturates the MAXI X-ray limits. However, the entire argument stands or falls on the assumption that the observed X-ray emission from TXS 0506+056 is jet-dominated, so that the true coronal LX is bounded by the disk-scaling relations. If instead the corona contributes significantly to the observed X-ray flux, the HL scenario's neutrino output could increase sufficiently to reach the IceCube 10-year flux, undermining the central claim. The paper explicitly acknowledges the factor-of-10 spread between the two empirical methods but does not propagate uncertainties or quantify the probability that the comet-like X-ray luminosity could be at the observed level. This is not an internal inconsistency but a genuine astrophysical uncertainty that warrants a conditional verdict, exactly as the reader recommended. Therefore no adjustment to the verdict is needed; the conditionality already captures the need for an uncertainty analysis on LX and a clearer justification for rejecting the observed X-rays as coronal.","tokens_in":17109,"tokens_out":7294,"duration_ms":71233,"concrete_test":"Recompute the HL scenario in AM3 with the same parameters (ηp=1, ηX=0.5, R=1.4rg, HGF proton spectrum) but set LX to the observed Swift-XRT lower bound 3.8e44 erg/s, and in a second run to the NuSTAR upper bound 3e45 erg/s. Compare the predicted 10-year neutrino flux to the IceCube point-source flux (Abbasi et al. 2022). If either run reaches or exceeds the IceCube flux, the central claim fails under the observer's own reported X-ray range, confirming that the conclusion depends on rejecting the observed X-rays as jet-dominated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's conclusion that coronal neutrinos from TXS 0506+056 are too low rests on the coronal X-ray luminosity range LX ≈ 4e43–4e44 erg/s derived in §3. This range is a factor of ~10–75 below the observed 0.3–10 keV Swift-XRT luminosity (3.8e44–3e45 erg/s) and the NuSTAR 3–78 keV luminosity (1.5e45–3e45 erg/s). The two empirical scaling methods used (Arcodia et al. 2019 vs. Runnoe et al. 2012) already disagree by a factor of 10, and no uncertainty propagation is provided. If the coronal LX were instead comparable to the observed X-ray luminosity, as assumed by Khatee Zathul et al. (2024) and Yang et al. (2024), then with ηp=1 (the HL upper limit) the proton luminosity Lp = (ηp/ηX) LX would rise by roughly the same factor, and the predicted neutrino flux in the HL scenario could plausibly reach the IceCube 10-year point-source flux. The EL scenario does not resolve this: it is designed to saturate the 2014/15 MAXI limits and is not a test of the steady-state LX assumption. Thus the load-bearing premise is the deconvolution of jet versus coronal X-ray components, not the reconnection model itself.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper tests whether the neutrinos from TXS 0506+056 can originate in the AGN corona rather than the blazar jet. Using the magnetic-reconnection corona model of Fiorillo et al. (2024b), the authors estimate the coronal X-ray luminosity of TXS 0506+056 from its masquerading-BL-Lac disk luminosity and empirical AGN scaling relations, obtaining L_X ≈ 4×10^43–4×10^44 erg/s. They then compute the predicted neutrino and cascade emission with the AM3 code in three scenarios: a low-luminosity (LL) case, a high-luminosity (HL) case with η_p=1, and an extreme-luminosity (EL) case saturating the 2014/15 MAXI X-ray limits. In all cases the predicted coronal neutrino flux falls below the IceCube 10-year point-source flux and below the 2014/15 and 2017 inferred fluxes. The authors conclude that the corona is disfavored as the dominant neutrino production site and that the blazar jet remains the preferred location. They also argue in an appendix that stochastic acceleration cannot energize protons to the required PeV energies.","tokens_in":17397,"tokens_out":4939,"duration_ms":51395,"significance":"If the central claim holds, the result is significant for multi-messenger astrophysics: it would provide a concrete counterexample to the suggestion that AGN-core/corona emission dominates neutrino production in all AGN, and it would reinforce the jet interpretation for TXS 0506+056. The paper's main strength is that the model is not fitted to the IceCube neutrino data; the neutrino flux is a genuine prediction compared against external measurements. The authors also give a transparent energy-budget argument and explicitly flag the unphysical nature of the HL scenario's η_p=1. The use of a public radiative code (AM3) and the decomposition of the cascade into individual components are additional positive features. The principal weakness is that the conclusion is contingent on the still-uncertain coronal X-ray luminosity, which is inferred through scaling relations rather than measured; this is the load-bearing premise that needs quantitative support.","major_comments":[{"comment":"The central conclusion rests on the assumed coronal X-ray luminosity range L_X ≈ 4×10^43–4×10^44 erg/s. The two empirical scaling methods used in this section differ by a factor of ~10 (L_2-10 keV ≈ 10^44 erg/s versus 8.2×10^42 erg/s), and no uncertainty propagation is provided. The observed Swift-XRT and NuSTAR luminosities are about 10–75 times higher. Since L_p = (η_p/η_X) L_X, a coronal X-ray luminosity comparable to the observed values would raise the HL prediction by roughly the same factor and could plausibly reach the IceCube 10-year flux. The paper dismisses the observed X-rays as jet-dominated, but this is an assumption rather than a demonstrated decomposition. A quantitative justification is needed—for example, a variability or spectral decomposition argument, or a plot of the coronal L_X required to match the IceCube flux against the observationally allowed range—before the 'too low' claim is robust.","section":"Section 3"},{"comment":"The HL scenario uses η_p=1 while η_X=0.5, so η_p+η_X=1.5, violating the energetic constraint η_p+η_X<1 that the authors themselves state. They acknowledge this, but the steady-state conclusion is therefore established only by comparing IceCube with an already unphysical upper limit. The numerical margin between the HL prediction and the IceCube flux is not stated; reporting the ratio of predicted to observed neutrino flux in each scenario would let the reader see exactly how much the coronal L_X would need to be revised upward to close the gap.","section":"Section 4, Table 1"},{"comment":"The EL case is constructed to saturate the 2014/15 MAXI 10 keV upper limit, but the resulting broadband L_X = 2×10^46 erg/s is not shown against the MAXI measurement in Fig. 2, and the paper does not quantify whether such a high L_X is consistent with the optical/UV SED (it notes a possible tension but does not evaluate it). Since the EL scenario already fails, the authors should demonstrate that the MAXI limit actually constrains the broadband coronal luminosity; otherwise the exclusion of a coronal origin for the 2014/15 flare is weaker than claimed.","section":"Section 4, EL scenario"}],"minor_comments":[{"comment":"The reference list contains duplicates: Abbasi et al. (2022) appears twice, Petropoulou et al. (2020) appears three times, Zhang et al. (2023) appears twice, and Wilkins & Gallo (2015) appears twice; please consolidate them.","section":"References"},{"comment":"There is a typo, 'magntization', which should be 'magnetization'.","section":"Appendix A"},{"comment":"The notation η_{p,-0.3} is used although the fiducial value adopted in Table 1 is η_p=0.1, not 0.3; please define the scaling variable explicitly.","section":"Eq. (2)"},{"comment":"The phrase 'AGN blazar' is redundant; consider using 'blazar' throughout.","section":"Section 6"}],"recommendation":"major_revision","confidential_remarks":"The paper's conclusion is plausible and well argued under the masquerading-BL-Lac interpretation, but the decisive L_X estimate is not yet robust enough for the strong claim made in the abstract. The authors should be asked to quantify the systematic uncertainty in Section 3 and to show the sensitivity of the neutrino flux to L_X, ideally including the observed Swift-XRT/NuSTAR luminosities as an upper limit. If even then the coronal flux falls short under physically allowed parameters, the paper would be very strong; as it stands, the central claim is contingent on an assumption that is disputed by other recent works."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one for the energy-budget argument, not for the corona model itself. The paper makes a clean, testable point: even under the most favorable parameters, the F24 coronal reconnection model cannot produce the IceCube flux from TXS 0506+056, and an extreme-luminosity version saturating the MAXI limits still misses the 2014/15 flare. That is a useful counter-result to the recent corona/accretion-flow proposals.\n\nWhat is actually new: they treat TXS as a masquerading BL Lac, use the inferred disk luminosity and AGN X-ray scaling relations to estimate a coronal X-ray luminosity of roughly 4e43-4e44 erg/s, one to two orders below the observed jet-dominated X-rays, and then run AM3 for LL, HL, and EL scenarios. The HL case uses eta_p=1, already unphysical when eta_X=0.5, and still falls short. That makes the null result fairly robust within the model.\n\nThe soft spot is the load-bearing premise: the coronal LX. The two scaling relations (Arcodia et al. vs. Runnoe et al.) already diverge by a factor of 10, and no uncertainty is propagated. The stress-test point is fair: if the corona's X-rays were at the Swift-XRT/NuSTAR level, as Yang et al. and Khatee Zathul et al. assume, the HL scenario could plausibly reach the IceCube flux. The paper argues those observed X-rays are jet-dominated, which is reasonable for a blazar, but it is an inference, not a measurement. The LL/HL bracketing is honest but does not quantify how likely each value is. A serious referee should ask for a Monte Carlo over the scaling-relation scatter and a sensitivity curve of the maximum neutrino flux to LX.\n\nThe model parameters come from prior PIC literature and the authors' own F24 framework, not fitted to TXS data, so the neutrino flux is a genuine prediction. Self-citation is appropriate here; the F24 model is the same machinery. AM3 is a public code; configuration files would help reproducibility, but the central argument is analytic enough that the numerics are not decisive.\n\nBottom line: a well-argued, significant within-subfield result that deserves peer review. The conclusion is contingent on the coronal LX estimate, but the paper states this clearly and brackets it. I would accept with a request for uncertainty propagation on LX and, ideally, AM3 configuration files. It is a solid paper with one soft spot, not a fatal one.","headline":"A clean energy-budget test that rules out a coronal origin for TXS 0506+056 neutrinos under the F24 reconnection model, with the main caveat being the factor-of-10 uncertainty in the inferred coronal X-ray luminosity.","tokens_in":18052,"tokens_out":2184,"would_cite":true,"duration_ms":21019,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper finds that the corona of the blazar TXS 0506+056 is too faint to produce the neutrinos IceCube detects, so the blazar jet remains the preferred origin.","keywords":["TXS 0506+056","AGN corona","neutrino astronomy","magnetic reconnection","blazar jet","IceCube","masquerading BL Lac","X-ray luminosity scaling"],"falsifier":"Measure the coronal X-ray luminosity of TXS 0506+056 with X-ray reflection spectroscopy or via a soft excess or reverberation lag during a quiescent state. If the 2-10 keV coronal luminosity is close to the source's observed X-ray level (about $10^{45}$ erg/s) rather than the scaling-relation value near $10^{44}$ erg/s, then the coronal model with $\\eta_p=1$ would predict neutrinos at the IceCube 10-year flux, directly contradicting the paper's central claim.","tokens_in":16820,"feed_emoji":"🔭","tokens_out":14046,"duration_ms":111543,"temperature":0.7,"pith_summary":"The paper asks whether the neutrinos observed from the blazar TXS 0506+056 by IceCube could be produced in the active galactic nucleus corona rather than in the relativistic jet. Using a magnetic reconnection model of proton acceleration and an observationally motivated estimate of the coronal X-ray luminosity, it predicts the neutrino flux the corona would emit and compares it with the observed steady signal, the 2017 flare event, and the 2014/15 flare. In all cases, even with the most optimistic assumptions about the proton energy budget and the coronal X-ray luminosity, the predicted flux falls short of what IceCube sees. The paper therefore concludes that the corona cannot be the dominant neutrino production site for this source, and the jet remains the preferred location.","feed_headline":"Corona can't make TXS 0506+056's neutrinos","feed_subtitle":"Even extreme coronal luminosity falls short of IceCube's flux, keeping the jet as the source.","key_machinery":"The central object is the magnetic reconnection layer model of the corona, denoted F24 in the paper, in which protons are accelerated by the reconnection electric field. This model fixes the coronal magnetic field from the X-ray luminosity and determines the injected proton luminosity as $L_p = (\\eta_p/\\eta_X)\\,L_X$, so that the predicted neutrino flux is directly proportional to the inferred coronal X-ray luminosity. The proton injection spectrum is a broken power law motivated by particle-in-cell simulations, with the break at $\\sigma_p m_p c^2$ and a high-energy cutoff set by the balance of acceleration and cooling; the paper considers both strong and weak guide-field cases. Secondary production and the electromagnetic cascade are simulated with a radiative code, producing neutrino and X-ray spectra that are compared with the multi-wavelength and IceCube observations. The energy-link between $L_X$ and $L_p$ is what makes a low coronal X-ray luminosity fatal for the coronal origin hypothesis.","core_discovery":"The paper's central claim is that the coronal neutrino emission from TXS 0506+056 is too low to account for the IceCube observations. The argument is energetics: the corona's X-ray luminosity, estimated from the disk luminosity of this masquerading BL Lac using scaling relations for non-jetted AGN, lies between about $4\\times10^{43}$ and $4\\times10^{44}$ erg/s, one to two orders of magnitude below the observed X-ray output of the source, which is attributed to the jet. Because the reconnection model ties the proton luminosity to the X-ray luminosity via $L_p=\\eta_p L_X/\\eta_X$, a fainter corona implies a fainter neutrino flux. Even in the high-luminosity scenario with $\\eta_p=1$, a value already in tension with the requirement $\\eta_p+\\eta_X<1$, the predicted neutrino flux cannot reach the 10-year IceCube flux; and in an extreme scenario that saturates the MAXI X-ray upper limits, the neutrinos still cannot match the 2014/15 flare.","pith_inferences":["If other masquerading BL Lacs have similarly faint coronae, their coronal neutrino fluxes would be proportionally suppressed, lowering the expected contribution of blazar cores to the diffuse neutrino background.","The same energetics test could be applied to other neutrino-bright AGN with known disk luminosities, converting this single-source study into a population-level discriminator between coronal and jet neutrino production.","A future detection of a soft X-ray excess or reflection component from TXS 0506+056 during a low-jet state would directly test the assumed coronal luminosity and could revive the coronal interpretation if it proved brighter than the scaling relations predict.","The paper's negative result strengthens the case that the jet (or possibly an external radiation field in the jet) is responsible for the 2017 flare, while leaving the 2014/15 flare as a separate puzzle."],"forward_implications":["The corona is ruled out as the dominant neutrino production site for TXS 0506+056, so the blazar jet remains the standard explanation for its IceCube neutrinos.","The 2014/15 neutrino flare cannot be explained by an extreme coronal outburst without violating X-ray upper limits; its origin remains an open problem also for jet models.","The X-ray emission from the corona and jet have sufficiently different spectra that multi-component fits to quiescent X-ray data can further constrain or confirm the coronal contribution.","The result provides a counter-example to the idea that all observed AGN neutrinos come from the core region, keeping non-jetted AGN like NGC 1068 as a separate class.","Population-level estimates of the blazar contribution to the diffuse neutrino flux should treat the coronal component of masquerading BL Lacs as subdominant."],"supporting_citations":[{"why":"Provides the 10-year IceCube point-source neutrino flux that the coronal models must reproduce.","marker":"Abbasi et al. 2022"},{"why":"Reports the 2017 flare-associated neutrino event and its inferred flux and upper limits.","marker":"IceCube Collaboration et al. 2018a"},{"why":"Provides the neutrino spectrum of the 2014/15 flare used to test the extreme coronal scenario.","marker":"IceCube Collaboration et al. 2018b"},{"why":"Classifies TXS 0506+056 as a masquerading BL Lac and gives the black hole mass, bolometric, and disk luminosities that seed the coronal X-ray estimate.","marker":"Padovani et al. 2019"},{"why":"Supplies the bolometric-to-UV and 2-10 keV luminosity relations used to infer the coronal X-ray luminosity.","marker":"Runnoe et al. 2012"},{"why":"Provides the UV-to-2 keV X-ray scaling used to derive the coronal X-ray luminosity for TXS 0506+056.","marker":"Arcodia et al. 2019"},{"why":"The magnetic reconnection model (F24) relating proton luminosity to X-ray luminosity and setting the acceleration timescale.","marker":"Fiorillo et al. 2024b"},{"why":"Motivates the broken power-law proton injection spectrum adopted in the coronal model.","marker":"Zhang et al. 2021"},{"why":"Provides the multi-epoch SED, Swift-XRT and NuSTAR X-ray measurements, and MAXI upper limits that bracket the coronal luminosity.","marker":"Petropoulou et al. 2020"},{"why":"Provides the radiative code used to simulate the multi-messenger spectra and cascades.","marker":"Klinger et al. 2024"}],"fun_headline_variants":["Corona can't explain TXS neutrino flux","TXS neutrinos not from corona, jet remains","Corona fails to make TXS neutrinos","TXS neutrino origin stays in the jet","Corona too weak for TXS 0506+056 neutrinos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the corona of TXS 0506+056 is only as bright in X-rays as ordinary AGN scaling relations predict, roughly $10^{43}$ to $10^{44}$ erg/s, rather than as bright as the observed, jet-dominated X-ray flux of about $3\\times10^{44}$ to $3\\times10^{45}$ erg/s; if the corona were as bright as the direct observations, the model's neutrino flux could reach the measured level.","fun_headline_variants_meta":{"raw":{"variants":["Corona can't explain TXS neutrino flux","TXS neutrinos not from corona, jet remains","Corona fails to make TXS neutrinos","TXS neutrino origin stays in the jet","Corona too weak for TXS 0506+056 neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000269,"raw_usage":{"total_tokens":1659,"prompt_tokens":1021,"completion_tokens":638,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":562}},"tokens_in":637,"tokens_out":638,"duration_ms":6005,"temperature":1.0,"reasoning_tokens":562,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T14:39:16.591200+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the coronal X-ray luminosity of TXS 0506+056 with X-ray reflection spectroscopy or via a soft excess or reverberation lag during a quiescent state. If the 2-10 keV coronal luminosity is close to the source's observed X-ray level (about $10^{45}$ erg/s) rather than the scaling-relation value near $10^{44}$ erg/s, then the coronal model with $\\eta_p=1$ would predict neutrinos at the IceCube 10-year flux, directly contradicting the paper's central claim.","supporting_citations":[],"review_version":1}