{"id":"f3f0784e-a9f2-48e3-8f7d-ff2aa074528f","arxiv_id":"1908.11849","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":14,"one_line_summary":"A proton blazar model in which cold jet protons supply pp targets reproduces the gamma-ray SED and the IceCube-170922A neutrino event rate from TXS 0506+056.","lead":"This paper models the 2017 gamma-ray flare and the co-located IceCube neutrino from the blazar TXS 0506+056 using proton-proton collisions inside the jet, with cold protons present under charge neutrality as the targets. It reports a consistent description of the multi-wavelength spectrum and the observed neutrino event rate with one set of parameters.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The pp target density is fixed by the assumed electron acceleration efficiency χ_e≈10^-3; for the Table I parameters the required gamma-ray luminosity has no charge-neutral solution for χ_e≳5×10^-3, so the central claim is not robust to the cited uncertainty.","rationale":"The reader's weakest-assumption analysis identified the same load-bearing element: the cold-proton target density is derived from the assumed electron acceleration efficiency χ_e ≈ 10^-3. I agree, and the stress test sharpens the objection from a mere energy-budget concern to an existence-of-solution constraint. For the numerical parameters in Table I, the pp luminosity needed to match the observed gamma-ray SED forces the product n_H n_p. But charge neutrality with a normal electron-proton plasma gives n_H = n_e - n_p, so this product is bounded above by n_e^2/4. Because n_e scales as 1/χ_e, the bound becomes (n'_e,h)^2/(4χ_e^2). Comparing with the required baseline product shows that no value of L'_p can reproduce the gamma-ray flux once χ_e exceeds about 5e-3. This is only a factor of about five above the assumed value, and the cited literature supports a general low efficiency rather than a precise source-specific number. Therefore the success of the model is conditional on χ_e being in a narrow range, and the paper should demonstrate robustness across that range. The proposed test is a focused parameter scan that would settle whether the concern lands: recompute the pp gamma-ray and neutrino fluxes with χ_e = 1e-2 and with the threshold value, allowing L'_p to vary, and verify whether a solution exists. This does not change the reader's conditional verdict; it strengthens the conditions that should be attached to acceptance.","tokens_in":13790,"tokens_out":12955,"duration_ms":121835,"concrete_test":"Using the published Table I parameters (R'_b = 2.2e16 cm, L'_e = 2.3e42 erg/s, n'_e,h = 1.7e3 cm^-3, L'_p = 1e46 erg/s), vary χ_e from 1e-4 to 1e-2. For each value, set n_e = n'_e,h/χ_e, express n_p through Eq. (2) as a function of L'_p, enforce n_H = n_e - n_p ≥ 0, and solve for L'_p that reproduces the hadronic gamma-ray component inferred from the Fermi-LAT/MAGIC excess over the fitted leptonic contribution. Record whether a solution exists, the required L'_p, and the resulting Nνµ with a Poisson 90% interval. If no solution exists for χ_e ≥ 5e-3, the paper should report this bound and justify the adopted χ_e for this source.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Within the assumed normal-plasma, charge-neutral jet, the pp target density is not an independent observable: n_H = n_e - n_p, with n_e = n'_e,h/χ_e. Using Table I, n'_e,h = 1.7e3 cm^-3, χ_e = 1e-3, and n_p ≈ 1.7e4 cm^-3 from L'_p = 1e46 erg/s via Eq. (2), one obtains n_H ≈ 1.68e6 cm^-3. The pp gamma-ray emissivity is proportional to n_H n_p = (n_e - n_p)n_p, which for fixed n_e has its maximum at n_p = n_e/2, with value n_e^2/4. To reproduce the observed Fermi-LAT/MAGIC gamma-ray component, this product must be at least the baseline value, roughly (1.68e6)(1.7e4). If χ_e were 1e-2, then n_e = 1.7e5 cm^-3 and the maximum possible product is only (8.5e4)^2 = 7.2e9, about a factor of four below the required value; no choice of L'_p can compensate. The threshold is χ_e ≈ 5e-3, only a factor of five above the assumed value. The cited references [25-27] provide a general electron acceleration efficiency for AGN jets, not a source-specific measurement for TXS 0506+056, so the allowed range plausibly spans this threshold. Thus the central claim rests on χ_e lying in a narrow window; an independent determination or a robustness scan is required.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a proton blazar model for the blazar TXS 0506+056, in which cold, non-relativistic protons present in the jet to satisfy charge neutrality provide the target matter for pp interactions with shock-accelerated protons. The model uses a one-zone blob with a broken power-law electron distribution for synchrotron and inverse-Compton emission, and a power-law proton distribution whose pp interactions produce gamma rays (via neutral-pion decay and cascades) and neutrinos. With parameters listed in Table I, the authors reproduce the multi-wavelength SED from radio to VHE gamma rays and compute an expected IceCube muon-neutrino event count of about 1.0 in 0.5 yr, consistent with the observed IceCube-170922A event. The paper also discusses the jet power budget and contrasts the scenario with cloud-in-jet or pγ models.","tokens_in":14224,"tokens_out":4109,"duration_ms":39120,"significance":"If the central claim holds, this is a useful consistency check showing that a standard electron-proton jet with charge neutrality can simultaneously account for the observed SED and the IceCube-170922A neutrino without invoking an external target cloud. The neutrino event rate is computed from the same parameters used for the gamma-ray fit (no additional parameter is tuned to the neutrino), which is a strength. The radiative calculations are standard and the numerical results are presented in a reproducible way. However, the model's significance is limited by the large number of free parameters in Table I and by the sensitivity of the pp target density to the assumed electron acceleration efficiency, as detailed below. The paper does not provide uncertainties or a robustness scan, so the consistency is indicative rather than demonstrative.","major_comments":[{"comment":"","section":"§II (after Eq. 2) and §III (n_H derivation)"},{"comment":"","section":"§III (jet power discussion)"}],"minor_comments":[{"comment":"","section":"§III and Table I (spectral index sign)"},{"comment":"","section":"Fig. 1 caption"},{"comment":"","section":"References [29] and [41]"},{"comment":"","section":"§III (goodness of fit)"},{"comment":"","section":"§IV (quiescent state)"}],"recommendation":"major_revision","confidential_remarks":""},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead Banik & Bhadra on TXS 0506+056. The new wrinkle is that they don't need external clouds or a dense photon field; cold protons required by charge neutrality in a normal electron-proton jet provide the pp targets. That's a genuinely different channel from Liu et al. (cloud-in-jet) and Ansoldi et al. (p-gamma). The model reproduces the SED's two bumps and gets ~1 IceCube event in 0.5 yr with the same parameters. The radiative calculations are standard, the machinery is transparent, and they are honest that the neutrino rate is a consistency check, not an independent prediction.\n\nThe soft spot is the target density. n_H = n_e - n_p, and n_e = n'_e,h / χ_e with χ_e = 1e-3 assumed from general AGN literature. The stress-test arithmetic holds up: with Table I, the required pp emissivity needs n_H n_p ≈ 2.9e10 cm^-6, and for χ_e = 1e-2 the maximum possible product (n_e^2/4) is a factor ~4 lower; the threshold is χ_e ≈ 5e-3, only five times the assumed value. Since χ_e is not measured for this source, the central conclusion rests on a narrow ledge. A robustness scan over χ_e (and over pair content) should be mandatory.\n\nOther soft spots: the model has 14 free parameters with no reported uncertainties or statistical fit quality, and the total jet power comes out near 1e49 erg/s, which they argue is acceptable for a 1e10 solar mass black hole but is nevertheless borderline. None of this kills the model, but it keeps the claim at 'viable' rather than 'demonstrated.'\n\nBottom line: this deserves a serious referee. The mechanism is worth taking seriously, the fragility is clearly diagnosable, and a revision that scans χ_e and adds error treatment would turn it into a solid paper. I would not cite it as established, but I would bring it to reading group to discuss the limits of charge-neutrality-based targets.","headline":"A viable but fragile hadronuclear channel for TXS 0506+056; the charge-neutrality target density hinges on χ_e within a factor of five of the threshold.","tokens_in":14769,"tokens_out":1972,"would_cite":false,"duration_ms":17299,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["96.50.S-","98.70.Rz","98.70.Sa"],"model":"deepseek-v4-flash","headline":"The paper claims that the same proton–proton collisions inside the jet of TXS 0506+056 produce both its gamma-ray flare and the IceCube-170922A neutrino, with the target protons supplied by charge neutrality.","keywords":["blazar","TXS 0506+056","IceCube-170922A","neutrino","gamma-ray flare","proton blazar model","pp interaction","charge neutrality"],"falsifier":"A decisive test would be a simultaneous gamma-ray and neutrino observation of a similar BL Lac flare: if the measured muon-neutrino event rate in IceCube deviates from the predicted $N_{\\nu_\\mu}\\approx 1$ per 0.5 yr by more than the uncertainties for a given Fermi-LAT/MAGIC spectrum, the charge-neutrality target density would be ruled out. A measurement of $\\chi_e$ for this source from the synchrotron-to-Compton luminosity ratio would settle the key assumption directly.","tokens_in":13578,"feed_emoji":"🔭","tokens_out":6623,"duration_ms":54277,"temperature":0.7,"pith_summary":"The paper argues that the gamma-ray flare and the IceCube neutrino event from the blazar TXS 0506+056 can both be produced by proton–proton collisions inside the jet, with no external cloud of gas needed. The target protons are the \"cold\" non-relativistic protons that must exist in the jet to keep it charge neutral while it carries many more electrons than relativistic protons. Using the proton blazar model, the authors fit the multi-wavelength emission with synchrotron and inverse-Compton radiation, and show that the same accelerated-proton population that generates the observed GeV–TeV gamma rays through pp collisions also yields about one IceCube muon-neutrino event in half a year, matching IceCube-170922A. The claim matters because it offers a parameter-free link between the two observed signals and removes the need for broad-line regions or dense clouds that BL Lac objects lack.","feed_headline":"A blazar's cold protons explain its neutrino and gamma-ray flare","feed_subtitle":"The same pp collisions that make the gamma rays also produce IceCube-170922A, with no extra parameters.","key_machinery":"The central object is the charge-neutrality condition for the jet plasma, $n_H = n'_e - n'_p$, with $n'_e = n'_{e,h}/\\chi_e$; it converts the observed synchrotron-emitting electron population into a dense sea of cold protons. The pp interaction emissivity of Kelner et al. (2006), together with electromagnetic cascades from $\\gamma\\gamma$ absorption and pair injection, then generates both the high-energy gamma-ray bump and the neutrino flux from the same accelerated proton spectrum.","core_discovery":"Within the proton blazar framework, the authors establish that charge neutrality of the jet fixes the number density of cold target protons to $n_H = 1.68\\times10^6$ cm$^{-3}$ once the relativistic electron population (with acceleration efficiency $\\chi_e\\approx 10^{-3}$) is known. Shock-accelerated protons with a $\\gamma_p^{-2.13}$ spectrum and maximum energy $E'_{p,\\max}=10$ PeV interacting with these targets reproduce the Fermi-LAT and MAGIC gamma-ray spectrum, and the same interaction, with no additional free parameters, produces a muon-neutrino flux corresponding to $N_{\\nu_\\mu}\\approx 1.0$ event in 0.5 yr, consistent with IceCube-170922A. The total jet power is $L_{\\rm jet}\\approx 1.3\\times10^{49}$ erg/s, mildly above the Eddington luminosity for a $\\sim10^{10}\\,M_\\odot$ black hole, which the authors argue is acceptable for an outburst.","pith_inferences":["A direct test: for a given gamma-ray flare amplitude, the model predicts the neutrino rate; stacking many BL Lac flares in IceCube and Fermi-LAT data would measure whether the pp target density inferred from charge neutrality is universal or source-specific.","The same charge-neutrality argument could be applied to other broad-line-free BL Lacs, predicting their neutrino rates from their measured SEDs without invoking external photon fields or clouds.","If future observations show a neutrino flare with no accompanying gamma-ray flare (as seen for TXS 0506+056 in 2014–2015), the pp picture would require the target density to persist while the accelerated proton luminosity varies independently of the electron component."],"forward_implications":["The same pp interaction that reproduces the gamma-ray spectrum yields the observed IceCube-170922A event rate without adding free parameters, so the two signals are tied one-to-one in this model.","No external photon fields or dense clouds near the black hole are needed; the jet's own charge-neutrality cold protons provide the target, which fits the lack of broad-line emission in BL Lacs.","The required maximum proton energy is about 10 PeV in the jet frame (one order below the cosmic-ray ankle, ~2e17 eV in the observer frame), so the source is a plausible but not extreme cosmic-ray accelerator.","Upcoming Cherenkov Telescope Array and LHAASO observations of similar flares should detect gamma rays up to ~100 TeV, providing a direct test of the pp + cascade spectrum.","The same hadronic mechanism can explain the quiescent gamma-ray state with a softer proton spectrum and lower maximum energy, producing a smaller but nonzero neutrino rate."],"supporting_citations":[{"why":"Reports the IceCube-170922A event and its spatial association with TXS 0506+056, the observation the model must explain.","marker":"[1]"},{"why":"Provides the variability timescale, source distance, and SED data used to fix the blob size and Doppler factor.","marker":"[5]"},{"why":"The cloud-in-jet pp model the paper contrasts with, which requires external dense gas clouds.","marker":"[9]"},{"why":"Defines the proton blazar model framework of relativistic protons, electrons, and cold protons used throughout.","marker":"[19, 20]"},{"why":"Supplies the secondary-particle spectra for pp collisions used to compute gamma-ray and neutrino emissivities.","marker":"[37]"},{"why":"Justifies the electron acceleration efficiency $\\chi_e\\approx 10^{-3}$ that sets the cold proton density through charge neutrality.","marker":"[25–27]"}],"fun_headline_variants":["Charge neutrality sets target protons for blazar neutrino and gamma flares","Cold protons from charge neutrality produce blazar's neutrino signal","Same protons make gamma rays and IceCube's neutrino from blazar","Neutrino and gamma from one cold proton population in blazar","No extra parameters: cold protons account for blazar's neutrino and gamma"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on assuming the electron acceleration efficiency is $\\chi_e\\approx 10^{-3}$; if electrons are accelerated more efficiently, or if the jet is pair dominated, the cold proton target density drops sharply and matching the gamma-ray flux would require an implausibly high proton power.","fun_headline_variants_meta":{"raw":{"variants":["Charge neutrality sets target protons for blazar neutrino and gamma flares","Cold protons from charge neutrality produce blazar's neutrino signal","Same protons make gamma rays and IceCube's neutrino from blazar","Neutrino and gamma from one cold proton population in blazar","No extra parameters: cold protons account for blazar's neutrino and gamma"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000958,"raw_usage":{"total_tokens":4102,"prompt_tokens":983,"completion_tokens":3119,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":3030}},"tokens_in":599,"tokens_out":3119,"duration_ms":19699,"temperature":1.0,"reasoning_tokens":3030,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:06:08.129389+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be a simultaneous gamma-ray and neutrino observation of a similar BL Lac flare: if the measured muon-neutrino event rate in IceCube deviates from the predicted $N_{\\nu_\\mu}\\approx 1$ per 0.5 yr by more than the uncertainties for a given Fermi-LAT/MAGIC spectrum, the charge-neutrality target density would be ruled out. A measurement of $\\chi_e$ for this source from the synchrotron-to-Compton luminosity ratio would settle the key assumption directly.","supporting_citations":[],"review_version":1}