{"id":"25ffc0e9-c47a-450b-a132-d2c0843f9218","arxiv_id":"2502.02738","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A second JWST epoch of WR 140 shows its dust shells expanding at nearly constant speeds near the stellar wind speed, with clumpy substructures persisting for over a year.","lead":"Using two JWST observations taken 14 months apart, astronomers measured the outward motion of the dust shells around the binary star WR 140, finding steady expansion speeds of about 2,500 to 2,700 kilometers per second. The work confirms the shells are real structures from the star system and that lumpy substructures persist over a year, supporting the idea that clumping in colliding stellar winds is what makes dust form.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PSF-subtraction systematics at 7.7 μm, not propagated into the quoted uncertainties, are the load-bearing risk; the abstract/body number mismatch (390±29 vs 332±20 mas/yr) is a concrete symptom that the peak-fitting pipeline may be sensitive to processing choices.","rationale":"The reader's CONDITIONAL verdict is appropriate. The two-epoch JWST dataset, publicly available code, and consistency across three filters are real strengths, and the qualitative expansion and clump persistence are visually compelling. However the quantitative central claim, steady ~2586 km/s expansion consistent with the WC7 wind, depends on small radial shifts (~0.3–0.4 arcsec) measured against a bright, complex PSF. The PSF subtraction is the least controlled step: a single scale factor for a radially varying spike pattern, acknowledged oversubtraction, and no propagated systematic uncertainty. The abstract/body number mismatch and the high-velocity outer shells are warning signs that the analysis pipeline is not stable under reasonable choices. The proposed test, varying the PSF scale and re-measuring, would directly determine whether the quoted velocities are astrophysical or artifacts. If the test shows median shifts below ~20 mas/yr, the central claim holds as stated and the paper could be accepted; if not, the conclusion would need substantial revision. Thus no verdict change from the reader's CONDITIONAL is needed yet.","tokens_in":13274,"tokens_out":5755,"duration_ms":58490,"concrete_test":"Re-run the 7.7 μm C1 analysis with the adopted WebbPSF scale factor varied by ±10% and, separately, with a radially dependent scale factor matched at several radii. Re-fit the shell peaks and recompute the per-shell and median proper motions. If the median shifts by more than the quoted 20 mas/yr uncertainty, or if shells 13/15/16 change by more than 50 mas/yr, the PSF-subtraction choice is a dominant systematic. Also independently recompute Table 1's mean/median and the mas/yr-to-km/s conversion to resolve whether 390±29, 332±20, or neither is the correct headline value.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central velocities come from Gaussian fits to peaks in radial profiles after subtracting a WebbPSF model scaled by a single factor. At 7.7 μm the diffraction-spike flux varies with radius (Sect. 3.2, Fig. 2), so one scalar cannot remove a radius-dependent spike pattern; the authors acknowledge oversubtraction in places. If residual PSF structure is radius-dependent and the two epochs are not at the same telescope roll angle, the residual pattern will shift the fitted peaks by different amounts in Cycle 1 and Cycle 2, biasing every proper motion. This systematic is absent from the quoted ±15–20 mas/yr errors. A concrete symptom is the internal inconsistency of the headline numbers: the abstract reports 390±29 mas/yr and 2714±188 km/s for C1, while the body states a median 332±20 mas/yr and 2586±152 km/s; at 1.64 kpc, 390 mas/yr corresponds to ≈3030 km/s, and 2714 km/s corresponds to ≈349 mas/yr. The discrepancy implies the result is sensitive to whether the mean or median is used, which shells are included, or the PSF-scaling choice, and that the systematic error may exceed the statistical error. The three outer 7.7 μm shells (Table 1, shells 13/15/16) have velocities of 3360–3770 km/s, well above the WC7 terminal speed; attributing these to background stars needs quantitative support rather than post-hoc exclusion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":13579,"tokens_out":10567,"duration_ms":98479,"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":[{"comment":"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.","section":"Abstract; Section 4.1"},{"comment":"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.","section":"Section 3.2, Figure 2"},{"comment":"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.","section":"Table 1, C1 rows 13/15/16"},{"comment":"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.","section":"Section 4.1, Equation (1), Table 1 (E feature)"}],"minor_comments":[{"comment":"Equation (1) is garbled in the printed text; rewrite it with clear parentheses and define all angles explicitly.","section":"Section 4.1, Equation (1)"},{"comment":"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.","section":"Figure 2 caption; Section 3.2"},{"comment":"The rows labeled 7a and 9a for the E feature are empty; add a footnote explaining that these shells were not measured.","section":"Table 1, E feature"},{"comment":"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.","section":"Table 1 note"},{"comment":"The horizontal axis is labeled only as 'shell'; label the shell numbers so that the points can be cross-referenced with Table 1.","section":"Figure 5"}],"recommendation":"major_revision","confidential_remarks":"This is a re-analysis of public JWST data with public code, so the internal inconsistencies I list can be checked and corrected by the authors without new observations. I recommend major revision rather than rejection, assuming the authors can supply the requested PSF-systematic tests and reconcile the headline numbers."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The first thing to know: this is a legitimate observational step forward—two epochs of JWST MIRI imaging 1.12 years apart, proper motions measured directly for 16 shells at three wavelengths, with code and data public. The central result, steady expansion at roughly 2500–2700 km/s with no measurable acceleration, holds up. The second thing to know: the abstract's headline numbers are internally inconsistent and don't match the body, and that mess is a symptom of an unquantified PSF-subtraction systematic the referee should pin down.\n\nWhat's genuinely new: Lau et al. 2022 had a single epoch and could not measure motion. This paper adds the second epoch and gives the first full-shell kinematic map of the WR 140 dust nebula. The clump-persistence result—identifiable substructures survive 14 months—is new and visually convincing, and it supports the idea that clumping in the wind collision region is required for dust production. The analysis is transparent: median of ten radial profiles, Gaussian fits, JHAT alignment to 0.1 pixels. The three-wavelength agreement (median projected velocities 2586, 2350, 2547 km/s at 7.7, 15, 21 μm) is real evidence that they are measuring dust motion rather than a single-wavelength artifact.\n\nThe soft spots, in order. First, the abstract quotes 390±29 mas/yr and 2714±188 km/s for C1, while the body reports a median 332±20 mas/yr and 2586±152 km/s at 7.7 μm. The abstract's two numbers don't even convert to each other at 1.64 kpc (390 mas/yr is about 3030 km/s; 2714 km/s is about 349 mas/yr). The 390 is likely the mean of all 16 shells, pulled up by three fast outer shells; the body's median is robust to those. But a headline that disagrees with the body and with itself means the result depends on aggregation choices, and that dependence is never discussed. Second, the 7.7 μm PSF subtraction uses a single scalar scale factor where the diffraction-spike flux varies with radius—the authors acknowledge oversubtraction in places. If residual spike structure is radius-dependent and differs between the two epochs, the fitted peaks shift; that systematic is absent from the quoted ±15–20 mas/yr errors. Third, shells 13/15/16 at 7.7 μm give 3360–3770 km/s, well above the WC7 wind speed of 2860 km/s. The paper blames background stars and low S/N, but the quoted errors on those shells are not larger than the inner ones, and no specific contaminants are identified. Those three need quantitative treatment.\n\nNone of this kills the central claim. The medians agree across three wavelengths, the velocities match prior estimates, and the clump persistence stands on the images themselves. The paper is for anyone working on WR binaries, episodic dust formation, or JWST MIRI PSF-subtraction methods—they get a new reference measurement and a reusable codebase. It deserves a serious referee. My recommendation: send it out; the referee's job is to fix the abstract, propagate the PSF systematic, and resolve the three fast shells.","headline":"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.","tokens_in":14272,"tokens_out":8250,"would_cite":true,"duration_ms":68156,"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":"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.","keywords":["Wolf-Rayet binaries","dust formation","colliding winds","WR 140","JWST MIRI","proper motions","dust shells","infrared imaging"],"falsifier":"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.","tokens_in":13061,"feed_emoji":"🔭","tokens_out":6680,"duration_ms":55115,"temperature":0.7,"pith_summary":"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.","feed_headline":"Dust shells around WR 140 expand at 2,600 km/s","feed_subtitle":"Two JWST images 14 months apart show the rings are real outflows moving at the star's wind speed, not artifacts.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Cycle 1 JWST imaging that first resolved 17 dust shells and provides the earlier epoch and shell numbering used here.","marker":"R. M. Lau et al. 2022"},{"why":"Geometric and orbital modeling that explains shell creation near periastron and supplies the wind-shock opening angle used to deproject the E feature.","marker":"Y. Han et al. 2022"},{"why":"Defines the C1, C0, and E dust features and reported earlier proper-motion differences that this paper compares with its JWST measurements.","marker":"P. M. Williams et al. 2009"},{"why":"Provides the WC7 terminal wind speed of 2860 km/s that the measured shell velocities are compared against.","marker":"P. R. J. Eenens & P. M. Williams 1994"},{"why":"Supplies the inclination angle used in the deprojection formula and an earlier dust stream velocity estimate.","marker":"R. Fahed et al. 2011"},{"why":"Provides the Gaia-based distance of 1.64 kpc used to convert proper motions to projected velocities.","marker":"C. A. L. Bailer-Jones et al. 2018"},{"why":"Describes the WebbPSF model used to subtract the diffraction spikes before measuring shell positions.","marker":"M. D. Perrin et al. 2012"},{"why":"Describes the JHAT alignment tool used to register Cycle 2 images to Cycle 1 to within 0.1 pixels.","marker":"A. Rest et al. 2023a, 2023b"}],"fun_headline_variants":["WR 140 dust shells expand at steady 2,700 km/s","JWST proves WR 140 dust rings are real outflows","Dust clumps around WR 140 persist across 14 months","WR 140's dust shells: constant speed, clumpy structure","Two JWST epochs confirm WR 140 shell motion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["WR 140 dust shells expand at steady 2,700 km/s","JWST proves WR 140 dust rings are real outflows","Dust clumps around WR 140 persist across 14 months","WR 140's dust shells: constant speed, clumpy structure","Two JWST epochs confirm WR 140 shell motion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000253,"raw_usage":{"total_tokens":1584,"prompt_tokens":986,"completion_tokens":598,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":510}},"tokens_in":602,"tokens_out":598,"duration_ms":5883,"temperature":1.0,"reasoning_tokens":510,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T11:18:14.123512+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}