{"id":"78633cf9-51e4-41fd-9ff2-53496195205b","arxiv_id":"2507.19579","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The phase spiral's winding and rotation phase are nearly constant across the Milky Way disk, suggesting a global perturbation and a winding time that increases steeply with Galactocentric radius.","lead":"Astronomers mapped the vertical 'phase spiral' of stars across the Milky Way disk using Gaia data, reaching 4 kpc from the Sun. They find the spiral's winding and orientation are surprisingly uniform across the disk, hinting that one or more large-scale events stirred the galaxy, and implying that the inferred winding time rises steeply with distance from the Galactic center.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Distant-bin selection effects could masquerade as a flat φ600 profile; the paper's own amplitude discontinuity at the v_l.o.s. cut shows spatially varying systematics that need a phase-bias test.","rationale":"The reader's weakest assumption identifies the same load-bearing point: the flat φ600 profile is central, and it relies on selection effects biasing only amplitude, not phase. I agree with that identification. The paper's own Figure 12/Appendix C showing an amplitude discontinuity at the v_l.o.s. cut is internal evidence that the two sample classes have different systematics; while amplitude bias does not prove phase bias, it raises the stakes for the phase-bias assumption. The reader's other concerns (lack of error bars, post hoc sample selection, model dependence of winding-time conversion) are real but secondary; the global-perturbation claim would fail specifically if the φ600 flatness is a selection artifact. The proposed synthetic-injection or overlap-region consistency test directly targets that failure mode. Hence I recommend keeping the CONDITIONAL verdict, with the condition being a dedicated selection-bias test for phase recovery in distant bins.","tokens_in":22259,"tokens_out":2373,"duration_ms":20342,"concrete_test":"Inject mock phase spirals with known radial slopes in φ600 and ω into synthetic Gaia-like catalogs with realistic dust extinction, crowding, parallax quality cuts, and BNN v_l.o.s. prediction errors; run the paper's exact fitting pipeline on both nearby (v_l.o.s.-required) and distant (BNN) spatial bins; if the recovered φ600 profiles are biased toward flatness in the distant bins by more than the intrinsic scatter, the global-perturbation conclusion is not secure. Alternatively, compare φ600 from the proper-motion sample against an independent v_l.o.s.-only sample in the overlap region (1.6–2.5 kpc), using only stars with measured v_l.o.s.; a systematic offset between the two would indicate selection-dependent phase bias.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that the phase spiral was sourced by one or many global perturbations, based on close-to-flat radial profiles of winding ω and especially rotation phase φ600 across a region spanning about 6 kpc in Galactocentric radius. This conclusion assumes that spatial selection effects in the distant proper-motion samples bias only the spiral amplitude, not its phase. Section 3 states that BNN v_l.o.s. uncertainties of 25–30 km/s were tested in Widmark et al. (2022a,b) to not significantly bias the inferred spiral shape, but the present analysis reaches larger distances and fits morphological parameters (ω and φ600) that those tests were not designed to validate. The paper itself shows indirect evidence of spatially varying systematics: in Appendix C and Figure 12, the inferred amplitude (α+β) has a discontinuity coinciding exactly with the v_l.o.s. cut at 1.6 kpc, demonstrating that the two sample classes do not share identical selection. If dust, crowding, or BNN systematics preferentially suppress stars in particular (z,w) phases, φ600 could be biased toward a common value, manufacturing the flat radial profile that drives the global-perturbation conclusion. The renormalization in Eq. (16) fixes the z-marginal density to sech² and does not remove phase-dependent selection. Since the paper provides no quantified uncertainty on the spatial variation of φ600, this unvalidated selection-bias assumption is the weakest link.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper maps the vertical phase spiral's properties across the Milky Way disk using Gaia DR3 proper motions with BNN-predicted line-of-sight velocities, together with a spatially local sample requiring measured line-of-sight velocities. For each spatial or phase-space bin, the authors fit a mixture-model background plus a spiral perturbation, extracting the potential scaling AΦ, winding ω, rotation phase φ600, and single/double-arm amplitudes. They report that ω and especially φ600 vary smoothly with close-to-flat radial profiles over roughly 6 kpc in Galactocentric radius, which they interpret as evidence that the phase spiral was sourced by one or many global perturbations. They further derive that the implied winding time increases strongly with radius, from about 150 Myr at 7 kpc to about 600 Myr at 9 kpc. The paper includes comparisons with a test-particle simulation and with prior observational work, and emphasizes the complementary nature of the spatial and phase-space binning schemes.","tokens_in":22708,"tokens_out":4916,"duration_ms":46757,"significance":"If the reported uniformity of the rotation phase is real, it is a strong and direct observational constraint on the origin of the phase spiral, favoring global perturbations (e.g., a satellite or large-scale disk response) over spatially local stochastic sources. The paper's use of two complementary binning schemes, its detailed morphological parametrization, and the explicit comparison with a test-particle simulation in Appendix B are notable strengths, and the derived winding-time profile offers a falsifiable benchmark for future self-consistent simulations. However, the central quantitative claim is currently supported without reported parameter uncertainties, and the distant-bin phase-bias assumption is validated only indirectly through earlier work. These issues are load-bearing for the flatness and global-perturbation conclusions.","major_comments":[{"comment":"The central claim that the winding ω and rotation phase φ600 have close-to-flat radial profiles is not supported by any reported uncertainty on the fitted parameters. Figures 4 and 6 show scatter points with no error bars, and Section 5 states that fits were selected by eye and 'dubious' fits omitted; Figure 9 additionally excludes seven high-winding-time outliers. Without a quantitative slope fit and a per-bin uncertainty budget (including the systematic contribution from selection), it is not possible to assess whether the data are actually inconsistent with the steep idealized-model profiles shown in Figure 8. Please provide parameter uncertainties and a statistical test of the flatness claim.","section":"§5.1, Figs. 4 and 6"},{"comment":"The global-perturbation conclusion rests on the assumption that spatial selection effects in the distant bins (dust, crowding, BNN v_l.o.s. uncertainties of 25–30 km/s) bias only the spiral amplitude and not its phase. The paper cites tests in Widmark et al. (2022a,b) for this claim, but those tests were designed for gravitational-potential inference, not for rotation-phase morphology at distances up to 4 kpc. Appendix C and Figure 12 show a discontinuity in the inferred amplitude α+β exactly at the 1.6 kpc v_l.o.s. cut, demonstrating that the two sample classes have different spatially varying systematics. The renormalization in Eq. (16) only enforces a fixed z-marginal density profile and cannot remove phase-dependent incompleteness. A quantitative phase-bias test (e.g., injecting a synthetic spiral into the distant samples and recovering φ600 as a function of distance and extinction) is needed before the flat φ600 profile can be used as evidence for a global perturbation.","section":"§3, §4.2, Appendix C, Eq. (16)"},{"comment":"The inferred strong slope of the winding time (150 Myr at R=7 kpc to 600 Myr at R=9 kpc) is not a direct observable but is derived from the fitted ω, the fixed shape of the vertical potential, and the paper's own fitted exponential disk scale length of 2.9 kpc. The text acknowledges that the anharmonicity parameter Γ can vary with radius, but only states that this affects ω by 'a few ten per cent' without evaluating the impact on the derived tω slope. Since the claimed slope spans a factor of roughly four, a 30% radial variation in Γ could absorb a substantial part of it. Please quantify the sensitivity of the winding-time profile to the assumed potential shape and disk scale length, and show the resulting range of tω(R).","section":"§5.3–5.4, Eq. (28), Fig. 9"}],"minor_comments":[{"comment":"The word 'complimentary' should be 'complementary' in the description of the two binning schemes.","section":"Section 1"},{"comment":"The amplitude threshold α+β≥0.14 applied to the phase-space-binned results is not justified in the text; please state how many data samples are removed and whether Figures 6 and 9 are robust to the chosen threshold.","section":"§5.2 and Figure 6"},{"comment":"The seven excluded high-winding-time outliers and the 'dubious' fits should be explicitly counted and, where possible, listed, so the reader can assess how much of the radial profile is determined by the retained samples.","section":"Figure 9 caption and Section 5"},{"comment":"The numerical constant 71.9 Myr/rad is introduced without derivation; it should be derived from the assumed vertical potential and anchor heights, or referenced to an appendix equation.","section":"Equation (28)"},{"comment":"The 'small Gaussian noise' added to Rg and vR should be described as a plotting jitter applied for visibility, not as measurement uncertainty, to avoid confusion.","section":"Figure 6 caption"},{"comment":"The statement that a uniform perturbation time cannot be reconciled with the data assumes self-gravity acts only in the vertical dimension; stating this assumption is useful, but a test against a three-dimensional simulation would strengthen the argument.","section":"Section 5.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for astro-ph.GA and presents a potentially important observational result. The main weaknesses are the absence of reported parameter uncertainties and the reliance on prior selection-effect tests beyond their validated distance range. The by-eye fit selection and explicit exclusion of outliers make it difficult to verify the flatness claim independently. I would encourage the editor to request a reproducibility table of all fitted parameters (including the excluded bins) and a quantitative selection-injection test for the rotation-phase bias."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The payoff here is the map: spatial maps of phase-spiral winding, rotation phase, and amplitude out to 4 kpc, plus a complementary phase-space binning that covers epicyclic motion. That is genuinely new and will be a useful benchmark. The paper also does a decent job of framing the physical question: the near-flat radial profiles of ω and φ600, and the contrast with the steep trend expected for a uniform perturbation time, are a sharp way to pose the global-versus-local origin problem. The authors are honest about the conundrums, and the comparison with prior work, especially Antoja et al. 2023, is balanced.\n\nThe soft spots are real but not fatal. The absence of any reported uncertainties on the fitted parameters is the biggest one. Without error bars, statements like “close-to-flat” are hard to evaluate, and the small-scale features in ω (the Local Arm correlation) could be noise. The by-eye rejection of fits is also uncomfortable; even if it is standard practice in this area, it makes the sample selection hard to audit. Same for the amplitude threshold and the exclusion of seven high-winding-time outliers. These choices are defensible, but they need to be stated as criteria, not post hoc curation.\n\nThe winding-time conversion is model-dependent: it uses the paper's own fitted disk scale length (2.9 kpc) and an assumed exponential profile, so the strong radial slope from ~150 Myr at 7 kpc to ~600 Myr at 9 kpc is as much a consequence of the model as of the data. The abstract overstates this as an implication. The selection-effect worry in the stress test is not demonstrated, but it is not dismissed either: the amplitude discontinuity right at the v_l.o.s. cut (Figure 12) shows that the two sample classes do differ systematically, and the assumption that this biases only amplitude, not phase, is borrowed from closer-range tests. That deserves a dedicated test.\n\nWho gets value: anyone working on phase spirals, disk perturbations, or the Milky Way's vertical potential. It deserves a serious referee. My recommendation: send it out, and require the authors to add error bars on all fitted parameters, make the sample selection algorithmic, and either validate or caveat the phase-bias assumption at 2–4 kpc.","headline":"Genuinely new maps of phase-spiral morphology across the disk, with a plausible but not airtight case for a global-perturbation origin; worth refereeing, but it needs error bars and a harder look at selection systematics.","tokens_in":23122,"tokens_out":1393,"would_cite":true,"duration_ms":17308,"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":"The Milky Way's phase spiral was likely set by a global perturbation, because its rotation phase is uniform across the disk.","keywords":["phase spiral","Milky Way disk","Gaia DR3","vertical phase space","disk perturbation","winding time","phase mixing"],"falsifier":"Restrict the analysis at distances beyond 1.6 kpc to stars with directly measured line-of-sight velocities, or redo the fit with a deeper future data release, and compare the inferred $\\varphi_{600}$ radial profile to the one obtained with predicted velocities; a flat profile that persists would support the global-perturbation claim, while a radial trend appearing in the velocity-confirmed sample would show the flatness was a selection artifact.","tokens_in":22056,"feed_emoji":"🌀","tokens_out":10061,"duration_ms":86535,"temperature":0.7,"pith_summary":"This paper maps the phase spiral—a spiral pattern in the distribution of stars' height and vertical velocity in the Milky Way disk—across a region several kiloparsecs wide. Using two complementary samples, one based on Gaia proper motions out to 4 kpc and one based on nearby line-of-sight velocities binned by orbit, the authors fit the spiral's shape in hundreds of small volumes. They find that two shape parameters, the winding and the rotation phase, vary smoothly and are nearly flat as functions of Galactocentric radius, even though the disk's vertical gravity changes considerably from the inner to the outer disk. The near-uniform rotation phase is presented as evidence that the spiral was produced by one or several global perturbations, rather than by many small local events. If that conclusion holds, the inferred winding time rises steeply with radius, from roughly 150 million years at 7 kpc to 600 million years at 9 kpc.","feed_headline":"One global disturbance shaped the Milky Way's phase spiral","feed_subtitle":"Uniform rotation phase across the disk points to one global source and a winding time that rises steeply with radius.","key_machinery":"The argument runs on a morphological parametrization of the spiral that avoids degeneracies between physical quantities. The model writes the spiral as a small relative overdensity on a smooth background, with a phase function $\\varphi(E_z)=\\varphi_{\\rm init.}+2\\pi t_\\omega/P(E_z)$ that advances with vertical period $P(E_z)$; the fitted parameters are the vertical potential scaling $A_\\Phi$, the winding $\\omega$ between anchor heights 300 and 800 pc, the rotation phase $\\varphi_{600}$ at 600 pc, and the one- and two-armed amplitudes $\\alpha$ and $\\beta$. The data side uses two complementary binnings: a hexagonal spatial grid with 400 pc spacing reaching 4 kpc, built from Gaia proper motions plus neural-network-predicted line-of-sight velocities where needed, and a nearby sample split in $v_R$ and $v_\\phi$ into 508 phase-space bins. Selection effects are handled by renormalizing each vertical phase-space histogram to a fixed density profile in $z$, under the assumption that incompleteness is mostly a function of height. This lets the authors extract the spiral's shape even where dust and crowding remove many stars, and then interpret the fitted morphology in terms of winding time using the vertical potential at each radius.","core_discovery":"The paper's central claim is that the phase spiral's present-day morphology is set by a globally sourced perturbation. Fitting the relative density perturbation $S(E_z,\\theta_z)=\\alpha\\cos(\\theta_z-\\varphi(E_z))+\\beta\\cos(2\\theta_z-2\\varphi(E_z))$, with vertical energy $E_z$ and phase angle $\\theta_z$, the authors measure the winding $\\omega$ between heights 300 and 800 pc and the rotation phase $\\varphi_{600}$ at height 600 pc across the disk. Both parameters are close to flat in Galactocentric radius over roughly 6 to 11 kpc, with only small-scale excursions, most notably a high-winding band near $R\\simeq 9$ kpc that lines up with the Local Arm. The uniformity of $\\varphi_{600}$ is the paper's key evidence that the spiral was sourced by one or many global perturbations: local stochastic sources would not produce a correlated rotation phase over such a large area. A direct corollary is that the winding time has a strong radial slope, increasing from roughly 150 Myr at 7 kpc to 600 Myr at 9 kpc, which the paper notes cannot be explained by a single uniform perturbation time in one-dimensional vertical dynamics.","pith_inferences":["If the global-perturbation interpretation is right, the smooth gradient of rotation phase across azimuth could encode the propagation of the perturbation through the disk; comparing it with back-propagated star positions might pin down the event time more precisely than current estimates.","A perturbation time that scales inversely with disk surface density would naturally produce flat $\\omega$ and $\\varphi_{600}$ profiles; this is a concrete, testable rule for self-consistent simulations of satellite encounters.","Because the winding responds to the vertical potential's anharmonicity, the measured winding map could be inverted, in combination with density maps, to constrain the cold-gas distribution across the disk.","The same analysis applied to deeper proper-motion catalogs could extend the flat-rotation-phase test beyond 4 kpc and separate the global source from outer-disk warps or tidal effects."],"forward_implications":["If the rotation phase is truly uniform, the phase spiral was seeded by a global event or several events acting coherently, rather than by many independent small-scale perturbations; this disfavors dark-matter subhalo or gas-turbulence sourcing as the dominant mechanism.","The inferred winding-time gradient means the same spiral is younger in the inner disk: about 150 Myr at 7 kpc versus 600 Myr at 9 kpc, so simple uniform-age models for the spiral must be revised.","The high-winding band near $R\\simeq 9$ kpc, aligned with the Local Arm, implies that local variations in the vertical potential's anharmonicity, possibly from cold gas, modulate how quickly the spiral winds.","The maps give quantitative benchmarks: any self-consistent simulation of a satellite or bar perturbation should reproduce the smooth rotation phase and the steep winding-time gradient, not just the spiral's existence."],"supporting_citations":[{"why":"Provides the Gaia astrometry that all stellar samples are built from.","marker":"Gaia Collaboration et al. 2016"},{"why":"Supplies the spectro-astrometric XP parallaxes used to replace Gaia parallaxes where more precise.","marker":"Zhang et al. 2023"},{"why":"Provides the Bayesian neural network line-of-sight velocity predictions used for distant stars lacking measured radial velocities.","marker":"Naik & Widmark 2024"},{"why":"Demonstrates that the Gaia proper motion sample can resolve phase spirals at distances of several kiloparsecs, the basis for the spatial binning.","marker":"Widmark et al. 2022b"},{"why":"Contains the tests the paper relies on to argue that selection effects bias spiral amplitude rather than phase.","marker":"Widmark et al. 2022a"},{"why":"Discovered the phase spiral and defined the phenomenon whose origin and evolution this paper constrains.","marker":"Antoja et al. 2018"},{"why":"Supplies the upper main sequence star map used to compare high-winding regions with Milky Way spiral arms.","marker":"Poggio et al. 2021"},{"why":"Provides the solar neighborhood matter density model from which the vertical potential $\\Phi_\\odot(z)$ is constructed.","marker":"Schutz et al. 2018"}],"fun_headline_variants":["Phase spiral hints at one global disturbance shaping the Milky Way disk","Uniform phase spiral reveals a single global trigger across the disk","Milky Way's phase spiral points to one global disturbance, not local jitters","Global origin for Milky Way's phase spiral: uniform rotation phase across the disk","Phase spiral's uniform spin across the Milky Way implies a single source"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The flat radial profile of the rotation phase, and with it the global-perturbation conclusion, assumes that dust extinction, stellar crowding, and the 25–30 km/s uncertainties of predicted line-of-sight velocities bias only the spiral's amplitude and not its measured phase.","fun_headline_variants_meta":{"raw":{"variants":["Phase spiral hints at one global disturbance shaping the Milky Way disk","Uniform phase spiral reveals a single global trigger across the disk","Milky Way's phase spiral points to one global disturbance, not local jitters","Global origin for Milky Way's phase spiral: uniform rotation phase across the disk","Phase spiral's uniform spin across the Milky Way implies a single source"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001109,"raw_usage":{"total_tokens":4647,"prompt_tokens":994,"completion_tokens":3653,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":3560}},"tokens_in":610,"tokens_out":3653,"duration_ms":19745,"temperature":1.0,"reasoning_tokens":3560,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:52:56.133188+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Restrict the analysis at distances beyond 1.6 kpc to stars with directly measured line-of-sight velocities, or redo the fit with a deeper future data release, and compare the inferred $\\varphi_{600}$ radial profile to the one obtained with predicted velocities; a flat profile that persists would support the global-perturbation claim, while a radial trend appearing in the velocity-confirmed sample would show the flatness was a selection artifact.","supporting_citations":[],"review_version":2}