{"id":"b4960199-e04c-462d-b784-25ca0f68c03a","arxiv_id":"1908.02563","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In N-body merger simulations, recoiling supermassive black holes escape major-merger remnants at velocities up to 25 percent lower than static analytical galaxy models predict, implying more spatially offset AGNs.","lead":"This paper compares how often recoiling supermassive black holes can escape their host galaxies, using static analytical galaxy models versus full N-body merger simulations. It finds that the static models overestimate escape velocities by up to 25 percent for major mergers, which implies more offset active galactic nuclei (AGNs) should exist than previously predicted.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative 25% escape-velocity gap rests on a single near-head-on merger orbit chosen for speed; more typical pericentric distances could substantially shrink it.","rationale":"The paper's strongest claim is explicitly quantitative: static analytical models overestimate SMBH escape velocities, with numerical major-merger escape velocities up to ~25 per cent lower, leading to more offset AGNs. The claim is supported internally by the pre-merger validation (Figs 2 and 3), the mass-profile evidence for violent relaxation (Figs 5 and 6), and the stability of remnant profiles after ejection. Those pieces make the qualitative direction credible. The load-bearing weak point is the single orbital configuration used for all mergers. Section 2.3.1 gives the choice away: the pericentre of 0.5% Rvir,1 is selected to make the encounter nearly head-on and computationally cheap. Violent relaxation, and therefore the mass loss that lowers escape velocities, is strongest for exactly such low-angular-momentum orbits. The cited cosmological motivation concerns eccentricity, not pericentre, so it does not rescue the particular value chosen. Because the paper does not vary orbital parameters or run multiple realizations, the 25 per cent figure is an upper limit for one geometrical extreme rather than a robust typical value. This does not require rejection: the paper consistently says \"can be up to\" and the qualitative conclusion is physically plausible. But it does justify the conditional verdict, since a realistic pericentre distribution could shrink the effect enough to change the predicted offset-AGN abundance. I therefore agree with the reader's identification of the weakest assumption and see no reason to change the conditional disposition.","tokens_in":19368,"tokens_out":6660,"duration_ms":82887,"concrete_test":"Repeat the Mgal = 10^12 Msun, 1:1, MBH = Mgal/10^5 major-merger case (open triangles in Fig. 4) with two additional orbital configurations, e.g. Rperi = 10% Rvir,1 and Rperi = 30% Rvir,1, keeping everything else identical. Recompute vesc,n/vesc,a from the SMBH trajectories. If the ratio moves from ~0.75 to above 0.90 in both cases, the headline percentage is an artifact of the head-on orbit; if the ratio remains below ~0.85, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.3.1 sets every merger to e = 1 and Rperi = 0.5% Rvir,1, explicitly described as leading to \"almost head-on collision and fast merging process that saves computational time.\" The central quantitative result — numerical major-merger remnants have escape velocities up to ~25 per cent lower than static analytical remnants (Abstract; Section 3.3; conclusion item i) — is attributed to violent-relaxation mass loss shown in Figs 5 and 6. Mass loss and the associated reduction of the potential at large radii are strong functions of the orbital angular momentum; low-angular-momentum, near-radial encounters are precisely the regime where violent relaxation is most efficient. The cited justification for e = 1 (Benson 2005; Khochfar & Burkert 2006) is that many cosmological orbits are near-parabolic, but e ~ 1 does not constrain the pericentric distance, so Rperi = 0.5% Rvir,1 remains an extreme choice. The paper runs one realization per configuration, with no variation in orbital parameters or numerical resolution, so the 25 per cent figure is not established as typical; it is an upper bound for one idealized geometry. If a more representative pericentre produces substantially less mass loss, the conclusion that numerical models predict a greater number of spatially offset AGNs weakens proportionally, even though the qualitative direction may survive.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares the trajectories of recoiling supermassive black holes (SMBHs) in static analytical galaxy potentials and in GADGET-2 N-body merger simulations, for remnant masses 10^11 and 10^12 Msun, major (1:1) and minor (1:10) mergers, and two SMBH mass scalings (Mgal/10^5 and Mgal/10^3). The numerical setup is first validated against analytical models for isolated progenitors by comparing escape velocities and maximal SMBH separations (Figs 2 and 3). The authors then simulate mergers with a single orbital configuration (e=1, pericentric distance 0.5% of the primary virial radius), construct analytical remnants from conserved progenitor masses, and compare escape velocities, mass profiles, and SMBH orbits. The central claim is that static analytical models overestimate SMBH escape velocities, with numerical major-merger remnants having escape velocities up to ~25 per cent lower, because violent relaxation during major mergers reduces the mass at large radii. The paper concludes that numerical models predict more spatially offset AGNs.","tokens_in":19589,"tokens_out":4223,"duration_ms":47549,"significance":"If the quantitative claim is robust, this is a useful and potentially important result: it identifies a systematic bias in static-potential models of post-merger galaxies and suggests that recoiling SMBHs are retained less often and found at larger offsets than those models predict. The paper has real strengths: the numerical method is validated against analytical models for progenitor galaxies (Figs 2 and 3), the stability of the initial galaxy models is checked, a softening-length test is reported, and the comparison with Illustris-based results of Blecha et al. (2016) grounds the work in the existing literature. The qualitative direction of the effect, driven by mass loss in major mergers, is plausible and physically well motivated. However, the headline magnitude of the effect rests on a narrow numerical setup: one orbital configuration, one realization per case, and an analytical remnant construction that needs to be stated more precisely. These limitations matter because the 25 per cent figure is the paper's main quantitative deliverable and is repeated in the abstract and conclusions.","major_comments":[{"comment":"The central quantitative claim that numerical major-merger escape velocities are up to ~25 per cent lower than analytical ones is established for only one orbital configuration: e=1 with Rperi=0.5% Rvir,1, which the manuscript itself describes as an almost head-on collision chosen to save computational time. The mass loss that drives the reduced escape velocity is expected to depend strongly on orbital angular momentum, so a near-radial encounter is likely to maximize the effect. I request either an additional merger simulation (or an energy-based estimate) with a more typical pericentric distance to show that the effect is not an artifact of this extreme orbit, or a revision of the abstract and conclusions to present the 25 per cent figure as an upper bound for near-radial mergers rather than as the typical difference.","section":"Section 2.3.1; Fig. 4; conclusion item (i)"},{"comment":"All merger outcomes are based on single realizations with no scatter estimates or multiple random seeds. The ratios rmax,n/rmax,a in Fig. 8 reach factors of ~10 for kicks near the escape velocity, and the claim that numerical models predict a greater number of spatially offset AGNs depends on these large ratios. Without a measure of realization-to-realization variance, the reader cannot assess whether the differences are statistically significant or partly due to the particular initial particle draw. At minimum, the authors should report the variance across realizations or soften the quantitative claims accordingly.","section":"Section 2.3.1; Figs 4 and 8"},{"comment":"The construction of the analytical merger remnants needs to be stated precisely because it is load-bearing for the main comparison. The sentence 'We make numerical models of isolated galaxies and then fit their mass profiles in order to produce analytical galaxies with the same properties' is ambiguous: if the analytical remnant parameters in Table 2 are fits to the numerical post-merger mass profiles, then the escape-velocity comparison would be partly circular. If instead the analytical remnants are constructed from conserved progenitor masses and the same density-profile shapes as the progenitors, with no fitting to the post-merger profiles shown in Fig. 5, that should be stated explicitly. The current wording does not rule out the circular reading.","section":"Section 2.3.2; Section 3.3"},{"comment":"The escape velocity is operationally defined as the kick velocity needed for the SMBH to return to the host halo after ~10 Gyr, rather than the classical escape speed of the instantaneous potential. This is a legitimate choice, but the comparison in Figs 4 and 5 is made at the time of SMBH merger/ejection, when the numerical remnant is still evolving (Fig. 6 shows continued central-profile evolution after ejection). The authors should justify that their conclusions are not sensitive to the chosen ejection time, for example by evaluating whether the early relaxation of the central regions is fully captured before the SMBH is kicked.","section":"Section 2.2.1; Fig. 6"}],"minor_comments":[{"comment":"The phrase 'so-cold final parsec problem' should read 'so-called final parsec problem'.","section":"Introduction, Section 1"},{"comment":"There are typographical repetitions: 'adopt the the similar approach' in Section 2.3.1 and 'the the galactic nucleus' in Section 3.5; these should be corrected.","section":"Section 2.3.1 and Section 3.5"},{"comment":"The caption says 'Panels show major and minor merger remnants' but the individual panels are not labeled in the caption text; identifying which panel corresponds to which case (for example, left versus right, upper versus lower) would improve readability.","section":"Figure 5 caption"},{"comment":"The text states that 'Escape velocities in weakened potential of major merger remnant are ≲ 25 per cent lower compared to minor merger remnants' before introducing the analytical-versus-numerical comparison in Section 3.3; the reader should be told explicitly that this is a statement about numerical models only, to avoid confusion with the 25 per cent figure presented later.","section":"Section 3.2"},{"comment":"The reference list contains a few formatting inconsistencies (for example 'Gonz´ alez' and the incomplete journal name in the 2007 reference), and the phrase 'decreasement' should be replaced by 'decrease' throughout the text.","section":"General"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid, clearly-written numerical study that makes a real point—live merger remnants have lower SMBH escape velocities than the static analytical potentials commonly used in recoil studies—but the headline 25 per cent is tied to a deliberately extreme merger orbit and should not be treated as a typical value.\n\nWhat's new: the quantitative comparison itself. Blecha et al. (2016), Smole (2015), and Choksi et al. (2017) all used static or pseudo-evolving potentials; this paper actually runs the merger and puts the kicked SMBH into the live remnant. The setup is careful: the same progenitor models are built analytically and numerically, and the agreement in isolated-galaxy escape velocities and maximum separations (Figs 2, 3) validates the integrator before the merger comparison. The mass profiles in Figs 5 and 6 show exactly the violent-relaxation signature you'd expect—central density increase, mass loss at large radii—so the mechanism for the lower escape velocities is well supported.\n\nThe soft spot is the one the stress-test flags, and it lands. Section 2.3.1 sets every merger to e=1 with pericenter 0.5% of the primary virial radius, explicitly \"almost head-on collision ... saves computational time.\" That is the regime where violent relaxation and mass loss are most efficient. A more typical pericenter would plausibly produce less mass loss and a smaller (though probably still nonzero) gap between numerical and analytical escape velocities. There is also one realization per configuration, no error bars, and no code or data release. So the \"up to 25 per cent\" is an upper bound for one idealized geometry, not a robust statistical statement. The qualitative direction—offset AGNs should be more common than static-potential models predict—is likely right, but the quantitative offset-AGN rates should not be taken from this paper alone.\n\nThe circularity worry is minor. Fitting analytical progenitors to numerical ones is a legitimate way to set up a controlled comparison; the post-merger analytical models conserve mass, so they aren't rigged to lose.\n\nVerdict: this deserves a serious referee. A good referee should ask for a more typical orbital configuration and a few realizations or a resolution test to see how the 25% moves. Without that, the paper is still worth publishing as a cautionary, proof-of-concept result. I'd bring it to reading group and probably cite it with a caveat.","headline":"A solid numerical check that static potentials overestimate recoil escape velocities in major mergers, but the headline 25 per cent rests on one extreme orbit and needs robustness testing before being used quantitatively.","tokens_in":20217,"tokens_out":3252,"would_cite":true,"duration_ms":35166,"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":"This paper establishes that evolving numerical merger remnant potentials give supermassive black hole escape velocities up to about 25 percent lower than static analytical potentials, making spatially offset active galactic nuclei more…","keywords":["supermassive black hole recoil","gravitational-wave kick","galaxy merger remnants","violent relaxation","escape velocity","offset active galactic nuclei","N-body simulations"],"falsifier":"Run a suite of equal-mass merger simulations of $10^{12}$ solar-mass progenitors with orbital eccentricities from 0.3 to 1.0 and pericentres from 0.5 to 10 percent of the virial radius at fixed resolution; if the mean numerical escape velocity stays within a few percent of the static analytical value for the non-parabolic orbits, the claim that numerical potentials reduce escape velocities by up to 25 percent would not generalise.","tokens_in":19119,"feed_emoji":"🕳️","tokens_out":7824,"duration_ms":77237,"temperature":0.7,"pith_summary":"The paper argues that the common shortcut of modelling a post-merger galaxy as a fixed, analytic potential makes recoiling supermassive black holes look harder to eject than they really are. In N-body simulations of major mergers, violent relaxation removes mass from the outer galaxy, so the escape velocity faced by a kicked black hole drops by up to about 25 percent relative to the analytic model. Since the same mass redistribution does not occur in minor mergers or in the static models, the difference is largest exactly for the major mergers that produce the biggest kicks. If the claim holds, offset active galactic nuclei—galaxies whose black hole is displaced from the centre—should be more numerous and more widely separated than static-potential estimates say.","feed_headline":"Merger simulations lower black hole escape speeds by 25%","feed_subtitle":"Static galaxy models overestimate escape speeds, so offset active galactic nuclei should be more common","key_machinery":"The engine of the comparison is the N-body merger simulation evolving the galactic potential while the recoiling black hole travels through it, set against a static analytical potential built from an NFW dark halo, a power-law bulge, and an exponential disc treated as spherical. Violent relaxation—the rapid, collisionless rearrangement of particle energies during a merger that pushes weakly bound stars beyond the virial radius—is the physical mechanism that reduces the outer mass profile and hence the escape velocity. The paper measures escape velocity as the kick speed needed for the black hole to return to the host halo after about 10 Gyr, and integrates the trajectory with Chandrasekhar dynamical friction included.","core_discovery":"The central claim is that escape velocities of recoiling supermassive black holes in numerical major merger remnant galaxies are up to about 25 percent lower than in analytical models, because violent relaxation during the merger lowers the galaxy's mass profile at large radii. The paper shows that the numerical and analytical models agree for isolated progenitor galaxies, isolating the merger itself as the source of the difference. In major merger remnants, the numerical potential has lower enclosed mass at large radii and a higher central density, and the outer mass loss wins: a black hole needs a smaller kick to leave and, for a given kick, reaches several to ten times larger maximum galactocentric distance over a Hubble time. The paper concludes that static analytical models overestimate escape velocities and underestimate the number of spatially offset active galactic nuclei.","pith_inferences":["I infer that the 25 percent headline is tied to the one orbital family tested (parabolic, pericentre 0.5 percent of the primary virial radius); real mergers with smaller eccentricity or larger pericentre would likely lose less mass, so the analytical–numerical gap should shrink.","A testable extension would be to rerun the same comparison varying orbital parameters and numerical resolution; if the gap persists across a realistic orbit distribution, observational campaigns for offset active galactic nuclei should weight the numerical predictions more heavily.","I infer that the enhanced separations at $v_{\\rm kick}/v_{\\rm esc} \\sim 0.2$ translate directly into discovery-rate predictions: the same spin-alignment assumptions that keep black holes in the central kiloparsec in analytic models place them several kiloparsecs out in numerical remnants."],"forward_implications":["Offset active galactic nuclei should be more common than static-potential estimates predict, because lower kick speeds can still displace the black hole and a displaced black hole spends more time on bound orbits outside the nucleus.","For a given kick-to-escape ratio, numerical major merger remnants put recoiling black holes up to about ten times farther from the centre than analytical models do, making them easier to detect.","Black hole retention is lower after major than after minor mergers, so merger-driven black hole growth should be somewhat suppressed in major-merger remnants.","Escape velocities of roughly $10^{11}$ and $10^{12}$ solar-mass remnants fall in ranges of about 170–350 and 500–700 km/s depending on black hole mass and merger type, so typical kicks below a few hundred km/s matter most for the smaller host."],"supporting_citations":[{"why":"Supplies the analytical merger-remnant construction, the kick-velocity distributions, and the offset-AGN predictions that this paper extends and compares with.","marker":"Blecha et al. 2016"},{"why":"Supplies the violent-relaxation mechanism invoked to explain the mass loss that lowers numerical escape velocities.","marker":"Lynden-Bell 1967"},{"why":"Supports the claim that violent relaxation is efficient mainly in major mergers, which drives the 25 percent gap.","marker":"Hilz et al. 2012"},{"why":"Supplies the N-body simulation code used to evolve the galaxy mergers whose remnant potentials produce the lower escape velocities.","marker":"Springel 2005"},{"why":"Supplies the NFW dark halo density profile used in the analytical models that are shown to overestimate escape velocities.","marker":"Navarro, Frenk & White 1997"},{"why":"Supplies the dynamical friction formula used in both the analytical and numerical trajectory integrations.","marker":"Chandrasekhar 1943"},{"why":"Supplies the spherical exponential-disc potential used as a component of the analytical galaxy models.","marker":"Geehan et al. 2006"}],"fun_headline_variants":["Numerical merger models cut escape speeds by 25%","Static potentials overestimate black hole escape speeds","Black holes roam farther in evolving merger galaxies","Merger simulations yield more offset AGNs","Escape velocity drop 25% in dynamic merger models"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 25 percent gap rests on one simulated merger geometry—a parabolic encounter with pericentre 0.5 percent of the primary virial radius and $10^6$ particles—being representative of the violent-relaxation mass loss in real major mergers; with less radial orbits or different resolution, the gap could shrink.","fun_headline_variants_meta":{"raw":{"variants":["Numerical merger models cut escape speeds by 25%","Static potentials overestimate black hole escape speeds","Black holes roam farther in evolving merger galaxies","Merger simulations yield more offset AGNs","Escape velocity drop 25% in dynamic merger models"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001364,"raw_usage":{"total_tokens":5513,"prompt_tokens":906,"completion_tokens":4607,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":4535}},"tokens_in":522,"tokens_out":4607,"duration_ms":35832,"temperature":1.0,"reasoning_tokens":4535,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:39:53.737415+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a suite of equal-mass merger simulations of $10^{12}$ solar-mass progenitors with orbital eccentricities from 0.3 to 1.0 and pericentres from 0.5 to 10 percent of the virial radius at fixed resolution; if the mean numerical escape velocity stays within a few percent of the static analytical value for the non-parabolic orbits, the claim that numerical potentials reduce escape velocities by up to 25 percent would not generalise.","supporting_citations":[{"cited_title":"Z., Torrey P., Vogelsberger M., Nelson D., Springel V., Snyder G., Hernquist L., 2016, MNRAS, 456, 961","cited_arxiv_id":null,"evidence_quote":"Supplies the analytical merger-remnant construction, the kick-velocity distributions, and the offset-AGN predictions that this paper extends and compares with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the violent-relaxation mechanism invoked to explain the mass loss that lowers numerical escape velocities."},{"cited_title":"P., Thomas J., Burkert A., Jesseit R., 2012, MNRAS, 3136, 311","cited_arxiv_id":null,"evidence_quote":"Supports the claim that violent relaxation is efficient mainly in major mergers, which drives the 25 percent gap."},{"cited_title":"J., Fardal M","cited_arxiv_id":null,"evidence_quote":"Supplies the spherical exponential-disc potential used as a component of the analytical galaxy models."}],"review_version":1}