{"id":"10bb5e39-a756-423b-ab21-dd1bae6e9fb0","arxiv_id":"2507.16929","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"During a neutrino-associated flare, the blazar PKS 0735+178 showed rotations in radio and optical polarization angles consistent with a propagating shock, indicating conditions for proton-photon neutrino emission.","lead":"Researchers studied the polarized radio and optical light of the blazar PKS 0735+178 around a December 2021 neutrino detection and found that the jet's magnetic field orientation rotated in a way consistent with a shock wave. The work suggests that polarization swings could serve as a new probe of neutrino production in blazar jets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Ad hoc two-component model selection is the linchpin; the three-component model favored by prior work would remove or alter the EVPA rotation underlying the shock claim.","rationale":"The reader identified the ad hoc chi2_tot<0.5 model-selection cutoff as the weakest assumption; my independent read agrees. The entire physical conclusion hinges on the two-component decomposition: the EVPA rotation of C and Q1 is what motivates the shock interpretation. If a third component is required, the EVPA time series changes, potentially eliminating the rotation or attributing it to the third component's evolution rather than a shock in the core. The paper's own citation of Kim & Kim (2025) as strengthening the case while rejecting their three-component model is internally inconsistent and underscores the fragility. A transparent model-comparison statistic would settle whether the rotation is real. Since the EVPA measurement itself is new and the authors hedge their neutrino-association assumption, conditional acceptance remains appropriate pending this model-selection check.","tokens_in":13401,"tokens_out":6510,"duration_ms":67862,"concrete_test":"Refit all six flaring-epoch VLBA datasets in the same eht-imaging framework with a three-component model and compare to the two-component model using AIC/BIC or leave-one-out cross-validation instead of the chi2_tot<0.5 cutoff. If the three-component model is preferred, or if the C/Q1 EVPA time series from the preferred model differs by more than the quoted 6 deg uncertainties or no longer shows the out-and-back rotation, the shock-front interpretation is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central shock-front inference rests on the two-component decomposition of the VLBI images. Section 3 rejects the three-component model solely because its chi2_tot falls below an uncalibrated cutoff (<0.5), while no statistically justified model comparison (AIC/BIC, Bayesian evidence, cross-validation) is reported. The adopted two-component solution has chi2_tot values 0.69-1.64, uncomfortably close to that cutoff, so the choice is not robustly established. This is not merely cosmetic: Kim & Kim (2025), using DIFMAP/CLEAN, find a third, superluminal component in the same data. If that component is real, the EVPAs assigned to 'C' and 'Q1' in Table 1 are contaminated by it, and the apparent ~70-90 deg rotation could be a consequence of that third component's flux and polarization evolution rather than a shock-induced rotation of the core. The paper even invokes the Kim & Kim component as support in Sec. 4.1 while rejecting the three-component model in Sec. 3, which creates an internal tension. Because the EVPA rotation is the only direct evidence for the shock, the shock claim is only as secure as the model selection.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes multi-epoch VLBA polarimetric imaging of the blazar PKS 0735+178 around the IceCube-211208A neutrino event of December 2021, together with Fermi-LAT gamma-ray data and optical polarimetry. Using eht-imaging geometric model fitting, the authors decompose the source into a core (C) and a jet component (Q1), and report a rotation of the radio EVPAs and an increase in fractional linear polarization during the multi-band flare, alongside an optical EVPA rotation of ~80 deg. They interpret these changes as evidence for a propagating shock front, argue that this shock could accelerate protons, and claim that the necessary conditions for proton-photon neutrino production are present. The paper explicitly acknowledges that the neutrino-blazar association is assumed, and it ends with a tentative 18-year prediction for the next neutrino association during a future quiescent phase.","tokens_in":13579,"tokens_out":4320,"duration_ms":48810,"significance":"If the central inference holds, the paper provides a rare direct structural probe of the jet conditions during a neutrino-associated blazar flare, with the novel twist that no new jet component is ejected; the shock interpretation would connect radio polarimetry to neutrino production in a concrete way. The work uses publicly available data, a less-supervised forward-modeling approach, and includes explicit caveats about the neutrino association. The 18-year prediction is falsifiable in principle. The main limitation is that the shock claim rests entirely on the two-component model selection, which is not statistically justified in the current manuscript; hence the significance is conditional on that choice being robust.","major_comments":[{"comment":"The choice of the two-component model over the three-component model is the linchpin of the shock interpretation, yet the only criterion given is an ad hoc chi2_tot < 0.5 cutoff for 'overfitting,' with no AIC/BIC, likelihood-ratio test, or cross-validation reported. The adopted two-component fits have chi2_tot = 0.69-1.64, uncomfortably close to the cutoff, so the rejection of the three-component model is not robustly demonstrated. Because the EVPAs assigned to C and Q1 would be contaminated if a third component is present, this model choice directly controls the ~70-90 deg rotation that underlies the shock claim.","section":"Section 3, Table 1"},{"comment":"There is an internal tension in the use of Kim & Kim (2025): Section 3 rejects their three-component model as overfitting in the eht-imaging framework, but Section 4.1 cites their superluminal component as 'further strengthening our case.' If the third component is real, the two-component EVPAs are not reliable; if it is not real, it cannot be used as supporting evidence. The authors should either present a statistically justified three-component analysis or remove the appeal to that component.","section":"Section 4.1"},{"comment":"The pre-flare baseline epoch MJD 59035 is separated by ~1.5 years from the neutrino-associated interval, and the authors themselves label its components Q01/Q02 because they cannot be certain of their correspondence to C and Q1. Yet the narrative that the EVPAs 'return to their initial quiescent state' treats this epoch as the baseline. In addition, Table 1 labels this epoch as '9035/2021.51' while the text states 'MJD 59035 (year 2020.35)'; this inconsistency should be corrected, and the stability of the component identification over the gap should be justified or the claim about returning to the quiescent state should be weakened.","section":"Section 3, Table 1"}],"minor_comments":[{"comment":"The first paragraph contains a typographical error: 'spurning the community' should likely be 'prompting the community.'","section":"Section 1"},{"comment":"The chi2_tot values in Table 1 would be more informative if reported as reduced chi-square with the number of degrees of freedom, and the statement that the one-component model yields 'chi2_tot >> 1' would benefit from a quantitative value.","section":"Section 3, Table 1"},{"comment":"In the bottom panel, the y-axis label 'Fractional linear polarisation' does not specify units; consider adding a percent sign or an explicit 'm' notation to match the text.","section":"Figure 2"},{"comment":"The 18-year prediction is explicitly conditional on the Britzen et al. (2010) periodicity and on the neutrino-quiescence connection, and the authors do note that further detections are not excluded; I suggest labeling this as a tentative expectation rather than a prediction to avoid overstatement.","section":"Section 4.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely multi-messenger question and the data are appropriate, but the central shock interpretation hinges on model selection that is currently defended by an uncalibrated cutoff. This is fixable within the manuscript's scope: the authors should add a proper model comparison (e.g., AIC/BIC or posterior evidence), quantify the impact of the third component on the EVPA rotation, and resolve the internal tension with Kim & Kim (2025). I also recommend checking the MJD/decimal-year inconsistency for the pre-flare epoch, as it directly affects the baseline comparison."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is a clean radio EVPA rotation in the core of PKS 0735+178 across the neutrino-associated 2021 flare, with no new component in their two-component RML fits. That is a real measurement, clearly presented with uncertainties. The paper does well to place it in the shock-in-jet framework and to report the optical rotation and fractional polarization rise. The 18-year prediction is explicitly labeled as based on Britzen's ~9-year cycle, so it is falsifiable, and the authors are honest about the working assumption that the neutrino event is associated with the blazar.\n\nThe load-bearing soft spot is exactly where the stress-test lands. Section 3 rejects the three-component model with an uncalibrated chi2_tot < 0.5 cutoff; the adopted two-component solutions sit at 0.69–1.64, so the margin is thin. No AIC/BIC, Bayesian evidence, or cross-validation is reported. The paper then cites Kim & Kim (2025) as supporting a superluminal component in Sec 4.1 while rejecting the three-component fit in Sec 3. That is an internal tension. If the third component is real, the assigned EVPAs for C and Q1 could be contaminated, and the apparent ~70–90 deg rotation could be a consequence of that component's flux and polarization evolution rather than a shock-induced core rotation. So the shock-front conclusion is only as secure as the model selection.\n\nAlso, the pre-flare baseline epoch is MJD 59035 (2020.35), far from the flare window; it is useful for illustration but not a matched baseline. The neutrino association is explicitly assumed, which is fair but limits the strength of the claim. The observation that the optical rotation likely continued beyond the last measurement is reasonable but speculative.\n\nThis paper is for the multi-messenger blazar community and VLBI polarimetry folks. It deserves a serious referee: a good referee can ask for a justified model comparison and a matched baseline without killing the core measurement. I would send it to review and would cite the EVPA rotation itself, while treating the shock interpretation as provisional.","headline":"Solid new EVPA rotation measurement, but the shock conclusion leans on a model-selection choice that needs better justification.","tokens_in":14170,"tokens_out":1221,"would_cite":true,"duration_ms":13923,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A propagating shock, not a new jet component, is the likely engine of the 2021 neutrino-associated flare in blazar PKS 0735+178.","keywords":["High energy astrophysics","Jets","Blazars","Cosmological neutrinos","Neutrino astronomy","Very long baseline interferometry","Polarisation","PKS 0735+178"],"falsifier":"A concrete observation that would settle the claim: high-cadence polarimetric VLBI of a future neutrino-coincident flare in PKS 0735+178 that resolves no EVPA rotation but instead shows a new superluminal component being ejected, or a three-component fit that matches the same epochs without overfitting, would falsify the shock-front explanation. Alternatively, if the reported neutrino event's association is shown to be due to chance—the event carries only a 30 percent probability of being astrophysical—the claimed connection would lose its anchor.","tokens_in":13175,"feed_emoji":"🔭","tokens_out":6989,"duration_ms":66259,"temperature":0.7,"pith_summary":"The paper claims that the December 2021 multi-band flare of the blazar PKS 0735+178, which coincided with a reported high-energy neutrino event, was caused by a shock front propagating through the jet rather than by the ejection of a new jet component. Using very-long-baseline interferometry images in polarised light, the authors track the electric vector position angle (EVPA) in the core and jet across six epochs; the radio EVPA rotates from parallel to the jet axis to perpendicular and back, while fractional linear polarisation rises and falls. This pattern matches a shock travelling through a stationary recollimation shock. The authors conclude that such a shock can accelerate protons and, with target photons supplied either by the surrounding medium or by proton synchrotron radiation, the necessary conditions for proton–photon neutrino production are present in this jet.","feed_headline":"Shock, not new blob, powered neutrino-linked blazar flare","feed_subtitle":"Radio polarisation maps trace the 2021 flare of PKS 0735+178 to a shock wave inside its jet.","key_machinery":"The central mechanism is the shock front propagating through the jet, diagnosed through its polarisation signature: a rotation of the electric vector position angle (EVPA, the orientation of the electric field of the radio emission) accompanied by a rise in fractional linear polarisation. That signature is read through a two-component Gaussian model (a stationary core 'C' and a quasi-stationary jet component 'Q1') fitted to VLBI images with regularised maximum-likelihood imaging; the model choice is what isolates the EVPA evolution. The shock-in-jet framework then connects the observed disturbance to proton acceleration and to the conditions for photomeson neutrino production.","core_discovery":"On the paper's own terms, the central discovery is that polarised radio structure reveals a shock front at work in PKS 0735+178 at the time of the neutrino-associated flare. The core component's EVPA, initially aligned with the bulk jet flow, turns perpendicular to the flow as the flare develops, then returns to alignment; the jet component Q1 shows the same behaviour, and fractional linear polarisation peaks at about 8 percent in the jet shortly after the neutrino detection. No new bright component appears in the images, so the flare is attributed to a propagating disturbance rather than an ejection. The simultaneous roughly 80-degree optical EVPA rotation, which likely continued further after observations stopped, reinforces the shock interpretation. The paper ends with the conclusion that the jet satisfied the energy and target-photon conditions for neutrino emission via proton–photon interactions.","pith_inferences":["Beyond the paper's data, the same EVPA-rotation diagnostic could be applied to other neutrino-associated blazars such as TXS 0506+056, where ejection and shock scenarios are still debated.","A decisive test of the model choice would be to re-fit these epochs with a three-component model under the same regularised likelihood framework; if it fits without overfitting, the inferred EVPA evolution and shock interpretation would need revision.","Continuous optical and radio polarisation coverage through the peak of gamma-ray flares would show whether such rotations are a systematic precursor of neutrino-coincident flaring rather than a single-source coincidence.","If the roughly 9-year activity cycle holds, quiescent phases of blazars—not just flaring peaks—may be the most promising windows for future neutrino searches."],"forward_implications":["Polarisation monitoring can flag neutrino-favourable states in blazars even when no new jet component is ejected.","Shock fronts interacting with the ambient medium or with their own synchrotron emission can supply the target photons needed for proton–photon neutrino production.","The paper predicts that the next neutrino association with PKS 0735+178 should occur during its next quiescent phase, in roughly 18 years, assuming the roughly 9-year activity cycle holds.","The absence of a new VLBI component during a neutrino-coincident flare is not evidence against the neutrino being real; a quasi-stationary shock can provide the acceleration.","Radio EVPA rotations may lag the optical ones by months, giving a longer-lived, later probe of the same physical disturbance."],"supporting_citations":[{"why":"Establishes that a shock front passing through a stationary recollimation shock produces the radio EVPA rotation and polarisation increase used here as the signature.","marker":"I. Liodakis et al. (2022a)"},{"why":"Supplies the shock-in-jet model connecting a travelling shock to particle acceleration and multi-band flaring.","marker":"A. P. Marscher et al. (2002)"},{"why":"Argues that mildly relativistic shocks near the jet base provide the proton energies needed for TeV–PeV neutrinos during quasi-stationary morphology.","marker":"A. V. Plavin et al. (2021)"},{"why":"Provides the long-term monitoring showing PKS 0735+178 alternates between stationary and twisted superluminal phases, used to identify the quiescent state.","marker":"S. Britzen et al. (2010)"},{"why":"Statistically connect gamma-ray flares to optical EVPA rotations, supporting the physical interpretation of the observed optical rotation.","marker":"D. Blinov et al. (2015, 2018)"},{"why":"Reports the multi-band flare observations used to time the flare relative to the neutrino event.","marker":"A. Acharyya et al. (2023)"},{"why":"The competing three-component model whose ejection interpretation the two-component shock scenario must be weighed against.","marker":"Y.-S. Kim & J.-Y. Kim (2025)"},{"why":"Provides the proton-synchrotron scenario in which protons interacting with their own synchrotron radiation create neutrinos via photomeson interactions.","marker":"A. Mastichiadis & M. Petropoulou (2021)"},{"why":"The neutrino event IceCube-211208A that anchors the multi-messenger coincidence.","marker":"IceCube Collaboration (2021)"}],"fun_headline_variants":["Polarisation pinpoints shock as source of blazar neutrino flare","Blazar neutrinos traced to shock wave in jet, not ejection","Polarised radio maps blame shock for neutrino-linked flare","Shock front, not new blob, explains blazar's neutrino flare","Neutrino blazar flare from shock, polarisation reveals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two-component Gaussian model accurately describes the VLBI jet structure; if the competing three-component model is right, the inferred EVPA rotation and the shock interpretation would change materially.","fun_headline_variants_meta":{"raw":{"variants":["Polarisation pinpoints shock as source of blazar neutrino flare","Blazar neutrinos traced to shock wave in jet, not ejection","Polarised radio maps blame shock for neutrino-linked flare","Shock front, not new blob, explains blazar's neutrino flare","Neutrino blazar flare from shock, polarisation reveals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000878,"raw_usage":{"total_tokens":3778,"prompt_tokens":906,"completion_tokens":2872,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":2783}},"tokens_in":522,"tokens_out":2872,"duration_ms":21969,"temperature":1.0,"reasoning_tokens":2783,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:59:57.196030+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete observation that would settle the claim: high-cadence polarimetric VLBI of a future neutrino-coincident flare in PKS 0735+178 that resolves no EVPA rotation but instead shows a new superluminal component being ejected, or a three-component fit that matches the same epochs without overfitting, would falsify the shock-front explanation. Alternatively, if the reported neutrino event's association is shown to be due to chance—the event carries only a 30 percent probability of being astrophysical—the claimed connection would lose its anchor.","supporting_citations":[{"cited_title":"P., Jorstad, S","cited_arxiv_id":null,"evidence_quote":"Supplies the shock-in-jet model connecting a travelling shock to particle acceleration and multi-band flaring."},{"cited_title":"P., et al","cited_arxiv_id":null,"evidence_quote":"Provides the long-term monitoring showing PKS 0735+178 alternates between stationary and twisted superluminal phases, used to identify the quiescent state."},{"cited_title":"The dynamics of the parsec-scale jet in the neutrino blazar PKS 0735+178","cited_arxiv_id":"2505.13876","evidence_quote":"The competing three-component model whose ejection interpretation the two-component shock scenario must be weighed against."},{"cited_title":"2021, ApJ, 906, 131, doi: 10.3847/1538-4357/abc952","cited_arxiv_id":null,"evidence_quote":"Provides the proton-synchrotron scenario in which protons interacting with their own synchrotron radiation create neutrinos via photomeson interactions."}],"review_version":1}