{"id":"96d70543-edc2-4952-aa3a-3ec5861c71eb","arxiv_id":"2608.10067","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The highest optical polarization flare ever seen from a blazar is explained by a sweeping, Doppler-boosted helical jet component, and the jet appears proton-dominated.","lead":"This paper argues that the record-breaking optical polarization flare of the blazar BL Lacertae in 2023 was caused by a sweeping, helical motion of its jet that temporarily aligned the jet with our line of sight. Using VLBI monitoring and relativistic magnetohydrodynamic simulations, the authors also infer that the jet is dominated by protons rather than positrons.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Q2–optical co-spatiality assumption is load-bearing, yet the EVPA offset and the absence of a VLBI epoch at the flare peak leave it unverified; the simulations cannot confirm it.","rationale":"The central claim requires that the radio component Q2 is the same physical region that produced the optical polarization flare. The paper explicitly assumes this in Sect. 2, but the optical data cannot resolve the emission region, and the VLBI epochs bracketing the flare peak do not include it. The only direct link offered is a claimed similarity between the Q2 EVPA rotation and the optical and radio EVPA rotation, but the absolute EVPAs are offset by about 50 degrees, and no Faraday-rotation correction or quantitative association test is provided. This matters because if the flare originated in a more compact region, neither the geometric and Doppler explanation nor the n>=100 baryon-loading inference would follow. The simulation confirmation is a secondary concern: even if the simulations were fully predictive, they illustrate how a rotating jet can produce similar polarization patterns but cannot identify the unresolved optical source. I therefore agree with the reader's weakest assumption. A Faraday-corrected EVPA comparison, or a contemporaneous VLBI observation at the flare peak, would settle the association. Without such a test, the paper should be read as presenting a plausible geometric scenario, not a confirmed explanation. The CONDITIONAL verdict is appropriate; I would not reject the paper, but the confirmation language in the abstract and Sect. 3 should be softened.","tokens_in":14583,"tokens_out":8372,"duration_ms":88162,"concrete_test":"Using the existing 43 GHz VLBA data, plus any contemporaneous 15 or 86 GHz observations from BEAM-ME or GMVA, compute the Faraday rotation measure toward Q2 and correct its EVPA at each epoch. If the corrected Q2 EVPA at the optical flare peak, or its interpolation, does not agree with the optical EVPA within the combined uncertainty, the assumed co-spatiality is falsified and the paper should be reframed as a consistency argument rather than a confirmation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sect. 2 states that the flare in P was caused by the sweeping motion of the Q2 component, and this assumption is the load-bearing step. The optical data have no spatial resolution, so the only evidence connecting the optical flare to Q2 is the claim that Q2's EVPA rotates smoothly with time, similar to the radio and optical measurements. Yet in Fig. 2 and Table B.1 the absolute EVPAs are offset by roughly 50 degrees: Q2 goes from 134 to 108 degrees, optical from about 70 to 30 degrees, and radio from about 86 to 53 degrees. The paper does not apply a Faraday-rotation correction to Q2 or present a quantitative association test. In addition, no VLBI epoch falls at the optical flare peak of MJD 60261-60262; Q2's behavior at the peak is interpolated between the epochs at MJD 60254 and 60273. If the optical flare came from a more compact, unresolved region, the geometric and Doppler explanation for the flare, and the derived baryon loading n>=100, would not follow. The RMHD simulations model a rotating jet and can be made to reproduce similar polarization patterns, but they cannot establish that the unresolved optical emission was co-spatial with Q2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper interprets multi-epoch 43 GHz VLBI observations of BL Lacertae taken during the record optical polarization flare of MJD 60261-60262.75. It argues that the sweeping/helical motion of jet component Q2 changed the viewing angle and Doppler factor, and that this geometric effect, rather than an intrinsic energy release, produced the flare and the smooth rotation of the polarization angle. An analytic estimate yields a viewing-angle change of about 15 degrees, and 3D RMHD simulations with a tuned sweep angle and proton-to-positron ratio are presented as confirmation. As a by-product, the simulations are used to infer a proton-dominated jet (n >= 100), and the paper discusses the possible neutrino implications.","tokens_in":14824,"tokens_out":4324,"duration_ms":41866,"significance":"If the interpretation holds, the result would establish that the most extreme optical polarization flare ever observed in a blazar was largely a geometric and relativistic-beaming phenomenon, and it would support substantial baryon loading in the jet of BL Lac. The paper's strengths are the coherent VLBI kinematics of Q2, the transparent analytic Doppler estimate, the use of publicly available multi-wavelength data, and the honest appendix-level treatment of neutrino inefficiency. However, the central conclusion rests on an unverified co-spatiality assumption and on simulation parameters that were tuned to match the observations, so the claim of confirmation is currently overstated.","major_comments":[{"comment":"The co-spatiality of the optical flare with VLBI component Q2 is the load-bearing assumption, but it is only asserted, not tested. The text states 'Assuming that the flare in P was caused by the sweeping motion of the Q2 component' and motivates this by the similarity of EVPA rotations. However, the absolute EVPA values in Fig. 2 and Table B.1 are offset by roughly 50 degrees (Q2 rotates from about 134 to 108 degrees, optical from about 70 to 30 degrees, radio from about 86 to 53 degrees), and no Faraday-rotation correction or quantitative association test is presented. In addition, no VLBI epoch falls at the optical peak on MJD 60261-60262.75; Q2's behavior is interpolated between MJD 60254 and 60273. If the optical flare arose in a more compact, unresolved region, the geometric explanation and the derived baryon loading would not follow. The authors should either provide a quantitative test of the Q2-optical association or clearly present the interpretation as assumption-dependent.","section":"Sect. 2, Fig. 2, Table B.1"},{"comment":"The simulation 'confirmation' is partly circular. The sweep angle Delta_theta ~ 7 degrees and the proton-to-positron ratio n >= 100 were selected by matching the simulated Stokes I and psi to the observed values (Figs. D.1, D.2 and Table D.1). The paper states that the simulations 'reproduced the observations most optimally' for these parameters, but this is a fit, not an independent confirmation. Because the same parameters also control the inferred baryon loading, the statement in the abstract that the flare is 'as also confirmed by our' RMHD simulations is too strong. The authors should either provide independent constraints on the sweep angle and composition, or explicitly reframe the simulation results as a consistency check rather than a confirmation.","section":"Sect. 3 and Appendix D"},{"comment":"The agreement between observations and simulations is presented qualitatively. In the right panel of Fig. 2, the VLBI and simulated polarized flux densities are rescaled by arbitrary factors (x5e-2 and x1e-3) to compare shapes, and the text calls the agreement 'excellent' without a quantitative metric. A reduced chi-square, a cross-correlation coefficient, or an equivalent statistic should be provided for both P and psi, or the claim of agreement should be explicitly stated as qualitative.","section":"Fig. 2 and Sect. 3"},{"comment":"The neutrino discussion mixes two separate inferences. The equipartition-based proton luminosity L_p in Eq. (E.1) is stated to be independent of n, yet the text says that the baryon loading inferred in Appendix D is 'a necessary but not sufficient condition' for neutrino emission. Since the non-thermal proton population is assumed to carry the entire kinetic energy independently of n, the value of n itself does not obviously enter the neutrino estimate. This logical connection should be clarified, and any statement that n >= 100 supports neutrino production should be justified more carefully.","section":"Appendix E"}],"minor_comments":[{"comment":"The symbol m_opt is used without prior definition, and the spacing in 'm opt =47.6%' should be corrected.","section":"Abstract"},{"comment":"The sentence 'the size of Θ FWHM;c is the speed of light' is garbled and should be rewritten; the light-crossing-time definition needs a proper formula.","section":"Sect. 3"},{"comment":"The phrase 'rationn' should be 'ratio n'.","section":"Appendix D"},{"comment":"The reference to 'Eq. C.2' for the dependence of the Faraday rotation coefficient on plasma composition is incorrect; Eq. (C.2) defines psi, while the Faraday rotation coefficient is not written out in the appendix.","section":"Appendix C"},{"comment":"MacDonald & Nishikawa 2021a and 2021b appear to be the same publication; if they are indeed a single work, the duplicate reference should be removed.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of A&A and uses public data plus reproducible simulation methods. The main concern is not the data quality but the strength of the interpretive claims: the word 'confirmed' in the abstract overstates what a parameter-tuned simulation can establish given the unresolved optical emission. I would ask the authors to substantially temper the confirmation language and to add a quantitative or at least explicitly assumption-dependent framing of the Q2-optical association before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me give you the short version. This paper explains the record 47.6% optical polarization flare of BL Lac as a geometric effect: the jet component Q2 sweeps across our line of sight, changing the Doppler factor and rotating the EVPA. The idea is not new—helical jets and Doppler boosting have been invoked for decades—but applying it to this specific historic flare, with VLBI kinematics and 3D RMHD synthetic polarization maps, is fresh. The VLBI tracking of Q2 is careful, and the simulations do reproduce the qualitative rise and fall of Stokes I and the smooth psi rotation.\n\nThe soft spots are real, and one is load-bearing. The identification of Q2 with the optical flare region is assumed, not shown: the optical flare is unresolved, no VLBI epoch falls at the flare peak (MJD 60261-60262, bracketed by 60254 and 60273), and the EVPA of Q2 is offset by roughly 50 degrees from the optical/radio EVPA. The paper never applies a Faraday rotation correction or a quantitative association test. If the flare came from a more compact, possibly different region, the geometric explanation and the baryon-loading inference fall apart. The simulations can produce similar polarization patterns for a rotating jet, but they cannot verify co-spatiality.\n\nSecond, the simulation \"confirmation\" is partly circular. The sweep angle (Δθ ≈ 7°) and the proton-to-positron ratio n>=100 were chosen to match the observed light curve and EVPA. That makes n a fit result, not a prediction. The abstract says \"confirmed by simulations\" and \"baryon loading required\"—that is overclaiming. Also, Eq. 1 (the analytic Doppler relation) is referenced but not displayed, no code or parameter files are released, and the Fig. 2 comparison is visual, not quantified.\n\nTo be fair, the paper is honest in places: it explicitly says \"assuming that the flare in P was caused by the sweeping motion of the Q2 component,\" and the neutrino appendix is appropriately cautious, giving upper limits and noting the photomeson efficiency is low. The analytic Δθ ≈ 15° is plausible and independent of the simulations to first order.\n\nVerdict: worth a serious referee. The subject is important, the VLBI data are real, and the modeling is sophisticated. But the authors need to fix the overclaim, release their pipeline, show Eq. 1, quantify the goodness of fit, and either measure or properly caveat the Q2-optical association. If those are addressed, this could be a solid contribution. If not, it remains a suggestive consistency argument.","headline":"A careful but overclaimed geometric interpretation of BL Lac's record optical flare; the n>=100 baryon loading is a fit parameter, not a prediction, and the Q2-optical co-spatiality is assumed, not shown.","tokens_in":15416,"tokens_out":3653,"would_cite":false,"duration_ms":31444,"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":"Simulations tie BL Lac's record optical flare to a jet swing","keywords":["BL Lacertae","blazars","relativistic jets","optical polarization flare","very long baseline interferometry","relativistic magnetohydrodynamics","Doppler beaming","baryon loading"],"falsifier":"Monitor BL Lac at 43 GHz with VLBI simultaneously with optical polarimetry through the next strong flare. If the optical electric-vector position angle rotates by tens of degrees while no parsec-scale component swings, or if a polarization-angle rotation occurs with no change in the component's position angle, the sweeping-motion interpretation is ruled out for that event; if a swing is always present, the geometric mechanism is supported.","tokens_in":14363,"feed_emoji":"🔭","tokens_out":8818,"duration_ms":75931,"temperature":0.7,"pith_summary":"The paper argues that the highest optical linear-polarization flare ever recorded from a blazar—BL Lacertae's $m_{\\rm opt}=47.6\\%$ peak—was a geometric event. The flare happened, the authors claim, because a component of the jet, Q2, swept in a helical path that briefly aligned the relativistic beam with our line of sight, raising the Doppler factor and with it the observed flux and polarization. The claim matters because it would mean that extreme blazar flares need not be powered by intrinsic energy releases, and because the simulations used to test it require a jet that is overwhelmingly proton-dominated. If confirmed, BL Lac becomes a stronger candidate source for astrophysical neutrinos, although the parsec-scale region studied here is not itself an efficient neutrino factory.","feed_headline":"Simulations tie BL Lac's record optical flare to a jet swing","feed_subtitle":"The flare was geometry and beaming, and BL Lac's jet is mostly protons.","key_machinery":"The load-bearing object is the helical sweep of the parsec-scale jet component Q2: a rotating, bending relativistic flow whose instantaneous viewing angle $\\theta$ changes with time, altering the Doppler factor $\\delta=[\\Gamma(1-\\beta\\cos\\theta)]^{-1}$ and therefore the beamed flux and polarization. The argument is carried by the analytical Hagen-Thorn relation connecting flux density to viewing angle; by 3D RMHD simulations of a cylindrical, non-axisymmetric jet computed with the PLUTO solver and post-processed with RADMC-3D to produce ray-traced full-Stokes synthetic maps; and by the proton-to-positron ratio $n$, which enters through the Faraday rotation coefficient and controls whether the simulated electric-vector position angles rotate with jet orientation. In short, the mechanism is geometry plus charge-asymmetric Faraday rotation acting on a Doppler-boosted, kink-unstable jet.","core_discovery":"On the authors' own terms, the central discovery is that the sweeping helical motion of the VLBI component Q2, tracked over five 43 GHz epochs around MJD 60238–60351, reproduces the record optical polarization flare of BL Lac. During the swing the component's position angle changes from about $175^\\circ$ to $185^\\circ$, its linear-polarization angle $\\psi$ rotates smoothly, and its polarized flux rises and falls in step with the optical and radio light curves. Using the analytical flux–viewing-angle relation of Hagen-Thorn et al., the authors find that at $\\Gamma=4.5$ a viewing-angle change of $\\Delta\\theta\\approx 15^\\circ$ (projected $\\le 10.7^\\circ$) is sufficient; their 3D RMHD simulations with a jet sweep of $\\Delta\\theta\\approx 7^\\circ$ match the observed Stokes $I$ and electric-vector position-angle evolution. The same simulations constrain the plasma composition: protons must outnumber positrons by $n\\ge 100$ for the simulated $\\psi$ to rotate as observed, because the Faraday rotation coefficient depends on the charge asymmetry of the flow. The authors infer a proton kinetic luminosity $L_p\\approx 5\\times 10^{-3}\\,L_{\\rm Edd}$, comfortably sub-Eddington but unable to make the parsec-scale region an efficient neutrino producer without additional, more compact target regions.","pith_inferences":["A testable extension the paper leaves implicit: if the geometric mechanism is generic, then every future optical flare in BL Lac should be accompanied by a measured swing of a parsec-scale VLBI component; a flare with no accompanying swing would falsify the mechanism for that event.","The paper's own finding that the parsec-scale region is neutrino-inefficient suggests that any detectable neutrino from BL Lac would have to be produced closer to the black hole or through a reconnection-type acceleration channel; this is an inference, since the paper only states that such conditions would be required.","The identification of the optical flare with Q2 could be tested with simultaneous optical and 43 GHz polarimetry at sub-day cadence during a future flare: correlated swings in the radio and optical electric-vector position angles would strengthen the identification, while decorrelation would point to a more compact, unresolved emission region."],"forward_implications":["Blazar flares that show a smooth rotation of the polarization angle and a rise-and-fall in flux can be produced by geometry alone, so a peak polarization like BL Lac's does not by itself require a new injection of energy or a disordered magnetic field.","The best-matching simulations require baryon loading $n\\ge 100$, making BL Lac a hadronic-jet candidate and a plausible target for neutrino searches, even though its parsec-scale emission region is not an efficient neutrino factory.","The geometric mechanism predicts a characteristic multi-epoch signature—component swing, flux increase while the jet points toward us, and decrease while it points away—that VLBI polarimetry can search for in other flaring blazars.","Because the simulated $\\psi$ rotation only appears when protons dominate the plasma, measurements of polarization-angle variability can serve as a remote diagnostic of jet composition."],"supporting_citations":[{"why":"Reports the record $m_{\\rm opt}=47.6\\%$ optical polarization flare and the mm-radio flare context that this paper sets out to explain.","marker":"Agudo et al. 2025"},{"why":"Supplies the analytical relation used to convert flux-density changes into a viewing-angle and Doppler-factor change for the Q2 component.","marker":"Hagen-Thorn et al. 2008"},{"why":"Provides the bulk Lorentz factor $\\Gamma\\approx 4.5$ and magnetosonic Mach number used in both the analytical estimate and the simulations.","marker":"Cohen et al. 2015"},{"why":"Establishes the core identification and large-scale viewing angle, and documents the jet collimation profile at the location of Q2.","marker":"Casadio et al. 2021"},{"why":"Provides the RMHD modelling framework, electron-scaling relations, and equipartition initial conditions used for the synthetic synchrotron emission.","marker":"Kramer & MacDonald 2021"},{"why":"Justifies the fluid-only simulations by showing that adding Lagrangian particles does not change the large-scale synchrotron intensity and polarization structure.","marker":"Kramer et al. 2024"},{"why":"Is the PLUTO code used to run the 3D RMHD simulations.","marker":"Mignone et al. 2007"},{"why":"Gives the Faraday rotation coefficient whose proton-to-positron dependence drives the $n\\ge 100$ constraint.","marker":"MacDonald & Nishikawa 2021a"}],"fun_headline_variants":["Jet swing, not particle acceleration, drove BL Lac's flare","BL Lac flare explained by jet bending, not shocks","Proton-loaded jet required to match BL Lac observations","Record blazar flare traced to jet's helical sweep","RMHD shows baryons dominate BL Lac's jet"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The optical flare is assumed to originate in the same physical region as the 43 GHz VLBI component Q2, even though the optical data cannot resolve that region; if the flare came from a more compact, unseen zone, the sweeping-motion explanation would not apply.","fun_headline_variants_meta":{"raw":{"variants":["Jet swing, not particle acceleration, drove BL Lac's flare","BL Lac flare explained by jet bending, not shocks","Proton-loaded jet required to match BL Lac observations","Record blazar flare traced to jet's helical sweep","RMHD shows baryons dominate BL Lac's jet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000641,"raw_usage":{"total_tokens":2960,"prompt_tokens":967,"completion_tokens":1993,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":1925}},"tokens_in":583,"tokens_out":1993,"duration_ms":14219,"temperature":1.0,"reasoning_tokens":1925,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:14:41.727997+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Monitor BL Lac at 43 GHz with VLBI simultaneously with optical polarimetry through the next strong flare. If the optical electric-vector position angle rotates by tens of degrees while no parsec-scale component swings, or if a polarization-angle rotation occurs with no change in the component's position angle, the sweeping-motion interpretation is ruled out for that event; if a swing is always present, the geometric mechanism is supported.","supporting_citations":[],"review_version":1}