{"id":"4cd31adc-a3fb-4bb3-90f2-33657fa31938","arxiv_id":"2411.18460","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"A three-component warp-and-ring model built from SITELLE ionized gas kinematics reproduces the perturbed velocity field in M31's central kiloparsec and supports a past head-on M32-like collision.","lead":"Using a new optical velocity map of ionized gas in the center of Andromeda, the authors build a geometric model with three gas structures: the main disk, a tilted ring, and a tilted warped nuclear disk. The model matches the complex velocity patterns and supports the earlier idea that a small galaxy crashed through Andromeda's center about 200 million years ago.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim depends on ignoring triaxial-bulge non-circular motions; the qualitative fit may let the imposed warp/ring geometry absorb bar-like flows.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the model assumes an almost-spherical potential and circular orbits, thereby excluding the triaxial box-peanut bulge and its associated non-circular motions that Section 6 discusses at length. The stress-test confirms this is the most serious threat to the central claim because the fitted warp and ring are inferred from a velocity field that could in principle be produced by bar-like flows. The reader's CONDITIONAL verdict already accounts for this uncertainty, so no verdict change is needed. The proposed test would settle the issue directly: include the published triaxial potential in the same static particle framework and check whether the observed velocity features require the extra geometrical components. The Table 2 orientation inconsistency is real but secondary; it affects reproducibility rather than the physical conclusion, and it is already noted in the reader's rationale. Overall, the paper presents a plausible qualitative scenario with new SITELLE data, but the central collision-and-reorientation narrative is not uniquely established until the circular-orbit assumption is tested against the triaxial-bulge interpretation.","tokens_in":21835,"tokens_out":3261,"duration_ms":35309,"concrete_test":"Re-run the Section 4 particle-forward model with the same Table 1 mass model, but replace the spherical bulge with the Blaña Díaz et al. (2017) triaxial box-peanut bulge potential at pattern speeds of 20 and 40 km/s/kpc, keeping only the main disk and no imposed nuclear warp or offset ring. Compare the predicted first- and second-moment maps pixel-by-pixel to the SITELLE velocity and dispersion maps in Figs. 5 and 8, restricted to the northeast region used for the fit, using a quantitative metric such as chi-square or median absolute deviation. If the triaxial model alone reproduces the central winding and the southeast crescent within the quoted 10-60 km/s velocity uncertainties, the warped nuclear disk and tilted ring are not supported; if it does not, the three-component geometry is genuinely required.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4.1 states that the central kpc potential 'can be considered almost spherical', and Section 4.2 places gas on circular orbits with only an asymmetric-drift correction. This assumption is load-bearing because the fitted inclinations, position angles, and offsets of the warped nuclear disk and tilted ring in Section 5.3 and Table 2 are all derived from this circular-orbit model. Section 6 itself summarizes the evidence that the central bulge is triaxial and box-peanut shaped, with pattern speeds of 20-40 km/s/kpc (Blaña Díaz et al. 2017, 2018; Feng et al. 2022, 2024), and that Opitsch et al. (2018) see non-circular, S-shaped features in the ionized gas velocity field. If non-circular motions contribute substantially to the observed central winding and the southeast velocity crescent, then the tilted nuclear warp and the offset ring are not uniquely required; they could be artifacts of fitting a spherical-potential, circular-orbit model to a barred or triaxial flow. This concern is compounded by the lack of a quantitative fit: Section 5.2 says the southwest region was excluded because of 'clear velocity perturbations, plausibly due to a large-scale shock', and the parameters are hand-fitted and compared only qualitatively. A secondary internal inconsistency in Table 2, where the main disk is listed with inclination 35 degrees and PA 77 degrees while the text and conclusion give 77 degrees and PA 37 degrees, further impedes reproduction, but the central issue remains the unmodelled triaxial potential.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a dynamical model for the central kiloparsec of M31, combining newly reduced SITELLE Hα/[NII] integral-field data with archival CO(2-1) and dust maps. The authors identify three gas components—a main disk, a tilted inner ring, and a warped nuclear disk—and construct a static, axisymmetric stellar potential in which gas particles are placed on circular orbits with an epicyclic/asymmetric-drift correction. The geometrical parameters of the three components are adjusted by hand to reproduce the observed velocity field, and the best qualitative match (scenario 3) is interpreted as evidence for a recent head-on collision with an M32-like galaxy. The paper includes a full description of the data reduction, velocity extraction, and the four geometric scenarios considered.","tokens_in":22197,"tokens_out":2932,"duration_ms":29950,"significance":"If the interpretation holds, the paper provides a coherent dynamical explanation for the complex central velocity field of M31, supporting the long-standing collision scenario of Block et al. (2006). The strengths are the fully sampled ionized-gas velocity field, the joint use of ionized and molecular gas tracers, and the explicit modeling of disk warping and tearing. However, the significance is limited by the qualitative nature of the fit: the parameters are not derived from a quantitative optimization, no uncertainties are given, and the comparison with observations is visual. The central kinematic features are real, but the paper does not yet demonstrate that the three-component geometry is uniquely required rather than one of many possible interpretations.","major_comments":[{"comment":"The fit is qualitative: Section 5.2 states that the geometrical parameters were 'essentially fitted on the northeast region' and Section 5.3 concludes with 'best qualitative model'. No residuals, chi-squared values, or parameter uncertainties are reported, and the southwest region is excluded from the fit because of a presumed shock. As a result, the central claim that the three-component model reproduces the observed velocity field is not quantitatively established. I ask the authors to provide a quantitative comparison (e.g., residual maps between the model and observed velocity fields, a goodness-of-fit statistic, and an exploration of parameter degeneracies) or to clearly state the model as a schematic interpretation rather than a fitted model.","section":"§5.2 and Tables 2–3"},{"comment":"The assumption that the central kpc potential is 'almost spherical' and that gas follows circular orbits (Eq. 4) is load-bearing, because the fitted warps, inclinations, and offsets are derived from this model. Yet Section 6 itself reviews evidence for a triaxial box-peanut bulge with pattern speeds of 20–40 km/s/kpc and for non-circular, S-shaped motions in the ionized gas (Opitsch et al. 2018; Feng et al. 2022, 2024). If non-circular motions contribute substantially to the observed velocity winding and the southeast crescent, the warped nuclear disk and tilted ring may be artifacts of fitting a circular-orbit model to bar-like flows. The authors should quantitively estimate the amplitude of non-circular motions within 1 kpc (e.g., from the Feng et al. models) and show that they are negligible compared to the 100–300 km/s features, or repeat the modeling in the published triaxial potential and demonstrate that the warps and ring are still required.","section":"§4.1 and §6"},{"comment":"The geometrical parameters (inclinations, position angles, offsets, warp extent) are fitted to the same velocity map that is subsequently compared with the model prediction. The agreement therefore mostly reflects the fitting procedure, not an independent prediction. This circularity should be acknowledged explicitly, and the authors should provide an independent test. For example, they could fit the model to only a subset of the velocity field (e.g., the NE region) and compare the prediction to the SE and SW regions, or predict the CO line profiles at positions not used in the fit and compare with the IRAM-30m data. Without such a test, the 'confrontation' of predicted and observed velocity fields does not validate the three-component scenario over alternative interpretations.","section":"§5.3 and Figs. 21–22"},{"comment":"There is an inconsistency in the orientation of the main disk: Table 2 lists inclination 35° and position angle 77°, while the text and the conclusion (Section 7) state inclination 77° and PA 37°. This is not a typographical nuance; it directly affects the reproduction of the model and the interpretation of the projected geometry. The authors must correct the table and confirm which values were actually used in the modeling.","section":"Table 2 and §7"}],"minor_comments":[{"comment":"The reference to 'According to ?' is an unresolved citation and should be completed.","section":"§1"},{"comment":"The notation '||− →v||(r) = r sqrt(r dΦ/dr)' is ambiguous; please clarify the square-root placement and the definition of r in cylindrical versus spherical coordinates.","section":"Eq. (4)"},{"comment":"The phrase 'cube moments 1 and 2' is used without definition; for clarity, specify that these are the intensity-weighted velocity and velocity-dispersion maps.","section":"§5.3"},{"comment":"The captions of the bottom panels of Fig. 5 and the color coding in Fig. 11 would benefit from a more explicit description of what the green triangles represent and how the flux-ratio threshold for the second component was chosen.","section":"Figs. 5 and 11"},{"comment":"The note that only varied parameters are displayed is helpful, but the table should explicitly state that unlisted parameters are identical to those in Table 2, to avoid confusion when comparing scenarios.","section":"Table 3"},{"comment":"The term 'teared apart' is used repeatedly; 'torn apart' is the standard English form. This is a wording issue only.","section":"§4.3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents interesting observational material and a plausible scenario, but the modeling falls short of current standards for quantitative comparison. The main issues are the hand-tuned parameters, the absence of uncertainties, and the treatment of the triaxial-bulge potential. These are addressable in revision, but they are load-bearing for the central claim. I recommend that the editor seek a revised version in which the authors provide a more rigorous fit and address the circularity concern explicitly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nBottom line: this paper's real asset is the fully sampled SITELLE Hα/[NII] velocity field of M31's central kiloparsec, plus a concrete three-component geometric model—warped nuclear disk, tilted offset ring, and main disk—built in the Raj et al. disk-tearing language. If that geometry holds, it's a useful step for M31 and minor-merger studies. The data cube is public, which is a point in its favor.\n\nWhat's actually new: the SITELLE map and the attempt to fit this three-component warped geometry to it. That is worth having. The paper also does a fair job arguing against the competing bar interpretation in Section 6, noting that hydrodynamic bar models produce a high central gas concentration that M31 does not show.\n\nThe soft spots are real and not manufactured. The fit is qualitative: parameters in Tables 2 and 3 are hand-adjusted, no uncertainties are quoted, and the southwest region is excluded because of a presumed shock. That alone would make me want a quantitative follow-up. More importantly, Section 4.1 assumes an almost spherical potential and puts gas on circular orbits with an asymmetric-drift correction. The paper itself reviews evidence that the central bulge is triaxial and box-peanut, with pattern speeds of 20–40 km/s/kpc, and that Opitsch et al. see S-shaped non-circular ionized-gas motions. If those non-circular flows are substantial, a spherical-potential model can absorb them into warp and ring geometry. The authors' rebuttal—that bar models fail to reproduce the gas morphology—is suggestive but does not prove the warped-ring geometry is unique. There is also a direct internal contradiction: Table 2 lists the main disk with inclination 35° and PA 77°, while the text and conclusion give 77° and PA 37°. The broken '?' citation placeholders should also be fixed. And the model-versus-data comparison is partly circular, since the geometry was fitted to the same velocity map it is then 'confronted' with; an independent tracer or a pre-registered model would be a cleaner test.\n\nWho is this for? M31 specialists, and anyone working on minor mergers, disk tearing, or nuclear gas kinematics. The new data alone justify a serious referee pass; the modeling needs more rigor.\n\nMy recommendation: send it to review, with the clear expectation that the quantitative fitting, the potential assumption, and the Table 2 inconsistency get addressed.\n\nBest,\n[You]","headline":"A genuinely useful SITELLE velocity field and a plausible three-component warp-and-ring model for M31's center, but the qualitative fit and the spherical-potential assumption keep it conditional rather than accepted.","tokens_in":22743,"tokens_out":1600,"would_cite":true,"duration_ms":16131,"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 paper claims that M31's chaotic central velocity field is a warped, torn nuclear disk reorienting after a recent head-on collision with an M32-like galaxy.","keywords":["galaxies: individual: M31","galaxies: kinematics and dynamics","nuclear disk","warped disk","disk tearing","integral field spectroscopy","molecular gas kinematics","galaxy collisions"],"falsifier":"A high-resolution stellar kinematic map of the central kiloparsec that measures the actual gravitational potential and its pattern speed would settle the claim: if the box-peanut bulge's non-axisymmetric torques are strong enough to drive the observed S-shape and velocity jumps, then a circular-orbit model in a spherical potential would fail, and a bar-plus-shock model would reproduce the same SITELLE cube. Concretely, measuring the full non-circular velocity field (radial and tangential streaming terms) from the SITELLE Hα/[NII] cube would show whether large radial streaming appears along the bar's major axis rather than the warp geometry.","tokens_in":21620,"feed_emoji":"🌀","tokens_out":5252,"duration_ms":45457,"temperature":0.7,"pith_summary":"This paper claims that the complex, asymmetric velocity field of ionized and molecular gas in the central kiloparsec of M31 is produced by three gas components moving in a stellar-dominated potential: the main disk (inclination 77°, position angle 37°), an off-centered tilted ring, and a warped nuclear disk. The authors build a static dynamical model in which 2.4 million particles follow near-circular orbits in an almost spherical, bulge-dominated potential, then rotate and offset each concentric ring to reproduce the observed twisting, double-peaked, and torn velocity structures. They find that the best configuration (scenario 3) has the three structures decoupled, with a warped zone where the nuclear disk tears into the 1 kpc ring, matching both the SITELLE Hα/[NII] velocity field and the IRAM CO(2-1) double components. If correct, this explains the long-standing puzzle of M31's central gas hole, double velocity components, and lopsided ring in a collision scenario rather than a bar-driven scenario.","feed_headline":"M31's chaotic core is a warped, torn disk settling after a collision","feed_subtitle":"New velocity maps fit M31's central gas to a tilted ring plus a warped nuclear disk left by an M32-like head-on hit.","key_machinery":"The central machinery is a static three-component dynamical model: a main disk, a tilted offset ring, and a nuclear warped disk, built as a set of misaligned concentric circular orbits in a gravitational potential composed of Plummer spheres for the dark halo, bulge, and nucleus, and Miyamoto-Nagai disks for the stellar and gaseous disks. Each gas particle is assigned the circular-orbit speed in this potential, corrected by an epicyclic asymmetric drift with Toomre parameter Q = 1.3, and then each radius is given its own inclination, position angle, and two sky-plane offsets to create the warp and ring geometry. The fitting procedure compares the predicted velocity field and second moment maps to the SITELLE Hα/[NII] kinematics and the IRAM-30m CO(2-1) data across four geometric scenarios, with scenario 3 (fully decoupled structures) giving the best qualitative match.","core_discovery":"The kinematical observations of the central kiloparsec of Andromeda correspond to a dynamical re-orientation of the perturbed nuclear disk: a series of warps tears the disk into an offset 1-kpc ring, following a recent head-on collision with an M32-like galaxy. The paper argues that the twisted velocity winding within about 200 pc is the signature of a warped nuclear disk that is almost face-on at the center, while the two widely separated velocity components seen along the minor axis in both Hα and CO are the projection of a tilted, offset inner ring superposed on the main disk. The best-fitting geometry has the main disk at inclination 77° and position angle 37°, an inner ring centered about 322 pc north of the nucleus with its own inclination and position angle, and a nuclear warped disk connecting to the ring through a torn transition zone. This structure is presented as the kinematic fingerprint of a collision that tilted and warped the nuclear disk, which is now settling back to equilibrium with the m=1 waves damping down.","pith_inferences":["Beyond the paper: if the collision scenario is correct, the fitted 322 pc north offset of the ring center should correlate with the wake of the M32-like encounter; a hydrodynamical simulation of the collision could predict the ring offset and warp amplitude that the static model leaves as free parameters.","Beyond the paper: the circular-orbit assumption could be tested with a stellar kinematic map of the central kiloparsec at matched resolution; if stars show no corresponding warp or ring structure, the gas-only warps would more plausibly be transient shock features rather than a settled disk reorientation.","Beyond the paper: a direct spectral fit of the SITELLE cubes with two-component line profiles at every pixel, rather than comparison to moment maps, would provide a quantitative likelihood comparison between scenario 3 and the bar-plus-shock interpretation."],"forward_implications":["The central gas hole and the double-velocity components along the minor axis are explained by projection of a tilted ring and a warped nuclear disk, rather than by bar-driven shocks.","The off-centered 1 kpc dust and gas ring is a torn remnant of the nuclear disk, not an inner Lindblad resonance ring of a bar.","The triaxial box-peanut bulge can remain as a pre-collision bar remnant while the thin bar in the disk was destroyed, reconciling stellar kinematics with the lack of a bar signature in the gas.","Scenario 3 predicts specific locations where two velocity components should appear along the line of sight, such as the crescent-shaped region in the southwest, which can be checked with higher-resolution integral-field data.","The nuclear warped disk connects naturally to the previously reported inner warp in HI, extending the known 'warp in the warp' down to the central kiloparsec."],"supporting_citations":[{"why":"Supplies the head-on collision scenario with an M32-like companion and the identification of the off-centered 1 kpc and 10 kpc rings that the model interprets as torn and warped gas.","marker":"Block et al. (2006)"},{"why":"Provides the disk-tearing and warping framework used to model misaligned precessing disks and justify the warped nuclear disk tearing into the ring.","marker":"Raj et al. (2021)"},{"why":"Supplies the molecular gas detection with large line splittings, the double-component interpretation, and parameters for the dark halo and bulge potentials plus Toomre Q = 1.3.","marker":"Melchior & Combes (2011)"},{"why":"Provides the IRAM-30m CO(2-1) observations with HERA used as the molecular gas kinematics constraint.","marker":"Melchior & Combes (2013)"},{"why":"Delivers the SITELLE SN3 data cube reduction and the ionized gas kinematics extraction method.","marker":"Martin et al. (2018)"},{"why":"Provides the previous [OIII] ionized gas velocity field with two components and the bar-interpretation baseline that the paper argues against.","marker":"Opitsch et al. (2018)"},{"why":"Supplies the triaxial box-peanut bulge density model and pattern speed used to discuss the bar alternative and the pre-collision bar remnant.","marker":"Blaña Díaz et al. (2017)"},{"why":"Relevant bar hydrodynamical simulations that require a tilted nuclear disk with no physical reason, which the paper contrasts with its collision-induced warp scenario.","marker":"Feng et al. (2024)"},{"why":"Defines the Miyamoto-Nagai potential used for the disk components in the mass model.","marker":"Miyamoto & Nagai (1975)"}],"fun_headline_variants":["M31's nuclear disk warps and tears after M32-like hit","Andromeda's core shows collision aftermath: warped, torn disk","M31's central gas reveals tilted ring from head-on collision","Warped nuclear disk in M31 points to recent galaxy collision","M31's inner disk is settling after a M32-style smash"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central kiloparsec potential is treated as almost spherical and bulge-dominated, so all gas is put on circular orbits with only a small asymmetric-drift correction; if non-circular motions from the triaxial box-peanut bulge contribute significantly to the velocity kinks, the fitted warped disk and tilted ring would be artifacts.","fun_headline_variants_meta":{"raw":{"variants":["M31's nuclear disk warps and tears after M32-like hit","Andromeda's core shows collision aftermath: warped, torn disk","M31's central gas reveals tilted ring from head-on collision","Warped nuclear disk in M31 points to recent galaxy collision","M31's inner disk is settling after a M32-style smash"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000239,"raw_usage":{"total_tokens":1594,"prompt_tokens":1106,"completion_tokens":488,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":722,"completion_tokens_details":{"reasoning_tokens":396}},"tokens_in":722,"tokens_out":488,"duration_ms":4626,"temperature":1.0,"reasoning_tokens":396,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:10:41.288434+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-resolution stellar kinematic map of the central kiloparsec that measures the actual gravitational potential and its pattern speed would settle the claim: if the box-peanut bulge's non-axisymmetric torques are strong enough to drive the observed S-shape and velocity jumps, then a circular-orbit model in a spherical potential would fail, and a bar-plus-shock model would reproduce the same SITELLE cube. Concretely, measuring the full non-circular velocity field (radial and tangential streaming terms) from the SITELLE Hα/[NII] cube would show whether large radial streaming appears along the bar's major axis rather than the warp geometry.","supporting_citations":[{"cited_title":"L., Bournaud, F., Combes, F., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the head-on collision scenario with an M32-like companion and the identification of the off-centered 1 kpc and 10 kpc rings that the model interprets as torn and warped gas."},{"cited_title":"J., & Do˘gan, S","cited_arxiv_id":null,"evidence_quote":"Provides the disk-tearing and warping framework used to model misaligned precessing disks and justify the warped nuclear disk tearing into the ring."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the molecular gas detection with large line splittings, the double-component interpretation, and parameters for the dark halo and bulge potentials plus Toomre Q = 1.3."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the IRAM-30m CO(2-1) observations with HERA used as the molecular gas kinematics constraint."},{"cited_title":"B., Drissen, L., & Melchior, A.-L","cited_arxiv_id":null,"evidence_quote":"Delivers the SITELLE SN3 data cube reduction and the ionized gas kinematics extraction method."},{"cited_title":"H., Saglia, R","cited_arxiv_id":null,"evidence_quote":"Provides the previous [OIII] ionized gas velocity field with two components and the bar-interpretation baseline that the paper argues against."}],"review_version":1}