{"id":"bcfc9c88-d32a-4e9d-ad2c-d309416c6f82","arxiv_id":"2601.14562","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Kinematics of 61,883 DESI stars indicate the Monoceros Ring is a Sagittarius-induced, corotating spiral arm with passage times 0.25 and 1.10 Gyr ago, and the Anticenter Stream is part of a broader vertical disk wave.","lead":"Using 61,883 DESI stellar spectra, this paper argues that the Monoceros Ring overdensity in the outer Milky Way is a spiral arm created by a past close encounter with the Sagittarius dwarf galaxy, and uses its star motions to estimate that Sagittarius last passed closest 0.25 and 1.10 billion years ago. It also finds the Anticenter Stream is not a separate stream but part of a larger vertical disk wave.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The FFT-based Sgr passage times are unsupported: the 0.25 Gyr peak has a period longer than the entire sampled L_Z^{-1} range, so it is unresolved and likely a window/trend artifact.","rationale":"The reader's CONDITIONAL verdict is appropriate: the kinematic maps and the qualitative association of the MRi with a tidally induced spiral arm are plausible and consistent with prior work, while the novel timing claim needs stronger support. My concern sharpens the reader's weakest assumption: it is not only that no null/mock test is shown, but that the lower FFT peak used for the 0.25 Gyr passage time is formally unresolved given the finite L_Z^{-1} span. This is a concrete, checkable issue rather than a generic appeal for more testing. The higher frequency peak (1.10 Gyr) is less obviously invalid, so I do not argue that the entire timing framework is false; rather, the onus is on the authors to demonstrate with synthetic recovery and detrended FFTs that their peaks correspond to physical winding frequencies. Since the reader already recommended CONDITIONAL, my read does not change the verdict, hence UNCHANGED. The condition should be sharpened to explicitly include resolving/detrending the low-frequency peak and a mock recovery test before the passage-time numbers are quoted as results.","tokens_in":33609,"tokens_out":6982,"duration_ms":72376,"concrete_test":"Reproduce the Section 3.4 FFT after subtracting a low-order polynomial (e.g., linear or quadratic) fit to the median V_R versus L_Z^{-1} before transforming, and repeat over a wider L_Z range (e.g., 1500–5500 km s^{-1} kpc) and over different l-slices (e.g., 170°<l<190°). If the f≈1313 peak disappears or shifts outside its quoted uncertainty, the 0.25 Gyr passage time is an artifact. Additionally, inject a synthetic tidally induced spiral-arm signal with known winding times into DESI-like mock data (using the Antoja et al. 2022 model or an N-body simulation) and run the full pipeline; require that both input passage times are recovered with the same binning, l-slice, and FFT procedure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative central claim — the two Sgr pericenter passage times — rests entirely on Section 3.4's FFT of median V_R versus L_Z^{-1} over L_Z ∈ [2000, 4500] km s^{-1} kpc. This span corresponds to a range in x = L_Z^{-1} of T = 1/2000 − 1/4500 ≈ 2.78×10^{-4} (km s^{-1} kpc)^{-1}. The two reported power-spectrum peaks are at f = 1313 and 5879 ((km s^{-1} kpc)^{-1})^{-1}. A peak at frequency f has period 1/f; for f=1313, the period is ≈ 7.62×10^{-4}, which is 2.7 times longer than the entire sampled span. The data therefore contain less than half a cycle of this putative oscillation. The Fourier resolution is ~1/T ≈ 3600, so the 1313 peak lies below the nominal resolution, and its quoted 1σ width (477) is smaller than the Rayleigh resolution — a signature of window/zero-padding artifacts or a smooth V_R gradient rather than a resolved winding frequency. This lowest peak is the one that yields the headline 0.25 Gyr passage time (Table 1). The higher peak (f≈5879, period ≈1.7×10^{-4}, ~1.6 cycles across the span) is at least in principle resolvable, but both peaks are quoted with formal widths that ignore the absence of a full cycle for the lower peak, the choice of l-slice (175°<l<185°), the completeness correction's restriction to RA/Dec bins, and the arbitrary impulse parameters (ΔV=10 km s^{-1}, D=20 kpc) in the Antoja et al. (2022) model. Without a synthetic recovery test or a null test on a smooth disk, the mapping from these FFT peaks to pericenter times is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses 61,883 main-sequence turn-off stars from DESI DR2 in the anticenter region (150° < l < 220°, 20° < b < 40°) to map the 6D kinematics of the Monoceros Ring (MRi) and Anticenter Stream (ACS). The authors report that the MRi overdensity has kinematics consistent with a tidally induced spiral arm driven by Sagittarius (Sgr), specifically a negative radial velocity region, a V_phi inflection line, and corotation at the overdensity. They then apply the Antoja et al. (2022) frequency-analysis method to the median V_R versus L_Z^{-1} curve and derive two recent Sgr pericenter passage times, 0.25 ± 0.09 Gyr and 1.10 ± 0.23 Gyr. For the ACS, they find positive V_R and V_Z kinematics that extend to lower Galactic latitudes, and they interpret the ACS as part of a broader vertical wave rather than a discrete kinematically distinct stream.","tokens_in":2346,"tokens_out":2675,"duration_ms":81695,"significance":"If the main claims hold, the paper would provide a direct observational identification of the MRi as a corotating, tidally induced spiral arm—a relatively rare classification in the Milky Way—and would demonstrate that outer-disk kinematic oscillations can be used to recover Sgr's recent pericenter timing. The DESI DR2 sample extends previous work to fainter magnitudes and larger distances, and the paper makes its data products publicly available, which is a strength. The MRi/ACS kinematic decoupling and the extension of ACS kinematics to lower latitudes are well supported by Figures 5, 12–14 and are consistent with APOGEE and Gaia-based studies. However, the quantitative passage-time claim rests on a Fourier analysis that, as presented, lacks validation against synthetic or null data and whose lowest-frequency peak is not resolved within the sampled baseline. The timing claim is therefore not yet on the same footing as the spatial/kinematic characterizations.","major_comments":[{"comment":"The lowest FFT peak is unresolved and is the basis for the headline 0.25 Gyr passage time. For L_Z ∈ [2000, 4500] km s^-1 kpc, the sampled range in x = L_Z^{-1} is T ≈ 2.78×10^-4 (km s^-1 kpc)^-1. The peak at f = 1313 ((km s^-1 kpc)^-1)^-1 has a period ≈ 7.62×10^-4, which is 2.7 times longer than the entire sampled span, so the data contain less than half a cycle of this putative oscillation. The nominal Fourier resolution ~1/T ≈ 3600 is larger than 1313, meaning this peak lies below the Rayleigh resolution; its quoted Gaussian width (477) is smaller than the resolution, consistent with a window/trend artifact rather than a resolved winding frequency. Because this peak directly yields the 0.25 ± 0.09 Gyr timing in Table 1, the central timing claim is not supported without a synthetic recovery test, a null test on a smooth disk model, or a demonstration that the peak is stable under chang","section":"§3.4, Fig. 10"},{"comment":"The quoted uncertainties on the passage times (0.09 and 0.23–0.28 Gyr) are only the formal widths from Gaussian fits to the FFT peaks after Monte-Carlo propagation of V_R errors. They do not include systematic errors from: the choice of the restricted l-slice 175° < l < 185°, the completeness correction being applied only in RA/Dec bins (not in distance or l,b), the number and width of L_Z bins, the subtraction of the background disk model, or the arbitrary impulsive-perturber parameters in the Antoja et al. (2022) model (ΔV = 10 km s^-1, D = 20 kpc). The mapping from FFT frequency to time via Eq. (1) assumes a single coherent, corotating tidally induced spiral arm winding in a flat rotation curve; no mock or N-body test is provided to show that a known input passage time is recovered within the claimed errors. The agreement with literature values in Table 1 is reassuring but is not a su","section":"§3.4, Eq. (1)"},{"comment":"The identification of the MRi as a tidally induced spiral arm rests on a visual / qualitative comparison: the overdensity, the −V_R region, and the fitted V_ϕ − ⟨V_ϕ(R)⟩ inflection line all appear spatially coincident. The inflection line itself is fit to the data in a way that is not fully specified (a cubic in each annulus, then a linear fit to the zero crossings), and no uncertainty is given for the magenta line. The paper does not quantitatively compare the observed V_R and V_ϕ patterns to the Antoja et al. (2022) models or the Stelea et al. (2024) simulations, nor does it test whether a simpler disk model (e.g., a warp, flare, or a smooth radial oscillation) could produce the same pattern. Since this classification is the foundation for the subsequent timing analysis, a quantitative model comparison or at least a null test on a smooth, unperturbed disk would materially strengthen th","section":"§3.2, Fig. 6"},{"comment":"The completeness correction is computed as the ratio of observed MAIN-BLUE targets to the full DESI target catalog in one-degree RA/Dec bins. This corrects for the survey footprint and tiling pattern but does not address completeness as a function of distance, magnitude, or Galactic latitude at fixed RA/Dec. If selection incompleteness correlates with kinematics—for example through distance-dependent MSTO sampling or dust-extinction-dependent proper-motion quality—the median V_R vs L_Z curve in Figure 9 could be biased. Given that the timing analysis uses the shape of this curve, the paper should quantify how robust the FFT peaks are to alternative completeness treatments (e.g., binning in l,b and distance, or weighting by the inverse completeness in each phase-space bin).","section":"§2.4, Fig. 4"}],"minor_comments":[{"comment":"The definition of L_Z is given as L_Z = R V_Z, which appears to be a typo: angular momentum about the Galactic center in cylindrical coordinates is L_Z = R V_ϕ (the azimuthal velocity). Please correct and ensure all subsequent uses of L_Z are consistent.","section":"§3.4"},{"comment":"The text states the median V_R is computed in 75 bins over L_Z ∈ [2000, 4500] km s^-1 kpc, while the Figure 9 caption says 70 bins over L_Z ∈ [1400, 45000] km s^-1 kpc. These numbers and ranges should be reconciled.","section":"§3.4 vs Fig. 9 caption"},{"comment":"The second passage time is reported as 1.10 ± 0.28 Gyr in the body text but as 1.10 ± 0.23 Gyr in the abstract and Table 1. Please make these consistent and specify which uncertainty (Gaussian width only, or including systematic terms) is being quoted.","section":"§3.4"},{"comment":"Equation (1) is reproduced from Antoja et al. (2022) but the typesetting is garbled and the units discussion is confusing ('time×length^-2' vs '(time/length)^2×time^-1'). Please provide a clean symbolic derivation or a clear statement of the units of each term so the conversion from frequency to Gyr is reproducible by a reader.","section":"Eq. (1)"},{"comment":"No significance threshold or false-alarm probability is assigned to the FFT peaks. The authors should state how 'significant' is defined (e.g., versus the noise floor or versus peaks in null simulations), or explicitly avoid the word 'significant'.","section":"Fig. 10"},{"comment":"Typo: 'We can assume the the completeness' should read 'We can assume that the completeness'.","section":"§2.4"},{"comment":"Several references appear duplicated with identical DOIs (e.g., Bernet et al. 2022 listed twice; Antoja et al. 2018a and 2018b share the same DOI/page numbers but are cited as separate works). Please check whether these are genuinely distinct papers or should be merged/corrected.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The core spatial/kinematic findings on the MRi and ACS are solid and well matched to prior work, and the paper has the potential to be a valuable contribution if the Sgr timing analysis is brought up to the same standard. However, the timing claim is the headline quantitative result and is currently supported only by an unresolved FFT peak and a model-mapping step that has not been validated against synthetic data. I would urge the editor to require either (a) a synthetic recovery test showing that the method recovers known input passage times from mock DESI-like catalogs, including a null test on a smooth disk, or (b) a substantial re-framing of the passage times as 'kinematically consistent with earlier pericenter estimates' rather than independently measured values. Without that change, the abstract and Section 3.4 overstate the certainty of the result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful part of this paper is the DESI DR2 kinematic mapping of the anticenter. The MRi/ACS decoupling (opposite V_R and V_Z) is clearly shown and consistent with Qiao et al. (2024), and the case that the MRi's -V_R region coincides with a corotating inflection in V_phi, following Antoja et al. (2022), is reasonable. Extending the ACS signature to lower latitudes as part of a broader vertical disturbance is a nice observation. They also post the data points on Zenodo and validate distances against SpecDis. This half deserves serious attention.\n\nThe Sgr timing claim is the weak half. Section 3.4 takes the FFT of median V_R vs L_Z^-1 over L_Z in [2000, 4500]. The peak that gives the headline 0.25 Gyr passage is at f ~ 1313, whose period is ~7.6e-4 - about 2.7 times the whole sampled span in L_Z^-1. You get less than half a cycle. That peak is not a resolved oscillation; it is likely a trend or window artifact. The quoted 1-sigma width (477) is smaller than the Rayleigh resolution (~3600), a red flag. The higher peak (f ~ 5879) is at least in principle resolvable, but without a null test on a smooth disk or a mock recovery test, neither peak is established as the winding signal. The mapping from FFT peaks to passage times via A22 Eq. 11 also carries systematic uncertainties (l-slice choice, per-RA/Dec completeness correction, A22's arbitrary impulse parameters) that are not in the quoted errors. The authors do test rotation-curve slope, which is good, but that is not the main fragility.\n\nThey are transparent about adopting A22's framework, and the times agree with literature - but that agreement is weak evidence when the model does the heavy lifting.\n\nBottom line: the kinematic maps and the MRi-as-corotating-arm interpretation are plausible and worth refereeing. The Sgr passage times should not be published as-is; they need a significance test, a synthetic recovery test, and an explicit treatment of the unresolved low-frequency peak. I'd send it to review with a request to fix or heavily caveat Section 3.4.","headline":"Solid observational kinematics for the MRi/ACS pair, but the FFT-based Sgr passage times rest on an unresolved peak and need synthetic-recovery work before they can be believed.","tokens_in":34845,"tokens_out":3756,"would_cite":false,"duration_ms":37391,"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 Monoceros Ring is a tidally induced spiral arm, and its winding pattern records the Sagittarius dwarf's last two disk passages.","keywords":["Monoceros Ring","Anticenter Stream","tidally induced spiral arms","Sagittarius dwarf galaxy","Galactic anticenter","disk kinematics","main-sequence turn-off stars","Milky Way outer disk"],"falsifier":"Construct a DESI-like mock catalog from a smooth, axisymmetric disk model (no satellite perturbation) with the same selection function and completeness; if the identical Fourier pipeline yields comparable power-spectrum peaks, the V_R oscillation is not uniquely attributable to a tidally induced spiral arm. Alternatively, if the two dominant peaks move significantly when the longitude window is expanded from 10° to 20° or shifted by 5°, the coherent single-arm assumption fails.","tokens_in":33590,"feed_emoji":"🌌","tokens_out":8808,"duration_ms":87224,"temperature":0.7,"pith_summary":"Analyzing full three-dimensional velocities of 61,883 main-sequence turn-off stars in the Galactic anticenter, this paper argues that the Monoceros Ring — a giant stellar overdensity in the outer disk — is not a separate accreted cloud of stars but a corotating spiral arm excited when the Sagittarius dwarf spheroidal galaxy last swept through the Milky Way. Because the arm corotates with the disk, its wind-up rate is a clock: the oscillation of radial velocity with angular momentum yields two recent Sagittarius pericenter passages, 0.25 ± 0.09 and 1.10 ± 0.23 billion years ago. The same data show the Anticenter Stream moves opposite the Ring in both radial and vertical velocity, and those stream-like motions continue below the stream's photometric boundary, suggesting the stream is part of a broader vertical wave in the outer disk rather than a discrete object. If the interpretation holds, the paper converts a known overdensity into a chronometer for the Milky Way's recent interactions with its satellites.","feed_headline":"Monoceros Ring is a tidal spiral arm carved by the Sagittarius dwarf","feed_subtitle":"The arm's wind-up rate dates the dwarf's last two disk passages: 0.25 and 1.10 billion years ago.","key_machinery":"The key instrument is the kinematic signature of a transient tidally induced spiral arm: a spatial coincidence between the minimum of radial velocity, the density peak, and the inflection point of Vφ − ⟨Vφ(R)⟩ (the corotation condition). To turn this into a clock, the paper computes the median V_R in bins of inverse angular momentum L_Z^{-1} (assuming a flat rotation curve) over one 10-degree longitude slice, takes a fast Fourier transform, and maps the two strongest frequency peaks to pericenter times using a published relation (Equation 1). Everything about the timing depends on that mapping's assumption that the observed oscillation is the coherent winding of a single arm.","core_discovery":"The central claim is that the Monoceros Ring (MRi) — a stellar overdensity at Galactocentric radii 14–18 kpc in the anticenter — shows the diagnostic kinematics of a tidally induced spiral arm: the most negative radial velocity V_R coincides with the arm's density peak, and the azimuthal velocity difference Vφ − ⟨Vφ(R)⟩ passes through zero there, meaning the arm corotates with the disk. This is the pattern predicted when a satellite galaxy delivers an impulsive gravitational kick to disk stars, causing their orbits to crowd into a winding arm. Because the arm winds at the disk's circular speed, the oscillation frequency of V_R as a function of inverse angular momentum gives the time since ea","pith_inferences":["If the identification is right, the paper implies that the outer disk's spiral pattern is predominantly transient and satellite-driven, rather than a long-lived density wave; a testable extension would be predicting a matching gas response at a slightly different arm phase.","The timing analysis could be stress-tested by repeating the Fourier fit in adjacent longitude windows (e.g., 165°–175° and 185°–195°); if the recovered peak frequencies do not stay within uncertainties, the single-slice coherent-winding assumption would be in question.","The vertical-wave interpretation of the ACS suggests the solar-neighborhood phase spiral and the anticenter vertical folds may share a common perturbing event; comparing their inferred perturbation ages (≈0.25 Gyr) could unify both phenomena.","A natural validation is to run the same pipeline on simulated DESI-like catalogs from an axisymmetric disk with no satellite; if power-spectrum peaks persist, the passage-time extraction would need an independent background subtraction."],"forward_implications":["The Monoceros Ring being a tidally induced spiral arm provides a direct observational link between the Sagittarius dwarf's orbit and the spiral structure of the outer Milky Way disk.","The two pericenter times (0.25 and 1.10 Gyr) add a kinematic constraint on Sgr's recent orbital history that is independent of stream modeling and star-formation fits.","Because the arm corotates, the same V_R wind-up method can in principle clock other satellite encounters, such as the LMC's first infall.","The Anticenter Stream's kinematic extension below its photometric boundary shifts its interpretation from a stripped stream to a vertical phase-space fold, meaning future studies should treat it as part of the disk's disequilibrium rather than an accreted remnant."],"fun_headline_variants":["Monoceros Ring is a tidal spiral arm from Sagittarius dwarf","Sagittarius dwarf carved a spiral arm in the Milky Way's disk","Tidal spiral arm reveals two recent Sagittarius passages","Monoceros Ring kinematics: a spiraling wake from Sgr impact","Sagittarius dwarf's last two passages dated from spiral arm"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The timing chain rests on treating the wiggles in the median radial velocity across one narrow 10-degree slice of sky as the clean imprint of a single corotating spiral arm winding up, so that the Fourier peak frequencies can be converted directly into Sagittarius passage times.","fun_headline_variants_meta":{"raw":{"variants":["Monoceros Ring is a tidal spiral arm from Sagittarius dwarf","Sagittarius dwarf carved a spiral arm in the Milky Way's disk","Tidal spiral arm reveals two recent Sagittarius passages","Monoceros Ring kinematics: a spiraling wake from Sgr impact","Sagittarius dwarf's last two passages dated from spiral arm"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1328,"prompt_tokens":818,"completion_tokens":510,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":417}},"tokens_in":562,"tokens_out":510,"duration_ms":6362,"temperature":1.0,"reasoning_tokens":417,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T09:08:55.257386+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct a DESI-like mock catalog from a smooth, axisymmetric disk model (no satellite perturbation) with the same selection function and completeness; if the identical Fourier pipeline yields comparable power-spectrum peaks, the V_R oscillation is not uniquely attributable to a tidally induced spiral arm. Alternatively, if the two dominant peaks move significantly when the longitude window is expanded from 10° to 20° or shifted by 5°, the coherent single-arm assumption fails.","supporting_citations":[],"review_version":1}