Pith. sign in

REVIEW 4 major objections 5 minor 2 references

Dynamic Imprints of Colliding-wind Dust Formation from WR140

T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

Pith's one-line read WR 140's infrared dust shells are genuine outflows expanding at a steady 2,600 km/s, roughly matching the Wolf-Rayet wind, with clumpy substructures that persist for at least 14 months.

desk verdict First two-epoch JWST measurement of WR 140's dust shell kinematics is a real advance with a robust central result, but the abstract's numbers are internally inconsistent and the PSF-subtraction systematic is unquantified. read the letter →

arxiv 2502.02738 v1 pith:S7RYVYWB submitted 2025-02-04 astro-ph.GA astro-ph.SR

classification astro-ph.GAastro-ph.SR
keywords Wolf-RayetbinariesdustformationcollidingwindsWR140JWSTMIRIpropermotionsshellsinfraredimaging
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

The paper aims to prove that the bright infrared rings around the Wolf-Rayet binary WR 140 are real dust shells ejected by the binary, not artifacts of JWST's optics. By comparing two JWST/MIRI images taken 14 months apart, it measures the shells' outward proper motion and finds they expand at a steady ~2,600 km/s, close to the WC7 star's wind speed, with no measurable acceleration. The authors also show that lumpy substructures within the shells survive from one epoch to the next, evidence that clumping in the colliding-wind region is what enables dust formation.

What carries the argument

The central objects are the dust shells themselves, traced by their thermal infrared emission at 7.7, 15, and 21 μm. The argument runs on a careful comparison of two epochs: a PSF-subtraction step removes the bright diffraction spikes using a detector-sampled model of the MIRI point-spread function; the Cycle 2 images are aligned to the Cycle 1 World Coordinate System to within 0.1 pixels; then ten radial rays per dust feature are median-combined into radial flux profiles, whose peaks are Gaussian-fitted to get shell positions in both cycles. The difference in peak positions yields proper motions and, with an assumed distance of 1.64 kpc, projected velocities.

What would settle it

A re-analysis of the same images with an alternative PSF model built from isolated field stars, or with a different scaling of the diffraction-spike model, that shifts the fitted shell peaks by more than the quoted 10–30 mas uncertainties at 7.7 μm would falsify the reported proper motions. Alternatively, a third JWST epoch whose shell positions do not continue the linear outward trend at ~390 mas yr⁻¹ would falsify the steady-expansion claim.

Watch

Extended reading notes

Core claim

The authors claim that the dust shells around WR 140 are astrophysical outflows originating from the binary, expanding at constant projected velocities of roughly 2,500 to 2,700 km/s (median 2586 ± 152 km/s at 7.7 μm along the C1 feature), consistent with the terminal wind speed of the WC7 star. The measured proper motions are the same across all visible shells, indicating no significant acceleration or deceleration, and identifiable clumps within the shells persist for at least 14 months. These observations confirm the shells' physical reality and support the hypothesis that dust formation in the wind-collision region requires clumping.

Load-bearing premise

The shell positions are measured only after subtracting a model PSF whose scaling is wavelength- and radius-dependent, so if residual diffraction-spike pattern shifts the apparent peaks with radius, the measured proper motions would be systematically biased.

Editorial extensions

If this is right

  • The dust shells are genuine physical structures, so future imaging can use them as tracers of the binary's wind history over the past century.
  • Steady expansion at the wind speed implies the local medium is tenuous; WR 140 may sit inside a cavity cleared by its own wind, with little drag on the shells.
  • Clumpy substructures surviving for over a year means clump lifetimes are long enough to be observable tracers of dust formation in the wind-collision region.
  • The PSF-subtraction and image-alignment pipeline can be applied to other bright, dusty point sources observed with JWST.
  • A third epoch of observations would turn the proper-motion measurement into a direct acceleration test and refine grain-size interpretations.

Reading between the lines

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

  • If the expansion velocities are truly constant across all shells, the spacing between shells records the binary's past periastron passages, so a detailed spacing model could date each shell and test the orbital history of WR 140.
  • The speed differences among the C1, C0, and E features may reflect grain-size sorting; a radial color analysis using the 7.7 and 21 μm images could test whether slower features host smaller grains.
  • The reported proper motions depend on the assumed distance of 1.64 kpc; adopting a 10% larger distance, as hinted by Gaia DR3, would raise the velocities by about 200 km/s, so an improved distance would sharpen the comparison with the wind speed.
  • The PSF-subtraction method relies on scaling a model to the observed diffraction spikes; an independent PSF reconstruction from isolated field stars would provide a cross-check on the absolute velocity scale.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. This Letter compares two JWST/MIRI epochs of WR 140 (2022 July and 2023 September) to measure the outward motion of the episodic dust shells. After subtracting a WebbPSF model, aligning the Cycle 2 images to Cycle 1 with JHAT, and extracting median radial flux profiles along the C1, C0, and E dust features, the authors derive shell-by-shell proper motions and projected velocities. They report median C1 velocities of 2586±152 km/s at 7.7 μm, 2350±109 km/s at 15 μm, and 2547±61 km/s at 21 μm, with lower median speeds along C0 and deprojected speeds along E, and conclude that the shells expand at a steady rate close to the WC7 wind speed and that clumpy substructures persist between the two epochs.

Significance. Two-epoch JWST imaging of the well-known episodic dust-making WC binary WR 140 is a valuable dataset, and the analysis relies on direct image differencing rather than model-dependent reconstruction. The agreement among the three wavelengths, the explicit use of previously published distance/inclination/cone-angle inputs, and the public release of data and PSF-subtraction code are notable strengths. If the velocities are correct, the result provides strong confirmation that the shells are astrophysical dust ejecta moving at roughly the terminal wind speed, and the persistence of clumps over 14 months supports clumping-based dust formation models. The main limitations are internal number inconsistencies and the need to quantify PSF-subtraction systematics; these are fixable and do not invalidate the core observational picture.

major comments (4)
  1. [Abstract; Section 4.1] The headline numbers are internally inconsistent. The abstract quotes an average C1 proper motion of 390±29 mas/yr and a projected velocity of 2714±188 km/s, while Section 4.1 reports median values of 332±20 mas/yr and 2586±152 km/s at 7.7 μm. The two proper motions differ by roughly 60 mas/yr, which is 2–3 times the stated statistical errors, and the abstract's velocity corresponds to about 350 mas/yr at 1.64 kpc, not 390 mas/yr. This suggests sensitivity to whether a mean or median is used, which shells are included, or the PSF-scaling choice. The abstract and body must be reconciled, and the exact shell sample and statistic must be stated.
  2. [Section 3.2, Figure 2] The PSF-subtraction systematic at 7.7 μm is not propagated into the quoted uncertainties. Section 3.2 and Figure 2 state that the 7.7 μm diffraction-spike flux varies with radius and that the scaling produced oversubtraction in some areas, yet the method uses a single average scale factor for this filter. Any residual PSF pattern is therefore radius-dependent, and if the telescope roll angle or PSF model errors differ between the two epochs, the residual pattern can shift Gaussian peak positions by different amounts in Cycle 1 and Cycle 2, biasing the proper motions. The quoted ±15–20 mas/yr errors contain no term for this. I request a quantitative test, for example varying the PSF scale factor or comparing independent subtraction recipes, and an estimate of the resulting systematic uncertainty on each velocity.
  3. [Table 1, C1 rows 13/15/16] The 7.7 μm velocities for C1 shells 13, 15, and 16 are 3773±166, 3365±210, and 3363±373 km/s, respectively, well above the adopted WC7 terminal speed of 2860 km/s. The text's statement that outlying measurements are 'likely due to low signal to noise in the outermost shells as well as the effect of background stars' is not supported by quantitative evidence for these specific shells, and the table does not flag them as excluded. If they are astrophysical, the claim that speeds are constant across all visible shells is not correct; if they are contaminated, a source check or explicit exclusion criterion is needed. This must be addressed directly.
  4. [Section 4.1, Equation (1), Table 1 (E feature)] The reported 'deprojected' E velocities are numerically identical to the projected velocities. For example, 289 mas/yr at 1.64 kpc corresponds to about 2250 km/s, and Table 1 reports 2252±490 km/s; similarly, 302 mas/yr gives about 2347 km/s and 271 mas/yr gives about 2107 km/s, matching the 15 and 21 μm entries. With i=119.6° and θ=40°, Equation (1) would multiply by sec(20.4°)≈1.07, not by unity. The E-feature deprojection therefore appears not to have been applied, or the equation and angle definitions are misstated, and the velocity comparison across C1/E/C0 is not valid as presented.
minor comments (5)
  1. [Section 4.1, Equation (1)] Equation (1) is garbled in the printed text; rewrite it with clear parentheses and define all angles explicitly.
  2. [Figure 2 caption; Section 3.2] The caption says the 7.7 μm PSF subtraction required 'complex scaling,' but Section 3.2 describes a single average scale factor for that filter; please clarify what was complex.
  3. [Table 1, E feature] The rows labeled 7a and 9a for the E feature are empty; add a footnote explaining that these shells were not measured.
  4. [Table 1 note] The table note contains the typo 'for for'; also state explicitly whether the E columns are deprojected, to avoid confusion with the C1 and C0 columns.
  5. [Figure 5] The horizontal axis is labeled only as 'shell'; label the shell numbers so that the points can be cross-referenced with Table 1.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: shell proper motions are direct two-epoch image-difference measurements, with external inputs stated.

full rationale

The kinematics measurement is self-contained: shell positions are Gaussian fits to radial flux profiles in each epoch, and proper motions are the direct differences divided by the 1.12 yr baseline. The PSF-subtraction scaling is a nuisance calibration, not a fit to the shell velocities. Adopted astrophysical inputs—distance 1.64 kpc (Bailer-Jones et al. 2018), inclination 119.6 deg (Fahed et al. 2011), cone half-angle 40 deg (Han et al. 2022)—are stated external parameters; although several of those papers share authors with this work, none of these inputs encodes the measured proper-motion result, and the paper's own deprojection formula (Eq. 1) would propagate any change in them linearly rather than being self-validating. The paper's self-citations (Lau et al. 2022 for shell numbering; Han et al. 2022 for geometry; Lau et al. 2023 for grain-size interpretation) are contextual, not load-bearing for the new velocity measurement. The acknowledged PSF oversubtraction at 7.7 microns (Sect. 3.2) and the discrepancy between the abstract's mean 390 +/- 29 mas/yr and the body's median 332 +/- 20 mas/yr are real statistical/systematic concerns, but they are accuracy issues, not circularity: no quantity in the derivation is defined in terms of the target result. No circular step can be exhibited.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

No new physical entities are introduced. The analysis relies on standard astronomical inputs (distance, inclination, cone angle) and on the fidelity of the PSF model and image alignment. The free parameters listed are adopted or chosen values that affect the derived velocities, not parameters fitted to the data in a model-dependent way.

free parameters (3)
  • Distance to WR 140 = 1.64 kpc (Gaia DR2, Bailer-Jones et al. 2018)
    Adopted distance used to convert proper motions to physical velocities; a 10 percent larger DR3 distance would increase velocities by ~200 km/s. This is an input, not fitted here, but the choice affects the reported speeds.
  • PSF scale fractions for 15 and 21 μm = 1/5 and 1/8 of peak flux ratio
    Chosen ad hoc to scale the WebbPSF model for the longer wavelengths; affects PSF subtraction quality but not directly the measured shell positions.
  • Deprojection angles for the E feature = i = 119.6 deg, theta = 40 deg
    Adopted from Fahed et al. 2011 and Han et al. 2022 to deproject E-feature velocities; uncertainty in these angles propagates into the deprojected speeds, though the correction factor is small (~1.07).
assumptions (5)
  • domain assumption The WebbPSF model accurately represents the MIRI PSF at the relevant wavelengths and the scaling method yields a clean subtraction without radial residuals.
    Section 3.2; if the PSF subtraction leaves radius-dependent residuals, the measured peak positions and proper motions could be biased.
  • domain assumption The C1 and C0 dust features move predominantly in the plane of the sky, so no deprojection is needed for their projected velocities.
    Section 4; based on geometrical modeling in Han et al. 2022 and Lau et al. 2022.
  • domain assumption The radial flux profile peaks correspond to the same physical shells in both epochs, and the Gaussian fits accurately locate the shell centroids.
    Section 4.1; if peak identification is inconsistent between epochs, the proper motions are wrong.
  • domain assumption The adopted distance of 1.64 kpc is correct.
    Section 4; the velocities scale linearly with distance.
  • domain assumption The nonstandard pipeline settings for saturated pixels produce reliable centroid positions for the central source.
    Section 3.1; the authors found empirically that the settings work for WR 140, but this is not independently validated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Dynamic Imprints of Colliding-wind Dust Formation from WR140." pith.science (2026). https://pith.science/paper/S7RYVYWB

@misc{pith2026250202738,
  author       = {Pith},
  title        = {Pith review of: Dynamic Imprints of Colliding-wind Dust Formation from WR140},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/S7RYVYWB}},
  note         = {Machine review of arXiv:2502.02738}
}
abstract

Carbon-rich Wolf-Rayet binaries are a prominent source of carbonaceous dust that contribute to the dust budget of galaxies. The "textbook" example of an episodic dust producing WR binary, WR140 (HD193793), provides us with an ideal laboratory for investigating the dust physics and kinematics in an extreme environment. This study is among the first to utilize two separate JWST observations, from Cycle 1 ERS (July 2022) and Cycle 2 (Sept. 2023), to measure WR140's dust kinematics and confirm its morphology. To measure the proper motions and projected velocities of the dust shells, we performed a novel PSF subtraction to reduce the effects of the bright diffraction spikes and carefully aligned the Cycle 2 to the Cycle 1 images. At 7.7 $\mu$m, through the bright feature common to 16 dust shells (C1), we find an average dust shell proper motion of $390\pm29$ mas yr$^{-1}$, which equates to a projected velocity of $2714\pm188$ km s$^{-1}$ at a distance of 1.64 kpc. Our measured speeds are constant across all visible shells and consistent with previously reported dust expansion velocities. Our observations not only prove that these dusty shells are astrophysical (i.e., not associated with any PSF artifact) and originate from WR140, but also confirm the "clumpy" morphology of the dust shells, in which identifiable substructures within certain shells persist for at least 14 months from one cycle to the next. These results support the hypothesis that clumping in the wind collision region is required for dust production in WR binaries.

Figures

Figures reproduced from arXiv: 2502.02738 by the authors.

Figure 1
Figure 1. False-color red, green, blue and orange–blue overlay. Point-spread function (PSF) subtracted Cycle 1 image (left), Cycle 2 image (middle), and comparison image of Cycle 1 (orange) to Cycle 2 (blue) at 7.7 μm (Section 3.2). The alignment of the two cycles to each other is highlighted in the background stars of the right panel (cycle overlay); their white color indicates that the orange Cycle 1 image is well aligned w… view at source ↗
Figure 2
Figure 2. PSF subtraction. Before (top) and after (bottom) PSF subtraction of Cycle 2 images at 7.7, 15, and 21 μm (Section 3.2). The varying flux of the diffraction spikes as a function of radius in the 7.7 μm image required a complex scaling of the PSF model, which resulted in some areas of oversubtraction. See the Appendix for the corresponding figure for the Cycle 1 images. 3 The Astrophysical Journal Letters, 979:L3 (10p… view at source ↗
Figure 3
Figure 3. Substructure persisting in WR 140. Left: 7.7 μm PSF-subtracted Cycle 2 image with red lines indicating the 10 radial flux profiles per dust structure (E, C0, and C1) that we used to calculate the proper motions (Section 4.1). The three colored boxes correspond to the zoomed-in regions shown on the right (yellow for E, blue for C1, and green for C0). Top right: zoomed-in 7.7 μm images for Cycle 1 (left) and Cycle 2 (… view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Radial flux profiles. Median of 10 radial flux profiles for Cycles 1 and 2 at 7.7, 15, and 21 μm through the major dust structures: C1, E, and C0 from left to right (Section 4). Gray vertical lines indicate intervals of 2.  67 (4380 au at our assumed distance of 1.64 …
Figure 5
Figure 5. Figure 5: shows that at each dust feature and for each wavelength, our calculated expansion velocities are broadly consistent across all shells. Outlying measurements exhibit large uncertainties, likely due to low signal to noise in the outermost shells as well as the effect of …
Figure 6
Figure 6. Figure 6: PSF subtraction. Before (top) and after (bottom) PSF subtraction of Cycle 1 images at 7.7, 15, and 21 μm. The varying flux of the diffraction spikes as a function of radius in the 7.7 μm image required a complex scaling of the PSF model which resulted in some areas of …
Figure 7
Figure 7. Figure 7: Animated videos. Top: video that blinks between the Cycle 1 and Cycle 2 images at 7.7 μm showing the entire system. Bottom: same as the top, but zoomed in on the C1 dust feature. All images have been PSF subtracted and the two cycles aligned. (An animation of this figu…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    WR 140 in the Infrared

    Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167 Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, AJ, 156, 58 Bradley, L., Sipo őcz, B., Robitaille, T., et al. 2024, astropy /photutils: v1.12.0, Zenodo, doi: 10.5281/zenodo.10967176 Dicken, D., Marn, M. G., Shivaei, I., et al. 2024, A&A,...

  2. [7]

    Top: video that blinks between the Cycle 1 and Cycle 2 images at 7.7 μm showing the entire system

    Animated videos. Top: video that blinks between the Cycle 1 and Cycle 2 images at 7.7 μm showing the entire system. Bottom: same as the top, but zoomed in on the C1 dust feature. All images have been PSF subtracted and the two cycles aligned. (An animation of this figure is available in the online article .) 10 The Astrophysical Journal Letters, 979:L3 (10...

Pith tools

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