{"id":"867bdcf3-71e0-4a53-b45f-ba2f90cd90ae","arxiv_id":"2507.19042","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"EDSER and Richardson-Lucy deconvolution of Chandra HETG zeroth-order images reveal a precessing two-knot structure in SS 433's jets, and PSF-subtracted spectra show no Fe lines in the outer regions, supporting non-thermal X-ray emission.","lead":"Using Chandra HETG data sharpened with subpixel repositioning and image deconvolution, this paper maps clumps in SS 433's jets about 1.7 arcseconds from the core and finds that the outer X-ray emission is consistent with the known jet precession pattern. A smart generalist should read it because it claims the outer jets shine mainly through non-thermal processes, which links jet power to particle acceleration in a famous Galactic microquasar.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The single 3.7 keV MARX PSF drives both the RL deconvolution and the core-leakage subtraction; without a spectral PSF variation check, the non-thermal outer-region claim is not yet settled.","rationale":"The reader's verdict is CONDITIONAL, and my independent review lands on the same condition. The paper is careful, externally anchored (precession model, radio comparison, HEG-calibrated spectral model), and the analysis steps are described in detail. However, the strongest claim in the abstract and Section 4.4—that outer-region Fe lines are absent after 'accounting for potential core contamination' and that a non-thermal power-law is required—rests almost entirely on the subtracted leakage fractions. Those fractions come from a single mono-energetic MARX PSF. The weakness is not that MARX is wrong in general; it is that the paper does not demonstrate that a 3.7 keV PSF adequately represents the energy-dependent leakage of the spectrally hard, line-rich core into the outer regions. Because the fit band is 2-7.5 keV and the Fe-K lines at 6.4-6.95 keV are the decisive diagnostics, a small error in the leakage fraction at those energies directly maps onto the existence or absence of Fe lines and onto the power-law normalization. The paper does show awareness of PSF limitations (Appendix 4 states the deconvolved image is 'a plausible approximation'), but no test of PSF-induced systematics is reported. I do not see this as fatal: the radio morphology agreement, the precession-model alignment, and the consistency with previous detections of knots by Migliari et al. (2002, 2005) give independent weight to the existence of extended arcsecond-scale structure. But the non-thermal interpretation of the outer spectra is exactly the kind of claim that should be conditional on an energy-resolved PSF robustness test. The VLA non-simultaneity is a genuine but secondary issue, affecting the contrast comparison (Section 4.3) rather than the core spectral conclusion. My verdict recommendation stays CONDITIONAL: the main claim is plausible and worth publishing, but the non-thermal spectral interpretation should be explicitly gated on the PSF-energy check and an injection test.","tokens_in":20094,"tokens_out":2291,"duration_ms":20822,"concrete_test":"Recompute the core-leakage fractions and the outer East/West spectra using energy-dependent MARX PSFs: simulate PSFs at 1.5, 3.7, and 6.4 keV (and ideally at the Fe XXV/XXVI line energies 6.7/6.95 keV), convolve the observed Center spectrum through each, and derive pcenter->east and pcenter->west per energy band before summing. Then re-run the Section 3.4 subtraction and fit. If the Fe-line normalization Agau in East/West changes from 0.0 to a value inconsistent with zero at 90% confidence, or if the power-law Apl becomes consistent with zero, the central non-thermal claim fails. As a second check, perform an injection test: place two simulated knot sources at 1.7 arcsec east/west with the claimed fluxes into the observed field, apply the same EDSER+RL pipeline, and verify the knots and their recovered fluxes are not artifacts of PSF mismatch or RL ringing.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central spectral conclusion—that outer East/West regions show no Fe lines and require a power-law component—depends on subtracting the core contribution with fixed leakage fractions pcenter->east ~0.130 and pcenter->west ~0.114 (Section 3.4). These fractions are computed from a MARX PSF simulated at a single energy, 3.7 keV, for the full 0.5-8 keV band. This is the paper's primary soft spot. The core is Fe-line-dominated and the PSF energy dependence is strong: the Chandra HETG zeroth-order PSF broadens and its core fraction changes with energy, so a monochromatic 3.7 keV PSF cannot accurately predict the fraction of 6.4-6.7 keV core line photons that leak into the outer regions. If the true Fe-K leakage is larger than the 3.7 keV estimate, the subtraction is incomplete and the remaining Fe-K residuals would mimic a power-law continuum in the limited 2-7.5 keV band. Conversely, if the 3.7 keV PSF overestimates leakage at line energies, the Fe lines could be over-subtracted and artificially erased. The same PSF is also used for the RL deconvolution that produces the two knotted structures (Section 2.2.1, Appendix 1), so a PSF error could affect both the morphology and the spectral conclusion. The paper also lacks an independent check: no injection/recovery test of simulated knots, no PSF energy-resolved verification, and no alternative PSF simulation (e.g., from flight calibration or at 6.4 keV). While the 3.7 keV choice is reasonable as a band-average, its adequacy is assumed rather than demonstrated. The non-simultaneous VLA comparison (2003 vs 2014) is a secondary concern about intensity contrast, not about the core spectral claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a spatial and spectral analysis of arcsecond-scale X-ray emission from SS 433 using zeroth-order Chandra HETG data, focusing on a 138 ks observation from 2014 (Obs. ID 15781). The authors apply energy-dependent subpixel event repositioning (EDSER) and Richardson-Lucy (RL) deconvolution with a MARX-simulated PSF at a monochromatic energy of 3.7 keV. They report two knotty east/west structures at ~1.7 arcsec from the core, consistent with the external kinematic precession model for an ejection age of ~200 days. Spatially resolved zeroth-order spectra of the East and West regions, after subtracting PSF-leaked core light using model-derived leakage fractions p_center->east ~ 0.130 and p_center->west ~ 0.114, show no significant Fe lines and require a power-law component, leading the authors to suggest a dominant non-thermal contribution to the outer X-ray emission. The analysis is complemented by a comparison with VLA radio data at a similar precession phase, a discussion of phase-dependent jet visibility, and an exponential decay timescale for the brightness profile.","tokens_in":20456,"tokens_out":3819,"duration_ms":40133,"significance":"If the central results hold, the paper provides an important step in resolving the long-standing ambiguity about the thermal versus non-thermal origin of arcsecond-scale X-ray emission in SS 433, with implications for jet physics and particle acceleration in microquasars. The work is also methodologically interesting as an application of EDSER and RL deconvolution to HETG zeroth-order data, and the comparison with an external kinematic model avoids circularity. The paper is clearly written, the data reduction is documented, and the statistical errors are carefully reported. The main spectral conclusion—absence of Fe lines in the outer regions and presence of a power-law component—is, however, contingent on the reliability of the PSF model used for both deconvolution and core-leakage subtraction; this is the principal weakness that needs to be addressed before the non-thermal interpretation can be considered robust.","major_comments":[{"comment":"The leakage fractions p_center->east ~0.130 and p_center->west ~0.114, which are subtracted from the outer-region spectra before concluding that Fe lines are absent, are computed from a MARX PSF simulated at a single energy of 3.7 keV, while the core spectrum is dominated by Fe K photons at 6.4–6.7 keV. The Chandra HETG zeroth-order PSF is known to have significant energy dependence, and the core fraction of the PSF changes with energy; hence the leakage fraction at Fe-line energies may differ substantially from the 3.7 keV value. If the true Fe K leakage is larger than the adopted value, residual Fe lines would remain in the East/West spectra and could imitate a power-law continuum over the fitted 2–7.5 keV band; if it is smaller, the subtraction would artificially erase real Fe lines. The paper does not provide an energy-resolved PSF check (e.g., a MARX simulation at 6.4 keV) or an injection/recovery test to validate the leakage subtraction. Because the non-thermal outer-spectrum claim directly rests on this subtraction, this issue is load-bearing for the central spectral conclusion.","section":"Section 3.4 and Section 2.2.1"},{"comment":"The RL deconvolution that produces the two knotty structures at ~1.7 arcsec uses the same monochromatic 3.7 keV MARX PSF. The stopping rule at 30 iterations is justified by a reduced chi-square convergence criterion and by a PSF deconvolution test in Appendix 4, but no independent validation shows that the knots are real rather than artifacts of deconvolution with a potentially mismatched PSF. An injection/recovery simulation—inserting simulated point-like or knot-like sources with known positions and fluxes into the observed event data and running the full EDSER+RL pipeline—would directly test whether the pipeline recovers such structures and whether the chosen iteration count suppresses or creates spurious features. Such a test is particularly important because the paper itself notes in Appendix 4 that the RL image is 'a plausible approximation' rather than a definitive reconstruction.","section":"Section 2.2.1 and Appendix 1"}],"minor_comments":[{"comment":"The X-ray/radio ratio map in Figure 8(c) is constructed from observations separated by ~6 years (2014 Chandra and 2003 VLA) with slightly different precession phases (0.43 vs 0.47) and different orbital phases. Although the authors acknowledge the non-simultaneity, a sentence in the text or caption quantifying the possible positional/phase misalignment and its effect on the ratio map would improve clarity.","section":"Section 4.3 and Figure 8"},{"comment":"The beaming correction in Equation (A4) adopts n=2 and alpha = Gamma-1 = 0.9, where Gamma = 1.9 is taken from the paper's own spectral fit. It would be helpful to state explicitly that the correction is therefore model-dependent and to note how the brightness profile would change if, for example, alpha were varied within the 90% confidence range.","section":"Appendix 3"},{"comment":"The lower limits on Abapec and Agau in the East and West regions are consistent with zero at 90% confidence, so the statement in the abstract that 'Fe lines are not evident' is appropriate, but the conclusion section's phrasing 'no Fe lines' could be slightly softened to 'no significant Fe lines' to match the non-detection nature of the measurement.","section":"Table 3 and Section 4.4"},{"comment":"The definition of the 1-sigma extent for the original and deconvolved PSF is clear, but it would be useful to report whether this measure is circular or elliptical and whether the same value 0.246 arcsec is used in the correction factor calculation in Appendix 3; the text implies this but does not explicitly restate it.","section":"Appendix 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely publishable in PASJ after a revision that addresses the PSF energy-dependence issue. The reader's stress-test concern is valid and central: the single-energy 3.7 keV MARX PSF drives both the deconvolution and the core-leakage subtraction, and without an energy-resolved PSF or injection/recovery test, the non-thermal outer-spectrum conclusion is not yet settled. I would suggest asking the authors to add a 6.4 keV MARX PSF leakage calculation and a simple injection/recovery test for the RL deconvolution. The VLA comparison, while non-simultaneous, is clearly flagged as such and is more of a secondary result; the main revision should focus on the spectral claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this paper reports a genuinely new spectral result for SS 433's arcsecond-scale jets—after subtracting core light with a PSF model, the outer east/west regions show no iron lines and need a power-law component, pointing to non-thermal emission rather than the in-situ reheating that earlier work inferred. The two-knot morphology is also new, though it depends on image deconvolution.\n\nWhat's good: the analysis is careful and well documented. The kinematic comparison uses an external precession model and radio ephemerides, so the match at ~1.7 arcsec is not circular. The VLA overlay, even though non-simultaneous, provides useful context, and the paper is transparent about its assumptions. The spectral fitting is detailed and the uncertainties are honestly reported.\n\nThe soft spots are real but do not kill the claim. The biggest one is the MARX PSF simulated at a single energy (3.7 keV) used both for the Richardson-Lucy deconvolution and for computing the core-leakage fractions that are subtracted from the outer spectra. The Chandra zeroth-order PSF does change with energy, and the core is Fe-line dominated, so the leaked fraction of 6.4-6.7 keV photons could differ from the 3.7 keV value. If leakage is underestimated, residual iron lines could masquerade as a power law; if overestimated, they could be artificially erased. The paper offers no injection/recovery tests or alternative PSF checks. This is addressable and I'd want to see it before treating the non-thermal interpretation as settled, but it is not obviously wrong.\n\nTwo minor issues: the VLA data from 2003 are compared to 2014 X-rays, so some brightness contrast differences could be variability. And the beaming correction uses the paper's own spectral index, which is mild self-reference but not a serious circularity.\n\nOverall: this is a solid observational paper that moves the debate. It deserves a serious referee; the main result should be published after robustness tests on the PSF. I'd bring it to our reading group.","headline":"Solid SS 433 study with a genuinely new non-thermal spectral result, but the mono-energetic PSF assumption needs robustness tests before the claim is settled.","tokens_in":21118,"tokens_out":2123,"would_cite":false,"duration_ms":22161,"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":"SS 433's arcsecond-scale X-ray emission is dominated by non-thermal processes.","keywords":["SS 433","X-ray jets","Chandra HETG","Richardson-Lucy deconvolution","EDSER subpixel repositioning","non-thermal X-ray emission","precessing jets","X-ray binaries"],"falsifier":"Recompute the deconvolution and the East/West leakage subtraction using energy-resolved PSFs evaluated at the iron-line energies (6.4, 6.7, and 6.95 keV); if the two knots disappear under the energy-resolved PSF, or if the leakage fractions at those energies exceed the assumed $\\sim0.11$--$0.13$ enough to hide the Fe lines, the non-thermal conclusion fails. A direct check would be a new Chandra observation at the same precession phase with the jets even closer to the sky plane, asking whether the two knots appear without any deconvolution.","tokens_in":19863,"feed_emoji":"🔭","tokens_out":13600,"duration_ms":117707,"temperature":0.7,"pith_summary":"SS 433 is a Galactic X-ray binary whose relativistic jets precess with a 162-day period, but at the arcsecond scale the bright core has made the jet's own X-ray emission hard to isolate. This paper analyzes the deepest available Chandra HETG observation, taken in 2014 during an orbital eclipse when the jets were nearly perpendicular to the line of sight. Using energy-dependent subpixel repositioning plus Richardson-Lucy deconvolution, it reveals two knot-like structures roughly 1.7 arcseconds east and west of the core, consistent with the kinematic precession model and an ejection age of about 200 days. After subtracting the core light that leaks into the outer regions through the telescope PSF, the outer X-ray spectra show no iron lines and require a power-law component, implying a dominant non-thermal contribution at arcsecond scales. The result matters because it turns a previously ambiguous, core-contaminated signal into a concrete, testable picture of where and how SS 433's jets emit X-rays.","feed_headline":"SS 433's arcsecond X-ray knots are non-thermal","feed_subtitle":"Sharpened 2014 Chandra data show the knots match the precession model and lack iron lines.","key_machinery":"The central object is an imaging chain, not a single identity: EDSER (energy-dependent subpixel event repositioning), which uses the detector's charge-diffusion model to assign each photon a subpixel impact position, followed by Richardson-Lucy deconvolution, an iterative Poisson-aware algorithm that sharpens the image by convolving a trial image with the PSF and comparing to the data. The PSF is a ray-trace simulation of the telescope optics evaluated at one energy, 3.7 keV, and the stopping point of 30 iterations is set by a reduced-chi-squared convergence rule. This same PSF is integrated over the East and West extraction regions to compute the core-leakage fractions $p_{\\rm center\\to east}\\sim0.130$ and $p_{\\rm center\\to west}\\sim0.114$, which are subtracted before the no-iron-line conclusion is drawn. The second load-bearing object is the kinematic precession model (jet speed $0.2602c$, precession period $162.15$ days, half-opening angle $19.85^\\circ$, inclination $78.83^\\circ$), whose sky projection, corrected by the light-travel-time relation $\\tau=t_{\\rm age}/(1-\\beta_\\ell)$, lets the knot positions be translated into an ejection age.","core_discovery":"On the paper's own terms, the central discovery is that the apparent arcsecond-scale X-ray emission of SS 433 is not simply scattered core light or purely thermal jet plasma. In the 2014 observation (Obs. ID 15781, $\\sim$140 ks, orbital eclipse, precession phase 0.43--0.44), EDSER repositioning and 30 Richardson-Lucy iterations produce two knot-like structures at $\\sim$1.7 arcsec ($\\sim 10^{17}$ cm) east and west of the core. Overlaid on the kinematic precession model with the light-travel-time correction $\\tau = t_{\\rm age}/(1-\\beta_\\ell)$, both knots fall on the predicted helical jet paths and share an ejection age $t_{\\rm age}\\sim 200$ days. Spectra extracted from the East and West regions, after subtracting PSF leakage fractions of $\\sim$0.130 and $\\sim$0.114 of the core spectrum, contain no Fe lines, while the core spectrum is well fitted by thermal plasma plus a fluorescent Fe I K$\\alpha$ line; a power-law component with photon index $\\Gamma\\simeq 1.9$--$2.0$ is required in the outer regions. The paper interprets this as thermal emission dominating the core and non-thermal, likely synchrotron, emission dominating the arcsecond-scale jet, and notes that the X-ray/radio spectral index from the ratio map, $\\alpha\\sim 0.6$--$0.9$, is consistent with that interpretation.","pith_inferences":["Our inference: if a single synchrotron electron population emits both bands, the faster X-ray decay can be converted into a magnetic-field estimate, since synchrotron cooling scales as $t_{\\rm cool}\\propto B^{-3/2}\\nu^{-1/2}$; a multi-epoch simultaneous radio/X-ray campaign would test this directly.","Our inference: the EDSER-plus-Richardson-Lucy pipeline should transfer to other X-ray binaries whose bright cores mask faint jets; applied with energy-resolved PSFs, it could resolve arcsecond knots in systems where pile-up previously blocked spatially resolved spectroscopy.","Our inference: the single-energy PSF assumption is the point most likely to affect the scientific conclusion, so redoing the leakage subtraction with energy-resolved PSFs at 6.4, 6.7, and 6.95 keV would either harden or overturn the non-thermal interpretation of the outer spectra.","Our inference: the 24-epoch dataset contains far more precession-phase information than the single epoch analyzed in detail here, so a systematic multi-epoch deconvolution could yield a tomographic map of SS 433's ejection history."],"forward_implications":["The two knots being symmetric in position and sharing one ejection age implies that a single ejection event around 200 days before the 2014 observation produced observable X-ray structure on both sides of the core.","The absence of Fe lines in the outer regions after PSF subtraction implies that the arcsecond-scale X-ray emission cannot be modeled as thermal plasma alone; any complete jet model must include a non-thermal power-law component.","The measured X-ray decay timescale, $\\tau'=35.3\\pm3.3$ days, being shorter than the radio value of $55.9\\pm1.7$ days, implies that X-ray-emitting electrons lose energy on a faster timescale, which constrains cooling in the jet if radio and X-rays share an electron population.","The visibility modulation seen across 24 observations implies that single-epoch images of SS 433's jets are phase-dependent: knots will be easiest to detect during eclipse and when Doppler beaming is minimized.","The spatial match between the X-ray knots and radio contours at similar precession phase implies that simultaneous future radio/X-ray observation can map the transition from core thermal to jet non-thermal emission directly."],"supporting_citations":[{"why":"Defines the EDSER energy-dependent subpixel event repositioning algorithm that produces the sub-pixel images.","marker":"Li et al. (2004)"},{"why":"Introduces the Richardson-Lucy deconvolution scheme used to sharpen the zeroth-order image.","marker":"Richardson (1972)"},{"why":"Adds the Poisson-noise-aware iterative scheme that the paper's Equation (3) applies.","marker":"Lucy (1974)"},{"why":"Supplies the ray-trace PSF simulation used for deconvolution and for the core-leakage fractions.","marker":"Davis et al. (2012)"},{"why":"Provides the reduced-chi-squared convergence criterion that selects 30 Richardson-Lucy iterations.","marker":"Marchenko et al. (2017)"},{"why":"Fixes the precession period and jet geometry in the kinematic model used to interpret knot positions.","marker":"Gies et al. (2002)"},{"why":"Gives the light-travel-time relation $\\tau=t_{\\rm age}/(1-\\beta_\\ell)$ and the comparable radio brightness plateau at 200-250 days.","marker":"Roberts et al. (2010)"},{"why":"Supplies the exponential radio decay timescale and the projection/beaming correction framework applied to the X-ray profile.","marker":"Bell et al. (2011)"},{"why":"First reported the arcsecond-scale east/west X-ray peaks that this paper resolves into knots.","marker":"Migliari et al. (2002)"},{"why":"Reported evidence for extended non-thermal X-ray emission that the paper's power-law result supports at arcsecond scales.","marker":"Khabibullin & Sazonov (2017)"}],"fun_headline_variants":["SS 433 jet knots are non-thermal, not core echoes","Non-thermal X-ray knots trace SS 433 jet ejection","Chandra sharpens SS 433: knots match precession model","SS 433 knots lack iron, pointing to non-thermal jets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the simulated point-spread function at a single energy, 3.7 keV, accurately matches the true Chandra HETG zeroth-order PSF after EDSER over the full 0.5--8 keV band, because that PSF is used both to deconvolve the image into knots and to compute the core-leakage fractions whose subtraction leaves the outer spectra with no iron lines.","fun_headline_variants_meta":{"raw":{"variants":["SS 433 jet knots are non-thermal, not core echoes","Non-thermal X-ray knots trace SS 433 jet ejection","Chandra sharpens SS 433: knots match precession model","SS 433 knots lack iron, pointing to non-thermal jets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000939,"raw_usage":{"total_tokens":4152,"prompt_tokens":1221,"completion_tokens":2931,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":837,"completion_tokens_details":{"reasoning_tokens":2857}},"tokens_in":837,"tokens_out":2931,"duration_ms":22555,"temperature":1.0,"reasoning_tokens":2857,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:02:12.343903+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the deconvolution and the East/West leakage subtraction using energy-resolved PSFs evaluated at the iron-line energies (6.4, 6.7, and 6.95 keV); if the two knots disappear under the energy-resolved PSF, or if the leakage fractions at those energies exceed the assumed $\\sim0.11$--$0.13$ enough to hide the Fe lines, the non-thermal conclusion fails. A direct check would be a new Chandra observation at the same precession phase with the jets even closer to the sky plane, asking whether the two knots appear without any deconvolution.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Richardson-Lucy deconvolution scheme used to sharpen the zeroth-order image."},{"cited_title":"E., Bautz, M","cited_arxiv_id":null,"evidence_quote":"Supplies the ray-trace PSF simulation used for deconvolution and for the core-leakage fractions."},{"cited_title":"E., Ostrowski, M., et al","cited_arxiv_id":null,"evidence_quote":"Provides the reduced-chi-squared convergence criterion that selects 30 Richardson-Lucy iterations."},{"cited_title":"R., McSwain, M","cited_arxiv_id":null,"evidence_quote":"Fixes the precession period and jet geometry in the kinematic model used to interpret knot positions."},{"cited_title":"H., Wardle, J","cited_arxiv_id":null,"evidence_quote":"Gives the light-travel-time relation $\\tau=t_{\\rm age}/(1-\\beta_\\ell)$ and the comparable radio brightness plateau at 200-250 days."},{"cited_title":"R., Roberts, D","cited_arxiv_id":null,"evidence_quote":"Supplies the exponential radio decay timescale and the projection/beaming correction framework applied to the X-ray profile."},{"cited_title":"2002, Science, 297, 1673","cited_arxiv_id":null,"evidence_quote":"First reported the arcsecond-scale east/west X-ray peaks that this paper resolves into knots."},{"cited_title":"I., & Sazonov, S","cited_arxiv_id":null,"evidence_quote":"Reported evidence for extended non-thermal X-ray emission that the paper's power-law result supports at arcsecond scales."}],"review_version":2}