{"id":"04322444-b6ff-43cb-b216-392ec667f7ba","arxiv_id":"2608.12696","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":22,"one_line_summary":"TXS 0506+056's delayed radio flare is modeled as the jet running into three parsec-scale electron clouds, reproducing frequency-dependent radio peak delays.","lead":"Astrophysicists modeled the blazar TXS 0506+056's 2017 gamma-ray flare and the radio flare that peaked three years later as one connected event: a jet disturbance travels outward and later hits dense electron clouds a few parsecs from the black hole. The model matches the delayed and frequency-dependent radio light curves, offering a concrete way to link neutrino candidate blazars to radio structure.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix A's nested-cloud luminosity assignment is internally inconsistent: Eqs. (A4)-(A6) credit the largest-cloud zone to the smallest cloud, so the 2 GHz flare interpretation may not follow from the stated model.","rationale":"The reader's CONDITIONAL verdict is reasonable, but the single most load-bearing point is not the VLBI sheath/spine tension, which the authors acknowledge and hedge. It is the internal consistency of the very calculation that produces the radio light curves. If Appendix A's printed equations are what was implemented, the assignment of cloud luminosities to the three zones is inverted: the zone nominally containing only the largest cloud is assigned L_cloud0, and cloud 0 is also weighted as though only a fraction of it lies in the inner zone. The volume factors in (A1)-(A3) also do not match the luminosity factors in (A4)-(A6). This matters because the paper's physical story—cloud 2 is the large, optically thin region producing the late 2 GHz flare while cloud 0 is compact and self-absorbed—is read directly off those zone assignments. The proposed check is decisive: compare the code to the equations and rerun the best-fit case. If the implementation is correct and only the equations are misprinted, the concern is a documentation issue; if the implementation follows the equations, the central conclusion is unsupported. I therefore keep the reader's CONDITIONAL verdict rather than moving to REJECT, because the flaw is checkable and possibly typographical. I disagree with the reader on the weakest assumption: the geometry tension is acknowledged and secondary, whereas the Appendix A inconsistency is unacknowledged and logically prior to the interpretation.","tokens_in":20580,"tokens_out":12746,"duration_ms":130766,"concrete_test":"Inspect the code that computes dN_cloud/dt and check whether the volume/luminosity weights match Eqs. (A4)–(A6) or the physically correct nested-sphere weights: V_a: L_cloud2(1−ξ1^3); V_b: L_cloud2(ξ1^3−ξ0^3)+L_cloud1(ξ1^3−ξ0^3)/ξ1^3; V_c: L_cloud2 ξ0^3 + L_cloud1 ξ0^3/ξ1^3 + L_cloud0. Then rerun the best-fit parameters (Table II) with the corrected weights. If the 2 GHz peak shifts by more than one 167-day bin, or the 5–22 GHz double-peak morphology changes materially, the stated 2–7 pc cloud configuration and the causal-delay interpretation are not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central conclusion—that the leading blob meets three cold-electron clouds at 2–7 pc and that the largest cloud produces the late, optically thin 2 GHz flare—rests on the numerical construction in Appendix A, and that construction is not internally consistent as written. Cloud 2 is defined as the largest and cloud 0 as the smallest/most upstream (Eq. 8, Table I). In the three-zone decomposition, V_a is described as the volume containing only the largest cloud, so its injected luminosity should be L_{e,cloud2}(t)(1−ξ1^3). Equation (A4) instead writes L_{e,cloud0}(t)(1−ξ1^3). Equation (A5) uses L_{e,cloud0} in the term that should carry the cloud-2 contribution to the middle shell, and Eq. (A6) multiplies the cloud-0 contribution by ξ0^3 even though cloud 0 lies entirely inside V_c and should enter with full weight. Moreover, Eqs. (A1)–(A3) express the shell volumes as (1−ξ1)^3, (ξ1−ξ0)^3, and ξ0^3, which do not sum to unity and disagree with the ξ^3 factors used in Eqs. (A4)–(A6). If the AM3 wrapper implements the printed equations, the radio light curves attributed to the 'largest' cloud are actually controlled by the smallest cloud, and the inferred cloud configuration and the claimed causal delay do not follow from the model as written.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a time-dependent leptonic jet model for the 2017 gamma-ray flare and the delayed multi-frequency radio flare of TXS 0506+056. A sequence of expanding blobs is launched near the BLR and propagated to parsec scales; the leading blob is assumed to accelerate electrons in three stationary cold-electron clouds, producing the delayed radio flare. The model is fitted to Fermi-LAT and RATAN-600 light curves using local MCMC sampling, with a constant steady-state baseline taken from Ref. [44] plus an additional by-eye one-zone component. The authors report best-fit cloud distances of 2–7 pc and interpret the frequency-dependent radio peak delays as signatures of the cloud configuration.","tokens_in":21115,"tokens_out":6771,"duration_ms":66856,"significance":"If the modeling is correct, it would be a meaningful advance: it offers a concrete, time-dependent mechanism connecting the high-energy flare to the late radio flare and demonstrates that multi-frequency radio light curves can constrain the parsec-scale environment of neutrino-candidate blazars. The work is computationally substantial, uses publicly available Fermi-LAT and RATAN data and the public AM3 framework, and is transparent about several limitations (local sampling, neglected hadronic processes, and tension with VLBI spine-sheath results). However, the central inference is not currently quantitatively established: the cloud configuration is fitted to the same light curves it is used to explain, no goodness-of-fit or model comparison is provided, and the printed equations in Appendix A contain internal inconsistencies that affect the attribution of the 2 GHz flare to the largest cloud.","major_comments":[{"comment":"The volume decomposition is internally inconsistent. For coaxial clouds, the swept volume fractions should be 1−ξ1^3, ξ1^3−ξ0^3, and ξ0^3; the printed (1−ξ1)^3 and (ξ1−ξ0)^3 do not sum to unity (with the best-fit ξ0=0.21, ξ1=0.34 they sum to 0.298). More importantly, Eq. (A4) assigns L_e,cloud0 to the zone that should contain only cloud 2, Eq. (A5) uses L_e,cloud0 in the cloud-2 term, and Eq. (A6) gives cloud 2 full weight in V_c while cloud 0 receives only ξ0^3. As written, the radio light curves attributed to the largest cloud are actually controlled by the smallest cloud, so the claimed causal attribution of the late 2 GHz flare (Sec. IV, Fig. 4) does not follow from the stated model. The authors must correct these equations and rerun the fit.","section":"Appendix A, Eqs. (A1)–(A6)"},{"comment":"No goodness-of-fit statistic or model comparison is reported. The statement that “the data are well described” (Sec. I C) is supported only by visual inspection of Fig. 3. Because the cloud luminosity profiles in Eq. (10) are free functions fitted to the very radio light curves they are used to reproduce, the paper should provide at least a reduced chi-square or a likelihood-ratio comparison against a simpler model (e.g., a single extended injection zone or a two-cloud configuration) to establish that three clouds and the resulting causal interpretation are required by the data.","section":"Secs. III C and IV"},{"comment":"The best-fit model under-predicts the 5–500 GeV Fermi band, as the authors acknowledge in Sec. V C. Since the abstract and Conclusions claim that the model describes the gamma-ray flare and establishes a causal link to the radio flare, this under-prediction means the gamma-ray flare is not fully accounted for by the leptonic component. The claim should either be restricted to the 0.1–5 GeV band or the model should be extended to include the additional component responsible for the highest-energy band.","section":"Sec. IV and Sec. V C, Fig. 3"},{"comment":"The steady-state baseline, which is added to the time-dependent model in Eq. (11), is partly constructed by eye: the additional synchrotron component in Appendix C has parameters obtained by a by-eye fit. The inferred cloud luminosities and distances depend on what is subtracted as baseline, and a by-eye component with no uncertainty can absorb part of the radio flare. The sensitivity of the cloud parameters to the baseline choice should be quantified, for example by allowing the baseline normalization to vary within a prior in the MCMC.","section":"Appendix C"}],"minor_comments":[{"comment":"In the first branch the exponent appears with a plus sign, exp[+1/2 (…)], whereas an asymmetric Gaussian should have a minus sign in both branches. Please correct the sign and confirm that the implementation uses the intended form.","section":"Eq. (10)"},{"comment":"The text says “the two lowest-frequency bands” and “the two highest-frequency bands”; please clarify the exact frequency pairs used for α_low and α_high.","section":"Sec. II A and Fig. 1"},{"comment":"The sentence “To reduce these computational resource requirements as much as possible, requirements, we used” contains a duplicated word; also “miliparsec” in Sec. I A should be “milliparsec.”","section":"Sec. V C"},{"comment":"The introduction refers to “PKS 1420+240”, while the rest of the paper uses PKS 1424+240; please make the name consistent.","section":"Sec. I C"},{"comment":"The parameters for which fixed values were obtained by eye (e.g., T_peak and the ω asymmetry parameters) should be marked explicitly in the table; the current note refers only to the lower section generally.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The Appendix A inconsistency is serious and must be fixed before publication; the corrected equations may change which cloud drives the late 2 GHz flare. I do not see this as a reason for rejection because the modeling framework is coherent and the error is local, but the authors should also address the absence of model comparison. If the corrected light-curve decomposition changes the inferred cloud configuration, the conclusions will need to be revised accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nFirst thing to know: this is the most serious attempt yet to reproduce the frequency-dependent delayed radio flare of TXS 0506+056 — 2 GHz peaking about 3.4 years after the 2017 gamma-ray flare, higher bands peaking earlier with a double-peaked structure. It is also honest about its own limits, which makes it worth engaging even where it overreaches.\n\nWhat is genuinely new is the three-cloud, multi-zone parameterization. The competing extended-jet model (ref. 67) predicts a simple simultaneous radio peak; this model gets the frequency-dependent peak delays by having the jet cross a set of embedded clouds of different transverse sizes, changing the self-absorption between zones. The Fermi-LAT reprocessing is standard and competent, and Sec. V C accurately admits the under-prediction in the 5-500 GeV band, the local nature of the MCMC, and the tension with the sheath/spine VLBI picture.\n\nThe soft spots, in proportion:\n\n1. Load-bearing, and real: Appendix A does not implement the stated model as written. Cloud 0 is defined as the smallest cloud, cloud 2 as the largest. V_a is \"the volume containing particles from only the largest cloud,\" yet Eq. (A4) places L_{e,cloud0} there; the cloud-2 contribution is missing from Eq. (A5); and in Eq. (A6) cloud 0 appears with an xi0^3 weight even though cloud 0 lies entirely inside V_c. The volume factors in Eqs. (A1)-(A3) also disagree with the shell factors in (A4)-(A6), and for the best-fit xi values they do not sum to unity. As printed, the equations credit the largest cloud's late, optically thin 2 GHz emission to the smallest cloud. If the code implements the printed equations, the central causal attribution does not follow from the model as written; if it implements something else, the equations still must be corrected. No code is shipped, so a reader cannot check.\n\n2. Real but partly acknowledged: the cloud complex is a fitting device. The Eq. (10) luminosity profiles are fit to the same radio light curves they are then used to explain. With no goodness-of-fit statistic, no model comparison against ref. 67, and a by-eye steady-state baseline, the abstract's \"causal relation\" overstates what the evidence supports.\n\n3. Minor: roughly 30 parameters versus ~10 for a single-zone model; the posteriors are local-likelihood contours, not Bayesian credible intervals.\n\nWho this is for: anyone working on blazar multi-messenger timing. It deserves a serious referee — the right outcome is conditional acceptance with a request to fix the appendix, release the code, and soften the causal language.","headline":"A serious multi-zone fit to TXS 0506+056's delayed radio flare, but the printed Appendix A equations do not implement the three-cloud picture, so the central causal claim needs fixing before the result can be trusted as stated.","tokens_in":21622,"tokens_out":9171,"would_cite":true,"duration_ms":83082,"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":"The delayed radio flare of TXS 0506+056 is produced when the jet launched in the 2017 gamma-ray flare reaches three cold-electron clouds at 2-7 pc and accelerates fresh electrons, causally linking the high-energy and radio emission.","keywords":["blazars","TXS 0506+056","delayed radio flares","neutrino multi-messenger astrophysics","jet-cloud interaction","synchrotron self-absorption","time-domain radiative modeling","Fermi-LAT"],"falsifier":"A direct VLBI measurement resolving the location of the radio flare region during 2019-2021: if the flaring radio knot is co-located with the inner jet at sub-parsec or milliparsec distances, or if its apparent motion implies a Lorentz factor or viewing angle substantially different from $\\Gamma = 4.54$ and $\\theta = 5^\\circ$, the 2-7 pc cloud configuration would be ruled out. A simpler test: high-cadence monitoring at 2 GHz, since the model predicts a single late peak about 3.4 years after the gamma-ray peak followed by a simultaneous decline across all bands, so observing a second peak or a frequency-independent rise would contradict the three-cloud geometry.","tokens_in":20323,"feed_emoji":"📡","tokens_out":8147,"duration_ms":68091,"temperature":0.7,"pith_summary":"The paper argues that the gamma-ray flare of TXS 0506+056 in 2017 and its delayed, frequency-dependent radio flare are two stages of one physical chain: plasma launched from the inner jet propagates outward, and when the leading blob reaches the parsec-scale jet it encounters three cold-electron clouds, accelerating fresh electrons whose synchrotron emission produces the later radio flare. The authors build a time-dependent leptonic jet model, couple it to the AM3 particle-interaction code, and fit Fermi-LAT and RATAN-600 light curves. The best fit places the clouds between 2 and 7 pc from the black hole and reproduces the frequency-dependent timing of the radio flare, including the roughly 3.4-year delay of the 2 GHz peak. What matters is that the model establishes a mechanism connecting high-energy and radio emission across spatial scales that the standard single-zone framework cannot describe.","feed_headline":"Model links blazar neutrino flare to radio peak three years later","feed_subtitle":"Jet blobs hitting electron clouds at 2-7 pc reproduce the delayed multi-frequency radio flare of TXS 0506+056.","key_machinery":"The central object is the leading blob of a sequence of expanding, coaxially propagating jet blobs: a sphere of radius $r'(t) = r'_0 + \\eta c t$ carrying a power-law electron population, launched near the broad-line region at $x_{0,\\mathrm{BH}}$ and evolved with the time-dependent radiative code AM3, which handles synchrotron, inverse Compton, pair production, synchrotron self-absorption, and adiabatic cooling. The load-bearing new element is the cloud complex: three stationary cold-electron clouds parameterized by asymmetric Gaussian density profiles with transverse radii $\\xi_0 r$, $\\xi_1 r$, and $r$ (cloud 2 spanning the full jet cross-section). As the leading blob sweeps through the overlapping volumes $V_a$, $V_b$, and $V_c$, it injects fresh non-thermal electrons according to Eq. (10), and the different optical depths of these volumes produce the frequency-dependent radio light curve. The machinery couples the inner-jet ejection profile of Eq. (2) to the cloud-injection profile, with 19 MCMC-searched parameters, to fit the joint gamma-ray and radio light curves.","core_discovery":"On the paper's own terms, the central result is a causal model rather than a statistical correlation: the perturbation that produced the 2017 gamma-ray flare in the inner jet of TXS 0506+056 continued to propagate outward as a chain of expanding blobs, and its leading front, upon reaching a complex of three stationary electron clouds at 2-7 pc, accelerated fresh electrons to a non-thermal spectrum. The superposition of three emitting regions with different transverse sizes, and therefore different self-absorption depths, is what generates the observed frequency-dependent radio light curve: a compact cloud produces the early rise at frequencies at and above 5 GHz while being self-absorbed at 2 GHz, and a larger cloud spanning the full jet cross-section produces the late 2 GHz peak. The model therefore claims that the delayed radio flare is the parsec-scale signature of the same energetic event that produced the high-energy flare, establishing a causal relation between the two.","pith_inferences":["Because the cloud configuration is inferred from a local MCMC search rather than a global posterior, a reader should treat the absolute 2-7 pc distances as provisional and the multi-zone structure with distinct self-absorption depths as the more robust result.","A testable extension: applying the same modeling pipeline to PKS 1424+240, the other neutrino blazar with a delayed radio flare mentioned in the paper, would either confirm that jet-cloud interactions generically produce such delays or reveal that the mechanism is source-specific.","If future VLBI proper-motion measurements track the leading blob as it crosses the cloud complex, the apparent speed should stay constant at the model's $\\Gamma\\beta c$; detecting acceleration or deceleration would directly constrain the mass loading of the clouds, a quantity the present model assumes rather than computes.","Adding proton acceleration to the same jet-cloud setup would allow a self-consistent estimate of the neutrino flux from the 2-7 pc region during the 2017 flare, since the electron clouds that supply the radio-emitting electrons could also provide target matter for hadronic interactions."],"forward_implications":["If the model is right, the 2017 gamma-ray flare and the delayed radio flare of TXS 0506+056 are causally connected stages of one propagation event, not independent activity.","The radio flare's frequency-dependent peak delays, with 2 GHz peaking about 3.4 years after the gamma-ray peak while higher frequencies show earlier double sub-peaks, are a direct consequence of the three-cloud geometry with different self-absorption depths.","The parsec-scale environment of neutrino blazar candidates can be constrained by fitting time-domain multi-frequency radio data, yielding cloud distances of 2-7 pc and individual cloud sizes as fractions of the jet cross-section.","The single-zone framework is insufficient for the radio band; any complete model of blazar multi-messenger emission must include spatially extended emission zones and multiple particle acceleration sites.","The same spatially resolved modeling approach can be applied to other neutrino-candidate blazars with delayed radio flares; the paper notes PKS 1424+240 as a candidate but limits its fit to TXS 0506+056."],"supporting_citations":[{"why":"Supplies the RATAN-600 multi-frequency radio light curves that the model is fitted to and that exhibit the frequency-dependent delayed flare.","marker":"[38]"},{"why":"Provides the multi-wavelength campaign that established the 2017 gamma-ray flare and the IceCube neutrino association, setting the temporal baseline.","marker":"[3]"},{"why":"Identified the delayed radio flare of TXS 0506+056, the phenomenon the model is built to explain.","marker":"[21]"},{"why":"Provides VLBI measurements placing the radio emission in a core of projected size several parsecs, motivating and constraining the 2-7 pc emission region and viewing angle.","marker":"[24]"},{"why":"Provides the previous extended-jet model whose steady-state spectrum is adopted as the baseline and which serves as the comparison for why continuous particle acceleration underproduces low-frequency radio emission.","marker":"[44]"},{"why":"Supplies the BLR and torus radii and photon-field prescriptions used to model the inner-jet environment where the gamma-ray flare originates.","marker":"[50]"},{"why":"Describes the AM3 time-dependent radiative code used to evolve the electron and photon populations in each jet blob.","marker":"[52]"},{"why":"Reports the VLBI evidence for a slow sheath and ultra-fast spine that conflicts with the model's one-dimensional coaxial jet geometry, the key tension the authors acknowledge.","marker":"[66]"},{"why":"Presents a competing model of the delayed radio flare based on jet deceleration; the frequency-dependent light-curve structure is the discriminator between the two explanations.","marker":"[67]"}],"fun_headline_variants":["Jet blobs explain delayed radio flare in neutrino blazar","Neutrino blazar's radio echo traced to parsec-scale electron clouds","Causal link found: blazar flare to delayed radio peak","Expanding blobs hit electron clouds, reproduce blazar's radio delay","Parsec-scale jetscope: delayed radio flare from blob-cloud collisions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the radio emission is generated at 2-7 pc from the black hole by the jet running into stationary cold-electron clouds, rather than at the much smaller distances suggested by the spine-sheath VLBI structure.","fun_headline_variants_meta":{"raw":{"variants":["Jet blobs explain delayed radio flare in neutrino blazar","Neutrino blazar's radio echo traced to parsec-scale electron clouds","Causal link found: blazar flare to delayed radio peak","Expanding blobs hit electron clouds, reproduce blazar's radio delay","Parsec-scale jetscope: delayed radio flare from blob-cloud collisions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000676,"raw_usage":{"total_tokens":3111,"prompt_tokens":1017,"completion_tokens":2094,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":633,"completion_tokens_details":{"reasoning_tokens":2003}},"tokens_in":633,"tokens_out":2094,"duration_ms":14528,"temperature":1.0,"reasoning_tokens":2003,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:22:23.170857+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct VLBI measurement resolving the location of the radio flare region during 2019-2021: if the flaring radio knot is co-located with the inner jet at sub-parsec or milliparsec distances, or if its apparent motion implies a Lorentz factor or viewing angle substantially different from $\\Gamma = 4.54$ and $\\theta = 5^\\circ$, the 2-7 pc cloud configuration would be ruled out. A simpler test: high-cadence monitoring at 2 GHz, since the model predicts a single late peak about 3.4 years after the gamma-ray peak followed by a simultaneous decline across all bands, so observing a second peak or a frequency-independent rise would contradict the three-cloud geometry.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identified the delayed radio flare of TXS 0506+056, the phenomenon the model is built to explain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides VLBI measurements placing the radio emission in a core of projected size several parsecs, motivating and constraining the 2-7 pc emission region and viewing angle."},{"cited_title":"Bhjet: a public multi-zone, steady state jet + thermal corona spectral model","cited_arxiv_id":"2108.12011","evidence_quote":"Provides the previous extended-jet model whose steady-state spectrum is adopted as the baseline and which serves as the comparison for why continuous particle acceleration underproduces low-frequency radio emission."},{"cited_title":"Padovani, F","cited_arxiv_id":null,"evidence_quote":"Supplies the BLR and torus radii and photon-field prescriptions used to model the inner-jet environment where the gamma-ray flare originates."},{"cited_title":"Ghisellini, F","cited_arxiv_id":null,"evidence_quote":"Describes the AM3 time-dependent radiative code used to evolve the electron and photon populations in each jet blob."}],"review_version":1}