{"id":"453dac73-f0d1-4009-bb6a-ab82ca5ed82b","arxiv_id":"2411.10391","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Using resonant angles rather than frequency ratios, the authors show that in a mature simulated bar, long-term vertical resonance trapping, not rapid passage, dominates the ongoing vertical growth of the boxy/peanut bulge.","lead":"This paper tracks individual star orbits in a simulated galaxy bar and shows that some orbits get stuck in a vertical resonance while others pass through it, and that the trapped orbits are the main ongoing source of the bar's boxy/peanut shape. It offers a new way to identify resonance in simulations using the resonance angle instead of frequency ratios, which helps settle a debate about how boxy/peanut bulges grow.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dominance claim is underdetermined: the paper shows more recently trapped orbits are librating and stay long, but it never measures the actual vertical-action growth attributable to trapped versus passage orbits; long residence in resonance is not shown to lift orbits.","rationale":"The paper introduces a genuinely useful technique, applying resonant-angle classification to self-consistent N-body orbits for the first time, and it correctly identifies coexistence of trapping and passage at the mature-bar stage. The central claim, however, is not supported by the presented metrics. Occupation fractions (30% trapped vs 21% passage) and residence times speak to how long orbits remain in a state, not to how much vertical action they gain; the paper's own Fig. 3 and the statement that maximum J_z is reached only after leaving the vILR indicate that vertical lifting is associated with leaving or passing through resonance. The reader's weakest assumption about the validity of unperturbed actions is legitimate, but it affects the classification itself; even granting the classification, the dominance inference does not follow without a per-state accounting of ΔJ_z. There are also secondary issues: the abstract's 'half of all bar orbits' overstates the body's statistic, which is restricted to banana-like orbits, and the conclusion rests on a single N-body realization. These are correctable through a reanalysis of the existing simulation data; hence the verdict remains CONDITIONAL, unchanged from the reader's assessment, pending the proposed test.","tokens_in":15521,"tokens_out":6492,"duration_ms":63079,"concrete_test":"Using the saved orbits from the t=300–400 trapped subset (N≈1.4e5), recompute secular J_z at t=400 and t=550. For each orbit, use the time-resolved resonant-angle classification from Appendix A to label every time segment as librating ('in resonance'), 'passage', or circulating. Compute the total growth ΔJ_z over t=400–550 and, for each orbit, partition this growth among the segments according to the time spent in each state (or, more directly, compute the mean secular J_z at entry and exit of each contiguous resonant/passage interval). Sum the growth over all orbits separately for trapped and passage segments. If the total ΔJ_z accumulated during passage segments is comparable to or larger than that accumulated during trapped (librating) segments, the claim that trapping dominates the ongoing B/PS growth fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sections 3.3 and 3.4 report that among orbits trapped into the bar during t=300–400, 30% are in the vILR and 21% in passage at t=500, and that trapped banana orbits reside in resonance for more than 1.9 Gyr. These statistics are used to conclude that long-term vertical trapping dominates the ongoing growth of the boxy/peanut bulge. But growth of a B/PS bulge is measured by the increase of the vertical action J_z, not by the instantaneous occupation fraction or residence time. The paper's own orbit examples (Fig. 3) show that an orbit librating in resonance maintains its average J_z, while an orbit that passes through the resonance increases J_z sharply. The text states that 'leaving the vILR, the orbit enters the circulation mode with the maximum value of J_z' and that the thickest part of the B/PS bulge is built from orbits that have already passed through the vILR and circulate with decreasing θ_res. Thus, long residence in resonance may simply delay the eventual J_z gain rather than enhance it. The conclusion that trapping dominates the ongoing process requires a quantitative attribution of the total ΔJ_z over t=400–550 to periods spent in the trapped (librating) state versus the passage state. Without such an accounting, the relative importance of the two mechanisms is not established. This is a gap in the argument, not a contradiction of known results, and it can be closed by reanalysis of the existing simulation output.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies vertical resonant trapping and resonant heating of orbits in a mature, self-consistent N-body galactic bar, using action-angle variables computed with AGAMA for the axisymmetrized potential. The authors track the resonant angle θ_z−θ_R and classify orbits as circulating, librating (in the vILR), or passing through resonance. They report that at t=500, among orbits trapped into the bar between t=300 and t=400, 30% librate in the vILR, 21% are in passage, and that trapped banana orbits can remain in resonance for more than 1.9 Gyr. From these statistics they conclude that in the later stages of bar evolution vertical trapping dominates the ongoing growth of the boxy/peanut bulge, contradicting some recent works.","tokens_in":15846,"tokens_out":5515,"duration_ms":49760,"significance":"The paper introduces a genuinely new diagnostic to the N-body study of B/PS bulges: the resonant angle as a function of time, as opposed to instantaneous frequency ratios. The qualitative demonstration that trapping and heating coexist at the mature-bar stage, and that individual orbits can librate for several Gyr and then escape, is valuable and goes beyond previous frequency-based analyses. The use of a fully self-consistent model without imposed symmetry is a strength. However, the central quantitative claim that trapping dominates the ongoing growth of the boxy/peanut shape is not supported by the evidence presented; the analysis shows that trapped orbits are numerous and long-lived, but it does not measure their actual contribution to the growth of the vertical action J_z. If the missing attribution analysis is provided, the result would be a significant contribution to the debate on B/PS formation.","major_comments":[{"comment":"The claim that vertical trapping dominates the ongoing growth of the B/PS bulge is not established because the paper never measures the actual growth of the vertical action J_z attributable to orbits in the trapped (librating) versus passage states. In Sec. 3.3 the authors state that while an orbit librates, J_z oscillates around the mean value, and that 'leaving the vILR, the orbit enters the circulation mode with the maximum value of J_z (Fig. 3)'; they also state that 'the thickest part of the B/PS bulge grows due to orbits that have passed through the vILR'. The residence-time statistics in Sec. 3.4 (half of trapped orbits in the vILR for >1.9 Gyr) do not imply a large ΔJ_z during trapping; the librating example in Fig. 3 maintains its average J_z. To support the dominance claim, the authors must compare the total ΔJ_z over t=400–550 contributed by each state, for example by integrating the time-dependent J_z over periods spent librating versus passing, or by comparing ensemble J_z growth for orbit subsets matched in J_v and initial J_z. Without such an accounting, the relative importance of the two mechanisms remains undetermined.","section":"Sec. 3.3 and Conclusions, second bullet"},{"comment":"The abstract states that 'Half of all bar orbits spend more than 2.5 Gyr in vertical resonance over an interval of 4 Gyr', but Sec. 3.4 reports this only for a restricted subset: 'Half of all bar orbits with a resonant angle ∼0(π) spend more than 180 time units in vertical resonance' (i.e., banana-type orbits identified at t=500). The fraction over all bar orbits is not stated and is likely much lower, so the abstract overstates the headline result. In addition, the abstract says 'Half of the orbits trapped into the bar over the last 3 Gyr of simulation remain captured in vertical resonance for more than 2 Gyr', whereas the body says 'half of trapped orbits stay in the vILR for more than 135 time units (1.9 Gyr)'. The abstract should be corrected to match the qualified statements in the body.","section":"Abstract and Sec. 3.4"},{"comment":"The entire classification into libration, circulation, and passage rests on 'unperturbed' actions and frequencies computed for the axisymmetrized potential, which the authors acknowledge in Sec. 1 are 'not really proper actions' in a strongly barred potential. If these quantities are not close to true invariants, the resonant angle θ_z−θ_R may show artificial libration or circulation unrelated to resonance. The paper does not test the robustness of the classification to the action definition—for example, by comparing with actions from a torus-fitting method in the barred potential, or by verifying that the resonant-angle behavior is consistent in a frozen barred potential where the Hamiltonian is time-independent. Because the central claim is a quantitative distinction between trapping and heating, this validation is load-bearing. The authors should add a test of the stability of the computed actions (e.g., conservation of J_z over several libration periods for representative orbits) or at least a discussion of the expected systematic error.","section":"Sec. 2.2 and Appendix A"}],"minor_comments":[{"comment":"The phrase 'contradicts the results of several recent works' is vague; please name the specific works (e.g., Sellwood & Gerhard 2020) and state precisely which of their conclusions are contradicted.","section":"Abstract"},{"comment":"The sentence 'In 100–200 time units, about 2·10^4 of 14·10^4 orbits were added to the B/PS bulge in this manner' is unclear; please specify the exact time interval and state the denominator (all orbits in the bar, or the trapped subset).","section":"Sec. 3.3"},{"comment":"The choices of 'flat orbit' (J_z < 0.05 at t=300) and the trapping interval t=300–400 are not justified or tested for sensitivity; please add a sentence explaining why these thresholds are robust or acknowledge the dependency of the statistics on them.","section":"Sec. 3.2"},{"comment":"In the description of the vertical action and frequency, the paper introduces separate θ_z^max and θ_z^min and then averages them; please clarify why this procedure differs from the standard single θ_z and how the initial phases are set consistently between the apocenter-based and z-maximum-based definitions.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The skeptic's central concern is valid and lands directly on the paper's headline claim: the residence-time statistics do not measure the contribution to vertical-action growth, and the paper's own examples show the opposite association (passage produces the largest ΔJ_z). This is repairable by reanalyzing the existing simulation outputs, so major_revision is appropriate rather than reject. The reliance on the authors' prior work (Zozulia et al. 2024) for the bar-trapping definition is a linkage rather than a circular argument, but it means the paper is not self-contained in that respect; the editor may wish to ensure that the definition is clearly understandable without that reference. The paper fits A&A's scope and the new diagnostic is potentially influential if the missing attribution analysis confirms the claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read it. First, the genuinely new piece is using the resonant angle θ_z−θ_R of individual orbits in a live N-body simulation, instead of frequency ratios, to tell vertical resonance trapping from resonant heating. That is a real methodological advance, and the orbit examples in Fig. 3 make the distinction concrete. Second, the central conclusion—that trapping dominates the ongoing growth of the boxy/peanut bulge in a mature bar—is not actually established by the analysis presented. The paper shows that among orbits trapped into the bar over t=300–400, about 30% are librating in the vILR at t=500, while 21% are in a fast passage. It also shows that some banana orbits stay in resonance for nearly 2 Gyr or longer. But those are occupation statistics, not a measurement of vertical-action growth. The authors never attribute the net increase in J_z over the mature phase to time spent in the trapped state versus the passage state. Their own example orbits show a librating orbit keeping its average J_z nearly constant, while a passing orbit gains J_z sharply. The text even says the thickest part of the bulge is built from orbits that have already left the resonance and are circulating. Long residence in the resonance might delay the eventual lift rather than amplify it. To support \"trapping dominates,\" you need a ΔJ_z budget split by orbital state; that is missing. It is a gap in the argument, not a contradiction, and it can be closed by re-analyzing the existing simulation snapshot data.\n\nTwo smaller points. The abstract's \"half of all bar orbits spend more than 2.5 Gyr in vertical resonance\" drops the qualifier in the body, where it's half of those with resonant angle near 0 or π (banana-type) identified at t=500. And the \"unperturbed\" actions in the barred potential are not proper actions, as the authors themselves note; the time-averaging helps, but it is worth remembering when interpreting quantitative results. Also, one N-body realization is a thin basis for the word \"typical.\"\n\nThe paper is clear, honest, and the method is worth publishing even if the dominance claim needs to be softened or quantified. I'd send it to peer review—it will benefit from a referee who asks for the J_z accounting. If that comes through, the paper could be a useful contribution. As is, I would not take the abstract's dominance statement at face value.","headline":"New method that separates trapping from heating in N-body bars, but the dominance claim is undercut by the lack of a measured ΔJ_z budget.","tokens_in":16395,"tokens_out":4446,"would_cite":true,"duration_ms":38168,"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":"In a mature galactic bar, vertical resonant trapping, not resonant heating, is the dominant ongoing process that grows the boxy/peanut bulge.","keywords":["galactic bar","boxy/peanut bulge","vertical inner Lindblad resonance","resonant trapping","resonant heating","action-angle variables","N-body simulations","orbital classification"],"falsifier":"Take orbits classified as trapped at t=500 and integrate them in a frozen potential built from the actual barred snapshot; if they do not librate around $\\theta_{\\rm res}=0$ or $\\pi$ in that potential, the trapping label is an artifact of using axisymmetric actions.","tokens_in":15326,"feed_emoji":"🌌","tokens_out":8713,"duration_ms":71712,"temperature":0.7,"pith_summary":"This paper asks why disk galaxies develop boxy or peanut-shaped bulges long after the initial buckling of the bar has finished. Using the actions and resonant angles of the actual orbits in a self-consistent N-body simulation, the authors distinguish two late-time mechanisms: resonant trapping, where orbits librate around the vertical inner Lindblad resonance, and resonant heating, where orbits simply pass through it. They report that trapping dominates: at t=500, 30% of recently trapped bar orbits are caught in the resonance versus 21% in passage, and half of all bar orbits spend more than 2.5 Gyr inside the resonance over a 4 Gyr interval. The implication is that real galaxies can keep growing their boxy/peanut bulges purely through resonant capture, without buckling or a frozen-potential artifact.","feed_headline":"Trapping, not heating, drives boxy-peanut bulge growth","feed_subtitle":"Resonant-angle tracking shows 30% of new bar orbits stay trapped versus 21% passing through, locking in for gigayears.","key_machinery":"The machinery is the resonant angle $\\theta_{\\rm res}=\\theta_z-\\theta_R$ built from the vertical and radial angles of each orbit, together with a four-state classification of its behavior: circulation with $\\theta_{\\rm res}$ increasing, circulation with $\\theta_{\\rm res}$ decreasing, libration around 0, $\\pi$ (banana orbits) or $\\pi/2$, $3\\pi/2$ (anti-banana orbits), and passage through the resonance. Unperturbed actions and frequencies are computed for the axisymmetrized potential, time-averaged over oscillation periods, and smoothed with mean-preserving splines to obtain secular values (Appendix A). The resonant angle, rather than the frequency ratio $\\omega_z/\\kappa$, is what lets the authors separate trapping from heating: a ratio near 1 can describe either a librating orbit or one that is merely passing through, while the angle's time history shows which is happening.","core_discovery":"On the paper's own terms, the central discovery is that in the mature, post-buckling phase of a typical bar, the vertical inner Lindblad resonance acts as a persistent trap rather than a transient heater. Tracking the resonant angle $\\theta_{\\rm res}=\\theta_z-\\theta_R$ over t=250–550 in a self-consistent N-body model, the authors classify every orbit into four behaviors: circulation with increasing or decreasing angle, libration in resonance (banana or anti-banana), and passage. Among orbits that joined the bar between t=300 and 400, 30% were librating at t=500 while only 21% were in passage; half of trapped orbits stayed in resonance for more than 1.9 Gyr, and half of all bar orbits logged more than 2.5 Gyr in resonance during the 4 Gyr window. The authors conclude that resonant trapping, not resonant heating, is the dominant ongoing mechanism that lifts disk orbits into the B/PS bulge, in direct contradiction of several recent works.","pith_inferences":["If trapping dominance is generic, the growth rate of a boxy/peanut bulge should track the bar's slowdown rate: a faster-decaying pattern speed sweeps the vILR outward more quickly, converting more passing orbits into trapped ones.","The same resonant-angle method could be turned on other bar resonances (corotation, outer Lindblad) to test whether trapping versus transit is a general dichotomy for secular bar-driven evolution, not just vertical thickening.","A natural stress test would be to recompute actions in a time-dependent barred potential fitted to the simulation snapshots; if the 30%/21% split and multi-gigayear librations survive, the conclusion is robust to the unperturbed-actions caveat, and if not, the distinction may be an artifact of the axisymmetric approximation."],"forward_implications":["After buckling has ended, resonant trapping alone can keep growing the boxy/peanut structure as the bar matures.","Frequency ratios near $\\omega_z/\\kappa = 1$ are ambiguous; the resonant angle's time history is required to tell trapped, passing, and circulating orbits apart.","The thickest part of the B/PS bulge formed early and is barely replenished, while the thinner, extended parts are fed by newly trapped orbits.","Some orbits remain locked in the vILR for more than 4 Gyr even while the bar slows down and the potential changes."],"supporting_citations":[{"why":"The recent study whose conclusion this paper directly contradicts: that buckling, resonant heating, and resonant trapping vary between models, with no single dominant late-stage process.","marker":"Sellwood & Gerhard 2020"},{"why":"Provides the Hamiltonian model of vILR orbital trapping that this paper tests in a self-consistent N-body context.","marker":"Quillen 2002"},{"why":"Formulates resonant heating as the resonance sweeping outward and lifting orbits; the alternative mechanism this paper argues is secondary at late times.","marker":"Quillen et al. 2014"},{"why":"The authors' previous work that defines the action and frequency calculation and the abnormal-orbit bar identification used here.","marker":"Zozulia et al. 2024"},{"why":"Introduces the concept of unperturbed actions for orbits in a barred potential, which the present analysis adopts.","marker":"Binney 2018"},{"why":"Develops resonance and trapped-orbit theory in action-angle space with unperturbed actions, underpinning the resonant-angle approach.","marker":"Binney 2020"},{"why":"Identifies anti-banana orbits as unstable, which the classification uses to expect shorter residence times for them.","marker":"Pfenniger & Friedli 1991"},{"why":"Shows banana orbits appear across the full range of Jacobi integral in self-consistent models, supporting the claim that resonant orbits persist in mature bars.","marker":"Parul et al. 2020"}],"fun_headline_variants":["Trapping trumps heating for peanut bulges","Resonant traps, not heating, form boxy bulges","Mature bar bulges grow by trapping orbits","Action-space view: trapping shapes peanut bulge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole classification rests on treating the time-averaged motions computed from the round, bar-free version of the galaxy's gravity as genuine underlying motions; if the strong bar breaks those approximations, an orbit's labeled trapping or passage could be a numerical artifact rather than real behavior.","fun_headline_variants_meta":{"raw":{"variants":["Trapping trumps heating for peanut bulges","Resonant traps, not heating, form boxy bulges","Mature bar bulges grow by trapping orbits","Action-space view: trapping shapes peanut bulge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000181,"raw_usage":{"total_tokens":1325,"prompt_tokens":983,"completion_tokens":342,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":280}},"tokens_in":599,"tokens_out":342,"duration_ms":3967,"temperature":1.0,"reasoning_tokens":280,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:40:49.003037+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take orbits classified as trapped at t=500 and integrate them in a frozen potential built from the actual barred snapshot; if they do not librate around $\\theta_{\\rm res}=0$ or $\\pi$ in that potential, the trapping label is an artifact of using axisymmetric actions.","supporting_citations":[{"cited_title":"C., Minchev , I., Sharma , S., Qin , Y.-J., & Di Matteo , P","cited_arxiv_id":null,"evidence_quote":"Formulates resonant heating as the resonance sweeping outward and lifting orbits; the alternative mechanism this paper argues is secondary at late times."},{"cited_title":"2018, , 474, 2706","cited_arxiv_id":null,"evidence_quote":"Introduces the concept of unperturbed actions for orbits in a barred potential, which the present analysis adopts."},{"cited_title":"2020, , 495, 886","cited_arxiv_id":null,"evidence_quote":"Develops resonance and trapped-orbit theory in action-angle space with unperturbed actions, underpinning the resonant-angle approach."},{"cited_title":"& Friedli , D","cited_arxiv_id":null,"evidence_quote":"Identifies anti-banana orbits as unstable, which the classification uses to expect shorter residence times for them."}],"review_version":1}