{"id":"c696c6d3-78b6-4a13-903b-536c1ecad4a2","arxiv_id":"2412.00752","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Orbital analysis of 1,879 RR Lyrae and 28,188 APOGEE giant stars separates a bar-aligned inner bulge from a slower-rotating central bulge, supporting a pseudo-bulge origin.","lead":"This paper sorts stars in the Milky Way's bulge by their orbits and finds two distinct populations: an inner bulge that rotates with the galactic bar and a central bulge that does not. It also compares boxy and X-shaped bulge models, reporting a very slight preference for the boxy shape.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The population split by apocentric distance is computed in one fixed barred potential, and the untested assertion that other potentials would not change the result is load-bearing for the central kinematic dichotomy.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the orbital labels that define the central bulge and inner bulge are computed in a fixed, idealized Milky Way potential, and the paper's assertion of insensitivity is untested. I agree with this assessment. The kinematic dichotomy is the paper's central claim, and the entire analysis depends on the classifier: \"central bulge\" and \"inner bulge\" are not independently defined populations but sets of stars selected by apocentric distance in one adopted potential. Because the bar potential is non-axisymmetric and the thresholds lie near dynamical resonances, modest changes in pattern speed or bar angle could move stars across the boundaries. The paper does include a useful robustness check in Section 5 using an alternative time-fraction classification, which mitigates concerns about threshold choice, but it does not address the potential dependence. A secondary issue is the boxy-versus-X-shape comparison: with reduced chi-square values of 0.99 versus 1.00 and probabilities near 50 percent, the reported \"preference for the boxy bulge\" is statistically inconclusive and should be softened. However, that sub-claim is not load-bearing for the central population-split argument, so the main verdict remains CONDITIONAL pending the recommended recomputation.","tokens_in":22096,"tokens_out":4300,"duration_ms":44575,"concrete_test":"Re-integrate the 1,879 RRab and 28,188 APOGEE stars with at least two alternative potentials, e.g., MWPotential2014 with a Dehnen bar pattern speed of 30 and 45 km/s/kpc and bar angle 20 and 30 degrees, plus an axisymmetric McMillan 2017 potential. Recompute r_apo, re-apply the central/inner/interloper classification, and regenerate the Figure 7 and Figure 11 kinematic maps. If the central-bulge versus inner-bulge differences in mean v_los, v_l*, and velocity dispersion persist in sign and approximate amplitude, the concern is settled; if population fractions shift by more than 10 percentage points or the kinematic dichotomy reverses, the central claim is conditional on the adopted potential.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 fixes the Milky Way potential to MWPotential2014 plus a Dehnen bar with pattern speed 52.25 km/s/kpc, bar angle 25 degrees, and bar radius 3.4 kpc. Section 4.1 then divides all stars into central bulge (r_apo < 1.8 kpc), inner bulge (1.8 <= r_apo < 3.5 kpc), and interlopers (r_apo >= 3.5 kpc) using apocenters from orbits integrated in this single potential. The paper states, \"different potentials will not significantly change our main physical results,\" but provides no test. This matters because the thresholds lie in a region where apocenters are sensitive to the adopted bar pattern speed and angle; a moderate change in either can move stars across the 1.8 and 3.5 kpc boundaries, directly reshaping the reported mean v_los and v_l* curves and the velocity dispersions attributed to the central versus inner bulge. The alternative time-fraction classification in Section 5 tests the choice of threshold criterion, not the choice of potential, so the robustness claim remains unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the structure, populations, and kinematics of the Milky Way bulge using 1,879 OGLE-IV RRab stars with Gaia EDR3 proper motions and Kunder et al. (2020) line-of-sight velocities, plus 28,188 APOGEE DR17 red giants and red clump stars with Gaia DR3 proper motions and StarHorse distances. Orbits are integrated in a fixed Milky Way potential (MWPotential2014 plus a Dehnen bar with pattern speed 52.25 km/s/kpc, angle 25 degrees, radius 3.4 kpc), and stars are classified by apocentric distance as central bulge (r_apo < 1.8 kpc), inner bulge (1.8 <= r_apo < 3.5 kpc), or halo/disk interlopers (r_apo >= 3.5 kpc). The central claim is that inner-bulge stars rotate and trace the bar, while central-bulge stars show lower rotation and lower velocity dispersion and do not trace the bar. The paper also argues that metallicity is not a good separator of bulge populations, that orbital classification is preferable, that chemical abundance maps show bimodal distributions in all three populations, and that the observed density is better fit by a boxy bulge model than by an X-shaped model. The conclusions support a secular-evolution, pseudo-bulge origin for the Milky Way bulge.","tokens_in":22327,"tokens_out":7150,"duration_ms":68595,"significance":"If the kinematic dichotomy between the central and inner bulge holds, the paper would strengthen the case that the Milky Way bulge is primarily a pseudo-bulge formed by disk/bar secular evolution, with a distinct old, pressure-supported central component. The use of two independent tracers (RR Lyrae and APOGEE giants) is a genuine strength, and the paper includes useful checks of distance/PLZ systematics (Fig. 2) and an alternative time-fraction classification (Section 5). However, the quantitative support currently lags behind the claims: there are no significance tests for the kinematic differences, no robustness analysis of the apocenter classification to the assumed potential, and no statistical basis for the claimed boxy-versus-X-shape preference. The paper is a reasonable confirmation and extension of Kunder (2022) and Olivares Carvajal et al. (2024), but the novel, load-bearing claims need stronger statistical backing before they can be regarded as established.","major_comments":[{"comment":"The central claim that the central bulge has slower rotation and lower velocity dispersion than the inner bulge is not backed by any significance test. The paper reports binned means and bootstrap error bars, but it never quantifies whether the differences between the blue and orange curves are statistically significant, accounting for the correlated bins and the sample sizes (451 vs 859 RRabs, and 5709 vs 11102 APOGEE stars). Please provide a quantitative test, for example a permutation or bootstrap p-value for the difference in mean v_los or v_l* between the central and inner samples, and an F-test or similar for the dispersion ratio.","section":"§4.1/§4.2, Figs. 7 and 11"},{"comment":"The robustness of the apocenter-based classification to the assumed potential is load-bearing and is not tested. The statement in Section 3 that 'different potentials will not significantly change our main physical results' is asserted without a test, and Section 5 itself attributes the 75% versus 70% bulge-fraction difference between this work and Kunder et al. (2020) to 'differences in the gravitational potential models used for orbital integration.' Because stars near the 1.8 and 3.5 kpc thresholds can be reclassified under moderate changes in bar pattern speed, angle, or strength, the kinematic dichotomy in Figures 7 and 11 could change. Please rerun the classification under at least a few plausible potential variants (or a Monte Carlo over bar parameters) and show that the population labels and the central/inner kinematic differences are stable.","section":"§3 and §5"},{"comment":"The claimed 'preference for the boxy bulge' is not supported by the reported statistics. The reduced chi-square values are 0.99 versus 1.00 for RRabs and 0.99 versus 1.00 for APOGEE stars, with p-values 53% versus 47% and 55% versus 45%, respectively. These are statistically equivalent fits; a 6–8 percentage point difference in p-values does not demonstrate a preference. The authors should either perform a proper model comparison (e.g., delta chi-square with the same bins, AIC/BIC, or bootstrap selection probabilities) or weaken the conclusion to state that both models are consistent with the data and no significant preference is found.","section":"§4.3, Eqs. (13)–(14) and Abstract/Conclusions"},{"comment":"The statement that classification based on orbital parameters rather than metallicity provides a 'more accurate population separation' is not quantified. The paper shows that the metallicity distributions of the three orbital populations overlap (Fig. 9) and that metallicity-selected bulge samples show similar kinematics (Fig. 8), but it never defines or measures classification accuracy. Please provide a quantitative comparison (e.g., contamination rates, separation metrics, or a formal comparison of kinematic homogeneity after each classification) or soften the claim to one of being 'more directly linked to dynamics' rather than 'more accurate.'","section":"§4.1–§4.2 and Conclusions"}],"minor_comments":[{"comment":"With one fitted amplitude, the number of degrees of freedom should be N - 2 rather than N - 1; although the effect is small for large N, the formula should be corrected.","section":"§4.3, Eq. (14)"},{"comment":"The caption says 'central bulge (left), inner bulge (right), and halo/disk interlopers (right)'; the three panels should be labeled '(left), (middle), (right).'","section":"Fig. 10 caption"},{"comment":"There is a typo: 'inner bugle' should be 'inner bulge.'","section":"Conclusions, first paragraph"},{"comment":"Please clarify the Monte Carlo distance error of 0.006 kpc: does it include photometric and reddening uncertainties or only magnitude errors? The later statement that the total distance error does not exceed 6% seems inconsistent with a purely statistical error of 0.006 kpc at roughly 8 kpc.","section":"§2.1"},{"comment":"The sentence beginning 'One such limitation is that RRab stars predominantly trace old and metal-poor stellar populations...' is a sentence fragment that should be joined properly to the following sentence.","section":"§5, first limitation paragraph"}],"recommendation":"major_revision","confidential_remarks":"I found no indication of fabrication or misconduct. The paper is honest about its limitations, but one of those limitations—the fixed gravitational potential used for apocenter classification—is load-bearing for the central kinematic claim and is only hand-waved away in Section 3. The boxy-versus-X-shape comparison should either be properly tested or removed from the abstract and conclusions. These issues are fixable within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The central result — that the inner bulge (1.8 ≤ r_apo < 3.5 kpc) rotates with the bar while the central bulge (r_apo < 1.8 kpc) shows slower rotation and lower dispersion — is visible in the figures and consistent with Kunder (2022) and earlier work. As a confirmation, it's solid. The new value is the joint chemo-dynamical comparison of two independent tracers (1,879 RR Lyrae and 28,188 APOGEE giants) with orbital classification, plus the spatial abundance maps and the alternate time-fraction classification. The distance tests with different PLZ relations are a careful touch.\n\nThe soft spots are real but not fatal. The boxy vs. X-shape 'preference' is a tie: reduced chi-square 0.99 vs. 1.00, probabilities around 53/47 and 55/45. Calling that a preference is an overstatement. Second, the apocenter classification is computed in one fixed potential (MWPotential2014 plus a Dehnen bar with pattern speed 52.25 km/s/kpc, angle 25°, radius 3.4 kpc). The paper asserts other potentials wouldn't change the results, but gives no test. Stars near the 1.8 and 3.5 kpc thresholds could cross boundaries under a different bar speed or angle, which would directly affect the mean velocity and dispersion curves. The time-fraction test in Section 5 tests the classification criterion, not the potential, so the robustness claim remains unsupported. Third, there's no significance test for the rotation/dispersion differences; the figures look convincing, but with overlapping populations a formal test would strengthen the claim. Finally, the pseudo-bulge 'primary origin' conclusion goes beyond the evidence: the inner bulge tracing the bar is good evidence for a disk origin of that component, but it doesn't rule out a classical bulge contribution in the central region.\n\nThis paper is for the Galactic bulge/archaeology community. It deserves a serious referee. The main dichotomy is likely right, but the manuscript needs revision: add significance tests, either test or soften the potential robustness claim, and fix the boxy/X-shape overstatement. I'd send it out.","headline":"Plausible, useful confirmation of the central/inner bulge split, but the boxy-vs-X-shape fit is a statistical tie and the potential robustness claim is untested.","tokens_in":22946,"tokens_out":4154,"would_cite":true,"duration_ms":38055,"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 bulge splits into a slowly rotating central core and a bar-following inner bulge.","keywords":["Galactic bulge","RR Lyrae stars","apocentric distance","orbital dynamics","Galactic bar","boxy bulge","stellar populations","APOGEE"],"falsifier":"Recompute the same orbits with a different bar pattern speed or angle (for example 40 km/s/kpc or 35 degrees), or with spiral arms added to the potential, and test whether stars near the apocenter thresholds swap populations; if the low-rotation central bulge signature vanishes under a plausible alternative potential, the division is an artifact of the assumed model.","tokens_in":21874,"feed_emoji":"🌌","tokens_out":7086,"duration_ms":60337,"temperature":0.7,"pith_summary":"The paper aims to establish that the Milky Way bulge is made of two dynamically distinct populations, separated not by metallicity but by how far their orbits reach. Using 1,879 ab-type RR Lyrae stars from OGLE-IV and 28,188 red giant stars from APOGEE, it classifies stars with apocenter below 1.8 kpc as central bulge, those between 1.8 and 3.5 kpc as inner bulge, and the rest as halo or disk interlopers. The inner bulge rotates with the Galactic bar while the central bulge shows slower rotation and lower velocity dispersion, and the same split appears in both old and young tracer populations. The paper further claims that orbital classification is more physically meaningful than metallicity cuts, and that the bulge density fits a boxy model better than an X-shaped one. These results matter because they say the bulge's dominant component is a bar-driven pseudo-bulge formed by secular disk evolution.","feed_headline":"Bulge splits into two stellar populations","feed_subtitle":"Orbits of 30,000 stars show the inner bulge follows the bar while the central bulge barely rotates.","key_machinery":"The load-bearing object is the apocentric distance $r_{\\rm apo}$ obtained from backward 5 Gyr orbital integrations in the MWPotential2014 gravitational potential, augmented by a Dehnen bar with pattern speed 52.25 km/s/kpc, angle 25 degrees, and radius 3.4 kpc. Apocenter thresholds at 1.8 kpc and 3.5 kpc divide stars into central bulge, inner bulge, and interlopers; this one-dimensional division is what creates the two populations whose rotation, dispersion, chemistry, and density are then compared. The thresholds come from prior work the paper adapts, and the alternative residence-time classification in the discussion gives the same kinematic result.","core_discovery":"On the paper's own terms, the central discovery is that the bulge should be divided into a central bulge ($r_{\\rm apo}<1.8$ kpc) and an inner bulge ($1.8 \\le r_{\\rm apo}<3.5$ kpc), with everything beyond classified as halo or disk contamination. Both the ancient RR Lyrae population and the APOGEE red giants show the same pattern: inner bulge stars rotate coherently with the Galactic bar and display the quadrupole velocity signature of barred orbits, while central bulge stars rotate slowly, have lower velocity dispersion, and show no bar alignment. Apocenter-based classification also reveals that metal-poor and metal-rich bulge stars have nearly identical kinematics, so metallicity is not the right divider of bulge populations. The paper concludes that the inner bulge's bar-like, disk-like orbits support secular evolution as the primary origin of the bulge, and that the observed density of both tracers is better matched by a boxy bulge model than by an X-shaped one.","pith_inferences":["Our inference: if the central bulge's low rotation and low dispersion are real, it may be an older spheroidal component that predates or coexists with the bar, not merely the bar's center; direct age dating of stars in each apocenter class could test this.","Our inference: the boxy-versus-X-shape comparison fixes all model parameters except amplitude, so the preference could shift if bar angle or axis ratios were also fitted; the boxy conclusion should be read as provisional.","Our inference: because interlopers dominate the high-dispersion tail, bulge samples lacking full 3D velocities will systematically overestimate velocity dispersion; applying the same cleaning to future catalogs should sharpen the kinematic split."],"forward_implications":["Bulge surveys that mix central and inner stars will dilute bar rotation signals, so apocenter-based cleaning should be applied before measuring bulge kinematics.","The inner bulge's majority share and bar-like, disk-like orbits support a pseudo-bulge formed by disk instability rather than a merger-built classical bulge.","Metallicity alone cannot define bulge populations; orbital parameters separate them more accurately.","Both ancient RR Lyrae stars and younger giants favor a boxy/peanut bulge over an X-shaped bulge in the adopted density fits."],"supporting_citations":[{"why":"Supplies the line-of-sight velocities for the RR Lyrae sample and the 3.5 kpc apocenter threshold used to separate bulge stars from halo interlopers.","marker":"Kunder et al. (2020)"},{"why":"Supplies the 1.8 kpc apocenter boundary separating central from inner bulge RR Lyrae stars.","marker":"Kunder (2022)"},{"why":"Provides the OGLE-IV RRab cleaning procedure, the period-luminosity-metallicity distances, and the metallicity calibration used here.","marker":"Pietrukowicz et al. (2015)"},{"why":"Provides the MWPotential2014 Milky Way potential and the coordinate and orbit machinery used for the integrations.","marker":"Bovy (2015)"},{"why":"Defines the bar potential whose pattern speed, angle, and radius are adopted for the orbit integrations.","marker":"Dehnen (2000)"},{"why":"Supplies the APOGEE sample's kinematic and chemical framework that this paper extends with apocenter classification.","marker":"Queiroz et al. (2021)"},{"why":"Provides the boxy bulge density model used in the model comparison.","marker":"López-Corredoira et al. (2005)"},{"why":"Provides the X-shaped bulge density model and its parameters used in the comparison.","marker":"Wegg & Gerhard (2013)"},{"why":"N-body simulation showing a bar and pseudo-bulge formed by disk instability, cited as the theoretical basis for the secular-evolution conclusion.","marker":"Shen et al. (2010)"},{"why":"Provides an alternative residence-time orbital classification used as a robustness check.","marker":"Olivares Carvajal et al. (2024)"}],"fun_headline_variants":["Orbits reveal inner bulge rotates with bar, center lags","Bulge classified by orbits, not metals, new study finds","Boxy bulge wins: new data on 30k stars","Three bulge populations, one bar-aligned core","Apocenter, not iron, splits bulge stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The two-population picture rests on apocentric distances computed in one fixed, idealized Milky Way potential with a rigid bar of assumed pattern speed, angle, and radius; if the true bar differs, stars near the 1.8 kpc and 3.5 kpc boundaries would be reassigned and the kinematic dichotomy could weaken or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Orbits reveal inner bulge rotates with bar, center lags","Bulge classified by orbits, not metals, new study finds","Boxy bulge wins: new data on 30k stars","Three bulge populations, one bar-aligned core","Apocenter, not iron, splits bulge stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000212,"raw_usage":{"total_tokens":1481,"prompt_tokens":1070,"completion_tokens":411,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":686,"completion_tokens_details":{"reasoning_tokens":330}},"tokens_in":686,"tokens_out":411,"duration_ms":4623,"temperature":1.0,"reasoning_tokens":330,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:02:29.426883+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the same orbits with a different bar pattern speed or angle (for example 40 km/s/kpc or 35 degrees), or with spiral arms added to the potential, and test whether stars near the apocenter thresholds swap populations; if the low-rotation central bulge signature vanishes under a plausible alternative potential, the division is an artifact of the assumed model.","supporting_citations":[{"cited_title":"M., et al.\\ 2020, , 159, 270","cited_arxiv_id":null,"evidence_quote":"Supplies the line-of-sight velocities for the RR Lyrae sample and the 3.5 kpc apocenter threshold used to separate bulge stars from halo interlopers."},{"cited_title":"M.\\ 2022, Universe, 8, 206","cited_arxiv_id":null,"evidence_quote":"Supplies the 1.8 kpc apocenter boundary separating central from inner bulge RR Lyrae stars."},{"cited_title":"Comparing bulge RR Lyrae stars with bulge giants -- Insight from 3D kinematics","cited_arxiv_id":"2405.08990","evidence_quote":"Provides an alternative residence-time orbital classification used as a robustness check."}],"review_version":1}