Pith. sign in

REVIEW 3 major objections 4 minor 20 references

Three-dimensional velocity fields in the silicon- and sulfur-reach ejecta in the remnant of Tycho supernova

T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Doppler shifts and proper motions reconstruct the 3D ejecta flow in Tycho's remnant, revealing a 20–30% front-to-back speed difference and distinct silicon versus sulfur distributions.

desk verdict First separate 3D velocity map for Si and S in Tycho, with careful Doppler maps, but the headline 20–30% asymmetry rests on an unnormalized sum that makes it grid-dependent. read the letter →

arxiv 2412.04096 v1 pith:O2ZKGZYJ submitted 2024-12-05 astro-ph.HE

O. Petruk , M. Patrii , T. Kuzyo , A. Baldyniuk , V. Marchenko , V. Beshley This is my paper · ORCID
classification astro-ph.HE
keywords supernovaremnantsTycho'sejectakinematicsDopplershiftspropermotionsilicon-richsulfur-richTypeIaasymmetry
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper reconstructs three-dimensional velocity fields for the silicon-rich and sulfur-rich ejecta in Tycho's supernova remnant by combining Doppler shifts of X-ray lines with proper motions of the remnant's edge. It reports that the speed of the ejecta on opposite sides of the remnant differs by 20–30%, and that the two elements have measurably different spatial distributions: silicon is more isotropic while sulfur has larger outward-directed components. If this is right, the Type Ia explosion that produced Tycho was only moderately asymmetric, and the progenitor's internal structure deviated from a perfectly layered, spherically symmetric configuration. Such a result gives a rare observational handle on the level of mixing and asymmetry in Type Ia explosions, which three-dimensional explosion models tend to predict as large.

What carries the argument

The load-bearing object is the two-zone line-of-sight asymmetry model described by Eqs. (4)–(6). In each cell of the projected remnant, the observed Doppler velocity $V_D$ is written as the sum of $\alpha v_x$ over the near half-space minus the sum of $(1/\alpha)|v_x|$ over the far half-space; solving the resulting quadratic for $\alpha$ gives the single multiplicative factor that turns a plane-symmetric expansion into one matching the observed blue- or redshift. The plane-of-sky components come from the proper motion of the remnant edge measured at every azimuth between two radio epochs, and the interior is filled by assuming homologous expansion, $v = (3V/4)\,r'/R'$, along ellipsoidal cross-sections whose line-of-sight size is set equal to the average projected radius.

What would settle it

Measure individual ejecta knots in three dimensions—matching proper-motion features in multi-epoch images with their Doppler shifts to get true space velocities—and count how many are moving toward versus away from the observer. A knot census that finds the front and back velocity distributions statistically identical would falsify the two-zone reconstruction and its 20–30% asymmetry.

Watch

Extended reading notes

Core claim

The central discovery, stated on the paper's own terms, is that the observed Doppler map and the measured expansion of Tycho's remnant imply a mildly asymmetric three-dimensional flow: the plasma speed on the near and far sides of the remnant differs by 20–30%, and the velocity fields of Si-rich and S-rich ejecta are not identical. In particular, the S-bearing material is preferentially shifted away from the observer, while the Si-bearing material is more isotropically distributed. The authors take this as evidence that the explosion asymmetry was moderate—enough to produce Doppler shifts of thousands of km/s but not an extreme global asymmetry—and that the progenitor's layered structure was disturbed, possibly by sloshing or mixing before the explosion.

Load-bearing premise

The reconstruction assumes that all ejecta on the near side share one velocity-scaling factor and all ejecta on the far side share its reciprocal; if the true line-of-sight flow is more complex, the 20–30% asymmetry is an artifact of that assumed shape.

Editorial extensions

If this is right

  • The 20–30% front-to-back speed difference becomes a benchmark that three-dimensional Type Ia explosion models must reproduce; models with much larger asymmetries can be ruled out for Tycho.
  • The different spatial distributions of Si and S imply the progenitor was not a perfectly layered sphere, so abundance and velocity maps like these can be compared with simulated pre-explosion mixing.
  • The reconstructed data cubes are directly usable as input for hydrodynamical models of Tycho's remnant evolution, connecting the explosion asymmetry to the present-day morphology.
  • The same Doppler-plus-proper-motion pipeline can be applied to other young supernova remnants with bright, spatially resolved ejecta lines, turning single remnants into a population of 3D kinematic probes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the asymmetry is real, the remnant's centroid of X-ray or radio emission should be slightly offset between approaching and receding sides; high-resolution broadband imaging could look for that surface-brightness asymmetry.
  • A natural testable extension is to apply the same two-zone decomposition to other ejecta lines such as Fe or Ar; if those show the same 20–30% ratio, the asymmetry is a global explosion property rather than a Si/S abundance effect.
  • Because the Doppler maps are brightness-weighted averages along the line of sight, the larger outward S velocity could partly reflect a brightness contrast between the front and back of the remnant rather than a pure kinematic difference; radiative-transfer modeling of the lines could separate the two.
  • The single-factor $\alpha$ could be generalized to a continuous line-of-sight velocity gradient; comparing that reconstruction with the two-zone one would quantify how much of the 20–30% asymmetry is driven by the model's simple form.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper reconstructs three-dimensional velocity fields for silicon- and sulfur-rich ejecta in Tycho's supernova remnant. Using Chandra ACIS spectra in a 20x20 grid, the authors map the Si XIII and S XV line centroids into line-of-sight velocities; using VLA radio images from 1994 and 2013, they measure proper motions of the remnant edge. They then combine these data with an assumed ellipsoidal remnant shape, a homologous radial flow profile, and a per-line-of-sight factor alpha that rescales the x-component of velocity in the near and far half-spaces. The central claims are that the ejecta speeds on opposite sides of the remnant differ by 20-30% and that the Si- and S-rich components have different three-dimensional spatial distributions.

Significance. The Doppler maps themselves appear carefully produced and agree with previously published maps, and the rest-energy analysis in Appendix 7.1 is a useful contribution. If the 3D reconstruction were reliable, the results would provide an important constraint on the degree of asymmetry in a Type Ia explosion and on mixing in the progenitor. However, as detailed below, the headline 20-30% asymmetry is not currently established because the central reconstruction equation lacks a proper normalization and the fitted parameter dominates the result.

major comments (3)
  1. [Sect. 4.2, Eqs. (4)-(6)] The observed Doppler velocity VD is a flux-weighted average of vx along the line of sight, but Eq. (4) writes it as a sum over grid elements. With S defined as a sum in Eq. (5), S scales with the number of resolution elements along x; consequently alpha from Eq. (6) and the inferred front/back asymmetry alpha - 1/alpha scale roughly as 1/N. Doubling the grid resolution would approximately halve the reported 20-30% asymmetry. The equations must be re-derived with S as a properly normalized mean (e.g., an integral divided by path length, or an emissivity-weighted average) and the reconstruction must be rerun to demonstrate resolution independence.
  2. [Sect. 4.2, Eq. (6), Figs. 9-10] Because alpha is solved separately for each line of sight from the observed VD, the projected Doppler map is reproduced exactly by construction. The 20-30% front/back speed difference is therefore an output of the assumed two-zone multiplicative model, not an independent measurement. The paper needs a validation step, such as forward-projecting the reconstructed cube through the same spectral pipeline and comparing predicted line widths or independent maps, and it should state explicitly what observational test could falsify the two-zone form.
  3. [Sect. 4.2, Eq. (2)] The unknown line-of-sight size of the remnant is set equal to the average projected radius, and this parameter enters R'(x, phi) and hence S and alpha. No sensitivity analysis is presented. The authors should vary the assumed ratio a/b over a plausible range (for example, 0.5 to 1.5) and show that the 20-30% conclusion and the Si/S spatial differences are stable; otherwise the quantitative claim is tied to an arbitrary geometric choice.
minor comments (4)
  1. [Abstract and title] The abstract, title, and text repeatedly use 'Si-reach' and 'S-reach'; these should be 'Si-rich' and 'S-rich'.
  2. [Sect. 3] The text states that 'the velocity for a cell is a sort of an average along the line of sight,' which is inconsistent with Eq. (4) where the velocity is written as a sum; this inconsistency should be resolved explicitly in the revised derivation.
  3. [Appendix 7.1] The statement that a 15 eV change in the rest energy is 'less than 1%' is misleading in velocity units: 15 eV at 1.86 keV corresponds to roughly 2400 km/s, comparable to the Doppler shifts reported here. The authors should quote the velocity equivalent and assess how the quoted rest-energy uncertainty affects the derived alpha values.
  4. [Sect. 3, Fig. 3] The figure caption refers to rstat values, and the text says values 'above 20' were corrected manually, but the caption for the 'worse' fit uses rstat = 51.4 while the text at that point says 'top' and 'bottom'; please make the labels and the fitting procedure unambiguous.

Circularity Check

2 steps flagged · score 4.0 of 10

Equation (4)-solved alpha reproduces the Doppler map by construction, so the 20-30% front/back asymmetry is a re-expression of the fitted LoS parameter rather than an independent measurement.

  1. fitted input called prediction [Section 4.2, Eqs. (4)-(6)]
    "we change the local values of vx along a given line of sight to v′x = αvx in the nearer (farther) half of SNR and to v′x = vx/α in the farther (nearer) half for VD > 0 (VD < 0), i.e. VD = Σ_{x>0} avx − Σ_{x<0} 1/a|vx| = (a − 1/a)S (4)... The value of α is given by a solution of a quadratic equation α = |VD|/(2S) + sqrt(V_D^2/(4S^2) + 1). (6)"

    Alpha is solved directly from the observed Doppler velocity VD in each cell, so the model's LoS velocity sum reproduces the input VD map by construction. The headline result that front and rear hemisphere speeds differ by 20-30% is a monotone function of alpha (alpha - 1/alpha roughly VD/S), i.e. it is the fitted parameter renamed as a physical finding, not an independent inference. The conversion from VD to alpha also depends on S, which is computed from the assumed plane-symmetric homologous model with the unknown LoS radius a set equal to the average projected radius, so the asymmetry value carries the modeling assumption rather than being a pure measurement.

  2. other [Section 4.2, Eq. (4) and preceding sentence]
    "The velocity VD in each cell on Fig. 6 is a residual of the sum of the vx components along the line of sight... we denoted the sum of vx along the LoS in one half-space S = Σ_{x>0} vx = Σ_{x<0} |vx|. (5)"

    Equation (4) equates the observed Doppler centroid VD, a velocity, to an unnormalized sum S of cell velocities over a half-space. If the grid has N cells along the LoS, S scales with N times the mean half-space velocity, so alpha from Eq. (6) satisfies alpha - 1/alpha ~ VD/S ~ VD/(N<vx>). Thus the inferred front/back asymmetry is inversely proportional to the discretization and to the assumed integration depth, making the 20-30% number an artifact of the chosen resolution and geometry unless S is defined as a properly normalized mean or integral.

full rationale

The paper is not globally circular: the Doppler maps and proper-motion speeds are external measurements, the Si vs S comparison uses two independently fitted lines, and the homologous expansion assumption is justified by standard Chevalier/Taylor results rather than by self-citation. However, the central 3D asymmetry claim is weakened by construction. In Eqs. (4)-(6), alpha is a per-cell fitting parameter chosen so that the model's LoS sum equals the observed Doppler velocity in every cell. The reported 20-30% difference between opposite sides is therefore an invertible function of this fitted alpha (and of the assumed S), not a prediction. Additionally, Eq. (4) is dimensionally inconsistent as written: an unnormalized sum of LoS velocities over cells cannot equal a Doppler centroid velocity, so the numerical value of alpha and hence the 20-30% asymmetry depend on grid resolution and on the arbitrary choice a = average projected radius. Because the independent Doppler maps and the Si/S differences still contain real observational content, the circularity is partial rather than total; score 4.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The reconstruction relies on four explicit assumptions: homologous expansion, an ellipsoidal shock surface with constant cross-section speed, a simplified two-zone scaling for line-of-sight asymmetry, and collisional ionization equilibrium for rest energies. Two free parameters, alpha per line of sight and the assumed LoS depth, control the shape of the derived 3D velocity field. No new physical entities are introduced.

free parameters (2)
  • alpha (per line-of-sight cell) = Per-cell values derived from observed Doppler velocity via Eq. (6)
    Multiplicative scaling of vx in the nearer and farther half-spaces, fitted so that the reconstructed Doppler velocity equals the observed one. This parameter encodes the entire front/back asymmetry of the 3D field.
  • Line-of-sight depth a = Set equal to the average projected radius in the yz plane
    The unknown extent of the remnant along the line of sight is chosen by hand as the average projected radius. This choice controls the integration depth and the velocity gradients in the reconstructed cube.
assumptions (4)
  • domain assumption Homologous expansion of the shocked ejecta, v proportional to r, with post-shock flow speed vs = 3V/4.
    Adopted in Sect. 4.2 based on Sedov and Chevalier solutions (refs. [18,19]); used in Eq. (3) to convert the shock speed V to interior flow velocities.
  • ad hoc to paper The shock speed along each radius within the half-ellipse cross-section equals the speed measured at the corresponding azimuth in the plane of the sky.
    Assumption i) in Sect. 4.2. This is a geometric modeling choice that directly determines the vy and vz components in the interior.
  • ad hoc to paper The line-of-sight asymmetry has the two-zone multiplicative form: vx scaled by alpha in one half-space and by 1/alpha in the other.
    Eqs. (4)-(6) in Sect. 4.2. This functional form is assumed, not derived from the physics of the remnant.
  • domain assumption Collisional ionization equilibrium for the line rest-energy calculation.
    Appendix 7.1 computes rest energies 1.856 keV (Si) and 2.450 keV (S) assuming CIE; the paper notes sensitivity to the ionization parameter of order 15 eV.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Three-dimensional velocity fields in the silicon- and sulfur-reach ejecta in the remnant of Tycho supernova." pith.science (2026). https://pith.science/paper/O2ZKGZYJ

@misc{pith2026241204096,
  author       = {Pith},
  title        = {Pith review of: Three-dimensional velocity fields in the silicon- and sulfur-reach ejecta in the remnant of Tycho supernova},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/O2ZKGZYJ}},
  note         = {Machine review of arXiv:2412.04096}
}
read the original abstract

The three-dimensional velocity structure of the shock-heated Si-reach and S-reach ejecta were reconstructed in Tycho supernova remnant from Doppler-shifted lines. The vector components along the line of sight were restored from the spatially resolved spectral analysis of the Doppler shifts of Si XIII and S XV lines. The components in the plane of the sky were derived from analysis of the proper motion of the remnant's edge at different azimuths. This has been done by using the data of X-ray observations from Chandra observatory as well as the radio data from the Very Large Array. Differences in Doppler velocities over the Tycho's SNR are of the order of thousands of km/s. The speed of the ejecta on the opposite sides of the remnant as a three-dimensional object differs on 20-30%. There are asymmetries and differences in the spatial distributions between the Si-reach and S-reach ejecta components. Namely, the level of isotropy is higher in Si while the vector components directed outward of the observer are larger in S. This puts limitations on the level of deviation of the internal structure of the progenitor star from the ideal layered structure as well as on the level of asymmetries in supernova explosion.

Figures

Figures reproduced from arXiv: 2412.04096 by the authors.

Figure 1
Figure 1. X-ray spectra of Tycho’s SNR from four epochs of observations marked [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. X-ray image of the remnant of Tycho’s supernova produced from the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Examples of shapes of Si lines in two cells, with a ‘good’ ( [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: The map of the parameter ampl from the model gauss1d (left) for the S line in comparison to the image of Tycho’s SNR in the photon energy range 2.3 − 2.6 keV (right). 5 [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: The map of the parameter pos for the Si line (left) and for the S line (right) [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: The map of the line-of-sight velocity for the Si line ( [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: North of Tycho’s SNR. Gray-scale shows the square root of the radio [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]
Figure 8
Figure 8. Figure 8: Geometry of the task. Axis x is directed toward the observer, z to the north, ϕ is an angle in the projection plane. 4.1 Expansion of SNR in the plane of the sky The shock speed V may be estimated by analyzing the proper motion of the SNR edge. It varies with azimuth i…
Figure 9
Figure 9. Figure 9: Cross-sections of the velocity data cubes by three planes passing through [PITH_FULL_IMAGE:figures/full_fig_p010_9.png]
Figure 10
Figure 10. Figure 10: The 3D velocity fields for Si (top panels) and S (bottom panels). [PITH_FULL_IMAGE:figures/full_fig_p010_10.png]
Figure 11
Figure 11. Figure 11: Intensity of emission for transitions (shown by colors) for Si [PITH_FULL_IMAGE:figures/full_fig_p015_11.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 19 canonical work pages

  1. [1]

    Woosley, T

    S. Woosley, T. Janka, Nature Physics 1, 147 (2005)

    work page 2005

  2. [2]

    Jordan et al., ApJ 681, 1448 (2008)

    I. Jordan et al., ApJ 681, 1448 (2008)

    work page 2008

  3. [3]

    Wongwathanarat, E

    A. Wongwathanarat, E. M¨ uller, H. T. Janka, A&A 577, A48 (2015)

    work page 2015

  4. [4]

    Reichert et al., MNRAS 518, 1557 (2023)

    M. Reichert et al., MNRAS 518, 1557 (2023)

    work page 2023

  5. [5]

    Supernova of 1572, Tycho’s Supernova

    A. Decourchelle, “Supernova of 1572, Tycho’s Supernova” in Handbook of Su- pernovae (Springer, 2017), p.117

    work page 2017

  6. [6]

    B. J. Williams et al. , ApJ 842, 28 (2017)

    work page 2017

  7. [7]

    T. Sato, J. P. Hughes, ApJ 840, 112 (2017) 13

    work page 2017

  8. [8]

    Godinaud, F

    L. Godinaud, F. Acero, A. Decourchelle, J. Ballet, A&A 680, A80 (2023)

    work page 2023

Show all 20 references
  1. [9]

    Godinaud, F

    L. Godinaud, F. Acero, A. Decourchelle, J. Ballet, arXiv e-prints, 2404.17296 (2024)

    2024 arXiv

  2. [10]

    Uchida et al., ApJ 962, 159 (2024)

    H. Uchida et al., ApJ 962, 159 (2024)

    2024

  3. [11]

    E. M. Reynoso et al. , ApJ 491, 816 (1997)

    1997

  4. [12]

    Katsuda et al

    S. Katsuda et al. , ApJ 709, 1387 (2010)

    2010

  5. [13]

    B. J. Williams et al. , ApJ 770, 129 (2013)

    2013

  6. [14]

    Tanaka et al

    T. Tanaka et al. , ApJL 906, L3 (2021)

    2021

  7. [15]

    B. J. Williams et al. , ApJL 823, L32 (2016)

    2016

  8. [16]

    Petruk et al

    O. Petruk et al. , ApJ 972, 63 (2024)

    2024

  9. [17]

    Z. Xue, B. E. Schaefer, ApJ 809, 183 (2015)

    2015

  10. [18]

    Taylor, Proc

    G. Taylor, Proc. Royal Soc. London Ser. A 201, 159 (1950)

    1950

  11. [19]

    R. A. Chevalier, ApJ 259, 302 (1982)

    1982

  12. [20]

    A. R. Foster, L. Ji, R. K. Smith, N. S. Brickhouse, ApJ 756, 128 (2012) 7 Appendix 7.1 Rest energies for the lines In order to translate the observed shifts in the energies of Si and S lines into the line-of-sight velocities, we need to know the reference (‘laboratory’ or ‘res...

    2012

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.