{"id":"2d8fb005-36ed-42a0-baa8-8734e847e4ed","arxiv_id":"2506.19074","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A 6D census of 8,456 bulge RR Lyrae stars shows that the Milky Way's spherical bulge component is made largely of stars orbiting backward relative to the bar, and this population can form secularly without a merger-built classical bulge.","lead":"Using 8,456 RR Lyrae stars with full 6D motions, this paper maps the orbits of the Milky Way's bulge population and shows that stars moving backward relative to the bar form a spherical, slowly rotating component. The authors argue, with a companion simulation, that this component may be a secular product of bar evolution rather than a merger-built classical bulge.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Simulation underpredicts observed retrograde fraction; paper's own Fig. 9 shows the gap, undermining the secular-origin claim without a quantitative classical-bulge comparison.","rationale":"The reader correctly identified reliance on a single isolated simulation as the weakest assumption. My stress-test refines this into a concrete, internal quantitative mismatch: the simulation's retrograde fraction is roughly a factor of two lower than the observed fraction in the central bins, a discrepancy the paper itself states in Section 5.3 and Figure 9. This is more specific than 'not representative' because the paper's own numbers show the secular mechanism underproduces the very population it is invoked to explain. The empirical core (6D orbits, interloper removal, prograde/retrograde separation, banana-orbit distance bimodality) is solid, and the paper is honest about the limitation, so the verdict should remain CONDITIONAL as the reader decided. However, the condition needs to be explicit: the authors should either demonstrate that the simulation can quantitatively reproduce the observed retrograde fraction under the correct selection, or present a two-component model that includes a classical bulge and show that the secular component is sufficient. My concern does not reject the paper; it tightens the required revision and prevents the interpretive claim from being read as stronger than the evidence supports.","tokens_in":49072,"tokens_out":8091,"duration_ms":81634,"concrete_test":"Recompute the HG1 simulation's retrograde fraction for the oldest stellar particles in the innermost radial bin (rGC<0.5 kpc), applying the exact observational selection of Eq. 14 plus the apocentric-distance and zmax cuts, and compare to the observed value of ~40% shown in Figure 9. If the simulated fraction falls below ~25%, the secular mechanism alone cannot account for the central spheroidal component, and a classical bulge contribution is required to match the data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central interpretive claim (Section 9) is that the observed spherical, slowly rotating bulge component can be attributed to secular bar evolution, minimizing the need for a classical bulge. The load-bearing quantity is the fraction of stars on retrograde (bar-frame) orbits, which the HG1 simulation should reproduce for the oldest, most centrally concentrated populations that correspond to metal-poor RR Lyrae stars. It does not. In Section 7.2 the oldest, bulge-confined stellar particles have a retrograde fraction of 20±1% (18±1% after footprint cuts), and Section 4.3 reports 14% for all selected particles. The observed RR Lyrae retrograde fraction is ≈28% overall (Section 4.3), ≈30% for [Fe/H]<-2 (Section 5.1), and nearly 40% at rGC<0.5 kpc (Section 5.3, Figure 9). The authors explicitly acknowledge this discrepancy in Section 5.3: the observed fraction is 'considerably higher at rGC<0.5 in the observations than in the simulation,' and 'an additional old, spheroidal classical bulge population at small Galactocentric radii' is also consistent with the data. This is not merely a single-simulation caveat: the paper's own quantitative comparison shows the secular mechanism produces too few retrograde stars by roughly a factor of two in the central region, exactly where the spheroidal component is most prominent. Without a quantitative model demonstrating that the secular population plus observational selection can reach the observed ~40% at rGC<0.5, the phrase 'secular spherical bulge' overstates the support; the data equally favor a mixed origin. The absence of any explicit classical-bulge baseline (e.g., an isotropic spheroid fit to the same volume) means the interpretive claim is not tested against the leading alternative, despite the paper's Section 8 statement that the connection to Rix et al. (2022) and Belokurov & Kravtsov (2022) is deferred to a forthcoming paper.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper assembles an 8,456-star sample of bulge RR Lyrae stars with full 6D phase-space information, computes orbits in a barred Milky Way potential, and classifies stars as prograde or retrograde in the bar frame. It reports that the retrograde fraction increases toward metal-poor and centrally concentrated stars, that retrograde orbits are comparable in regularity to prograde orbits, and that RR Lyrae stars on banana orbits show a bimodal distance distribution similar to the X-shaped bulge traced by red clump stars. Comparing these observations to an isolated N-body+SPH simulation (HG1), the authors argue that the previously identified spheroidal, slowly rotating bulge component can be produced secularly by bar evolution and angular-momentum exchange, reducing the need for a classical bulge formed by mergers.","tokens_in":49340,"tokens_out":8145,"duration_ms":83385,"significance":"The observational dataset is a substantial advance: 8,456 bulge RR Lyrae stars with 6D kinematics, 100 Monte Carlo orbit realizations per star, a publicly available barred potential, pattern-speed robustness tests in Appendix C, and an explicit discussion of photometric-metallicity systematics. The empirical detection of a bar-frame retrograde population with a centrally concentrated, rounder spatial distribution, and the recovery of the X-shape signature through banana-orbit selection, are valuable results that stand independently of the interpretive claims. The data products and derived velocities are made available on Zenodo. If the secular-origin interpretation survives quantitative scrutiny, it would have broad implications for the role of classical bulges in the Milky Way; as it stands, however, the paper's central conclusion in Section 9 is stronger than the simulation comparison supports.","major_comments":[{"comment":"The quantitative comparison between observations and the HG1 simulation does not support the strength of the secular-origin conclusion. For stars at r_GC < 0.5 kpc the observed retrograde fraction is almost 40% for RR Lyrae stars and about 35% for red giants (Fig. 9), while the oldest bulge-confined simulation particles reach only 20±1%, or 18±1% after the Eq. 14 footprint cuts (Section 7.2), and the full selected particle sample gives 14% (Section 4.3). The authors themselves state in Section 5.3 that the observed fraction is \"considerably higher at rGC<0.5 in the observations than in the simulation\" and that \"an additional old, spheroidal classical bulge population at small Galactocentric radii\" is also consistent with the data. Because the central claim of Section 9 is that the spheroidal component \"can be related to the secular evolution\" on the basis of this simulation, the factor-of-two shortfall in the very region where the spheroidal component is most prominent is load-bearing. A quantitative model that adds a classical-bulge component, or a revised claim limited to a secular contribution to part of the spheroidal population, is needed before the interpretive conclusion can be drawn.","section":"Section 5.3, Fig. 9; Section 7.2; Section 4.3"},{"comment":"The secular-origin interpretation rests on a single isolated N-body+SPH simulation (HG1), and the paper's own caveat in Section 8 concedes this: \"additional investigations with a diverse set of N-body and cosmological simulations would be useful.\" No independent model is presented to show that the retrograde fraction, its age/metallicity trend, and its orbital stability are generic consequences of bar slowdown rather than specific to the initial conditions of HG1. The claim that the Milky Way's spherical component \"can be related to the secular evolution\" is therefore under-supported even setting aside the quantitative mismatch in Fig. 9. I would like to see either an additional simulation (for example, a cosmological run or a model with different initial angular-momentum distribution) or a more restricted statement that the simulation demonstrates a mechanism by which some retrograde stars can arise secularly.","section":"Section 8.1 and Section 9; Section 4.1"}],"minor_comments":[{"comment":"The percentages for the [Fe/H] < -2.0 dex sample are inconsistent between the two sections: Section 5.1 gives 40% interlopers and 41% prograde stars, while Section 9 gives 38% interlopers and 43% prograde stars. These should be harmonized to the same sample definition and rounding.","section":"Section 5.1 vs. Section 9"},{"comment":"The phrase \"while more than44 remains prograde\" is missing a percent sign and should read \"more than 44% remain prograde.\"","section":"Section 7.2"},{"comment":"The text says the figure shows \"the dependence of Ωrot_phi on stellar age and rGC,\" but the y-axis of Fig. 9 is the percentage of stars with Ωrot_phi < 0. Please clarify the wording to match the plotted quantity.","section":"Section 5.3, Fig. 9"},{"comment":"The sentence \"which makes distraction and clear separation difficult\" should read \"which makes distinction and clear separation difficult.\"","section":"Section 3.1"},{"comment":"The description of the simulation binning (\"a box with a size of 1000 and a step equal to 350\") does not state the units or what quantity is being binned; please specify, for example the number of particles per boxcar window.","section":"Section 5.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of Astronomy & Astrophysics, and the observational data products are clearly valuable. The main issue is that the title and Section 9 overstate the support for a purely secular origin of the spheroidal bulge component, given the factor-of-two discrepancy acknowledged in Section 5.3 and the reliance on a single isolated simulation. The authors may wish to reframe the conclusion as 'secular contribution to the spherical component' or add a quantitative composite-model comparison."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one for the data product, not for the headline. The 6D sample of 8,456 bulge RR Lyrae stars is the largest so far, the RRc systemic velocities are genuinely first, and the orbit classification is careful. The empirical findings—banana-orbit bimodal distance distribution, linear decline of prograde fraction with metallicity, retrograde stars forming a centrally concentrated round component—are new and probably durable. The authors run 100 Monte Carlo orbit realizations per star, test pattern-speed choices in Appendix C, and are upfront about the limits of the photometric metallicities. That part earns credit.\n\nThe soft spot is exactly where the stress-test note lands. Section 9 concludes the spheroidal component can be related to secular evolution based on one isolated N-body+SPH simulation, HG1. But the simulation's retrograde fraction is ~18–20% for the oldest bulge-confined particles, while the observed RR Lyrae fraction is ~28–30% overall and ~40% at rGC<0.5. The authors see this and say in Section 5.3 that the fraction is 'considerably higher' in the observations and that an additional old classical-bulge population is also consistent. That is not a minor caveat; it is the central interpretive claim failing its own quantitative test. A secular model that underpredicts the retrograde population by a factor of two in exactly the region where the spherical component is strongest cannot, on its own, retire the classical-bulge hypothesis. There is also no explicit classical-bulge baseline fit to the same volume. Without that comparison, 'secular spherical bulge' overstates what the data show; 'a secular component is present, but a classical-bulge contribution is still allowed' would be accurate.\n\nOne smaller point: the prograde/retrograde split is defined by the sign of the azimuthal frequency in the adopted barred potential, so the roundness of the retrograde subset is partly a consequence of which orbits get labeled retrograde. The paper does not use that roundness to validate the classifier, and it should not lean on the roundness as independent evidence.\n\nThe empirical machinery is sound, and the paper is honest about most of these gaps. It deserves a serious referee, but the referee should push hard on Section 9. I would cite the dataset and the orbital classifications, not the secular-bulge conclusion, without qualification. A revision that either adds a quantitative classical-bulge comparison or lowers the claim to match the evidence would be citable on both counts.","headline":"Read it for the 6D RR Lyrae data and orbit classifications; the 'secular spherical bulge' conclusion is not supported by the paper's own simulation comparison.","tokens_in":50113,"tokens_out":2114,"would_cite":true,"duration_ms":22034,"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 paper argues that the Milky Way's spheroidal, slowly rotating bulge component—long read as a merger-built classical bulge—can instead be produced secularly by the bar's angular-momentum exchange, supported by 8,456 RR Lyrae orbits and…","keywords":["Galactic bulge","RR Lyrae stars","stellar orbits","bar-driven secular evolution","retrograde orbits","Milky Way bar","N-body simulations","bulge formation"],"falsifier":"Measure [α/Fe] abundances for the retrograde, centrally concentrated bulge RR Lyrae stars. If the secular scenario is right, these stars should be an old, relatively α-rich population present before bar formation; if a large fraction turn out to be young or α-poor, the claim that the spheroidal component is made by bar-driven angular-momentum exchange would be falsified.","tokens_in":48815,"feed_emoji":"🌌","tokens_out":9929,"duration_ms":97125,"temperature":0.7,"pith_summary":"This paper tries to establish that the Milky Way's spheroidal, slowly rotating bulge component—usually read as the imprint of a merger-built classical bulge—can instead be produced entirely by secular bar evolution. Using full six-dimensional phase-space information for 8,456 bulge RR Lyrae stars, it classifies orbits in the bar's rotating frame and finds that the retrograde (counter-rotating) fraction rises linearly toward metal-poor stars and toward the Galactic center. Those retrograde orbits are as regular as prograde ones and together form a centrally concentrated, nearly spherical structure. An isolated N-body+SPH simulation that forms a bar without any merger reproduces the same trends, with stars oscillating between prograde and retrograde states as the bar exchanges angular momentum and slows down. The paper concludes that the spheroidal element is a real secondary component of the bulge built by the bar itself, so a massive classical bulge is not required.","feed_headline":"The bar alone can build the Milky Way's spherical bulge","feed_subtitle":"A slowly rotating spherical bulge can grow from angular-momentum exchange with the bar, no merger required.","key_machinery":"The load-bearing object is the sign of the azimuthal frequency in the bar's rotating frame, $\\Omega^{\\rm rot}_{\\varphi} = \\Omega^{\\rm ine}_{\\varphi} - \\Omega_{\\rm P}$, which separates bar-supporting prograde stars from retrograde stars (negative values). Around this, the paper builds a classification of orbital families (banana, brezel, fish, and x- and z-tube orbits) from frequency maps, plus a chaoticity measure, the frequency drift $\\log_{10}\\Delta\\Omega$, showing that retrograde and prograde orbits have comparable regularity. The mechanism is completed by the isolated N-body+SPH simulation: tracking old stellar particles across snapshots shows them migrating through the vertical inner Lindblad resonance and inner Lindblad resonance, repeatedly exchanging angular momentum and flipping between prograde and retrograde states, with the retrograde fraction stabilizing within a few gigayears after bar formation.","core_discovery":"On the paper's own terms, the discovery is that the spheroidal component of the Galactic bulge is a genuine secondary structure whose origin is secular rather than accretional. The evidence is orbital: in the bar's rotating frame, roughly 72% of bulge RR Lyrae stars are prograde while the rest are retrograde, and the retrograde fraction grows monotonically from about 10–15% at the metal-rich end to about 30% at the metal-poor end, and toward the center reaches about 40% at Galactocentric radius near 0.5 kpc. Retrograde orbits are just as regular (low frequency drift and Lyapunov exponents) as prograde orbits, so they form a stable, nearly spherical central concentration. In the isolated HG1 simulation, only about seven percent of the retrograde particles stay retrograde over five gigayears; most oscillate across the inner Lindblad resonance, and the prograde/retrograde ratio stabilizes a few gigayears after bar formation. The paper concludes that angular-momentum exchange with the slowing bar can manufacture a classical-bulge-like spheroid from old, low-angular-momentum stars already present before the bar, making a massive merger-built classical bulge unnecessary.","pith_inferences":["If the paper is right, estimates of the Milky Way's accreted stellar mass that count the inner spheroid as a classical bulge would need to be scaled down; separating the secular retrograde component first is the obvious correction.","The paper's own comparison leaves a gap at small Galactocentric radii, where the observed retrograde fraction exceeds the isolated simulation; that region is the natural place to look for a residual merger-built component.","The same orbit-separation test could be applied to other old tracer populations or to barred galaxies in cosmological simulations; a universal linear retrograde-metallicity trend would strengthen the case for bar-driven angular-momentum exchange.","A clean test would be to measure alpha-element abundances for the retrograde bulge stars: the secular scenario requires them to be old and relatively alpha-rich, whereas a large young or alpha-poor retrograde population would point to an accreted origin."],"forward_implications":["Removing retrograde stars from the metal-poor RR Lyrae sample reveals a weak rotation signal that is otherwise invisible, explaining the apparent absence of rotation in the metal-poor bulge.","RR Lyrae stars on banana orbits show a double-peaked distance distribution, so the X-shaped bulge is imprinted in the old stellar population, not only in younger red clump giants.","The retrograde fraction has been roughly stable for several gigayears after bar formation, so today's observations reflect a settled equilibrium rather than a transient state.","A high fraction of retrograde stars near the center naturally explains why centrally concentrated surveys find slower rotation and a rounder bulge.","The simulation predicts that most retrograde stars are not permanently retrograde but oscillate across the inner Lindblad resonance, so a star's current retrograde status is a snapshot rather than a fixed identity."],"supporting_citations":[{"why":"Supplies the photometric metallicities, distances, and the earlier demonstration that metal-rich RR Lyrae stars follow the bar; the dynamical sample is built on this.","marker":"Prudil et al. (2025)"},{"why":"Provides the BRAVA-RR systemic velocities and the rotation-lag measurement that this paper confirms and seeks to explain.","marker":"Kunder et al. (2020)"},{"why":"Supplies the frequency-estimation and chaoticity framework (frequency drift, vILR/ILR maps) used to classify orbits.","marker":"Beraldo e Silva et al. (2023)"},{"why":"Provides the barred analytic Milky Way potential within which all stellar orbits are integrated.","marker":"Portail et al. (2017)"},{"why":"Supplies the bar-model rotation and dispersion curves used as the comparison baseline for the observed kinematics.","marker":"Shen et al. (2010)"},{"why":"Earlier simulation-based work linking kinematic fractionation and bar-built bulges; underpins the secular interpretation.","marker":"Debattista et al. (2017)"},{"why":"Supplies the spatial scaling and alignment procedure used to compare the HG1 simulation with the Milky Way.","marker":"Gough-Kelly et al. (2022)"},{"why":"Previously identified retrograde stellar particles inside 4 kpc of this simulation and tied them to bar evolution.","marker":"Fiteni et al. (2021)"},{"why":"Establishes the bar angle and the X-shaped red clump distance bimodality that the banana-orbit RR Lyrae distribution is compared with.","marker":"Wegg & Gerhard (2013)"}],"fun_headline_variants":["Bar alone builds bulge sphere, no merger needed","Milky Way's round bulge grows from bar, not crash","Spherical bulge traced to bar, not ancient merge","Bar-driven bulge sphere: no classical merger required","Bulge's spherical core from bar, say RR Lyrae orbits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole secular interpretation rests on a single isolated simulation in which a bar forms without any merger; if that simulation's retrograde fraction, orbital stability, or age trends are not representative of the real Milky Way, the central claim weakens even though the observed orbit classifications remain valid.","fun_headline_variants_meta":{"raw":{"variants":["Bar alone builds bulge sphere, no merger needed","Milky Way's round bulge grows from bar, not crash","Spherical bulge traced to bar, not ancient merge","Bar-driven bulge sphere: no classical merger required","Bulge's spherical core from bar, say RR Lyrae orbits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000217,"raw_usage":{"total_tokens":1536,"prompt_tokens":1148,"completion_tokens":388,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":764,"completion_tokens_details":{"reasoning_tokens":309}},"tokens_in":764,"tokens_out":388,"duration_ms":4563,"temperature":1.0,"reasoning_tokens":309,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:38:15.329149+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure [α/Fe] abundances for the retrograde, centrally concentrated bulge RR Lyrae stars. If the secular scenario is right, these stars should be an old, relatively α-rich population present before bar formation; if a large fraction turn out to be young or α-poor, the claim that the spheroidal component is made by bar-driven angular-momentum exchange would be falsified.","supporting_citations":[{"cited_title":"2025, , 695, A211","cited_arxiv_id":null,"evidence_quote":"Supplies the photometric metallicities, distances, and the earlier demonstration that metal-rich RR Lyrae stars follow the bar; the dynamical sample is built on this."},{"cited_title":"2017, , 465, 1621","cited_arxiv_id":null,"evidence_quote":"Provides the barred analytic Milky Way potential within which all stellar orbits are integrated."},{"cited_title":"M., Kormendy , J., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the bar-model rotation and dispersion curves used as the comparison baseline for the observed kinematics."}],"review_version":2}