{"id":"9216f6f6-fdda-4d40-8a00-c7a435cabfe4","arxiv_id":"2412.04096","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using Doppler maps and proper motions, the authors build 3D velocity fields showing 20-30% front-to-back speed differences and distinct Si versus S spatial distributions in Tycho's remnant.","lead":"This paper reconstructs three-dimensional velocity fields for silicon- and sulfur-rich ejecta in Tycho's supernova remnant using Chandra X-ray line shifts and VLA radio proper motions. The result constrains how asymmetric the explosion and the progenitor's internal structure were, and reveals that sulfur moves preferentially away from the observer while silicon is more isotropic.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (4) equates the observed Doppler velocity to an unnormalized sum of line-of-sight velocities, making the asymmetry parameter alpha grid-dependent and the 20–30% result an artifact unless S is properly normalized.","rationale":"The reader's verdict is CONDITIONAL based on the model dependence of the two-zone alpha scaling. I agree that the two-zone multiplicative model is a strong assumption, but I find a more fundamental, internal inconsistency in the same equations. The observed Doppler velocity is an average along the line of sight, while Eq. (4) defines it as a sum. This dimensional mismatch makes the derived alpha sensitive to the arbitrary number of cells used in the numerical cube. Because S enters the denominator of Eq. (6), increasing resolution reduces the inferred asymmetry; the 20–30% value therefore cannot be separated from the chosen discretization. The paper even notes (Sect. 4.2) that a is taken equal to the average projected radius, which affects S, but the sum/mean issue is never addressed. A simple resolution test would settle this. If the result is resolution-independent, the concern would be resolved; but as written, the central quantitative claim is not supported. The qualitative differences between Si and S Doppler maps (Sect. 3, Figs. 5–6) remain valuable, so the paper has merit, but the 3D reconstruction should not be accepted as a measurement until Eqs. (4)–(6) are corrected and validated.","tokens_in":9153,"tokens_out":7617,"duration_ms":78127,"concrete_test":"Repeat the reconstruction of Sect. 4.2 on the same Chandra/VLA data but with two different cell sizes along the line of sight, e.g., 20 and 60 equal steps instead of the original grid. If the derived alpha map and the front/back speed ratio (the 20–30% figure) change by more than ~5% between resolutions, the unnormalized sum in Eq. (4) is the cause and the 3D velocity field is resolution-dependent.","verdict_should_be":"REJECT","load_bearing_attack":"In Sect. 4.2, the reconstruction hinges on Eqs. (4)–(6). The observed Doppler velocity VD in each spatial cell is the centroid shift of the Si/S line, i.e., a flux-weighted average of vx along the line of sight. The paper instead writes VD = Σ_{x>0} v'_x − Σ_{x<0} |v'_x|, with S = Σ vx over one half-space. Summing N resolution elements changes the units: S is proportional to N times the average half-space velocity, not to the observed velocity. As a result, alpha from Eq. (6) satisfies alpha ≈ 1 + VD/(2S), and the inferred front/back asymmetry (alpha − 1/alpha) scales inversely with N. Doubling the grid resolution roughly halves the asymmetry. No normalization by path length or cell count is provided. The choice a = average projected radius (Sect. 4.2) also enters R'(x, phi) and hence S, so the 20–30% value mixes an arbitrary discretization with a geometric assumption. Unless Eqs. (4)–(6) are re-derived with S as the mean (or the integral properly normalized), the headline 3D asymmetry is not a measurement.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":9380,"tokens_out":5400,"duration_ms":56384,"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":[{"comment":"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.","section":"Sect. 4.2, Eqs. (4)-(6)"},{"comment":"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.","section":"Sect. 4.2, Eq. (6), Figs. 9-10"},{"comment":"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.","section":"Sect. 4.2, Eq. (2)"}],"minor_comments":[{"comment":"The abstract, title, and text repeatedly use 'Si-reach' and 'S-reach'; these should be 'Si-rich' and 'S-rich'.","section":"Abstract and title"},{"comment":"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.","section":"Sect. 3"},{"comment":"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.","section":"Appendix 7.1"},{"comment":"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.","section":"Sect. 3, Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The normalization problem in Eq. (4) is serious and directly affects the headline result; I would not recommend acceptance until the authors recompute with a properly normalized S, demonstrate grid independence, and add sensitivity tests for the assumed ellipsoid depth. The observational Doppler maps and the cross-checks with earlier work are valuable and should be preserved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Paul,\n\nThe headline: this is the first paper to reconstruct separate 3D velocity cubes for Si and S ejecta in Tycho from Doppler maps plus proper motions, and the Doppler maps themselves look solid. The quantitative claim that front and back speeds differ by 20–30% is not yet supported: the reconstruction in Sect. 4.2 is built on an unnormalized sum that makes the asymmetry grid-dependent, plus a two-zone line-of-sight model and an assumed ellipsoidal depth.\n\nWhat is genuinely new: combining the existing Doppler maps [7,8,10] with proper-motion measurements into separate 3D fields for Si and S, and the finding that S is more redshifted and less isotropic than Si. That difference is directly visible in their Doppler maps and Table 2, so it does not depend on the 3D model. The rest-energy analysis in the appendix is careful, and the cross-checks against earlier independent maps are the right thing to do.\n\nThe soft spots are real. The stress-test note about Eq. (4) lands. The observed Doppler velocity in a cell is a flux-weighted centroid of vx along the line of sight, not a sum of per-cell vx values. Eq. (4) writes VD as (alpha - 1/alpha) S, with S a sum over N resolution elements. If you refine the grid, S roughly doubles, so alpha from Eq. (6) and the resulting 20–30% asymmetry change. No normalization by path length or cell count appears anywhere. The choice a = average projected radius also feeds into R'(x, phi) and hence S, so the asymmetry number mixes geometry with discretization. The paper's own admission that reconstructed speeds exceed the post-shock limit (3500 vs 2900 km/s) is honest, but it reinforces that the 3D field is a model output, not a measurement.\n\nThe two-zone multiplicative alpha is also arbitrary: alpha is solved per cell to reproduce the observed Doppler velocity, so the Doppler map is incorporated by construction rather than predicted. The qualitative conclusion that the asymmetry is moderate may survive, but the quoted percentage needs a proper derivation and error bars.\n\nWho benefits: SNR kinematicists and anyone testing Type Ia explosion asymmetry against Tycho. The Doppler maps deserve publication; the 3D reconstruction is a promising but unproven method.\n\nRecommendation: send to peer review, not desk reject. A referee should push for reworking Eqs. (4)–(6) with correct normalization and for uncertainty propagation on alpha. If the authors fix that, the 20–30% result could become a real measurement.","headline":"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.","tokens_in":9935,"tokens_out":3333,"would_cite":true,"duration_ms":33877,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["supernova remnants","Tycho's supernova","ejecta kinematics","Doppler shifts","proper motion","silicon-rich ejecta","sulfur-rich ejecta","Type Ia asymmetry"],"falsifier":"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.","tokens_in":8891,"feed_emoji":"🔭","tokens_out":7856,"duration_ms":71270,"temperature":0.7,"pith_summary":"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.","feed_headline":"Tycho's silicon and sulfur flows differ 20–30% front-to-back","feed_subtitle":"Doppler shifts plus proper motions map the 3D silicon and sulfur flows, testing how asymmetric the Type Ia blast was","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Provides the azimuth-dependent shock speeds and the 2.3 kpc distance used for the plane-of-sky expansion.","marker":"[13]"},{"why":"Supplies the two-epoch 1.4 GHz radio maps whose edge contours define the proper motion of the remnant's boundary.","marker":"[15]"},{"why":"Gives the Sedov solution with $v \\propto r$ used for the interior radial velocity profile.","marker":"[18]"},{"why":"Gives the ejecta-dominated homologous expansion $v \\propto r$ used for the same interior profile.","marker":"[19]"},{"why":"Atom database used to fix the rest energies of the Si XIII and S XV lines.","marker":"[20]"},{"why":"Prior Doppler map of the Si line used to check the new map and to estimate the uncertainty in the Doppler velocity.","marker":"[8]"},{"why":"Prior spatially resolved spectral analysis used for comparison and for the dense-wall interpretation of the fastest vectors.","marker":"[10]"},{"why":"Prior Si line energy map whose agreement with the new map supports the reliability of the fitting method.","marker":"[7]"}],"fun_headline_variants":["Tycho supernova: Si and S ejecta show 3D velocity asymmetry","Tycho's Si and S flows differ 20–30% front-to-back in 3D","3D map of Tycho's ejecta reveals S flowing outward, Si isotropic","Tycho's near and far ejecta speeds differ by 20–30%","Tycho's Si and S ejecta move at different 3D speeds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Tycho supernova: Si and S ejecta show 3D velocity asymmetry","Tycho's Si and S flows differ 20–30% front-to-back in 3D","3D map of Tycho's ejecta reveals S flowing outward, Si isotropic","Tycho's near and far ejecta speeds differ by 20–30%","Tycho's Si and S ejecta move at different 3D speeds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000833,"raw_usage":{"total_tokens":3626,"prompt_tokens":925,"completion_tokens":2701,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":541,"completion_tokens_details":{"reasoning_tokens":2587}},"tokens_in":541,"tokens_out":2701,"duration_ms":18778,"temperature":1.0,"reasoning_tokens":2587,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:45:20.966523+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the azimuth-dependent shock speeds and the 2.3 kpc distance used for the plane-of-sky expansion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-epoch 1.4 GHz radio maps whose edge contours define the proper motion of the remnant's boundary."},{"cited_title":"Taylor, Proc","cited_arxiv_id":null,"evidence_quote":"Gives the Sedov solution with $v \\propto r$ used for the interior radial velocity profile."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the ejecta-dominated homologous expansion $v \\propto r$ used for the same interior profile."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Atom database used to fix the rest energies of the Si XIII and S XV lines."},{"cited_title":"Godinaud, F","cited_arxiv_id":null,"evidence_quote":"Prior Doppler map of the Si line used to check the new map and to estimate the uncertainty in the Doppler velocity."},{"cited_title":"Uchida et al., ApJ 962, 159 (2024)","cited_arxiv_id":null,"evidence_quote":"Prior spatially resolved spectral analysis used for comparison and for the dense-wall interpretation of the fastest vectors."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior Si line energy map whose agreement with the new map supports the reliability of the fitting method."}],"review_version":1}