{"id":"59b80c42-d9b0-461f-9339-356237829054","arxiv_id":"2507.20790","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A bulgeless galaxy simulation with a dense Sérsic disc and 7% gas forms a bar and box/peanut bulge whose size, pattern speed, and line-of-sight structure match the Milky Way.","lead":"This paper shows that a Milky Way-like bar and boxy bulge can form in a computer galaxy that starts with no central bulge, if the initial disc is dense in the center and contains a small amount of gas. The best model matches the bar's size, speed, and shape, supporting the idea that the Milky Way's bulge grew from its disc rather than from mergers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The TG07 match to the Milky Way rests on a metastable bar state that the authors' own Appendix C flags as likely fragile to substructure and gas accretion; the claim that the model reproduces the MW bar and BP bulge therefore requires a perturbation test the paper defers.","rationale":"Reading the paper in good faith, the central claim has two layers. The existence layer — a bulgeless Sérsic disc with a 7% gas disc can evolve a bar and BP bulge whose properties resemble the MW — is well supported: three independent TG07 variants with different random seeds and star-formation recipes converge on the same steady bar from ~3.5 to 13 Gyr (§4.2.1), and the D-series shows a monotonic gas-fraction trend consistent with Athanassoula et al. (2013) and Beane et al. (2023). The representational layer — TG07 is a plausible model of the present-day MW — is the least secure. The match exists only while the bar sits in the Sellwood & Debattista (2006) metastable state, and the authors' Appendix C concedes that this state is 'quite fragile, and liable to be broken by external perturbations, such as would arise in a fully cosmological setting,' and that 'confirmation needs further simulations with substructure taken into account.' The actual MW is perturbed (Sagittarius, Magellanic Clouds, accretion); on the authors' own account those perturbations would eject the bar from the metastable state and resume the secular growth/slowdown that produces bars inconsistent with the MW (the D00/TG00 runaway regime). Separately, at the selected best snapshot (Table 2) the match is genuine but partial: Rbar and Ωp agree at ~1–1.4σ, the slow-down rate is ~2.5σ high, and the BP line-of-sight distributions miss the central dips at l = −8.5° (§4.2.2). I agree with the reader's direction but sharpen it: the key test is not snapshot selection per se (Fig. 6 shows Rbar and Ωp are steady from 3.5–13 Gyr, so the reader's 'about 1 Gyr' understates the persistence of the individual bar properties) but survival of the metastable state under realistic perturbations. If the proposed perturbation run breaks the state within ~1 Gyr, the abstract's 'reproduces' overclaims; if not, the concern is retired. Because the existence claim is solid and the paper is transparent, the CONDITIONAL verdict stands.","tokens_in":22746,"tokens_out":22268,"duration_ms":239442,"concrete_test":"Evolve TG07 from its t = 3 Gyr matched snapshot (and from the original ICs) with a realistic perturbation history: a Sagittarius-like satellite (≈10^9 M☉ on a decaying orbit) and/or ongoing gas accretion, tracking Rbar, Ωp, and the BP line-of-sight bimodality for 5–10 Gyr. If the bar resumes secular growth and slowdown within ~1 Gyr — as Sellwood & Debattista (2006) observed for their metastable bars — the MW match is an artifact of the isolated-galaxy idealization and the abstract's 'reproduces' overclaims. If the matched properties (Rbar within 4.6±0.3 kpc, Ωp within 39±3.5 km/s/kpc, BP bimodality) persist for several Gyr under perturbation, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (§5) is that TG07 shows it is 'possible to build a plausible MW bar and BP bulge using bulgeless ICs, provided that a gas disc is present.' The least secure condition is not bar formation but the persistence of the match: TG07 matches the MW's bar radius, pattern speed, and BP shape only while its bar occupies the metastable state identified with Sellwood & Debattista (2006). The authors' own Appendix C states that this state is 'quite fragile, and liable to be broken by external perturbations, such as would arise in a fully cosmological setting,' and that 'confirmation needs further simulations with substructure taken into account.' The real Milky Way is not isolated: the Sagittarius dwarf has crossed the disc repeatedly, the Magellanic Clouds are massive perturbers, and gas inflow is ongoing. On the authors' account, such perturbations would eject the bar from the metastable state, resuming the secular growth and slowdown seen by 13 Gyr in Fig. C1 — the 'too long, too slow' regime the paper argues is incompatible with the MW. The three TG07 variants (§4.2.1) establish numerical robustness to seeds and feedback, not physical robustness to environment. A secondary aggravator: the matching snapshot is selected (Table 2, §4.2.2), and at that snapshot the slow-down rate is −8 km/s/kpc/Gyr versus the MW's −4.5±1.4; the dismissal of Chiba et al. (2021) as assuming linear Ωp(t) is an argument about the model, not a refutation of the MW constraint. The match is genuine but partial, and its relevance to the present-day MW is untested.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents isolated galaxy simulations built with a modified GalactICS that generates bulgeless initial conditions with Sérsic stellar discs and optional gas discs. A D-series varies gas fraction from 0% to 30% and a TG-series is based on the Tepper-Garcia et al. (2021) Milky Way model. The authors report that gas suppresses the secular growth of the bar and that their best model, TG07, at a selected t=3 Gyr snapshot roughly matches the Milky Way's bar length, pattern speed, and box/peanut bulge line-of-sight density distributions. The central claim in Section 5 is that it is possible to build a plausible Milky Way bar and BP bulge from bulgeless initial conditions provided a gas disc is present.","tokens_in":23098,"tokens_out":6457,"duration_ms":78202,"significance":"If the central claim holds, this is a useful proof-of-possibility result: it demonstrates that a classical bulge is not required to reproduce several Milky Way bar/bulge observables. The paper has clear strengths: the modified GalactICS code is publicly available, the three TG07 variants test robustness to random seeds and feedback implementation, and the mock red clump/RGBB comparisons against Gonzalez et al. (2015) are more independent than the bar-length and pattern-speed comparisons. The work also includes unusually explicit caveats about snapshot matching and metastability in Appendices B and C. The main value is as a proof-of-concept and as a tool release rather than as a unique or fully constrained model of the Milky Way.","major_comments":[{"comment":"The stress-test concern about the metastable state lands. The TG07 match to the Milky Way is evaluated at t=3 Gyr, a snapshot selected because it best satisfies the three bar measurements in Table 2, and at that time the bar is in the metastable state discussed in Appendix C. Appendix C explicitly states that this state is 'quite fragile, and liable to be broken by external perturbations, such as would arise in a fully cosmological setting,' and Figure C1 shows that by 13 Gyr all three TG07 models leave the state, with TG07v2 and TG07v3 growing rapidly and TG07 becoming a slow rotator. Because the real Milky Way has experienced satellite impacts and ongoing gas inflow, the Section 5 claim requires either a direct perturbation test (for instance, a delayed satellite flyby or a subhalo encounter in the same numerical setup) or a clear scope restriction to isolated evolution. Without one of these, the conclusion is supported only by an unperturbed, selected snapshot.","section":"§4.2.2, Table 2, Appendix C, Figure C1"},{"comment":"The treatment of the Chiba et al. (2021) slowdown constraint is not fully convincing. At the selected t=3 Gyr snapshot, Table 2 lists Omega_dot_p = -8 km/s/kpc/Gyr for TG07, which is about 2.5 sigma from the Milky Way value of -4.5 +/- 1.4. Section 5 dismisses this constraint because Chiba et al. assume a linear Omega_p(t), but the observable is a local deceleration estimate, and the paper does not demonstrate that the TG07 nonlinearity makes the comparison invalid. The authors should either apply the Chiba et al. measurement methodology to the simulation to compute what that method would infer from TG07, or explicitly list the slow-down rate as an unmatched property in Section 5.","section":"Table 2, §4.2, §5"},{"comment":"The rotational-curve and surface-density comparisons are shown only at t=0 and t=10 Gyr, even though the Milky Way match is claimed at t=3 Gyr. At t=10 Gyr the TG07 rotation curve is 15-20 km/s lower than the Eilers et al. (2019) data between 5 and 10 kpc and is still rising where the observed curve falls, and the surface density profile in Figure B1 shows a double inflection not present in the Bovy & Rix (2013) data. Because the t=3 Gyr snapshot is the one claimed to represent the Milky Way, these diagnostics should be shown at t=3 Gyr as well; otherwise the 'plausible Milky Way' statement should be explicitly limited to the bar and BP-bulge observables.","section":"Appendix B, Figure B1"},{"comment":"There is a degree of fit-driven agreement that should be quantified. The Sersic index n=1.75 and the 7% gas fraction were chosen after inspecting the D-series outcomes, and the snapshot in Table 2 is selected using the same bar measurements against which the model is judged. This is not fatal for a proof-of-possibility, but the paper should state how sensitive the match is to these choices, for example by indicating how far n and the gas fraction can be varied before the bar length or pattern speed leaves the observed ranges. The three TG07 variants establish robustness to seeds and feedback but not to the initial-condition parameters themselves.","section":"§3.2, §4.1, Table 2"}],"minor_comments":[{"comment":"The caption says the TG sequence is shown 'at t = 10 Gyr,' but the figure and the associated text in Section 4.2.2 compare snapshots at t=3, 5, and 10 Gyr; the caption should be corrected.","section":"Figure 8 caption"},{"comment":"The text refers to 'model T00' where the context indicates model TG00; please correct the typo.","section":"§4.2.1"},{"comment":"The column header 'TG07 - 5 Gy' is truncated; it should read 'TG07 - 5 Gyr'.","section":"Table B1"},{"comment":"The alternate normalization used in Figure D4 is important for interpreting the shapes of the line-of-sight distributions, but readers meet it only in an appendix; a brief explanation in Section 4.2.2 would make the main comparison easier to follow.","section":"Appendix D, Figure D4"}],"recommendation":"major_revision","confidential_remarks":"The paper's own appendices are unusually candid about the metastable-state fragility and the snapshot-selection issue, so the required revision path is feasible. I would ask for either a perturbation experiment or a clearly scoped claim, rather than treating the isolated, selected-snapshot result as a complete Milky Way model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read this one carefully. The core result matches the abstract: a bulgeless Sérsic disc plus a 7% gas disc forms a bar and BP bulge whose size, pattern speed, and line-of-sight stellar density distributions are close to the Milky Way's at a 3 Gyr snapshot. That result is new in this specific combination, and the Sérsic disc implementation in GalactICS is the real hard artifact—it's public on GitHub and straightforwardly reusable. The simulation suite is well designed: the D-series varies gas fraction 0–30% and shows gas suppresses bar growth, and the three TG07 variants show the outcome is robust to random seeds and feedback details. The mock observations against Gonzalez et al. (2015) are also done carefully, with convolution by red-clump/RGBB magnitude distributions and shape comparisons rather than just surface density overlays.\n\nThe soft spots are mostly ones the paper itself flags. The matched snapshot is selected, and the bar only sits near the MW constraints for about 1 Gyr; by 13 Gyr it has left the metastable state and resumed slowdown and growth. Appendix C states plainly that the metastable state is fragile and likely to be broken by perturbations 'such as would arise in a fully cosmological setting.' That matters because the Milky Way is not isolated: Sagittarius, the Magellanic Clouds, and ongoing gas inflow would, on the authors' own account, knock the bar out of the state that produces the match. So this is a demonstration of a plausible formation path in isolation, not a demonstrated match to the present-day Galaxy. The slowdown rate at the matched snapshot is about -8 km/s/kpc/Gyr versus the MW's -4.5±1.4; the paper's dismissal of the Chiba et al. (2021) measurement as assuming a linear pattern speed is a model-side argument, not a refutation of the data point. Some of the agreement is also fit-driven: the Sérsic index and gas fraction are tuned to the target observables. The less tunable comparison, the BP shape along the line-of-sight distributions, is the most convincing piece.\n\nNone of this sinks the paper. The authors are transparent about the tuning and the fragility, and the central claim is modest: it is possible to build a plausible MW bar and BP bulge from bulgeless ICs with gas. I think that claim holds. What is untested is whether such a state survives in the real, clumpy, accreting Milky Way.\n\nThis paper is for people working on bar formation, secular bulge evolution, and Galactic archaeology. It deserves a proper referee. The referee should push for a perturbation test—add a satellite or clumpy halo to TG07 and see if the metastable state breaks—or at least a stronger caveat about the present-day applicability, and a more careful handling of the Chiba slowdown constraint. Those are revision-level asks, not reasons to desk reject.","headline":"A bulgeless Sérsic disc with 7% gas gives a plausible MW bar and BP bulge, but the match sits on a metastable state the paper itself expects perturbers to break.","tokens_in":23670,"tokens_out":3751,"would_cite":true,"duration_ms":40564,"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":"An initially bulgeless galaxy with a Sérsic stellar disc and 7% gas can grow a bar and box/peanut bulge matching the Milky Way's observed bulge+bar.","keywords":["Milky Way bar","box/peanut bulge","Sérsic disc","bulgeless galaxy models","bar pattern speed","gas-regulated bar growth","N-body simulation","Galactic dynamics"],"falsifier":"Measure the Milky Way bar's corotation ratio and pattern-speed history with Gaia-era data: the TG07 scenario requires a fast bar ($R_{\\rm CR}/R_{\\rm bar} < 1.4$) whose speed is roughly flat over the last ~1 Gyr, so finding a slow bar or steady multi-gigayear deceleration rules the model out. The numerical complement is to rerun the TG07 initial conditions with a satellite or subhalo perturber, which the paper's Appendix C identifies as the likely route out of metastability; if the bar then resumes lengthening and spindown within a few gigayears, the isolated simulation overstates the durability of the matched state.","tokens_in":22550,"feed_emoji":"🌌","tokens_out":15532,"duration_ms":148288,"temperature":0.7,"pith_summary":"The paper claims that a classical bulge is not needed to explain the Milky Way's boxy bulge: a galaxy that starts bulgeless, with a centrally concentrated Sérsic stellar disc and a modest gas disc, can evolve a bar and box/peanut (BP) bulge whose length, pattern speed, and shape match the observed bulge+bar. The authors build equilibrium initial conditions with a modified version of the GalactICS code that replaces the usual exponential disc with a Sérsic profile (best fit $n = 1.75$) and run isolated N-body plus hydrodynamics simulations with gas fractions from 0% to 30%. Gasless models develop runaway bars that grow to 8-10 kpc and spin down far below the Milky Way value, while adding a gas disc of only 7% of the stellar mass produces a metastable bar that stays fast and short for several gigayears; three independent 7%-gas runs reach the same configuration, suggesting the result does not depend on the numerical recipe. If the claim is right, the Milky Way's bulge formed secularly from the disc under the bar's influence, and gas, rather than a classical bulge or a specially tuned halo, is the main regulator of bar growth.","feed_headline":"A bulgeless disc builds the Milky Way's bar and boxy bulge","feed_subtitle":"The match implies the Galaxy's boxy bulge formed from the disc alone, with no classical bulge.","key_machinery":"The central object is the Sérsic stellar disc, with surface density $\\Sigma(R) = \\Sigma_0 e^{-(R/R_d)^{1/n}}$, which for $n > 1$ is denser in the centre than an exponential disc and thereby mimics the excess central light that other Milky Way models attribute to a classical bulge. The mechanism that carries the argument is gas-regulated metastability: the 7% gas disc is funnelled inward by the bar, briefly raising the pattern speed and trapping the bar in a metastable state (in the sense of Sellwood & Debattista 2006) in which the resonances face a rising phase-space density, suppressing angular-momentum loss to the halo. This keeps the bar short and fast for several gigayears, while gasless discs shed angular momentum steadily and grow bars of 8-10 kpc. The enabling tool is a modified GalactICS code that constructs equilibrium Sérsic discs with a double-exponential radial velocity dispersion, which is what allows the bulgeless initial conditions to exist at all.","core_discovery":"The paper's central claim, stated in its conclusions, is that it is indeed possible to build a plausible Milky Way bar and box/peanut bulge from bulgeless initial conditions, provided a gas disc is present. The best-fit model, TG07, starts from a Sérsic stellar disc with index $n = 1.75$ and no classical bulge, a gas disc with 7% of the stellar mass, and a dark halo; at $t = 3$ Gyr it reproduces the Milky Way's bar length ($\\sim 4.25$ kpc), pattern speed ($\\sim 44$ km s$^{-1}$ kpc$^{-1}$), and the double-peaked line-of-sight number counts that trace the X-shaped BP bulge in the Gonzalez et al. (2015) data. The gas is the load-bearing ingredient: the gasless equivalents grow bars of 8-10 kpc with pattern speeds of roughly 20-35 km s$^{-1}$ kpc$^{-1}$, far outside the observed range, while three different 7%-gas realizations converge on the same steady configuration. The authors interpret this as a metastable state in which gas funneled to the centre briefly raises the pattern speed, trapping the bar against secular slowdown, and they emphasise that the match to the Milky Way lasts only about 1 Gyr before the model evolves away.","pith_inferences":["An implication the paper leaves implicit is that real galaxies with satellite encounters would likely be knocked out of the metastable state, so for the Milky Way the match would need either recent re-entry into metastability or an additional stabilising agent; re-running these initial conditions with live substructure would settle this.","The same construction predicts a population-level trend the paper does not state: gas-rich, bulgeless barred galaxies should preferentially host short, fast bars, which is testable with edge-on samples of barred galaxies.","The paper's unconvolved line-of-sight distributions show a sharper double-peaked bulge signature than the magnitude-convolved data reveal, so surveys with individual stellar distances could recover that intrinsic bimodality and test the BP shape more strictly than the histograms used here."],"forward_implications":["The Milky Way's bulge can be entirely secular: a classical bulge is not required to produce the observed box/peanut shape, so the Galaxy's central light excess may come from the disc itself.","A gas disc of only a few percent of the stellar mass can hold a bar in a fast, metastable configuration for several gigayears, so reproducing the Milky Way's bar length and pattern speed does not require a specially tuned halo.","Because a simulation can match the Milky Way's bar for only about 1 Gyr before evolving away, models should be compared through snapshot matching, and the time a model takes to reach an observed configuration could be used to date Galactic structures.","The 7%-gas models host fast bars ($R_{\\rm CR}/R_{\\rm bar} < 1.4$) for nearly their entire evolution, implying that if this scenario describes the Galaxy, the Milky Way's bar should likewise be fast."],"supporting_citations":[{"why":"Supplies the metastability mechanism invoked to explain why the 7%-gas bars stay fast and short: a brief rise in pattern speed traps the bar against slowdown.","marker":"Sellwood & Debattista 2006"},{"why":"The Milky Way-tuned model whose bulge-plus-disc parameters the TG-series initial conditions replace with a Sérsic disc; provides the rotation-curve and surface-density baseline.","marker":"Tepper-Garcia et al. 2021"},{"why":"The Milky Way model forming the basis of the D-series initial conditions, whose classical bulge is swapped for a Sérsic disc.","marker":"D'Onghia & Aguerri 2020"},{"why":"The line-of-sight red-clump and RGBB number-count data used as the observational shape test for the simulated BP bulge.","marker":"Gonzalez et al. 2015"},{"why":"Provides the Milky Way bar pattern speed ($39.0 \\pm 3.5$ km s$^{-1}$ kpc$^{-1}$) that the models must reproduce.","marker":"Portail et al. 2017"},{"why":"Provides the Milky Way bar length ($4.6 \\pm 0.3$ kpc) used as the size target for the simulations.","marker":"Wegg et al. 2015"},{"why":"Provides the bar slowdown-rate estimate that defines the ~1 Gyr window in which any simulation can match the Milky Way bar.","marker":"Chiba et al. 2021"},{"why":"The version of GalactICS with gas discs and double-exponential velocity dispersions that this paper extends to Sérsic discs.","marker":"Deg et al. 2019"},{"why":"Observational evidence that gas-rich galaxies host bars that do not slow down, cited as motivation for gas-regulated bar evolution.","marker":"Beane et al. 2023"}],"fun_headline_variants":["Gas-rich disc builds the Milky Way's bar and boxy bulge","Bulgeless disc with gas matches Milky Way's bulge","Gas is the load-bearing ingredient for the galaxy's bar","Milky Way's boxy bulge from a gas-rich disc alone","Gas turns a plain disc into the Milky Way's bulge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Milky Way's bar is right now sitting in the same short-lived metastable state that the simulated bar occupies at $t = 3$ Gyr, and that the real bar has not been perturbed out of that state; the paper's own runs show the state breaks down by 13 Gyr and that the model drifts from the observed rotation curve and surface density at late times, so if the real bar has been significantly disturbed, the resemblance is a snapshot coincidence rather than an explanation.","fun_headline_variants_meta":{"raw":{"variants":["Gas-rich disc builds the Milky Way's bar and boxy bulge","Bulgeless disc with gas matches Milky Way's bulge","Gas is the load-bearing ingredient for the galaxy's bar","Milky Way's boxy bulge from a gas-rich disc alone","Gas turns a plain disc into the Milky Way's bulge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000565,"raw_usage":{"total_tokens":2703,"prompt_tokens":991,"completion_tokens":1712,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":607,"completion_tokens_details":{"reasoning_tokens":1637}},"tokens_in":607,"tokens_out":1712,"duration_ms":19040,"temperature":1.0,"reasoning_tokens":1637,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T13:15:03.941939+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Milky Way bar's corotation ratio and pattern-speed history with Gaia-era data: the TG07 scenario requires a fast bar ($R_{\\rm CR}/R_{\\rm bar} < 1.4$) whose speed is roughly flat over the last ~1 Gyr, so finding a slow bar or steady multi-gigayear deceleration rules the model out. The numerical complement is to rerun the TG07 initial conditions with a satellite or subhalo perturber, which the paper's Appendix C identifies as the likely route out of metastability; if the bar then resumes lengthening and spindown within a few gigayears, the isolated simulation overstates the durability of the matched state.","supporting_citations":[{"cited_title":"A., Zoccali M., Debattista V","cited_arxiv_id":null,"evidence_quote":"The line-of-sight red-clump and RGBB number-count data used as the observational shape test for the simulated BP bulge."}],"review_version":1}