{"id":"b446889b-cd9d-4c16-a86c-8f1d17560fdc","arxiv_id":"2504.20962","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":15,"one_line_summary":"A new analytical fit of jet emission profiles reveals that the jets of Swift J1727.8-1613 were intrinsically asymmetric, with the approaching jet stronger and the magnetic field increasing with time.","lead":"Astronomers introduce a new method for fitting the brightness profiles of black hole jets and apply it to the X-ray binary Swift J1727.8-1613. They find the approaching jet is intrinsically brighter than the receding one, which suggests the jet's magnetic field was strengthening with time.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The beam convolution in Eq. (8) treats each jet as an even two-sided source, injecting a mirror-image term into the fitted profiles; this can bias the jet/counterjet asymmetry on which the central claim rests.","rationale":"The reader's stated weakest assumption, the steady-state plus time-delay interpretation, is explicitly acknowledged by the authors in Section 6 and is a legitimate concern about interpretation. My stress-test found a more concrete, internal issue that precedes that interpretation: the asymmetric fits use Eq. (8) to convolve a one-sided jet with the beam as though the jet were even on both sides of the centroid. The physical image of a jet/counterjet system is not even unless the two sides have identical intrinsic parameters, which is exactly what the paper argues against. The mirror term G(r+x) f_app(x) that Eq. (8) adds to the approaching-side profile is not a contribution of the approaching jet to its own side; it is the approaching jet's leakage onto the counterjet side, which should instead be represented by the counterjet profile f_rec(x) when modeling the total image. Because the paper's central claim (intrinsically asymmetric jets, slow velocity, magnetic field increasing with time) is inferred from p_app, p_rec, A1, and A2, a reanalysis with the correct convolution is required before those source-specific claims can be accepted. I therefore keep the conditional verdict but for a different and sharper reason: the current numerical implementation of the convolution may be biasing the exact asymmetry it is used to measure. If the recomputation changes the best-fit parameters, the source-specific conclusions should be rejected or substantially revised; if it leaves them unchanged, the reader's original caveat about time dependence remains the main residual limitation.","tokens_in":20060,"tokens_out":15722,"duration_ms":182336,"concrete_test":"Re-run the MCMC fits of Section 5 with a signed, physically correct convolution: for the approaching-side data at r>0, use C(r) = integral_0^inf [G(r-x) f_app(x) + G(r+x) f_rec(x)] dx, and for the receding side use the mirrored expression, with the same priors, data cuts, and Delta_xi values (1.5 and 2.0 mas). Then compare the best-fit p_app, p_rec, A1, A2, beta, and the log-likelihood against the published Table 1, and also test whether an intrinsically symmetric jet model (same p, A1, A2 for both sides, with Doppler factors from Eq. A11) is rejected at the same significance. If the asymmetry and beta ~ 0.3-0.4c survive, the central claim holds; if not, the mirror term in Eq. (8) is the cause.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The load-bearing step is the convolution of the intrinsic profile with the 2.5 mas restoring beam, Eq. (8). The model profile dF_nu/d|xi| is defined as a function of absolute offset, and Eq. (8) integrates it over x from -infinity to +infinity. For the approaching jet this is the convolution of an even function, i.e., it is equivalent to assuming the approaching jet also emits on the counterjet side. The observed approaching-side profile should be the one-sided convolution I_app(r) = integral_0^inf G(r-x) f_app(x) dx plus the counterjet leakage integral_0^inf G(r+x) f_rec(x) dx, not the self-mirror integral_0^inf G(r+x) f_app(x) dx. The authors state explicitly that they 'calculate dF_nu,c/d|xi| for the jet and counterjet using their respective parameters for both positive and negative xi', so the same mirror contamination is applied on both sides. With sigma_t = 1.06 mas and data fitted from |xi| ~ 0.5-1 mas, the spurious G(r+x) term is comparable to the true G(r-x) term near the core, so the fitted p, A1, and A2 for the two jets can shift materially. The large-scale flatter approaching jet may survive, but the quantitative asymmetry and the derived time evolution of the magnetic field are not secure until the convolution is done with the correct one-sided support.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an analytical model, based on Blandford & Königl (1979), for the spatial profile dFν/d|ξ| of a compact jet (Eq. 6, Appendix A), convolves it with the VLBI restoring beam (Eq. 8), and fits the result to the 8.37 GHz jet and counterjet profiles of Swift J1727.8–1613 from Wood et al. (2024) using MCMC. Symmetric-jet fits fail clearly (Figs. 3–4), so the authors fit the two jets with independent parameters, finding the approaching jet intrinsically stronger and flatter. They attribute the asymmetry to the observed rise in radio flux combined with a light-travel time lag, derive β≈0.35 at i=45°, and infer B1, ṃeff, jet powers, and magnetic flux, concluding that the magnetic field increased with time and that the magnetic flux is well below the MAD limit.","tokens_in":20554,"tokens_out":12834,"duration_ms":127282,"significance":"If the quantitative results survive scrutiny, the method is a valuable complement to core-shift measurements: it uses the full one-dimensional brightness distribution, gives closed-form expressions for physical quantities (Appendix A), and is generally applicable to XRB and AGN jets. The derivation in Appendix A is complete, the forward-modeling logic is transparent, and the failure of symmetric fits is a robust qualitative result. However, the convolution step in Eq. (8) mis-treats the support of the jet emissivity for asymmetric jets, which can bias exactly the fitted parameter differences on which the central quantitative claim rests; in addition, the time-evolution interpretation is an assumption rather than a fitted result. These issues make the quantitative conclusions tentative pending revision.","major_comments":[{"comment":"The convolution is performed on dFν/d|x|, i.e., an even function of x, so the model adds to each side a mirror image of that jet's own emission. For a one-sided jet, the observed approaching-side profile should be ∫_0^∞ G(ξ−x) f_app(x) dx + ∫_0^∞ G(ξ+x) f_rec(x) dx, i.e., direct approaching emission plus physical beam leakage from the counterjet, with the analogous expression for the counterjet side. The current implementation instead adds ∫ G(ξ+x) f_app(x) dx to the approaching side and ∫ G(ξ+x) f_rec(x) dx to the receding side, replacing true cross-contamination with self-mirror terms. With σ_t = 1.06 mas and the fit starting at |ξ′|≈0.5–1 mas, the spurious term is a large fraction of the direct term, so the fitted p, A1, and A2 for the two jets can shift substantially. Since the intrinsic-asymmetry claim and the derived time evolution of B rest on the difference between the two fitted parameter sets, the fits must be redone with the correct one-sided convolution and cross-contamination before the quantitative conclusions can be accepted.","section":"Section 4, Eq. (8)"},{"comment":"The inference that the magnetic field increased with time is not directly fitted. The model treats the jet and counterjet as two independent steady-state solutions with different A1, A2, and p, and the time evolution is imposed post-hoc through the light-travel-time argument. As the authors themselves note in Section 6, a more realistic model would need to couple the time dependence and the time lags; an intrinsically asymmetric steady jet would produce the same type of fit. The abstract and conclusions should therefore state the time-evolution result as an interpretation with substantial systematic uncertainty, not as a direct finding.","section":"Sections 5–7"},{"comment":"The systematic uncertainty in the assumed core offset Δξ is not propagated into the derived quantities. The two adopted values, Δξ = 1.5 and 2.0 mas, yield B1_app = 3.3 G versus 9.0 G and ϕ_BH = 0.044 versus 0.12, among other changes. Quoting single numbers in Table 2 understates the model dependence; please present the derived physical quantities with the Δξ systematic included, or state explicitly that the values are conditional on Δξ.","section":"Section 5, Table 2"}],"minor_comments":[{"comment":"The caption's second sentence says 'Figure 7(a) compares the emitted spectra (before smoothing) with the data'; this should refer to Figure 7(b).","section":"Figure 7 caption"},{"comment":"Convergence is judged by manual inspection that the walkers are 'no longer significantly evolving'; a quantitative convergence criterion (e.g., the Gelman–Rubin statistic) would be more reproducible.","section":"Section 5, MCMC"},{"comment":"In Eq. (6) the optical-depth exponent is written −bp/2−b−1, which corresponds to a=2 in Eq. (A12); please state a=2 explicitly in the main-text equation to avoid confusion.","section":"Equations (6) and (A12)"},{"comment":"The spectral index α = 0.19±0.07 used to fix b = 1.17 was measured four days before the VLBA observation; the paper should note the possible effect of spectral variability on the fixed value of b.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The convolution problem in Eq. (8) is the main technical blocker. If the authors re-fit the data with the correct one-sided convolution and cross-contamination, and the qualitative asymmetry persists, the paper is likely publishable. Please also ask them to temper the time-evolution language in the abstract and conclusions, and to propagate the Δξ systematic into Table 2. No concerns about novelty or scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper to know about: Zdziarski et al. propose a new way to extract jet physics from resolved radio images of compact jets, fitting the spatial profile with an analytically normalized version of the Blandford-Konigl model. Applied to Swift J1727.8-1613, they find the approaching jet is intrinsically flatter and stronger than the receding one, attribute the asymmetry to light-travel delay during a radio flare, and derive beta~0.35c, a rising magnetic field, and sub-MAD magnetic flux. The analytic formula (Eqs A12-A16) with physical constants is genuinely new; Paragi et al. (2013) had similar functional forms without the physical normalization. The derivation is complete and the symmetric-jet fits clearly fail, so the empirical asymmetry is solid.\n\nThe soft spots are real. The steady-state plus time-delay interpretation is explicitly tentative and rests on the rising radio light curve; the authors say a more realistic time-dependent model is needed. That's fine as a first cut. The more serious issue is the beam convolution. Equation 8 convolves dFnu/d|xi|, which is an even function, over both positive and negative xi. When fitting the approaching jet, they use the approaching jet parameters on both sides, effectively mirroring the approaching emission onto the counterjet side rather than including the actual counterjet leakage. With a 2.5 mas beam and data points starting near the core, this mirror term is not negligible; it can shift p, A1, and A2 enough to change the quantitative asymmetry and the derived magnetic field evolution. The qualitative flatter approaching jet may survive, but the quoted numbers are not secure until the convolution is done with one-sided support and the opposite jet's contamination included. This is fixable, but it needs to be redone.\n\nThere are also smaller concerns: the fits truncate at 4 and 6 mas with post-hoc exclusion of wiggles, which is understandable given ISM interactions but adds systematic uncertainty. The derived magnetic flux depends on the assumed opening angle and other parameters; the authors are upfront about this.\n\nWho is this for? Anyone modeling resolved XRB or AGN jets. The method is a real step forward, but the application as presented needs major revision. I'd send it to a serious referee, and I'd expect them to require a corrected convolution and a sensitivity check. The underlying physics and data are good enough that the paper should not be desk-rejected.","headline":"A genuinely useful new fitting method for resolved compact jets, but the beam convolution in the Swift J1727.8-1613 application is flawed and needs to be redone.","tokens_in":21053,"tokens_out":5181,"would_cite":true,"duration_ms":53036,"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":"This paper shows that the apparent jet/counterjet brightness asymmetry in Swift J1727.8-1613 is a light-travel effect: the receding jet is seen at an earlier, weaker epoch, and the jets move at about 0.3–0.4 times the speed of light.","keywords":["compact jets","X-ray binaries","synchrotron self-absorption","Blandford-Königl model","jet/counterjet asymmetry","core shift","Swift J1727.8-1613","magnetic flux"],"falsifier":"A decisive test would be to observe Swift J1727.8-1613 with VLBA at two epochs separated by about the light-travel lag, roughly one day, while the radio flux is changing: the steady-state-plus-delay model predicts that the counterjet profile at the later epoch should match the approaching-jet profile at the earlier epoch after beam convolution, and that the asymmetry should track the derivative of the radio light curve. A persistent asymmetry in the opposite sense to the flux derivative, or an asymmetry that persists when the flux is flat, would falsify the time-delay interpretation.","tokens_in":19889,"feed_emoji":"📡","tokens_out":6827,"duration_ms":63919,"temperature":0.7,"pith_summary":"This paper develops a way to read a compact jet's physics directly from its resolved radio image: instead of using only the frequency-dependent core position, it fits the full spatial profile of flux density along the jet and counterjet with analytical formulae derived from the standard self-absorbed synchrotron model. Applied to Swift J1727.8-1613, the most resolved continuous jet from a stellar-mass black hole, the method yields jet velocities around 0.3–0.4c and implies that the magnetic field strength grew with time. The key finding is that the approaching jet is intrinsically stronger than the receding one even after Doppler factors are included, which the authors attribute to a rise in the radio flux combined with the light-travel delay that makes the receding jet appear at an earlier epoch. A sympathetic reader would care because the technique uses far more of the available information than core-shift methods and is proposed as a general tool for jets from both stellar-mass and supermassive black holes.","feed_headline":"Swift J1727.8-1613's jet asymmetry is a time-delay effect","feed_subtitle":"Fitting the full radio brightness profile of a black-hole jet reveals slow jets, a brighter near side, and a rising magnetic field.","key_machinery":"The load-bearing object is the analytical brightness profile $$\\frac{dF_\\nu}{d|\\xi|} = A_1 \\$delta^{{1/2}}$ |\\xi|^{1 + b/2} \\left[1 - \\exp\\left(-A_2 \\$delta^{{1+p/2}}$ |\\xi|^{-bp/2 - b - 1}\\right)\\right],$$ which gives the flux density per unit angular separation along the jet in terms of the Doppler factor $\\delta$, the magnetic-field spatial index $b$, the electron power-law index $p$, and two amplitudes $A_1$ and $A_2$ tied to the field strength and mass-flow rate. This profile is convolved with the telescope restoring beam before fitting, and the two jets are allowed to have different parameters so that a time-delay asymmetry can be represented. Fitting $A_1$ and $A_2$ and then inverting the model formulae yields the magnetic field, the effective mass-flow rate, the jet power components, equipartition and magnetization parameters, and the magnetic flux.","core_discovery":"On the paper's own terms, the central claim is that the spatial distribution of synchrotron flux density along a compact jet can be modelled analytically and fitted to resolved images, and that doing so for Swift J1727.8-1613 reveals an intrinsic jet/counterjet asymmetry. The paper shows that no symmetric steady-state model can reproduce both the approaching and receding profiles: the approaching jet is flatter and stronger at large separations. Because the source's radio flux was rising on roughly a day timescale, matching the light-travel lag between the two jets, the asymmetry is interpreted as evolution of the same underlying jet seen at different epochs. The best-fit solutions give a bulk velocity of about 0.35–0.36c at the adopted inclination of 45 degrees, with the counterjet parameters corresponding to an earlier, weaker state. The inferred magnetic field at the jet base increased by a factor of tens between the counterjet and jet epochs, while the conserved magnetic flux threading the black hole stays well below the magnetically arrested disk limit.","pith_inferences":["A natural extension would be to observe Swift J1727.8-1613 again with very long baseline interferometry while the radio flux is changing: the time-delay model predicts that the asymmetry should reverse when the light curve turns from rising to falling.","With a dense series of images through one radio flare, one could attempt tomographic reconstruction of the evolving magnetic-field profile rather than treating each jet side as a separate steady state.","The method's analytical profile assumes self-similar electron reacceleration; the paper itself notes wiggles at large separations that indicate clumpy interstellar-medium interactions, so adding a clump-scattering term seems a testable refinement."],"forward_implications":["For Swift J1727.8-1613, the jets are slow, with bulk velocity roughly $0.3$–$0.4c$, slower than is often assumed for compact X-ray binary jets.","The magnetic field in the jet grew over the roughly one-day light-travel lag between counterjet and jet, while the magnetic flux stayed low; this hard-state jet is not magnetically arrested.","The counterjet profile is effectively a time-delayed image of the approaching jet, so a rising or falling radio light curve should control which side appears intrinsically brighter.","The method affords an independent route to magnetic field and power estimates that uses the whole spatial profile rather than a single core displacement, and it can be applied to extragalactic jets on parsec scales."],"supporting_citations":[{"why":"Establishes the flat-spectrum self-absorbed synchrotron jet model whose spatial flux distribution the new formulae are derived from.","marker":"Blandford & Königl 1979"},{"why":"Generalizes the model to power-law spatial dependencies of electron density and magnetic field, providing the $b$ and $p$ parametrisation used here.","marker":"Königl 1981"},{"why":"Supplies the VLBA observation and the 8.37 GHz jet and counterjet emission profiles that the method is fitted to.","marker":"W24"},{"why":"Provides the core-shift measurement and fast knot velocities used to set the core offset and to constrain distance and inclination.","marker":"Wood et al. 2025"},{"why":"Gives the Appendix-A formulation of the Blandford–Königl model from which the analytical $dF_\\nu/d|\\xi|$ formula is obtained.","marker":"Zdziarski et al. 2019"},{"why":"Supplies the MCMC sampling algorithm used to fit the model parameters and their uncertainties.","marker":"Foreman-Mackey et al. 2013"},{"why":"Provides the radio spectrum measurement with $\\alpha = 0.19 \\pm 0.07$ that fixes the magnetic-field index $b$ in the model.","marker":"Miller-Jones et al. 2023"}],"fun_headline_variants":["Light travel time explains Swift J1727's jet asymmetry","New analytical model reveals slow jets in Swift J1727","Rising magnetic field inferred from Swift J1727 jet profile","Swift J1727's jets: slow, asymmetric, and evolving"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The jets are treated as steady-state flows during the observation, with the entire jet/counterjet difference caused by a time lag between epochs; if the asymmetry instead comes from intrinsic differences between the two sides, the inferred time evolution of the magnetic field would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Light travel time explains Swift J1727's jet asymmetry","New analytical model reveals slow jets in Swift J1727","Rising magnetic field inferred from Swift J1727 jet profile","Swift J1727's jets: slow, asymmetric, and evolving"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00133,"raw_usage":{"total_tokens":5471,"prompt_tokens":1067,"completion_tokens":4404,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":683,"completion_tokens_details":{"reasoning_tokens":4334}},"tokens_in":683,"tokens_out":4404,"duration_ms":33489,"temperature":1.0,"reasoning_tokens":4334,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:15:06.621903+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be to observe Swift J1727.8-1613 with VLBA at two epochs separated by about the light-travel lag, roughly one day, while the radio flux is changing: the steady-state-plus-delay model predicts that the counterjet profile at the later epoch should match the approaching-jet profile at the earlier epoch after beam convolution, and that the asymmetry should track the derivative of the radio light curve. A persistent asymmetry in the opposite sense to the flux derivative, or an asymmetry that persists when the flux is flat, would falsify the time-delay interpretation.","supporting_citations":[{"cited_title":"A., Stawarz , ., & Sikora , M","cited_arxiv_id":null,"evidence_quote":"Gives the Appendix-A formulation of the Blandford–Königl model from which the analytical $dF_\\nu/d|\\xi|$ formula is obtained."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the radio spectrum measurement with $\\alpha = 0.19 \\pm 0.07$ that fixes the magnetic-field index $b$ in the model."}],"review_version":1}