{"id":"6bc1fd4e-bc68-4892-869f-a2b02dab735e","arxiv_id":"2505.19784","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A proton-synchrotron model with a magnetized, parabolic jet naturally reproduces the low X-ray polarization and high optical polarization of BL Lac.","lead":"This paper shows that the large optical polarization and weak X-ray polarization of the blazar BL Lac can be explained by a model in which X-rays come from proton synchrotron radiation, not from electron inverse Compton scattering. If correct, the IXPE polarimetric data do not rule out a hadronic origin for the X-ray emission.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Adiabatic cooling length (~2 z0) is shorter than the Δz=3–4 z0 needed: at physical Δz, constant-density cases exceed the IXPE 7% limit.","rationale":"The paper's strongest claim is that the IXPE polarization limit does not rule out a hadronic proton-synchrotron origin for BL Lac's X-rays. I looked for the least secure condition needed for that claim to hold. The mechanism that lowers ΠX below 7% is the combination of a hard proton spectrum (pp=2.3) and a spatially extended proton emission region. The extended region is the more fragile ingredient: Appendix B gives l_adiab=3z/(2q), which is ~2.1–2.5 z0 for the q=0.6–0.7 values that are otherwise favored, while Table 1 and Fig. 2 only reach ΠX<0.07 for constant-density cases at Δz=3–4 z0. The paper itself states lp,cool=Δz, so using Δz=4 z0 is inconsistent with its own cooling estimate unless the proton energy distribution is evolved along z, which is not done in the Stokes calculation. At Δz=2 z0, which is close to the adiabatic length, the constant-density cases violate the IXPE bound. I also confirmed the reader's observation that the 'both cases' sentence overstates the tabulated results: q=0.6 with constant density gives ΠX=0.118 at Δz=4 z0, not below 7%. A self-consistent constant-magnetization case with q=0.7 and Δz≈2 z0 may still satisfy the data, so I do not think the central claim is destroyed; it is, however, supported only by a much narrower corner of parameter space than the text claims, and the physical consistency of that corner must be demonstrated. The reader's weakest assumption pointed at the extended Δz and the post-hoc selection of pp; my concern is more specific, namely that the extended Δz is not actually the adiabatic cooling length the model invokes. Hence partial agreement. The appropriate verdict remains CONDITIONAL, matching the reader's assessment, with the added condition that the authors verify the model with a self-consistent cooling profile and correct the 'both cases' statement.","tokens_in":10387,"tokens_out":8847,"duration_ms":86282,"concrete_test":"Recompute ΠX in the Table 1 setup with a self-consistent proton cooling profile: either truncate the z-integration at l_adiab=3z/(2q) or evolve the proton distribution along the jet using adiabatic losses, then recompute the Stokes parameters in Appendix D. As a direct check of Table 1, rerun the q=0.6, constant-density case at Δz=4 z0; if ΠX remains above 0.07, the 'both cases' claim in Sect. 2.3 is internally contradicted. If the self-consistent computation still gives ΠX<0.07 for the constant-magnetization q=0.7 case at Δz≈2 z0, the hadronic scenario survives in a narrower parameter space; otherwise the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The low X-ray polarization that carries the central claim is obtained for proton emission regions Δz≥3–4 z0 (Table 1, Fig. 2). However, the adiabatic cooling length derived in Appendix B is l_adiab=3z/(2q), which for q=0.6–0.7 is only about 2.1–2.5 z0, and Sect. 2.1 explicitly identifies lp,cool=Δz with this length; Fig. 1’s caption even states Δz∼z0. At the physically allowed Δz=2 z0, the constant-density entries in Table 1 give ΠX=0.164 (q=0.6) and 0.101 (q=0.7), both above the IXPE upper limit of 0.07; only the constant-magnetization q=0.7 case passes (ΠX=0.032). The Stokes integration in Appendix D integrates a fixed power-law proton population over the whole volume without depleting protons beyond l_adiab, so the outermost ~1.5–2 z0 of the claimed 4 z0 region contributes spurious low-polarization emission. Additionally, the statement that the observed polarization can be reproduced 'in both cases ... for q=0.6−0.7 and Δz≳3' is contradicted by Table 1, where q=0.6 with constant density gives ΠX=0.118 even at Δz=4 z0. The central claim therefore rests on an internally inconsistent volume, not merely on a post-hoc choice of parameters.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits IXPE multifrequency polarimetry of BL Lac and argues that the observed low X-ray polarization (ΠX ≲ 7%) with high optical polarization (ΠO ~ 25–30%) does not rule out a hadronic, proton-synchrotron origin of the X-ray emission. The authors construct a stationary, axisymmetric, magnetically dominated MHD jet with a helical magnetic field, assume that electrons emit optical synchrotron radiation from a thin shell while protons emit X-rays from an extended region, and compute Stokes parameters. They claim that for jet shape parameter q = 0.6–0.7 and proton emission length Δz ≳ 3 z0, the model reproduces the observed polarization in both a constant-density and a constant-magnetization prescription. The paper concludes that a hadronic model remains viable.","tokens_in":10664,"tokens_out":8070,"duration_ms":85300,"significance":"If the calculation were self-consistent, this paper would be an important counterpoint to the widely cited conclusion that IXPE observations of BL Lac favor a leptonic origin of its X-ray emission. The MHD field setup, the treatment of synchrotron polarization, and the resulting Table 1 are transparent and allow the reader to check the parameter dependence; the appendices provide useful derivations of the adiabatic cooling length, jet power, and Stokes integrals. However, the central claim is presently carried by an internally inconsistent choice of the proton emission volume and by a misstatement of the results in Table 1. After correcting these issues, only a narrow parameter subset survives, so the paper is a useful contribution in need of substantial revision rather than a definitive demonstration that hadronic models are not ruled out.","major_comments":[{"comment":"The proton emission lengths Δz = 3–4 z0 adopted to obtain the low X-ray polarization are longer than the adiabatic cooling length derived in the paper. Appendix B gives l_adiab = 3z/(2q), which is 2.5 z0 for q = 0.6 and 2.1 z0 for q = 0.7, and Sect. 2.1 identifies the proton cooling length with this scale; Fig. 1's caption even states Δz ∼ z0. At the physically allowed Δz = 2 z0, the constant-density entries in Table 1 give ΠX = 0.164 (q = 0.6) and 0.101 (q = 0.7), both above the IXPE upper limit of 0.07; among the q = 0.6–0.7 cases highlighted in the text, only the constant-magnetization q = 0.7 entry passes (ΠX = 0.032). The claim that the low X-ray polarization is naturally produced by the larger proton cooling volume is therefore not supported for the parameter range emphasized in the text.","section":"Sect. 2.1, Appendix B, Table 1"},{"comment":"The statement that the observed polarization can be reproduced 'in both cases ... for q = 0.6−0.7 and Δz ≳ 3' is contradicted by Table 1. For q = 0.6 with constant density, Δz = 4 z0 gives ΠX = 0.118, which exceeds the IXPE upper limit; for q = 0.6 with constant magnetization, Δz = 2 z0 gives ΠX = 0.112, also above 0.07. Only a subset of the quoted parameter combinations actually satisfies ΠX < 0.07. The text, Table 1, and Fig. 2 must be brought into agreement, and the claimed parameter range must be corrected.","section":"Sect. 2.3, Table 1"},{"comment":"The Stokes integrals in Appendix D integrate a fixed power-law proton distribution over the entire volume z0 < z < z0 + Δz, with no depletion or spectral evolution of the protons due to adiabatic losses. Since the adiabatic cooling length is only about 2 z0, the outer part of a 4 z0 emission region should have a lower normalization and a high-energy cutoff; its contribution to the computed low polarization is therefore spurious. The calculation should either restrict Δz to the cooling length or include a distance-dependent proton distribution that accounts for adiabatic cooling.","section":"Appendix D"},{"comment":"The proton spectral slope pp = 2.3 is adopted from the observed X-ray spectrum (Agudo et al. 2025), and the emission region size Δz is chosen after the fact to reach ΠX < 7%. Consequently, the low X-ray polarization is a consistency check rather than an independent prediction of the hadronic model. The paper would be strengthened by a sensitivity study showing how ΠX varies with pp and by an explicit statement that pp and Δz are tuned to the observations.","section":"Sect. 2.3"}],"minor_comments":[{"comment":"The caption says protons are advected 'up to a distance Δz ∼ z0', which is inconsistent with the values Δz = 3–4 z0 used in the calculations; please reconcile the sketch and caption with the quantitative model.","section":"Fig. 1 caption"},{"comment":"The header is confusing: the first numerical column (labeled Δz → 0) contains ΠO, while the following columns contain ΠX, but this is not stated explicitly. Please label the columns as ΠO and ΠX to avoid ambiguity.","section":"Table 1"},{"comment":"The function A(q) in Eq. (B.5) is described as 'a factor of order unity' but no explicit expression is given; please provide the expression or a precise citation to the derivation in Bolis et al. (2024b).","section":"Appendix B"},{"comment":"The sentence 'lp,cool = Δz ≃ 3z0/2q' is at odds with both Fig. 1's caption (Δz ∼ z0) and the later use of Δz = 3–4 z0; this should be clarified in the text.","section":"Sect. 2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an timely question, but the central claim currently rests on an internally inconsistent emission volume and on a text/table mismatch. A revision that restricts Δz to the adiabatic cooling length, includes proton cooling in the Stokes integrals, and re-derives the allowed parameter region could make the hadronic interpretation viable; the current version overstates the case."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. It's a genuine counterargument to the recent IXPE-based conclusion that BL Lac's X-rays must be leptonic, and it names two physical effects, a hard proton spectrum and a larger proton cooling volume, that push X-ray polarization down. The MHD field setup is standard but clean, and the extension of Bolis et al.'s polarization machinery to proton emitters is a real step forward. The paper correctly notes that Liodakis et al. and Agudo et al. only considered radiative losses and uniform fields, so this is not a strawman.\n\nBut the quantitative match is less solid than the text suggests. The adiabatic cooling length derived in Appendix B is l_adiab = 3z/2q, about 2.1-2.5 z0 for the q values used. The paper then sets Δz = 3-4 z0 and integrates proton emission over the whole volume without depleting the population beyond the cooling length. At the physically allowed Δz = 2 z0, the constant-density cases in Table 1 give ΠX = 0.164 (q=0.6) and 0.101 (q=0.7), both above the IXPE 7% limit. Only the constant-magnetization, q=0.7 case passes. So the claim that 'both cases' work for q=0.6-0.7 and Δz≳3 is contradicted by their own Table 1: the q=0.6 constant-density row gives 0.118 even at 4 z0. The parameter selection is also post-hoc, with pp=2.3 taken from the observed X-ray slope and Δz basically chosen to land below 7%. That is not a fatal flaw in a modeling paper, but it should be labeled as a demonstration of consistency, not a prediction.\n\nNone of this kills the central point. A proton-synchrotron model with a hard spectrum and a large cooling volume can, in a corner of parameter space, reproduce the observed polarization pattern. That is enough to show the IXPE data alone do not rule out hadronic emission. The paper deserves a serious referee, but the referee should ask for three things: fixing the emission-volume inconsistency so protons are not counted beyond their cooling length; correcting the overstated 'both cases' claim; and an explicit statement of which parameters are free and how Δz is constrained. With those changes, this would be a solid contribution. For now, treat the quantitative match with caution.\n\nWho should read: people working on IXPE blazar polarimetry and hadronic jet models. I would bring it to a reading group and cite it as the strongest current hadronic counterexample.","headline":"A serious but over-claimed attempt to revive a hadronic explanation for BL Lac's X-ray polarization; the physical mechanism is real, but the quantitative match depends on a cooling-volume inconsistency and post-hoc parameters.","tokens_in":11293,"tokens_out":2846,"would_cite":true,"duration_ms":27831,"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":"A hadronic jet model with proton-synchrotron X-rays can reproduce BL Lac's high optical and low X-ray polarization, so the IXPE measurements do not rule out a hadronic origin.","keywords":["BL Lac","blazar jets","hadronic models","proton synchrotron","X-ray polarization","IXPE","magnetically dominated jets","polarization degree"],"falsifier":"Measure the X-ray polarization of BL Lac during a similar state to the IXPE campaign and see whether it stays below about $7\\%$; a value above roughly $15\\%$ would contradict the model. Independently, determine the proton injection slope from the gamma-ray SED and the shape of the high-energy hump; a slope steeper than $p_p \\approx 2.3$ would remove the parameter choice that keeps $\\Pi_X$ low.","tokens_in":10112,"feed_emoji":"🔭","tokens_out":5900,"duration_ms":58399,"temperature":0.7,"pith_summary":"The paper argues that the multifrequency polarization data for BL Lac do not force a leptonic interpretation of its X-ray emission. It constructs a hadronic model in which relativistic electrons produce the optical synchrotron component while protons, radiating through proton synchrotron, produce the X-rays and gamma-rays. In this model the X-ray polarization is naturally small, mostly because the proton spectrum is hard (slope $p_p = 2.3$) and because slow-cooling protons fill a much larger jet volume than fast-cooling electrons. The authors conclude that a proton-synchrotron origin for BL Lac's X-rays remains viable, with jet powers below the Eddington luminosity.","feed_headline":"Proton-synchrotron model survives BL Lac's polarization test","feed_subtitle":"With a hard proton spectrum and a large emission zone, a hadronic jet matches IXPE's 7% X-ray limit.","key_machinery":"The central mechanism is the polarization of synchrotron radiation computed from the MHD model of magnetically dominated, stationary, axisymmetric outflows of Lyubarsky (2009), in which the jet is parabolic ($R_0 \\propto z_0^q$) and carries a helical magnetic field. For such fields, the polarization degree depends strongly on the slope of the emitting particle energy distribution, more so than in the textbook uniform-field case, and also on the size and shape of the emission region. The optical versus X-ray contrast comes from two different cooling regimes: electrons emitting in the optical band are fast-cooling and radiate from an extremely thin shell ($\\Delta z \\to 0$), while protons are slow-cooling and are advected outward until adiabatic losses set a scale $\\Delta z \\simeq 3z_0/(2q)$, a few times $z_0$. The hard proton distribution and the larger proton-emitting volume jointly push $\\Pi_X$ below the IXPE limit, while the steeper electron distribution keeps $\\Pi_O$ high.","core_discovery":"On the paper's own terms, the low X-ray polarization upper limit of about $7\\%$ measured by IXPE does not rule out hadronic X-ray production. Using a magnetically dominated, stationary, axisymmetric jet with a helical magnetic field, the authors compute synchrotron polarization from two particle populations: electrons with spectral slope $p_e = 4.6$ radiating from a very thin shell at the acceleration site, and protons with slope $p_p = 2.3$ radiating from a region extending to $\\Delta z \\simeq 3\\text{--}4\\,z_0$ set by adiabatic cooling. The calculation yields optical polarization $\\Pi_O \\sim 25\\text{--}45\\%$ and X-ray polarization $\\Pi_X \\lesssim 7\\%$, matching the observed averages, with the electric vector position angle aligned with the jet axis in both bands as observed. The earlier conclusion that hadronic models are ruled out is attributed to single-zone models with uniform fields and radiative-only cooling, both of which this paper drops.","pith_inferences":["The parameter combination that reproduces the observations ($p_p = 2.3$ and $\\Delta z \\simeq 3\\text{--}4\\,z_0$) is selected after the IXPE results were known; an independent test would be to derive the proton injection spectrum from the gamma-ray SED and check whether it is actually as hard as 2.3.","A testable extension of the jet geometry is that an optical flaring episode produced by a small moving region on a helical path should show a large increase in optical polarization with little change in X-ray polarization; simultaneous IXPE and optical monitoring during a flare would distinguish this picture from single-zone models.","If a future MeV polarimeter measures the band below the high-energy peak and finds polarization significantly above the X-ray value, that would challenge the model's prediction that adiabatic cooling keeps the polarization uniformly low across the whole high-energy hump."],"forward_implications":["If this model is right, the IXPE polarization measurement no longer uniquely favors leptonic X-rays; proton-synchrotron X-rays remain observationally viable for BL Lac.","The standard objections to hadronic models, namely long variability timescales and too-high predicted X-ray polarization, are resolved once adiabatic cooling and an extended proton emission region are included.","Because the proton cooling length is essentially energy-independent below the SED peak, the MeV band should show a similarly low degree of polarization, while above the roughly 100 MeV peak the polarization should rise, reaching about $15\\text{--}20\\%$ in the constant-density case.","The required proton energy flux, about $3\\times 10^{45}\\,\\mathrm{erg\\,s^{-1}}$, keeps the total jet power below the Eddington luminosity of the estimated black hole mass, so the hadronic scenario is not energetically extreme.","The model predicts a very low neutrino flux, well below current detector sensitivities, because photomeson production is negligible in the adopted parameter regime."],"supporting_citations":[{"why":"Supplies the multifrequency optical and X-ray polarization measurements of BL Lac and the leptonic interpretation that this paper re-examines.","marker":"Agudo et al. (2025)"},{"why":"Provides the polarization formalism for synchrotron radiation from magnetically dominated jets on which the present calculation is built.","marker":"Bolis et al. (2024a)"},{"why":"Gives the MHD solution for magnetically dominated stationary axisymmetric outflows used to set the jet shape and helical magnetic field.","marker":"Lyubarsky (2009)"},{"why":"Supplies the proton-synchrotron frequency-energy relation used to fix the proton Lorentz factors needed for the X-ray and gamma-ray bands.","marker":"Aharonian (2000)"},{"why":"Is the companion modeling study concluding that inverse Compton emission dominates the X-rays; its hadronic scenarios are the ones this paper shows are not uniquely disfavored.","marker":"Liodakis et al. (2025)"},{"why":"Establishes the earlier single-zone expectation that hadronic X-rays should be highly polarized, which the extended emission region of this paper revises.","marker":"Zhang & Böttcher (2013)"},{"why":"Used to estimate proton-photon cooling and the expected neutrino flux, finding both negligible in the adopted regime.","marker":"Murase et al. (2012)"}],"fun_headline_variants":["BL Lac's low X-ray polarization does not rule out protons","Proton synchrotron matches BL Lac's polarization data","Hadronic jet survives BL Lac's IXPE polarization test","Proton-synchrotron model fits BL Lac's polarization limits","Low X-ray polarization still allows hadronic BL Lac jet"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The match to the IXPE upper limit depends on adopting a hard proton spectrum with slope 2.3 and letting the protons radiate from a region roughly three to four times the distance from the black hole; if the real proton spectrum were steeper or the proton emission region smaller, the predicted X-ray polarization would rise above 7%.","fun_headline_variants_meta":{"raw":{"variants":["BL Lac's low X-ray polarization does not rule out protons","Proton synchrotron matches BL Lac's polarization data","Hadronic jet survives BL Lac's IXPE polarization test","Proton-synchrotron model fits BL Lac's polarization limits","Low X-ray polarization still allows hadronic BL Lac jet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000425,"raw_usage":{"total_tokens":2178,"prompt_tokens":941,"completion_tokens":1237,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":1152}},"tokens_in":557,"tokens_out":1237,"duration_ms":9441,"temperature":1.0,"reasoning_tokens":1152,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:06:10.845520+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the X-ray polarization of BL Lac during a similar state to the IXPE campaign and see whether it stays below about $7\\%$; a value above roughly $15\\%$ would contradict the model. Independently, determine the proton injection slope from the gamma-ray SED and the shape of the high-energy hump; a slope steeper than $p_p \\approx 2.3$ would remove the parameter choice that keeps $\\Pi_X$ low.","supporting_citations":[{"cited_title":"2009, , 698, 1570","cited_arxiv_id":null,"evidence_quote":"Gives the MHD solution for magnetically dominated stationary axisymmetric outflows used to set the jet shape and helical magnetic field."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the proton-synchrotron frequency-energy relation used to fix the proton Lorentz factors needed for the X-ray and gamma-ray bands."},{"cited_title":"& B \\\"o ttcher , M","cited_arxiv_id":null,"evidence_quote":"Establishes the earlier single-zone expectation that hadronic X-rays should be highly polarized, which the extended emission region of this paper revises."},{"cited_title":"D., Takami , H., & Migliori , G","cited_arxiv_id":null,"evidence_quote":"Used to estimate proton-photon cooling and the expected neutrino flux, finding both negligible in the adopted regime."}],"review_version":1}