{"id":"d5d2b961-c590-4ab2-9c08-65169dc1e42c","arxiv_id":"2606.27435","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"VLBI observations of TXS 0506+056 show a spine-sheath jet with aligned inner and perpendicular outer EVPAs plus a new superluminal component near the IceCube neutrino detection.","lead":"This paper uses 15 GHz VLBI polarimetry from 2009-2025 to map the jet in blazar TXS 0506+056, identifying a spine-sheath structure and a new component ejected near the 2017 neutrino event. A smart generalist might read it to see how layered jet geometry could connect radio observations to high-energy neutrino production.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"EVPA-to-B-field mapping in the stacked map assumes negligible differential Faraday rotation at 15 GHz","rationale":"The reader's weakest assumption (EVPA patterns arising from geometry rather than propagation) is exactly the load-bearing step. Because the observations are single-frequency, that assumption remains untested even in the full manuscript. The concern is therefore unchanged from the abstract-only review, but now has a concrete observational test that would decide it.","tokens_in":1861,"tokens_out":363,"duration_ms":15793,"concrete_test":"Re-observe the source at 8 GHz and 22 GHz with the same VLBI array and epochs; construct an RM map across the jet. If the RM-corrected EVPAs retain the parallel/perpendicular stratification reported in the 15 GHz stack, the geometric claim is supported; if the corrected pattern collapses to a single orientation, the original interpretation is propagation-dominated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim rests on the stacked 15 GHz polarization map showing an inner spine (EVPA parallel to jet) and outer sheath (EVPA perpendicular). At a single frequency, any line-of-sight RM gradient or internal Faraday rotation can rotate the observed EVPA by tens of degrees without altering the underlying B-field geometry. The paper reports no multi-frequency data, no RM synthesis, and no explicit test that the perpendicular sheath component survives after correcting for plausible RM values (10^3–10^4 rad m^{-2} are common in blazar cores). If the observed 90° flip is produced by differential rotation rather than a toroidal sheath field, the spine-sheath interpretation is no longer required by the data. The kinematic and variability arguments are secondary and do not constrain the propagation contribution.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper analyzes 15 GHz VLBA full-polarization observations of TXS 0506+056 from 2009–2025. Kinematic modeling identifies moderate superluminal speeds (∼1–2 c), two quasi-stationary components, and a new component ejected near the 2017 neutrino event. The stacked polarization map is presented as evidence for a stratified spine-sheath structure (inner spine with EVPA parallel to the jet axis, outer sheath with perpendicular EVPA). Variability in total intensity and polarization is interpreted as indicating interactions between the layers, consistent with prior spine-sheath models invoked to explain the neutrino emission; the source is proposed as archetypal for similar signatures in other neutrino blazars.","tokens_in":2007,"tokens_out":577,"duration_ms":31434,"significance":"If the polarization interpretation survives scrutiny of propagation effects, the multi-epoch dataset supplies direct observational constraints on jet stratification in a neutrino-associated blazar, reinforcing the spine-sheath scenario previously used to address the Doppler crisis in TeV sources. The timing coincidence between component ejection and the IceCube event adds a temporal dimension to the structural claim.","major_comments":[{"comment":"The central claim that the stacked polarization map demonstrates a spine-sheath geometry rests on mapping the observed EVPA orientations directly to the underlying B-field (parallel in spine, toroidal in sheath). At a single frequency of 15 GHz, differential Faraday rotation (RM gradients of 10^3–10^4 rad m^{-2} typical in blazar cores) can rotate the observed EVPA by tens of degrees without any change in the intrinsic field geometry, potentially producing the reported 90° flip. No multi-frequency observations, RM synthesis, or explicit test against plausible RM values are reported to exclude this alternative.","section":"stacked polarization map and its interpretation"},{"comment":"The kinematic analysis reports moderate superluminal speeds and a new component ejected contemporaneously with the 2017 neutrino event, but the paper does not quantify the uncertainty in component identification criteria or proper-motion fitting that would be required to establish the ejection timing as statistically robust rather than coincidental.","section":"kinematic analysis"}],"minor_comments":[{"comment":"The abstract refers to 'characteristic linear polarization flares' associated with EVPA rotations; the main text should define the quantitative criteria used to identify these flares and report their measured amplitudes and timescales.","section":"variability discussion"},{"comment":"Notation for EVPA and jet position angle should be stated explicitly once in the methods or results section to avoid ambiguity when comparing spine and sheath regions.","section":"polarization results"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive and detailed report. We address each major comment below and indicate the revisions we will make to the manuscript.","responses":[{"response":"We acknowledge that single-frequency 15 GHz data alone cannot exclude differential Faraday rotation as a contributor to the observed EVPA pattern. The stacked map is formed from multiple epochs spanning 2009–2025 and shows a spatially consistent orientation, which we interpret as favoring an intrinsic spine-sheath geometry over transient propagation effects. Nevertheless, we agree that this alternative cannot be ruled out without additional data. We will revise the manuscript to add an explicit discussion of possible RM effects, cite typical RM values in blazar cores, and state that multi-frequency observations would be required to confirm the magnetic-field interpretation.","revision_made":"yes","referee_comment":"[stacked polarization map and its interpretation] The central claim that the stacked polarization map demonstrates a spine-sheath geometry rests on mapping the observed EVPA orientations directly to the underlying B-field (parallel in spine, toroidal in sheath). At a single frequency of 15 GHz, differential Faraday rotation (RM gradients of 10^3–10^4 rad m^{-2} typical in blazar cores) can rotate the observed EVPA by tens of degrees without any change in the intrinsic field geometry, potentially producing the reported 90° flip. No multi-frequency observations, RM synthesis, or explicit test against plausible RM values are reported to exclude this alternative."},{"response":"We agree that a quantitative assessment of uncertainties is necessary to support the claimed ejection timing. In the revised version we will expand the kinematic section to describe the component identification criteria in detail, report the formal uncertainties from the model fits, and include a statistical evaluation of the ejection epoch relative to the 2017 neutrino event.","revision_made":"yes","referee_comment":"[kinematic analysis] The kinematic analysis reports moderate superluminal speeds and a new component ejected contemporaneously with the 2017 neutrino event, but the paper does not quantify the uncertainty in component identification criteria or proper-motion fitting that would be required to establish the ejection timing as statistically robust rather than coincidental."}],"tokens_in":1589,"tokens_out":470,"duration_ms":22800,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main thing here is new 15 GHz VLBI monitoring of TXS 0506+056 through 2025. The paper identifies a new jet component ejected near the 2017 IceCube neutrino detection and presents the first stacked polarization map for this source.\n\nThe kinematics section is straightforward and useful. It reports moderate superluminal speeds of 1-2c, two quasi-stationary components, and the timing match between the new component and the neutrino flare. That timing stands on its own as an observational fact worth recording. The description of total intensity and polarization variability, including flares tied to EVPA rotations in the stationary features, adds detail on how the jet behaves over time.\n\nThe stacked map is the clearest addition. It shows an inner region with EVPA aligned to the jet direction and an outer layer with perpendicular EVPA, which the authors link to a spine-sheath geometry. This fits with earlier suggestions that such layering could help explain neutrino production and the Doppler issues in TeV blazars.\n\nThe soft spot is the single-frequency limitation. At 15 GHz alone, differential Faraday rotation or line-of-sight RM gradients can rotate the observed EVPA by tens of degrees without any change in the underlying B-field. The paper does not report multi-frequency observations, RM synthesis, or tests against typical blazar RM values, so the perpendicular sheath signature could arise from propagation rather than a toroidal field. The kinematic and timing results do not resolve that ambiguity.\n\nThis work is for people following blazar jets and high-energy neutrino associations. The new measurements on this specific source give it value even with the polarization caveat. It deserves peer review because it supplies concrete extended data on a source already under scrutiny.","headline":"New VLBI data on TXS 0506+056 includes a component timed with the 2017 neutrino event and a stacked polarization map, but the spine-sheath claim at single frequency leaves the EVPA pattern open to Faraday rotation effects.","tokens_in":2602,"tokens_out":445,"would_cite":false,"duration_ms":32008,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"VLBI observations reveal a spine-sheath structure in the jet of the neutrino blazar TXS 0506+056.","keywords":["TXS 0506+056","spine-sheath jet","VLBI polarimetry","neutrino association","EVPA variability","jet kinematics","blazar jet structure"],"falsifier":"If higher-frequency VLBI observations show EVPA patterns that do not match the expected perpendicular sheath or if Faraday rotation measures vary strongly across the jet in a way inconsistent with the layers.","tokens_in":2760,"feed_emoji":"🔭","tokens_out":633,"duration_ms":31372,"temperature":0.7,"pith_summary":"The paper analyzes 15 GHz full-polarization VLBI data spanning 2009 to 2025 to study the parsec-scale jet of TXS 0506+056, the blazar linked to a 2017 neutrino detection. Kinematic analysis shows moderate superluminal speeds of 1-2c, two quasi-stationary components, and a new component ejected around the neutrino event. The key finding is that the stacked polarization map matches a stratified jet with an inner spine having EVPA parallel to the jet and a sheath layer with perpendicular EVPA. This configuration aligns with models attributing the neutrino emission to interactions in such a layered jet. Readers would care because it provides direct radio evidence for the jet structure thought to produce high-energy neutrinos in blazars.","feed_headline":"Spine-sheath jet mapped in neutrino blazar","feed_subtitle":"Polarization data from 15 GHz VLBI show inner aligned and outer perpendicular layers consistent with neutrino production models.","key_machinery":"The spine-sheath jet structure identified through the stacked polarization map, where the inner spine shows EVPA aligned with the jet direction and the outer sheath shows perpendicular EVPA, which accounts for the observed polarization patterns and flares.","core_discovery":"The stacked polarization map is consistent with a stratified spine-sheath jet structure: an inner spine with EVPA aligned with the jet, surrounded by a sheath layer with perpendicular EVPA. Variability in total intensity and polarization indicates interaction between these layers, and the multi-layered jet is consistent with models explaining the neutrino emission.","pith_inferences":["If the spine-sheath model holds, targeted polarization monitoring of other IceCube-associated blazars could identify candidates with similar layered jets.","The quasi-stationary components likely mark sites of layer interaction where particle acceleration for neutrinos occurs.","Higher-resolution or multi-frequency observations could map the radial extent of the sheath layer to test the stratification depth."],"forward_implications":["The interaction between spine and sheath produces linear polarization flares associated with EVPA rotations in quasi-stationary components.","A new jet component was ejected contemporaneously with the 2017 neutrino event.","The configuration explains the neutrino emission via the spine-sheath structure.","TXS 0506+056 represents an archetypal case where similar polarization signatures may appear in other neutrino-emitting sources.","This may offer a solution to the Doppler-crisis phenomenon in TeV-emitting blazars."],"fun_headline_variants":["Spine-sheath structure in TXS 0506+056 jet","Layered EVPA in neutrino blazar jet","VLBI polarimetry of TXS 0506+056 spine-sheath","Aligned spine with perpendicular sheath"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The EVPA orientations and variability are caused by the spine-sheath geometry instead of projection effects, Faraday rotation, or other propagation effects.","fun_headline_variants_meta":{"raw":{"variants":["Spine-sheath structure in TXS 0506+056 jet","Layered EVPA in neutrino blazar jet","VLBI polarimetry of TXS 0506+056 spine-sheath","Aligned spine with perpendicular sheath"]},"model":"grok-4.3","cost_usd":0.009031,"raw_usage":{"total_tokens":4106,"prompt_tokens":772,"num_sources_used":0,"completion_tokens":64,"cost_in_usd_ticks":90312000,"prompt_tokens_details":{"text_tokens":772,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3270,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":772,"tokens_out":64,"duration_ms":36920,"temperature":1.0,"reasoning_tokens":3270,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-29T01:22:58.801453+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"If higher-frequency VLBI observations show EVPA patterns that do not match the expected perpendicular sheath or if Faraday rotation measures vary strongly across the jet in a way inconsistent with the layers.","supporting_citations":[],"review_version":1}