{"id":"334dd4b6-9523-4d5c-80d7-5e018257961c","arxiv_id":"2507.01084","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Chandra observations 15 years apart yield a 24.7 ± 6.8 mas/yr proper motion for the neutron star in G18.9-1.1, corresponding to transverse velocities of about 264 or 474 km/s depending on distance, directed opposite the ejecta.","lead":"By comparing X-ray images taken 15 years apart, astronomers measured the neutron star in supernova remnant G18.9-1.1 moving across the sky. The speed is roughly 260 to 470 kilometers per second depending on the remnant's distance, and the motion points away from the bulk of the supernova debris, supporting the idea that neutron star kicks come from momentum conservation during the explosion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The four-parameter Gaia registration cannot exclude position-dependent astrometric distortion; the reported post-fit residual floor (Eq. 3) does not bound this systematic, so the 24.7 mas/yr proper motion may be biased at the ~1-sigma level.","rationale":"The reader's weakest-assumption analysis correctly identifies the four-parameter epoch tie as the load-bearing step. The paper otherwise contains the right ingredients: a 15-year baseline, PSF-aware centroiding via ChaRT and Sherpa, reference stars tied to Gaia DR3 with proper-motion propagation, and a reported post-fit residual floor near 0.015 arcsec. However, that residual floor is an output of the same four-parameter model and therefore cannot certify the absence of position-dependent errors. Because the two epochs use largely disjoint reference-source sets, a spatially structured distortion can masquerade as proper motion. This does not make the measurement wrong, but it means the central 24.7 +/- 6.8 mas/yr value should be regarded as provisional until the distortion check is performed. The conservation-of-momentum and distance conclusions downstream of the proper motion are secondary and already presented with caveats, so they do not change the conditional nature of the verdict. I therefore keep the reader's CONDITIONAL verdict unchanged, with the distortion test as the condition for lifting it.","tokens_in":12013,"tokens_out":7612,"duration_ms":180753,"concrete_test":"For each epoch, after applying the best-fit four-parameter correction, regress the RA and Dec residuals of the Gaia registration sources against off-axis angle (and, as a secondary check, against source brightness). Then refit Equation (1) with an additional radial distortion term Delta-r = a (r - r0)^2 and recompute the NS displacement. If the NS position shifts by more than ~0.05 arcsec (~3 mas/yr) or the residual-versus-off-axis correlation is significant at 2-sigma, the four-parameter registration is inadequate and the quoted proper motion uncertainty must be inflated accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on Equation (1), a four-parameter similarity transform (scale, rotation, translation) that registers each Chandra epoch to Gaia DR3. This model absorbs global shifts, rotations, and plate scale, but not position-dependent distortions or PSF-dependent centroid biases that differ between the 2009 (ObsID 10098) and 2024 (merged 26656+29478) observations. The two epochs use different reference-source sets (13 vs 9 Gaia matches, only 4 in common), so any distortion term not common to both configurations projects directly onto the NS displacement. The stated ~0.015 arcsec systematic floor in Section 2.4 is the weighted mean residual of the fit (Equation 3), not a bound on unmodeled higher-order terms; a four-parameter fit to 13/9 points can leave small residuals even when such terms are present. The measured NS displacement is ~0.37 arcsec over 15 years, so a residual position-dependent error of ~0.10 arcsec would shift the proper motion by ~6.7 mas/yr, comparable to the quoted 6.8 mas/yr uncertainty. The paper reports no residual-versus-off-axis-angle analysis, no leave-one-out test over reference sources, and no check for a radial distortion term, leaving this systematic unquantified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports the first direct measurement of the proper motion of the neutron star CXOU J182913.1-125113 in the supernova remnant G18.9-1.1, using Chandra observations separated by 15 years (2009 and 2024). The authors register each epoch's astrometry to the Gaia DR3 frame via a four-parameter similarity transform using faint X-ray sources matched to Gaia, accounting for the reference stars' proper motions and using PSF-fitted centroids. They measure a total proper motion of 24.7 ± 6.8 mas yr^{-1} at position angle 336° ± 16° east of north, which translates to Galactocentric transverse velocities of 264 km/s (at 2.1 kpc) or 474 km/s (at 3.8 kpc). They compare the NS motion to the SNR's bulk ejecta asymmetry, finding the NS moves roughly opposite the X-ray-emitting ejecta, and use a power-ratio multipole analysis to argue that the larger distance and velocity are slightly favored. The paper also estimates the progenitor mass as 15-20 M_sun and discusses implications for NS kick mechanisms.","tokens_in":12274,"tokens_out":4671,"duration_ms":58137,"significance":"If the astrometric measurement is robust, this is a valuable addition to the small sample of directly measured neutron-star proper motions in young, thermally dominated supernova remnants. The measurement would strengthen the observational case that NS kicks are anticorrelated with bulk ejecta motion, and the derived velocities would help constrain the distance and age of G18.9-1.1. The paper's methodology is careful in several respects: it uses Gaia DR3 reference sources with propagated proper motions, PSF-based centroiding, and explicit quadrature combination of centroid, registration, and transformation uncertainties. These strengths are real and deserve credit. However, the central claim depends on the validity of the four-parameter registration model, and the paper does not fully demonstrate that position-dependent astrometric distortions are negligible at the claimed precision.","major_comments":[{"comment":"The absolute astrometric correction uses a four-parameter similarity transform (scale, rotation, translation) fitted to only 13 Gaia-matched sources in the 2009 epoch and 9 in the 2024 epoch, with only 4 sources in common between the two epochs. The stated residual, σRes,tot in Eq. (3), is a weighted mean offset of the registration-source residuals, not a measure of unmodeled position-dependent distortion. A residual radial or higher-order plate distortion of order 0.1 arcsec—which cannot be excluded from a four-parameter fit to sparse source samples—would shift the NS displacement (~0.37 arcsec over the 15-year baseline) by ~0.1 arcsec and bias the proper motion by ~6.7 mas/yr, comparable to the quoted 1σ uncertainty of 6.8 mas/yr. Please add a residual-versus-off-axis-angle analysis, a leave-one-out test over the registration sources, and explicit tests for radial or quadratic distortion terms. If such terms are present at the claimed systematic floor, a systematic error term must be added; if they are not, the paper should demonstrate this quantitatively.","section":"§2.3–2.4, Eqs. (1)–(3)"},{"comment":"The power-ratio distance preference is not a strong constraint and is presented in a slightly misleading way. The multipole moments are centered on the back-evolved NS birth site, whose angular distance from the current NS position depends on the assumed distance and age (and on the measured proper motion). The sample of comparison remnants is small (six objects), and the text admits the conclusion is 'tenuous.' In the abstract, however, the power-ratio analysis is stated as favoring the higher velocity without noting how weak this preference is. Please provide a quantitative measure of the significance of the PRM difference (e.g., the scatter among comparison remnants and the shift in PRM when moving the center within the birth-site uncertainty) or explicitly relegate this to an inconclusive consistency check.","section":"§3 and §3.2"},{"comment":"There is an inconsistency in the velocity numbers as printed: the text first quotes raw transverse velocities of 246d2.1 ± 78 km/s and 445d3.8 ± 129 km/s, then after the Galactic-rotation correction quotes 264d2.1 ± 78 km/s and 474d3.8 ± 129 km/s. The abstract quotes 264d2.1 ± 79 km/s and 474d3.8 ± 129 km/s. Please label the raw and rest-frame values explicitly and harmonize the numerical uncertainties (78 vs. 79 km/s) so a reader can follow which quantity is being cited.","section":"§3, first paragraph"}],"minor_comments":[{"comment":"The typesetting of the remnant name is inconsistent (e.g., 'G18.9 −1.1' vs. 'G18.9-1.1'); please standardize the hyphenation and spacing.","section":"Throughout"},{"comment":"Equations (2) and (3) are rendered with broken formatting in the manuscript (the fractions appear as '1PN' and 'PN'), which makes them difficult to read. Please ensure the inverse-variance weighted sums are typeset correctly.","section":"§2.4, Eqs. (2)–(3)"},{"comment":"The conversion from 2D to 3D velocity using a multiplicative factor of sqrt(3/2) assumes that the unseen radial velocity component has the same magnitude as each tangential component. This assumption should be stated explicitly when the 3D velocity estimates of ~323 and ~580 km/s are introduced.","section":"§3.3"},{"comment":"There is a punctuation error at 'Holland-Ashford et al. 2024. Specifically, young SNRs...'; a period and comma should be separated properly, and the sentence continues correctly.","section":"§1"},{"comment":"Several unicode/TeX artifacts appear in the reference list (e.g., 'T¨ ullmann' with a combining diaeresis). Please check all author names for correct rendering of diacritics.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The central astrometric measurement appears to be carefully done and is likely correct, but the manuscript currently does not rule out a position-dependent astrometric distortion at the level that would undermine the claimed 1σ precision. Because the proper-motion value is the main result, the added systematic tests should be required before publication. The power-ratio distance preference is overstated in the abstract and should be toned down. I recommend major revision rather than rejection, because the concerns are addressable with additional analysis, not fundamental to the method."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a careful first direct proper motion measurement for the NS in G18.9-1.1, and I think the central number, 24.7 ± 6.8 mas/yr, is probably right. The two auxiliary claims—the slight preference for the 3.8 kpc distance and the anti-aligned kick—are weaker than the abstract suggests, and the authors mostly admit as much.\n\nThe astrometry is done well: PSF-based centroiding, Gaia DR3 reference frame with proper motions propagated, quadrature combination of centroid, registration, and transformation uncertainties. The method is standard, but it's applied cleanly. The paper is also honest that the power-ratio distance conclusion is tenuous.\n\nThe main soft spot is the one the stress-test flags. The four-parameter transform (scale, rotation, translation) doesn't model position-dependent distortion. The two epochs use different reference-source sets (13 vs 9, with only 4 in common), so any distortion term that isn't common to both configurations projects directly onto the NS displacement. The ~0.015 arcsec residual floor is the mean residual of the fit, not a bound on higher-order terms. The NS moved ~0.37 arcsec over 15 years; a 0.1 arcsec position-dependent error would shift the proper motion by ~6.7 mas/yr, comparable to the quoted 6.8 mas/yr uncertainty. So the systematic could be at the 1-sigma level. The paper doesn't test this with a leave-one-out analysis or residual vs off-axis angle plot. That's a legitimate gap, but not fatal: even a 1-sigma systematic leaves the detection at roughly 2.5–3 sigma.\n\nThe power-ratio distance preference has a mild circularity: the PRM origin is the back-evolved NS position, which depends on the assumed distance and age. The comparison to other SNRs is suggestive, not conclusive, and no PRM uncertainties are reported. The authors call it tenuous, which is fair.\n\nThe 'directly opposite the bulk ejecta' statement is based on morphology (X-ray center vs NS position), not a measured ejecta velocity vector. It's consistent with the PWN elongation, but it's not a new kinematic measurement of the ejecta.\n\nWho is this for? People working on NS kicks and SNR-embedded NS velocities. It adds one more object to a small sample. It deserves peer review—the astrometry should be scrutinized, and the PRM usage could be tightened. I'd engage with it.","headline":"Careful first proper motion measurement for this neutron star, probably right but with an unquantified systematic that could bite at the 1-sigma level; the distance and kick claims are softer than the abstract suggests.","tokens_in":12804,"tokens_out":3306,"would_cite":true,"duration_ms":35950,"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":"The first direct proper-motion measurement of the neutron star in SNR G18.9-1.1 finds it moving northwest at 24.7 mas/yr, opposite the remnant's bulk ejecta.","keywords":["neutron star proper motion","supernova remnant","G18.9-1.1","Chandra astrometry","Gaia reference frame","neutron star kick","ejecta asymmetry","power ratio method"],"falsifier":"Re-observing the field with a third Chandra epoch in the 2030s and checking whether the neutron star lands on the $24.7$ mas yr$^{-1}$ track extrapolated from 2009–2024, within the quoted uncertainties, would settle whether the two-epoch transformation absorbed a systematic astrometric error.","tokens_in":11808,"feed_emoji":"⭐","tokens_out":10315,"duration_ms":96831,"temperature":0.7,"pith_summary":"This paper claims to have measured, for the first time, the motion of the neutron star in the supernova remnant G18.9-1.1 across the plane of the sky. Using two Chandra observations separated by 15 years and tying both epochs to the Gaia reference frame, the authors find a total proper motion of $24.7 \\pm 6.8$ mas yr$^{-1}$ at $336^\\circ \\pm 16^\\circ$ east of north. For the two published distances of 2.1 and 3.8 kpc, this corresponds to Galactic-rotation-corrected transverse velocities of about 264 and 474 km s$^{-1}$. The motion points nearly opposite the bulk of the X-ray ejecta, which the authors interpret as support for a momentum-conservation kick mechanism, and a multipole analysis of the ejecta slightly favors the larger distance and higher velocity.","feed_headline":"Neutron star in G18.9-1.1 flies at 264–474 km/s","feed_subtitle":"A 15-year Chandra baseline places its motion opposite the ejecta, supporting momentum-conservation kicks.","key_machinery":"The load-bearing object is a four-parameter astrometric alignment: a scale factor $r$, a rotation $\\theta$, and translations $\\Delta x,\\Delta y$ that map each Chandra image onto the absolute Gaia frame, solved by least squares using point sources matched to Gaia DR3 and weighted by their centroid uncertainties. Source positions themselves come from PSF-aware centroid fits, so the neutron star's position at each epoch is measured relative to the same external frame rather than relative to the other epoch. The other main element is the power ratio method, a multipole decomposition of the 0.5–2.1 keV ROSAT image whose normalized moments ($P_1/P_0$, $P_2/P_0$, $P_3/P_0$) quantify dipole, quadrupole, and octupole asymmetry of the ejecta around the back-evolved birth site. Uncertainties are formed by adding the neutron star centroid error, the inverse-variance-weighted registration source errors, and the weighted residual scatter of the transformation in quadrature.","core_discovery":"The central claim is that the neutron star CXOU J182913.1-125113 in G18.9-1.1 is moving at $24.7 \\pm 6.8$ mas yr$^{-1}$ toward position angle $336^\\circ$, measured from a 2009 Chandra epoch and a merged 2024 epoch after correcting both to Gaia DR3 stellar positions. Back-evolving this motion for the age implied by each distance places the neutron star's birth site well away from the remnant's geometric center and its X-ray brightness center. The paper also claims that the kick direction is nearly opposite the bulk ejecta motion, consistent with the gravitational tugboat picture in which the neutron star is accelerated by asymmetric ejecta, and that the power ratios of the ejecta asymmetry are closer to those of other core-collapse remnants when the 3.8 kpc distance is adopted, mildly favoring the $\\sim474$ km s$^{-1}$ velocity.","pith_inferences":["This measurement opens a cheap archival route: any young SNR with two Chandra epochs separated by more than a decade and a handful of Gaia-registered sources can now yield a NS proper motion without new observing time.","The 3.8 kpc preference rests on a six-object comparison and is not decisive; an independent distance (for example from HI absorption or a future astrometric counterpart) would sharpen both the velocity and the SNR age.","A direct test of the method's systematic floor would be to register a field containing an extragalactic point source at both epochs and confirm that its apparent proper motion is consistent with zero."],"forward_implications":["The neutron star's birth site lies several arcminutes from the SNR's geometric and X-ray centers, so velocity estimates based on the geometric center (700–960 km/s) are rejected by this measurement.","The nearly opposite alignment of NS motion and bulk ejecta motion adds a new object to the small sample supporting momentum-conservation (gravitational tugboat) kick mechanisms.","If the 3.8 kpc distance is correct, the SNR is about 10 kyr old and its thermal plasma is close to collisional ionization equilibrium, matching the lower-temperature spectral fits preferred by earlier analysis.","The implied 3D velocities (about 323 and 580 km/s) match 3D core-collapse simulations with 15–20 solar-mass progenitors.","Longer X-ray baselines (20–40 years) will make NS proper-motion measurements routine and improve astrometric registration."],"supporting_citations":[{"why":"Defines the absolute reference frame (Gaia DR3 positions and proper motions) to which both Chandra epochs are corrected.","marker":"Gaia Collaboration et al. 2023"},{"why":"Supplies the astrometric-correction procedure, including the soft-l1 least-squares fit used for the transformation.","marker":"Long et al. 2022"},{"why":"Provides the comparison sample of SNRs with NS velocities and the power-ratio analysis this paper extends.","marker":"Holland-Ashford et al. 2017"},{"why":"Introduces the power ratio method used to quantify ejecta asymmetry.","marker":"Lopez et al. 2009a"},{"why":"Models the gravitational tugboat mechanism that the opposite-direction kick is taken to support.","marker":"Wongwathanarat et al. 2013"},{"why":"Identified the neutron star and reported the elongated pulsar wind nebula and the earlier geometric velocity estimate.","marker":"Tüllmann et al. 2010"},{"why":"Supplies the SNR age and spectral modeling that set the back-evolution timescale for each distance.","marker":"Harrus et al. 2004"},{"why":"One of the two distance estimates (2.1 kpc) used to convert proper motion to velocity.","marker":"Ranasinghe et al. 2019"},{"why":"The other distance estimate (3.8 kpc) used to convert proper motion to velocity.","marker":"Zhou et al. 2023"}],"fun_headline_variants":["G18.9-1.1's neutron star: measured 264–474 km/s kick","Kick velocity 264–474 km/s for NS in G18.9-1.1","Chandra data pin down NS kick in G18.9-1.1","Neutron star in G18.9-1.1 moves opposite its ejecta","15-year Chandra baseline reveals NS kick in G18.9-1.1"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes that a single overall shift, rotation, and stretch of each Chandra image fully removes the difference between the telescope's astrometric solution and the Gaia frame, so that any leftover position-dependent distortion or centroid bias cannot masquerade as neutron-star motion.","fun_headline_variants_meta":{"raw":{"variants":["G18.9-1.1's neutron star: measured 264–474 km/s kick","Kick velocity 264–474 km/s for NS in G18.9-1.1","Chandra data pin down NS kick in G18.9-1.1","Neutron star in G18.9-1.1 moves opposite its ejecta","15-year Chandra baseline reveals NS kick in G18.9-1.1"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001482,"raw_usage":{"total_tokens":5985,"prompt_tokens":1005,"completion_tokens":4980,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":621,"completion_tokens_details":{"reasoning_tokens":4864}},"tokens_in":621,"tokens_out":4980,"duration_ms":39652,"temperature":1.0,"reasoning_tokens":4864,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:00:50.777864+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-observing the field with a third Chandra epoch in the 2030s and checking whether the neutron star lands on the $24.7$ mas yr$^{-1}$ track extrapolated from 2009–2024, within the quoted uncertainties, would settle whether the two-epoch transformation absorbed a systematic astrometric error.","supporting_citations":[{"cited_title":"J., Plucinsky, P","cited_arxiv_id":null,"evidence_quote":"Supplies the astrometric-correction procedure, including the soft-l1 least-squares fit used for the transformation."},{"cited_title":"A., Auchettl, K., Temim, T., & Ramirez-Ruiz, E","cited_arxiv_id":null,"evidence_quote":"Provides the comparison sample of SNRs with NS velocities and the power-ratio analysis this paper extends."},{"cited_title":"2013, A&A, 552, A126","cited_arxiv_id":null,"evidence_quote":"Models the gravitational tugboat mechanism that the opposite-direction kick is taken to support."},{"cited_title":"M., Slane, P","cited_arxiv_id":null,"evidence_quote":"Supplies the SNR age and spectral modeling that set the back-evolution timescale for each distance."},{"cited_title":"Distance and Evolutionary States of Supernova Remnant G18.9-1.1 and Candidate G28.6+0.0","cited_arxiv_id":"1910.05407","evidence_quote":"One of the two distance estimates (2.1 kpc) used to convert proper motion to velocity."}],"review_version":1}