Post-Newtonian dynamical modeling of supermassive black holes in galactic-scale simulations

Abstract

We present KETJU, a new extension of the widely-used smoothed particle hydrodynamics simulation code GADGET-3. The key feature of the code is the inclusion of algorithmically regularized regions around every supermassive black hole (SMBH). This allows for simultaneously following global galactic-scale dynamical and astrophysical processes, while solving the dynamics of SMBHs, SMBH binaries and surrounding stellar systems at sub-parsec scales. The KETJU code includes Post-Newtonian terms in the equations of motions of the SMBHs which enables a new SMBH merger criterion based on the gravitational wave coalescence timescale pushing the merger separation of SMBHs down to $\sim 0.005$ pc. We test the performance of our code by comparison to NBODY7 and rVINE. We set up dynamically stable multi-component merger progenitor galaxies to study the SMBH binary evolution during galaxy mergers. In our simulation sample the SMBH binaries do not suffer from the final-parsec problem, which we attribute to the non-spherical shape of the merger remnants. For bulge-only models, the hardening rate decreases with increasing resolution, whereas for models which in addition include massive dark matter halos the SMBH binary hardening rate becomes practically independent of the mass resolution of the stellar bulge. The SMBHs coalesce on average 200 Myr after the formation of the SMBH binary. However, small differences in the initial SMBH binary eccentricities can result in large differences in the SMBH coalescence times. Finally, we discuss the future prospects of KETJU, which allows for a straightforward inclusion of gas physics in the simulations.

28 pages, 18 figures. Accepted for publication in ApJ

References in corpus (27)

• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
• 3D-HST+CANDELS: The Evolution of the Galaxy Size-Mass Distribution since $z=3$
• A unified model for AGN feedback in cosmological simulations of structure formation
• The Illustris simulation: the evolving population of black holes across cosmic time
• Intermediate and Extreme Mass-Ratio Inspirals -- Astrophysics, Science Applications and Detection using LISA
• A massive binary black-hole system in OJ287 and a test of general relativity
• Rapid Formation of Supermassive Black Hole Binaries in Galaxy Mergers with Gas
• The Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation
• Equal- and unequal-mass mergers of disk and elliptical galaxies with black holes
• Massive Black Hole Binary Mergers in Dynamical Galactic Environments
• Interaction of massive black hole binaries with their stellar environment: I. Ejection of hypervelocity stars
• Evidence for Wide-Spread AGN Driven Outflows in the Most Massive z~1-2 Star Forming Galaxies
• A Candidate Sub-Parsec Supermassive Binary Black Hole System
• Feedback from Active Galactic Nuclei: Energy- versus momentum-driving
• Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
• Off the Beaten Path: A New Approach to Realistically Model The Orbital Decay of Supermassive Black Holes in Galaxy Formation Simulations
• Performance Analysis of Direct N-Body Algorithms on Special-Purpose Supercomputers
• A 17-billion-solar-mass black hole in a group galaxy with a diffuse core
• A close-pair binary in a distant triple supermassive black-hole system
• Algorithmic regularization with velocity-dependent forces
• Long-Term Evolution of Massive Black Hole Binaries. III. Binary Evolution in Collisional Nuclei
• A Displaced Supermassive Black Hole in M87
• Dynamical friction of massive objects in galactic centres
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods
• VINE -- A numerical code for simulating astrophysical systems using particles II: Implementation and performance characteristics
• A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces

Cited by in corpus (52)

• Simulating cosmic structure formation with the GADGET-4 code
• The Gravitational Wave Background from Massive Black Hole Binaries in Illustris: spectral features and time to detection with pulsar timing arrays
• Stellar-mass black holes in young massive and open stellar clusters and their role in gravitational-wave generation II
• The Role of Black Hole Feedback on Size and Structural Evolution in Massive Galaxies
• Post-Newtonian evolution of massive black hole triplets in galactic nuclei -- II. Survey of the parameter space
• Origins and Demographics of Wandering Black Holes
• Galaxy mergers in EAGLE do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN
• The formation of extremely diffuse galaxy cores by merging supermassive black holes
• Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems
• Eccentricity evolution of massive black hole binaries from formation to coalescence
• The simultaneous formation of cored, tangentially biased, and kinematically decoupled centers in massive early-type galaxies
• MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates
• Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe
• Black Hole Growth and Feedback in Isolated Romulus25 Dwarf Galaxies
• Dynamical Evolution and Merger Time-scales of LISA Massive Black Hole Binaries in Disk Galaxy Mergers
• Signatures of the Many Supermassive Black Hole Mergers in a Cosmologically Forming Massive Early-Type Galaxy
• Resolving the Complex Evolution of a Supermassive Black Hole Triplet in a Cosmological Simulation
• Defeating stochasticity: coalescence timescales of massive black holes in galaxy mergers
• Inward Bound: The incredible journey of massive black holes as they pair and merge; I. The effect of mass ratio in flattened rotating galactic nuclei
• Gravitational Waves from the Inspiral of Supermassive Black Holes in Galactic-scale Simulations
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• SMART: A new implementation of Schwarzschild's Orbit Superposition technique for triaxial galaxies and its application to an N-body merger simulation
• Gravitational Wave Driven Mergers and Coalescence Time of Supermassive Black Holes
• The two phases of core formation -- orbital evolution in the centres of ellipticals with supermassive black hole binaries
• Close Encounters of Wide Binaries Induced by the Galactic Tide: Implications for Stellar Mergers and Gravitational-Wave Sources
• KETJU -- resolving small-scale supermassive black hole dynamics in GADGET-4
• A New Discrete Dynamical Friction Estimator Based on $N$-body Simulations
• Extremely efficient mergers of intermediate mass black hole binaries in nucleated dwarf galaxies
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers
• Morphological evolution of supermassive black hole merger hosts and multimessenger signatures
• RABBITS -- I. The crucial role of nuclear star formation in driving the coalescence of supermassive black hole binaries
• Accuracy and precision of triaxial orbit models I: SMBH mass, stellar mass and dark-matter halo
• Evolution of supermassive black hole binaries and tidal disruption rates in nonspherical galactic nuclei
• Dynamics and spin alignment in massive, gravito-turbulent circumbinary discs around supermassive black hole binaries
• The isotropic center of NGC 5419 -- A core in formation?
• NGC 6240 Supermassive Black Hole Binary dynamical evolution based on Chandra data
• Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
• The supermassive black hole merger driven evolution of high-redshift red nuggets into present-day cored early-type galaxies
• RABBITS -- II. The impact of AGN feedback on coalescing supermassive black holes in disc and elliptical galaxy mergers
• Mergers all the way down: stellar collisions and kinematics of a dense hierarchically forming massive star cluster in a dwarf starburst
• Identifying supermassive black hole recoil in elliptical galaxies
• On the detectability of massive black hole merger eventsby LISA
• Wandering and escaping: recoiling massive black holes in cosmological simulations
• Astrophysics with the Laser Interferometer Space Antenna
• Properties of loss cone stars in a cosmological galaxy merger remnant
• Premature supermassive black hole mergers in cosmological simulations of structure formation
• Tracking on-the-fly massive black hole binary evolution and coalescence in galaxy simulations: RAMCOAL
• FROST-CLUSTERS -- III. Metallicity-dependent intermediate mass black hole formation by runaway collisions in dense star clusters
• Dynamical evolution timescales for the supermassive black hole system in the galaxy NGC 7727 (Arp 222)
• Recoiling supermassive black holes in analytical and numerical galaxy potential
• Rapid sinking and efficient mergers of supermassive black holes in compact high-redshift galaxies