{"id":"04ba3dbb-530f-47fb-be9d-1698af9737f9","arxiv_id":"2411.17389","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Simulations show that decreasing the magnetic pitch suppresses the Rayleigh-Taylor instability in over-pressured jets and instead excites a current-driven kink instability, offering a qualitative explanation for BL Lac quasi-periodic oscillations.","lead":"This paper runs 3D computer simulations of magnetized relativistic jets and finds that a Rayleigh-Taylor instability grows where the jet meets the surrounding gas. It also argues that a helical kink instability can develop after the jet's recollimation shocks, which may explain periodic brightness changes seen in the blazar BL Lac.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The kink identification in §3.6 rests on morphology and advection speed, not on mode analysis or a growth rate, and §3.4/Fig. 9 contain an internal reversal of the pitch ordering; the headline pitch-to-kink transition therefore needs a quantitative mode test before the BL Lac QPO link is accepted.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing point: the helical structure is labelled a CD kink without a quantitative mode analysis. My review confirms that this is the step on which the abstract's central claim and the BL Lac QPO interpretation depend. The additional internal inconsistency in §3.4 (pitch ordering of MHD2/MHD3 reversed relative to the figure captions and Table 2) is a presentation/correctness issue that needs fixing, but it can be resolved in the paper's favour by taking the captions and Table 2 as authoritative. I therefore do not move the verdict: the paper remains CONDITIONAL on a mode-based demonstration of the kink and on correcting the pitch labels. I disagree with no part of the reader's assessment, and I do not see a basis for rejecting the paper's simulations or for accepting the kink claim as currently supported. A Fourier-mode re-analysis of the stored simulation outputs is a well-posed, single check that would settle the question one way or the other.","tokens_in":13787,"tokens_out":9160,"duration_ms":86652,"concrete_test":"Perform an azimuthal-mode analysis on the MHD3 run (P=0.16): on cylindrical shells 1Rj≤R≤2Rj and z=30–50Rj, decompose the transverse velocity and magnetic-field perturbations into exp(i m phi) harmonics and fit the m=1 amplitude in time over t≈350–400. If m=1 grows exponentially with a rate consistent with the CDI dispersion relation for the local helical field and has a phase speed near 0.8c, the kink identification is supported. If the m=1 amplitude does not grow, or if the twisted structure is equally well described by an m=0 pinch or by nonlinear RTI harmonics, the claimed pitch-to-kink transition and the BL Lac QPO link are not established. Separately, re-plot Fig. 9 with line colors keyed to the P values in Table 2 to determine which of the text and captions gives the correct pitch ordering.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that decreasing magnetic pitch suppresses RTI and instead excites the CD kink instability. The kink is identified in §3.6 from the helically twisted density structure in Fig. 11 and from a propagation speed of about 0.8c inferred from four azimuthally averaged density profiles (Fig. 12), with qualitative comparison to Mizuno et al. (2014) and Singh et al. (2016). No azimuthal Fourier decomposition, no m=1 growth-rate fit, and no comparison with the linear CDI dispersion relation is given. A nonlinearly evolved RTI finger, a recollimation-shock remnant, or a numerical artifact could produce a morphologically similar twisted boundary; the smooth-boundary test in Appendix B checks only the initial sharpness, not the mode content. Independent of this, the manuscript's own labelling is internally inconsistent: §3.4 calls MHD2 (P=0.50 in Table 2) the lower-pitch case and MHD3 (P=0.16) the higher-pitch case, while the Fig. 8/9 captions and the abstract use the opposite ordering. If the captions are authoritative, the pitch dependence survives, but the kink identification is still not established. Since the QPO interpretation depends on the kink being real, this is the load-bearing weak point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents three-dimensional relativistic magnetohydrodynamic (RMHD) simulations of over-pressured relativistic jets with helical magnetic fields, extending earlier two-dimensional work. The authors study how the magnetic pitch (ratio of poloidal to toroidal field) and field strength affect the growth of instabilities at the jet-ambient interface. In 2D axisymmetric simulations the recollimation shock structure is stable; in 3D, the pure hydrodynamic case develops Rayleigh-Taylor instability (RTI), while magnetized cases show that increasing toroidal field (lower pitch) suppresses RTI. The paper claims that in the lowest-pitch cases the jet instead develops a current-driven (CD) kink instability, and it connects this qualitatively to quasi-periodic oscillations observed in BL Lac.","tokens_in":14119,"tokens_out":6923,"duration_ms":60267,"significance":"If confirmed, the claimed pitch-controlled transition from RTI to CD kink would be a valuable result for interpreting jet stability and variability in AGN jets. The paper includes a useful smooth-boundary test (Appendix B), a higher-Lorentz-factor test (Appendix C), and an analytic radial-motion equation (Appendix A) that is checked against the simulations. However, the central kink identification is based on morphology and advection speed rather than a quantitative mode analysis, and the reported pitch ordering is internally inconsistent. The result is therefore not yet established, but the underlying numerical campaign is appropriate and the question is timely.","major_comments":[{"comment":"The pitch assignments are contradictory. In the first paragraph of §3.4 the text refers to \"lower magnetic pitch (MHD2) with P = 0.16 and higher magnetic pitch (MHD3) with P = 0.5\", but Table 2 lists MHD2 with P = 0.50 and MHD3 with P = 0.16. The Fig. 8 captions are consistent with Table 2, but the Fig. 9 caption says the green line is MHD2 (P = 0.5) and the purple line is MHD3 (P = 0.16), while the text immediately below states that the \"lower pitch case (green line)\" has larger radial velocity and the \"higher pitch\" (purple line) has reduced velocity. Depending on which statement is taken as authoritative, the central trend either supports or contradicts the abstract's claim that lower pitch suppresses RTI. The paper must be corrected so that the text, tables, and figure captions all use the same pitch ordering.","section":"§3.4 and Fig. 9"},{"comment":"The identification of the CD kink instability is not quantitative. The claim rests on the helically twisted density structure in Fig. 11 and a propagation speed of about 0.8c inferred from the azimuthally averaged density profiles in Fig. 12, compared morphologically with earlier simulations by Mizuno et al. (2014) and Singh et al. (2016). Because the initial jet already contains a helical magnetic field, a passively advected helical pattern can produce the same appearance without any instability growing. The paper does not provide an azimuthal Fourier decomposition, a time trace of the m = 1 mode amplitude, a growth rate, or a comparison with the linear CDI dispersion relation; the smooth-boundary test in Appendix B checks only the initial sharpness, not the mode content. Please add a quantitative mode analysis (e.g., evolution of the m = 1 Fourier amplitude at fixed radius and axial position) and a growth-rate estimate, or rephrase the claim as a morphological similarity rather than an established instability.","section":"§3.6"},{"comment":"The RTI strength is characterized by the azimuthally averaged radial velocity at a single axial position, z = 60 R_j, with one realization per setup and no error bars. The text in §3.3 notes that \"axial variation by the development of instability depends a little on the spatial resolutions\" but no resolution test is shown. The reader therefore cannot assess whether the pitch-dependent differences in the radial velocity traces are robust or within numerical noise. Please provide a convergence study for at least one magnetized case and state the time-averaging or error-bar methodology used for the diagnostics in Figs. 7 and 9.","section":"Figs. 7 and 9"},{"comment":"The radial-motion equation (A.7) is shown to reproduce the recollimation kinematics in Fig. A.1, but the paper does not connect this equation to the observed pitch dependence of the RTI. If the magnetic tension term b_φ^2/R or the centrifugal term γ^2 ω_t v_φ^2/R in Eq. (A.7) is responsible for the suppression at lower pitch, this should be stated and tested; otherwise the claim that magnetic pitch controls RTI remains purely empirical. Please clarify whether Eq. (A.7) can be used to interpret the pitch trend, or explicitly defer that connection to future work.","section":"Appendix A and §4"}],"minor_comments":[{"comment":"The word \"excitement\" should be \"excitation\" throughout the manuscript (e.g., abstract, §3.6).","section":"Abstract and §3.6"},{"comment":"The caption contains a typo: \"higher magnetic pith case\" should be \"higher magnetic pitch case\".","section":"Fig. 9 caption"},{"comment":"The labels \"MHD1-2D\" and \"MHD-2D\" are used interchangeably in the text; please unify the notation.","section":"§3.1 and Table 1"},{"comment":"The paper structure described in the Introduction says \"We analyze ... in Section 4 and summarize and discuss our findings in Section 4\", but the results for the kink instability are in §3.6 and the summary is in §4; please correct the section references.","section":"Introduction and §4"},{"comment":"The MHD1 case is shown at ts = 365 while the other cases are at ts = 400; please state explicitly that the comparison is made at slightly different times and justify why this does not affect the conclusions.","section":"§3.3"},{"comment":"The maximum magnetization values are quoted without specifying the time or radius at which they are evaluated; please provide these details.","section":"§3.5"},{"comment":"Quantitative comparison with the growth rates of Mizuno et al. (2014) and Singh et al. (2016) would strengthen the kink discussion; currently the comparison is only qualitative.","section":"§3.6 and references"},{"comment":"The scaling argument l ∝ 1/γ is correct, but the numerical comparison \"49/14 ≃ 10/3\" should be given with the actual values and a brief explanation of how the recollimation length was measured.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of A&A and addresses a topic of interest to the jet physics community. The two issues that require major revision are the inconsistent pitch labeling (which must be fixed before the paper can be evaluated) and the lack of a quantitative kink-mode analysis. The kink identification is the main technical risk: a skeptical reader could argue that the helically twisted structure is simply the advected initial field, especially since the propagation speed is close to the jet speed. A Fourier mode analysis and growth-rate measurement would settle this. The RTI suppression by lower pitch appears to be a genuine trend, but the reversed statements in §3.4 and Fig. 9 need to be resolved. The paper's observational connection to BL Lac is qualitative and appropriately framed as such."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the systematic 3D RMHD scan of magnetic pitch in over-pressured jets: higher pitch (more poloidal field) promotes RTI, lower pitch suppresses it and leaves a helically twisted structure that the paper interprets as the CD kink. The trend itself is credible—the density slices and radial velocity traces point the same way, and the comparison runs (2D, HD, weak-field) give a useful control set. The paper also does several things right: it rules out KHI, RMI, and CFI using published criteria, tests sensitivity to a smooth boundary, checks a higher Lorentz factor, and derives a radial acceleration formula that matches the recollimation motion in both 2D and 3D. That is real work, honestly presented.\n\nThe soft spot is exactly where the stress-test puts it. The kink identification in §3.6 is morphological: a helically twisted density structure plus a propagation speed of about 0.8c inferred from four density peaks. There is no azimuthal Fourier decomposition, no m=1 growth-rate fit, no comparison to the linear CDI dispersion relation. A nonlinearly evolved RTI finger or a shock-driven deformation could look similar, so the headline pitch-to-kink transition is not yet demonstrated. The authors are careful to call the BL Lac QPO link \"qualitative,\" which is fair, but they still lean on the kink interpretation as if it were established. Also, §3.4 has an internal labeling reversal: the text calls MHD2 the lower-pitch case and MHD3 the higher-pitch case, while Table 2 and the Figure 8/9 captions say the opposite. That is likely a typo, but it makes the already subtle pitch argument harder to follow and must be fixed. Minor: the axial resolution is modest (10 cells per jet radius), and the fiducial run stops at t=365; the paper justifies this, but it does limit confidence in the kink morphology.\n\nWho is this for? People working on jet stability, recollimation shocks, and blazar variability interpretation. It is a useful parameter study, not a paradigm shift. I would cite it for the pitch dependence of RTI in 3D RMHD, and I would send it to review—the core question is legitimate and the simulations are reproducible within the PLUTO framework. But I would press the authors to either add a mode analysis or soften the kink claim, and to clean up the labeling.\n\nRecommendation: send to peer review, with a request for quantitative kink diagnostics or a more cautious interpretation.","headline":"Solid 3D RMHD parameter study of magnetic pitch in over-pressured jets; the pitch trend for RTI is credible, but the claimed kink transition needs a quantitative mode analysis before it can carry the BL Lac QPO interpretation.","tokens_in":14578,"tokens_out":2851,"would_cite":true,"duration_ms":27662,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76W05","76E17","85A30","85-08"],"pacs":[],"model":"deepseek-v4-flash","headline":"Magnetic pitch decides which instability shreds over-pressured relativistic jets.","keywords":["galaxies: jets","magnetohydrodynamics (MHD)","instabilities","relativistic jets","Rayleigh-Taylor instability","kink instability","BL Lacertae objects","numerical simulations"],"falsifier":"A Fourier decomposition of the jet density or magnetic field in the twisted region showing a dominant azimuthal mode other than m=1, or a measured growth rate of the helical perturbation inconsistent with kink instability, would show that the structure is not the CD kink.","tokens_in":74,"feed_emoji":"🔭","tokens_out":4274,"duration_ms":84832,"temperature":0.7,"pith_summary":"This paper uses three-dimensional relativistic magnetohydrodynamic simulations to study how the magnetic field structure of an over-pressured jet controls which instability destroys its collimation. The authors find that when the magnetic pitch is high (poloidal field dominant), Rayleigh-Taylor fingers grow at the jet–ambient interface inside the recollimation shock structure. When the pitch is low (toroidal field dominant), those fingers are suppressed, but a current-driven kink instability twists the jet into a helical shape; the twist propagates at about 0.8c, close to the jet flow speed. The paper argues this kink-excitation route matches the quasi-periodic oscillations observed in the blazar BL Lac after the passage of recollimation shocks. The result matters because it gives a single physical parameter—magnetic pitch—that selects between two very different disruption mechanisms in relativistic jets.","feed_headline":"Magnetic pitch sets which instability disrupts relativistic jets","feed_subtitle":"Simulations show lower pitch suppresses Rayleigh-Taylor fingers and excites a kink linked to BL Lac oscillations.","key_machinery":"The controlling object is the magnetic pitch P = a/(k Rj), the ratio of poloidal to toroidal field components on the jet axis, with lower P meaning a more toroidally dominated field. The simulations also use the effective inertia I = γ²ρh + Bz² + Bφ², whose radial variation traces the recollimation shock structure and locates the low-inertia interface where RTI develops. The authors derive a radial momentum equation for the jet boundary and use it to compute the effective gravity that drives RTI and to rule out centrifugal instability via the magnetization criterion σ/(1+σ) > (θ0γ)²/16. The kink identification relies on comparing the helically twisted morphology and propagation speed with those in earlier CD-kink simulations of relativistic jets.","core_discovery":"The paper's central claim is that the magnetic pitch of a helical field, not the magnetization alone, sets the dominant instability in over-pressured relativistic jets. Comparing 3D RMHD runs with pitch values P = 0.50, 0.22, and 0.16, they find that decreasing the pitch weakens the Rayleigh-Taylor instability at the jet–external-medium interface, and instead excites the CD kink instability, producing a helically twisted density structure downstream of the recollimation shock. They identify the kink by the helical morphology and a measured propagation speed of about 0.8c, consistent with earlier kink-jet simulations, and match the interpretation of BL Lac's quasi-periodic oscillations. The paper is an extension of the earlier 2D axisymmetric simulations of over-pressured jets, which showed no instability; the 3D runs reveal that non-axisymmetric modes are essential for RTI and kink growth.","pith_inferences":["A testable next step would seed simulations with BL Lac parameters and synthesize light curves and polarization maps to compare directly against the observed 2020 outburst quasi-periodic oscillations.","If the kink identification were confirmed by a Fourier mode analysis measuring the growth of the azimuthal m=1 amplitude, the pitch-to-kink transition could be turned into a quantitative stability diagram in jet parameter space.","The stabilizing effect of lower pitch suggests that the presence of Rayleigh-Taylor fingers versus helical twisting could serve as an observational diagnostic of the magnetic configuration of relativistic jets in very long baseline interferometry images."],"forward_implications":["If magnetic pitch controls the RTI-to-kink transition, observations of jet morphology alone could constrain the magnetic field structure of AGN jets.","The result provides a physical mechanism linking recollimation shocks to BL Lac quasi-periodic oscillations: kink waves launched after the shock travel along the jet at about 0.8c.","Two-dimensional axisymmetric simulations miss both instabilities; instability studies of jets need full three-dimensional treatment.","The pitch dependence suggests jets with stronger toroidal fields (lower pitch) may remain collimated longer before being twisted, while higher-pitch jets are disrupted by Rayleigh-Taylor fingers closer to the core.","The stability criteria against centrifugal instability used in the paper can be applied to future jet simulations to distinguish instability drivers without expensive mode decomposition."],"supporting_citations":[{"why":"Provides the 2D RMHD over-pressured jet setup and recollimation shock structure that this work extends to 3D.","marker":"Mizuno et al. 2015"},{"why":"Supplies the effective inertia criterion for Rayleigh-Taylor instability in the interface between jet and ambient medium.","marker":"Matsumoto & Masada 2013"},{"why":"Gives the updated analytical stability criterion for RTI used to interpret the hydrodynamic comparison run.","marker":"Matsumoto et al. 2017"},{"why":"Establishes that magnetization stabilizes RTI, the trend the paper confirms and then refines in terms of pitch.","marker":"Millas et al. 2017"},{"why":"Provides the reference CD kink instability simulations whose helical morphology is compared to identify the kink in the present runs.","marker":"Mizuno et al. 2014"},{"why":"Offers the earlier kink-jet simulation with similar propagation speed, used to support the 0.8c measurement.","marker":"Singh et al. 2016"},{"why":"Supplies the centrifugal instability stability criterion used to rule out CFI in the magnetized cases.","marker":"Matsumoto et al. 2021"},{"why":"Reports the BL Lac quasi-periodic oscillations that the simulated post-recollimation kink is invoked to explain.","marker":"Jorstad et al. 2022"}],"fun_headline_variants":["Magnetic pitch decides which instability disrupts jets","Lower magnetic pitch suppresses RTI and excites kink in jets","Pitch controls jet instability: low pitch swaps RTI for kink","Magnetic pitch selects jet instability mode between RTI and kink","Jet instability type flips with magnetic pitch in simulations"],"cache_read_input_tokens":16768,"weakest_assumption_plain":"The helically twisted structure observed in the simulations is assumed to be the CD kink instability, inferred from morphology and propagation speed rather than from a mode decomposition or a measured growth rate.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic pitch decides which instability disrupts jets","Lower magnetic pitch suppresses RTI and excites kink in jets","Pitch controls jet instability: low pitch swaps RTI for kink","Magnetic pitch selects jet instability mode between RTI and kink","Jet instability type flips with magnetic pitch in simulations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000406,"raw_usage":{"total_tokens":2113,"prompt_tokens":950,"completion_tokens":1163,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":1078}},"tokens_in":566,"tokens_out":1163,"duration_ms":9300,"temperature":1.0,"reasoning_tokens":1078,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:09:10.216938+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A Fourier decomposition of the jet density or magnetic field in the twisted region showing a dominant azimuthal mode other than m=1, or a measured growth rate of the helical perturbation inconsistent with kink instability, would show that the structure is not the CD kink.","supporting_citations":[],"review_version":1}