{"id":"90796739-e03b-41c7-9393-bd55db8177f6","arxiv_id":"2505.08499","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Following binaries through up to three encounters with SgrA* boosts disruption fractions by roughly 20 percent and, for an example system, makes mergers about 31 percent of outcomes.","lead":"This paper tracks what happens to binary stars that pass close to the Milky Way's central black hole, following them through up to three encounters instead of just one. It finds that later encounters significantly increase the number of binaries torn apart and that a large share of systems merge.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eccentric-binary phase sampling, not CM-frame neglect, is the soft spot in the multi-passage boost claim.","rationale":"The paper's central quantitative claim is that multiple pericentre passages boost disruptions by 20% or more, and the mechanism is the evolution of eccentric, more prograde CBs into deeper encounters. The most delicate step is therefore how the initial conditions for the second and third passages are drawn. The reader identified the CM-frame orientation update (footnote 4) as the weakest assumption. That is a real acknowledged gap, but a magnitude estimate shows it is almost certainly subdominant: the change in the CM angular momentum vector is at most of order L_b/L_cm ~ Q^{-2/3} ~ 10^{-4}, so the missing rotation of the CM plane and periapse direction changes the next-passage geometry at the 10^{-4} level. By contrast, the phase distribution for returning eccentric binaries is a first-order statistical choice. The paper explicitly randomizes the binary phase for CBs, and the natural reading of 'random phases' is a uniform draw in true anomaly, the variable used in Eq. (8). For an eccentric Keplerian orbit, a random time corresponds to uniform mean anomaly, not uniform true anomaly; the two differ by factors of several near pericentre for e_b ~ 0.6-0.8. Since the tidal outcome is sensitive to the binary phase at encounter, this can bias the very fractions that support the boost claim. The proposed test isolates this effect cleanly by re-running the same recursion with the physically correct phase distribution while leaving everything else fixed. If the boost is robust to this change, the central claim stands; if not, the paper needs to specify and justify the phase sampling. I therefore keep the reader's CONDITIONAL verdict unchanged rather than escalating it, because the concern is concrete, testable, and already within the scope of the reader's call for documentation of the subsequent-passage initial conditions.","tokens_in":30027,"tokens_out":16923,"duration_ms":185274,"concrete_test":"Re-run the Sec. 4 three-passage Monte Carlo for the example system, drawing the initial binary phase of every CB from the physical time-of-flight distribution: sample M uniformly in [0,2pi), solve Kepler's equation M = E - e_b sin E, then phi = 2 atan2(sqrt(1+e_b) sin(E/2), sqrt(1-e_b) cos(E/2)), with e_b taken from the end of the previous passage. Keep the CM update (Eqs. 37-44) and all cuts unchanged. Compare the cumulative fractions D3, F3, CB3, M3 and the boost D3/D1 - 1 in Fig. 10 to the uniform-phi run. If the boost changes by more than ~3 percentage points, or if the '20% or more' statement fails for beta_0 > 1, the phase sampling is a load-bearing flaw.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The multi-passage iteration in Sec. 4.1 carries forward every binary element except the phase, which is redrawn 'randomly' for each CB. Because Eq. (8) uses the true anomaly phi, a uniform draw over [0,2pi) is not the phase distribution of an eccentric binary at a random time: the physical distribution is uniform in mean anomaly, giving dN/dphi proportional to (1-e^2)^{3/2}/(1+e cos phi)^2. For the eccentric CBs produced by the first passage (e_b up to ~0.8), uniform phi overweights pericentre phases by a factor ~5 and underweights apocentre phases relative to the time-weighted ensemble. The second- and third-passage D/F/C/M fractions, and hence the headline '20% or more' disruption boost, are averaged over these phases. The paper nowhere states that the redrawn phase is sampled from the time-weighted distribution, and the same uniform-phi convention is used throughout. This is an O(1) bias in the initial conditions of the very passages that constitute the paper's central new result, not a small correction. The CM-frame omission flagged in footnote 4 is likely much less important, since Delta L_b/L_cm ~ Q^{-2/3} ~ 10^{-4}, so the missing reorientation changes the geometry at the 10^{-4} level.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the restricted three-body treatment of the Hills mechanism to multiple pericentre passages around SgrA*. Starting from initially circular binaries on parabolic CM orbits, it integrates the linearized equations of motion to classify outcomes as disruptions (Ds), fly-away binaries (FAs), coming-back binaries (CBs), and, for a chosen example system, mergers (Ms). The central claims are that multiple encounters boost the disruption fraction by 20% or more relative to a single passage, that the example system produces mergers 31% of the time, and that these effects depend strongly on inclination and diving factor. The paper presents phase-space maps of outcomes, distributions of ejected and captured star properties, and the effect of finite stellar lifetimes and sizes.","tokens_in":30463,"tokens_out":7174,"duration_ms":78056,"significance":"If the central claim holds, single-passage treatments systematically underestimate the disruption and merger yields of binaries interacting with SgrA*, with consequences for HVSs, S-stars, EMRI progenitors, and Galactic-centre transients. The paper's main strength is that the restricted formalism yields predictions that are claimed to be independent of binary physical properties, and the results are obtained by direct numerical integration rather than by fitting to data. The fraction table (Table 1) is internally consistent, and the paper gives useful, checkable definitions of the outcome channels. The significance is currently conditional, however, because the multi-passage results rest on a phase-sampling convention that is not stated or justified, on a validation claim that is not shown quantitatively, and on an admitted inconsistency in the CM-frame update between passages.","major_comments":[{"comment":"The re-initialization of binary phases for CBs is not specified as time-weighted. Since Eq. (8) parameterizes the binary orbit by the true anomaly φ, drawing phases uniformly on [0,2π) samples a non-uniform distribution in time: for an eccentric binary the physical phase density is dN/dφ ∝ (1−e^2)^{3/2}/(1+e cosφ)^2. CBs after the first passage reach e_b up to about 0.8, so a uniform-φ draw overweights pericentre phases by a factor of several and underweights apocentre phases relative to the time-weighted ensemble. The second- and third-passage D/F/C/M fractions, and hence the headline '20% or more' disruption boost, are averaged over this initial-condition distribution. Please state the sampling distribution explicitly, recompute the multi-passage fractions using a mean-anomaly (or otherwise time-uniform) phase draw, and at minimum quantify the sensitivity of the reported boosts to this choice.","section":"§4, phase randomization; Eq. (8)"},{"comment":"The validation against full three-body integration is asserted but never shown quantitatively. The text says that comparison simulations with REBOUND show 'excellent agreement' and footnote 1 mentions a test, but no comparison figure, table, parameter ranges, or error statistics are provided. Because the restricted EOM and, especially, the iterative multi-passage update depend on this approximation, please supply a quantitative validation, for example relative differences in final a_b, e_b, and outcome classifications as a function of β0 and inclination, including at least one repeated passage.","section":"§3.1; footnote 1; §6"},{"comment":"Footnote 4 acknowledges that L_cm is inconsistent with the simulated r_cm and v_cm and that the CM orbital frame is not reoriented between passages, with ω_cm and Ω_cm left undefined. The multi-passage geometry therefore assumes that this omission is negligible. Given Q ≈ 10^6 the effect is plausibly tiny, but the paper should demonstrate this quantitatively—for example by propagating the full L_cm direction and recomputing the fractions, or by estimating the resulting change in the binary orientation relative to the updated CM orbital plane—before the multi-passage boosts can be taken at face value.","section":"§4.1, footnote 4"},{"comment":"The abstract states that, for the example system with finite stellar sizes and lifetimes, mergers occur 31% of the time, but Table 1's 'Mergers and lifetime' row at passage 3 gives 26.54%, with 30.70% corresponding to the merger cut only. The §6 bullet similarly reports ≈31% for mergers before the lifetime cut and 26.54% after combining both cuts. Please harmonize the abstract, the text, and Table 1, and specify explicitly which combination of cuts produces each number.","section":"Abstract; Table 1; §6 bullets"},{"comment":"Reclassifying CBs whose CM period exceeds the primary's lifetime as FAs conflates 'does not return for another passage' with 'unbound from the MBH'. These binaries remain on bound CM orbits but their stars expire before the next pericentre passage; they are not physically flying away. This reclassification changes the reported FA fractions and the physical interpretation of FAs as surviving binaries on unbound trajectories. Please either introduce a separate category for such systems or clearly state in all fraction tables and discussion that the FA fraction in the lifetime-cut case includes binaries that are not actually unbound.","section":"§3.3.1 and §4.2"}],"minor_comments":[{"comment":"The caption contains a typo: 'calcualtions' should be 'calculations'.","section":"Fig. 1 caption"},{"comment":"There are several typos in the introduction, including 'SgRA∗' for 'SgrA*' and 'torwards' for 'towards'.","section":"§1"},{"comment":"In the paragraph beginning 'In Fig. 2 we provide some examples', 'dirsupted' should be 'disrupted'.","section":"§3.2"},{"comment":"The first bullet of the Discussion contains 'thid' instead of 'third'; please proofread throughout.","section":"§6"},{"comment":"The text says initial orientations and phases are chosen 'uniformly and randomly', but it is not explicit whether the inclination i0 is drawn uniformly in cos(i0) or in i0. Since Fig. 3 uses cos(i0) as the ordinate, please state the sampling convention explicitly.","section":"§3.1 and Fig. 3"},{"comment":"The lower panels of Fig. 15 show |1−e_cm|/δ_B and v_ej/ν_B, but the text and caption do not state whether the medians shown refer to the first or third passage in the bottom two panels; please clarify.","section":"§5.3, Fig. 15"}],"recommendation":"major_revision","confidential_remarks":"The phase-sampling issue is the main correctness risk in the multi-passage boost claim; if the authors can show that a time-weighted phase draw does not change the conclusions, the paper would be close to acceptable. The REBOUND validation and CM-frame footnote are secondary but should be addressed in the revision. The self-citation pattern is natural given the formalism builds directly on Sari et al. (2010) and Kobayashi et al. (2012), and I do not see a circularity problem. The data availability statement is weak; recommending code or data release would strengthen reproducibility, but this is not a blocker."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the restricted three-body treatment of multiple pericentre passages, which no one has done systematically, and the detailed description of fly-away binaries, which Mandel & Levin only mentioned. The framework is transparent, the fractions in Table 1 sum consistently, and the rescaled units let the results scale to arbitrary binary masses and separations. The headline claim that multiple encounters boost disruptions by 20% or more is plausible and worth taking seriously.\n\nThe soft spot I would focus on is not the one the authors flag. Footnote 4 admits the CM orbital frame is not reoriented between passages, but the stress-test note is right: that inconsistency changes the encounter geometry only at the 1e-4 level because Delta L_b / L_cm ~ Q^{-2/3}. The real problem is the redrawn binary phase. Section 4.1 says each comeback binary is simulated with random phases, but it never states the distribution. If the phase is uniform in true anomaly phi, as Eq. (8) suggests, then for the eccentric binaries that come back (e_b up to ~0.8) the draw overweights pericentre by a factor of several relative to a time-weighted distribution. The physical phase at a random time is uniform in mean anomaly, and dN/dphi is proportional to (1-e^2)^{3/2}/(1+e cos phi)^2. This is an O(1) bias in the initial conditions of the very passages that produce the paper's central boost claim. The authors need to sample uniform mean anomaly, or demonstrate that the choice does not matter.\n\nTwo smaller issues: the REBOUND validation is claimed but never shown quantitatively, and no code or data are released, which makes it hard to check the phase convention retroactively. The merger fraction (31%) and the period/merger cuts are sensitive to the example system, which is fine, but the phase sampling affects the general claim independently of that example.\n\nIf the phase issue is fixed, the paper becomes a solid contribution to GC rate estimates for HVSs, S-stars, mergers, and possibly EMRI progenitors. As it stands, the qualitative conclusion (multi-encounters matter) is likely correct, but the quantitative boost should not be taken at face value. This deserves peer review, with a required revision on the phase sampling and a shown REBOUND comparison.","headline":"A useful multi-passage extension of Hills-mechanism calculations held back by an unsampled binary phase distribution and a few unshown validations.","tokens_in":30873,"tokens_out":1317,"would_cite":true,"duration_ms":15813,"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":"Multiple periapsis passages increase the Hills-mechanism disruption rate by at least 20 percent, and can make mergers a 31 percent outcome for example binaries.","keywords":["Hills mechanism","hypervelocity stars","SgrA*","tidal disruption","stellar binaries","multiple encounters","restricted three-body problem","galactic centre"],"falsifier":"A full three-body (or N-body) integration that propagates the centre-of-mass orbital plane and orientation between encounters, without the fixed-frame assumption, would settle the claim: if the resulting fractions of disruptions, fly-aways, and mergers after three passages differ by more than the quoted boosts (roughly 20 percent for disruptions), the paper's central conclusion fails.","tokens_in":29851,"feed_emoji":"🕳️","tokens_out":6146,"duration_ms":55445,"temperature":0.7,"pith_summary":"This paper argues that a stellar binary's encounter with the supermassive black hole SgrA* cannot be judged from a single close passage. Using a restricted three-body formalism, the authors follow initially circular binaries through up to three periapsis passages and show that binaries that survive the first passage often return, become more eccentric and more prograde, and then resolve into disruptions or ejections. The cumulative disruption fraction rises from about 46 percent after one passage to about 60 percent after three, a boost of 20 percent or more, and for an example massive binary, mergers occur 31 percent of the time when finite stellar sizes and lifetimes are included. If correct, single-passage treatments systematically underestimate how many binaries are disrupted or merged by the Hills mechanism.","feed_headline":"Multiple close passes boost black-hole binary disruptions by 20%","feed_subtitle":"A three-passage model shows surviving binaries often come back to be disrupted or to merge.","key_machinery":"The load-bearing object is the restricted three-body formalism, in which the binary's centre of mass follows a fixed Keplerian orbit around the stationary massive black hole and the internal binary motion is integrated under the linearized tidal force. Between passages the binary's final internal energy and angular-momentum changes are transferred to the centre-of-mass orbit through $E_{\\rm cm}^{n+1}=E_{\\rm cm}^n-\\Delta E_b$ and $\\mathbf{L}_{\\rm cm}^{n+1}=\\mathbf{L}_{\\rm cm}^n-\\Delta \\mathbf{L}_b$, yielding a new periapsis distance and eccentricity for the next encounter. The encounter is parameterized by the diving factor $\\beta=r_t/r_p$ (the ratio of tidal radius to periapsis distance) and the binary inclination; because the integration is rescaled by $\\lambda=(m/M)^{1/3}r_p$ and $\\tau=\\sqrt{r_p^3/GM}$, the results are independent of the binary's physical properties.","core_discovery":"The central claim is that multiple periapsis passages materially change the outcome distribution of the Hills mechanism. A binary whose centre of mass remains bound to SgrA* after one encounter comes back on a slightly altered, now elliptical orbit; because the binary itself has gained eccentricity and often a more prograde orientation, the next encounter is statistically more destructive. Tracked across three passages, the overall fraction of disrupted binaries increases from 45.6 percent to 59.6 percent (a boost of 20 percent or more for $\\beta_0 > 1$), the fly-away binary fraction grows from 28.4 to 32.3 percent, and, for an example 0.1 AU, 4 solar mass binary with radii following mass, mergers make up 31 percent of outcomes after three passages, reducing the disruption boost to roughly 10 percent once stellar lifetimes and mergers are included.","pith_inferences":["If multi-passage boosts hold generally, rate estimates for transient phenomena tied to the Hills mechanism—tidal disruption events, extreme-mass-ratio inspirals, and quasi-periodic eruptions in galactic nuclei—should be revised upward by tens of percent.","The formalism's scale-independence means the same qualitative boost should apply to compact-object binaries (white dwarfs, neutron stars, black holes) around SgrA*, for which stellar lifetime cuts vanish; the merger fraction there could be larger than the 31 percent quoted for main-sequence stars.","A direct test would compare the predicted low-velocity tail of ejected stars and the predicted population of eccentric fly-away binaries against future survey data in the Galactic Centre.","The three-passage truncation likely captures most of the effect, but the persistent retrograde come-back population suggests a slow tail: extending to more passages would test how quickly the cumulative fractions converge."],"forward_implications":["Single-passage estimates of Hills-mechanism outcomes undercount disruptions: including three passages raises the disruption fraction by 20 percent or more for deep encounters, and extends disruptions to shallower encounters below the single-passage threshold $\\beta_{\\rm lim}$.","Fly-away binaries form a distinct population ejected at speeds about two orders of magnitude lower than individual ejected stars, on marginally hyperbolic orbits, and their numbers are boosted by about 10 percent by later passages.","For stellar binaries with finite sizes, mergers are a major channel: about 31 percent of example systems merge by the third passage, mostly systems that would otherwise disrupt or fly away.","Hypervelocity stars are produced predominantly on the first passage; later passages add lower-velocity ejecta and matter most for shallow encounters near $\\beta_{\\rm lim}$.","Inclination and eccentricity control fate: non-retrograde, highly eccentric, or deep encounters preferentially disrupt, while retrograde binaries resist disruption and dominate the remaining come-back population."],"supporting_citations":[{"why":"Defines the tidal separation mechanism that the paper's multiple-encounter analysis extends.","marker":"Hills 1988"},{"why":"Supplies the single-passage restricted three-body formalism and outcome fractions that this paper generalizes to multiple passages.","marker":"Sari et al. 2010"},{"why":"Provides the energy/angular-momentum transfer framework and characteristic units used to update centre-of-mass orbits between passages.","marker":"Kobayashi et al. 2012"},{"why":"Earlier single-passage distributions of ejection velocities and captured orbits that the multi-passage results are compared against.","marker":"Brown et al. 2018a"},{"why":"Previous N-body follow-up of binaries surviving several revolutions around the SMBH; the paper compares its merger and HVS trends with this result.","marker":"Antonini et al. 2010"},{"why":"Prior estimate of merger fractions for a population of eccentric binaries; the paper's fly-away analysis is contrasted with its scope.","marker":"Mandel & Levin 2015"},{"why":"Provides alternative merger-fraction estimates for radial and shallow encounters used as comparison for the merger channel.","marker":"Bradnick et al. 2017"}],"fun_headline_variants":["Multiple SgrA* passes boost binary disruptions by 20%","Three encounters raise Hills-mechanism disruption odds 20%","Repeated black-hole encounters boost stellar binary disruptions","Multiple close passes increase binary disruption rate by 20%","Hills mechanism: extra passes boost disruptions by a fifth"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a returning binary's centre-of-mass orbit keeps the same spatial orientation between passages, so only the energy and angular-momentum magnitudes need updating; if the orbital plane or orientation rotates significantly, the second- and third-passage geometries—and the reported disruption, fly-away, and merger fractions—would change.","fun_headline_variants_meta":{"raw":{"variants":["Multiple SgrA* passes boost binary disruptions by 20%","Three encounters raise Hills-mechanism disruption odds 20%","Repeated black-hole encounters boost stellar binary disruptions","Multiple close passes increase binary disruption rate by 20%","Hills mechanism: extra passes boost disruptions by a fifth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000174,"raw_usage":{"total_tokens":1305,"prompt_tokens":991,"completion_tokens":314,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":607,"completion_tokens_details":{"reasoning_tokens":232}},"tokens_in":607,"tokens_out":314,"duration_ms":3697,"temperature":1.0,"reasoning_tokens":232,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:53:39.506224+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A full three-body (or N-body) integration that propagates the centre-of-mass orbital plane and orientation between encounters, without the fixed-frame assumption, would settle the claim: if the resulting fractions of disruptions, fly-aways, and mergers after three passages differ by more than the quoted boosts (roughly 20 percent for disruptions), the paper's central conclusion fails.","supporting_citations":[],"review_version":1}