-
RABBITS -- I. The crucial role of nuclear star formation in driving the coalescence of supermassive black hole binaries
Authors:
Shihong Liao,
Dimitrios Irodotou,
Peter H. Johansson,
Thorsten Naab,
Francesco Paolo Rizzuto,
Jessica M. Hislop,
Alexander Rawlings,
Ruby J. Wright
Abstract:
In this study of the `Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations' (RABBITS) series, we focus on the hardening and coalescing process of supermassive black hole (SMBH) binaries in galaxy mergers. For simulations including different galaxy formation processes (i.e. gas cooling, star formation, SMBH accretion, stellar and AGN feedback), we systematically control…
▽ More
In this study of the `Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations' (RABBITS) series, we focus on the hardening and coalescing process of supermassive black hole (SMBH) binaries in galaxy mergers. For simulations including different galaxy formation processes (i.e. gas cooling, star formation, SMBH accretion, stellar and AGN feedback), we systematically control the effect of stochastic eccentricity by fixing it to similar values during the SMBH hardening phase. We find a strong correlation between the SMBH merger time-scales and the presence of nuclear star formation. Throughout the galaxy merging process, gas condenses at the centre due to cooling and tidal torques, leading to nuclear star formation. These recently formed stars, which inherit low angular momenta from the gas, contribute to the loss cone and assist in the SMBH hardening via three-body interactions. Compared to non-radiative hydrodynamical runs, the SMBH merger time-scales measured from the runs including cooling, stellar and SMBH physical processes tend to be shortened by a factor of ${\sim}1.7$. After fixing the eccentricity to the range of $e \sim 0.6$--$0.8$ during the hardening phase, the simulations with AGN feedback reveal merger time-scales of ${\sim} 100$--$500$ Myr for disc mergers and ${\sim} 1$--$2$ Gyr for elliptical mergers. With a semi-analytical approach, we find that the torque interaction between the binary and its circumbinary disc has minimal impact on the shrinking of the binary orbit in our retrograde galaxy merger. Our results are useful in improving the modelling of SMBH merger time-scales and gravitational wave event rates.
△ Less
Submitted 1 February, 2024; v1 submitted 2 November, 2023;
originally announced November 2023.
-
RABBITS -- II. The impact of AGN feedback on coalescing supermassive black holes in disc and elliptical galaxy mergers
Authors:
Shihong Liao,
Dimitrios Irodotou,
Peter H. Johansson,
Thorsten Naab,
Francesco Paolo Rizzuto,
Jessica M. Hislop,
Ruby J. Wright,
Alexander Rawlings
Abstract:
In this study of the `Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations' (RABBITS) series, we investigate the orbital evolution of supermassive black holes (SMBHs) during galaxy mergers. We simulate both disc and elliptical galaxy mergers using the KETJU code, which can simultaneously follow galaxy (hydro-)dynamics and small-scale SMBH dynamics with post-Newtonian c…
▽ More
In this study of the `Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations' (RABBITS) series, we investigate the orbital evolution of supermassive black holes (SMBHs) during galaxy mergers. We simulate both disc and elliptical galaxy mergers using the KETJU code, which can simultaneously follow galaxy (hydro-)dynamics and small-scale SMBH dynamics with post-Newtonian corrections. With our SMBH binary subgrid model, we show how active galactic nuclei (AGNs) feedback affects galaxy properties and SMBH coalescence. We find that simulations without AGN feedback exhibit excessive star formation, resulting in merger remnants that deviate from observed properties. Kinetic AGN feedback proves more effective than thermal AGN feedback in expelling gas from the centre and quenching star formation. The different central galaxy properties, which are a result of distinct AGN feedback models, lead to varying rates of SMBH orbital decay. In the dynamical friction phase, galaxies with higher star formation and higher SMBH masses possess denser centres, become more resistant to tidal stripping, experience greater dynamical friction, and consequently form SMBH binaries earlier. As AGN feedback reduces gas densities in the centres, dynamical friction by stars dominates over gas. In the SMBH hardening phase, compared to elliptical mergers, disc mergers exhibit higher central densities of newly formed stars, resulting in accelerated SMBH hardening and shorter merger time-scales (i.e. $\lesssim 500$ Myr versus $\gtrsim 1$ Gyr). Our findings highlight the importance of AGN feedback and its numerical implementation in understanding the SMBH coalescing process, a key focus for low-frequency gravitational wave observatories.
△ Less
Submitted 24 April, 2024; v1 submitted 2 November, 2023;
originally announced November 2023.
-
Reviving stochasticity: uncertainty in SMBH binary eccentricity is unavoidable
Authors:
Alexander Rawlings,
Matias Mannerkoski,
Peter H. Johansson,
Thorsten Naab
Abstract:
We study supermassive black hole (SMBH) binary eccentricity of equal-mass galaxy mergers in $N$-body simulations with the KETJU code, which combines the GAGDET-4 fast multipole gravity solver with accurate regularised integration and post-Newtonian corrections around SMBHs. In simulations with realistic, high eccentricity galactic merger orbits, the hard binary eccentricity is found to be a non-li…
▽ More
We study supermassive black hole (SMBH) binary eccentricity of equal-mass galaxy mergers in $N$-body simulations with the KETJU code, which combines the GAGDET-4 fast multipole gravity solver with accurate regularised integration and post-Newtonian corrections around SMBHs. In simulations with realistic, high eccentricity galactic merger orbits, the hard binary eccentricity is found to be a non-linear function of the deflection angle in the SMBH orbit during the final, nearly radial close encounter between the SMBHs before they form a bound binary. This mapping between the deflection angle and the binary eccentricity has no apparent resolution dependence in our simulations spanning the resolution range of $1\times10^5$-$8\times10^6$ particles per galaxy. The mapping is also captured using a simple model with an analytic potential, indicating that it is driven by the interplay between a smooth asymmetric stellar background potential and dynamical friction acting on the SMBHs. Due to the non-linearity of this mapping, in eccentric major merger configurations small, parsec-scale variations in the merger orbit can result in binary eccentricities varying in nearly the full possible range between $e=0$ and $e=1$. In idealised simulations, such variations are caused by finite resolution effects, and convergence of the binary eccentricity can be achieved with increasing resolution. However, in real galaxies, other mechanisms such as nuclear gas and substructure that perturb the merger orbit are likely to be significant enough for the binary eccentricity to be effectively random. Our results indicate that the distribution of these effectively random eccentricities can be studied using even moderate resolution simulations.
△ Less
Submitted 21 September, 2023; v1 submitted 17 July, 2023;
originally announced July 2023.
-
KETJU -- resolving small-scale supermassive black hole dynamics in GADGET-4
Authors:
Matias Mannerkoski,
Alexander Rawlings,
Peter H. Johansson,
Thorsten Naab,
Antti Rantala,
Volker Springel,
Dimitrios Irodotou,
Shihong Liao
Abstract:
We present the new public version of the KETJU supermassive black hole (SMBH) dynamics module, as implemented into GADGET-4. KETJU adds a small region around each SMBH where the dynamics of the SMBHs and stellar particles are integrated using an algorithmically regularised integrator instead of the leapfrog integrator with gravitational softening used by GADGET-4. This enables modelling SMBHs as p…
▽ More
We present the new public version of the KETJU supermassive black hole (SMBH) dynamics module, as implemented into GADGET-4. KETJU adds a small region around each SMBH where the dynamics of the SMBHs and stellar particles are integrated using an algorithmically regularised integrator instead of the leapfrog integrator with gravitational softening used by GADGET-4. This enables modelling SMBHs as point particles even during close interactions with stellar particles or other SMBHs, effectively removing the spatial resolution limitation caused by gravitational softening. KETJU also includes post-Newtonian corrections, which allows following the dynamics of SMBH binaries to sub-parsec scales and down to tens of Schwarzschild radii. Systems with multiple SMBHs are also supported, with the code also including the leading non-linear cross terms that appear in the post-Newtonian equations for such systems. We present tests of the code showing that it correctly captures, at sufficient mass resolution, the sinking driven by dynamical friction and binary hardening driven by stellar scattering. We also present an example application demonstrating how the code can be applied to study the dynamics of SMBHs in mergers of multiple galaxies and the effect they have on the properties of the surrounding galaxy. We expect that the presented KETJU SMBH dynamics module can also be straightforwardly incorporated into other codes similar to GADGET-4, which would allow coupling small-scale SMBH dynamics to the rich variety of galactic physics models that exist in the literature.
△ Less
Submitted 8 June, 2023;
originally announced June 2023.
-
Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers
Authors:
Shihong Liao,
Peter H. Johansson,
Matias Mannerkoski,
Dimitrios Irodotou,
Francesco Paolo Rizzuto,
Stuart McAlpine,
Antti Rantala,
Alexander Rawlings,
Till Sawala
Abstract:
We introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the KETJU code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi--Hoyle--Lyttleton accretion into the binary phase and incorporates preferential…
▽ More
We introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the KETJU code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi--Hoyle--Lyttleton accretion into the binary phase and incorporates preferential mass accretion onto the secondary SMBH, which is motivated by results from small-scale hydrodynamical circumbinary disc simulations. We perform idealised gas-rich disc galaxy merger simulations using pure thermal or pure kinetic active galactic nuclei (AGN) feedback. Our binary accretion model provides more physically motivated SMBH mass ratios, which are one of the key parameters for computing gravitational wave (GW) induced recoil velocities. The merger time-scales of our simulated SMBH binaries are in the range $t_{\rm merge}{\sim} 10$--$400$ Myr. Prograde in-plane equal-mass galaxy mergers lead to the shortest merger time-scales, as they experience the strongest starbursts, with the ensuing high stellar density resulting in a rapid SMBH coalescence. Compared to the thermal AGN feedback, the kinetic AGN feedback predicts longer merger time-scales and results in more core-like stellar profiles, as it is more effective in removing gas from the galaxy centre and quenching star formation. This suggests that the AGN feedback implementation plays a critical role in modelling SMBH coalescences. Our model will be useful for improving the modelling of SMBH mergers in gas-rich galaxies, the prime targets for the upcoming LISA GW observatory.
△ Less
Submitted 18 February, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
The complex evolution of supermassive black holes in cosmological simulations
Authors:
Peter H. Johansson,
Matias Mannerkoski,
Antti Rantala,
Shihong Liao,
Alexander Rawlings,
Dimitrios Irodotou,
Francesco Rizzuto
Abstract:
We present here self-consistent zoom-in simulations of massive galaxies forming in a full cosmological setting. The simulations are run with an updated version of the KETJU code, which is able to resolve the gravitational dynamics of their supermassive black holes, while simultaneously modelling the large-scale astrophysical processes in the surrounding galaxies, such as gas cooling, star formatio…
▽ More
We present here self-consistent zoom-in simulations of massive galaxies forming in a full cosmological setting. The simulations are run with an updated version of the KETJU code, which is able to resolve the gravitational dynamics of their supermassive black holes, while simultaneously modelling the large-scale astrophysical processes in the surrounding galaxies, such as gas cooling, star formation and stellar and AGN feedback. The KETJU code is able to accurately model the complex behaviour of multiple SMBHs, including dynamical friction, stellar scattering and gravitational wave emission, and also to resolve Lidov-Kozai oscillations that naturally occur in hierarchical triplet SMBH systems. In general most of the SMBH binaries form at moderately high eccentricities, with typical values in the range of e =0.6-0.95, meaning that the circular binary models that are commonly used in the literature are insufficient for capturing the typical binary evolution.
△ Less
Submitted 23 March, 2022;
originally announced March 2022.
-
Signatures of the Many Supermassive Black Hole Mergers in a Cosmologically Forming Massive Early-Type Galaxy
Authors:
Matias Mannerkoski,
Peter H. Johansson,
Antti Rantala,
Thorsten Naab,
Shihong Liao,
Alexander Rawlings
Abstract:
We model here the merger histories of the supermassive black hole (SMBH) population in the late stages of a cosmological simulation of a $\sim 2 \times 10^{13} M_\odot$ galaxy group. The gravitational dynamics around the several tens of SMBHs ($M_{\bullet} \gtrsim 7.5\times 10^7 M_\odot$) hosted by the galaxies in the group is computed at high accuracy using regularized integration with the KETJU…
▽ More
We model here the merger histories of the supermassive black hole (SMBH) population in the late stages of a cosmological simulation of a $\sim 2 \times 10^{13} M_\odot$ galaxy group. The gravitational dynamics around the several tens of SMBHs ($M_{\bullet} \gtrsim 7.5\times 10^7 M_\odot$) hosted by the galaxies in the group is computed at high accuracy using regularized integration with the KETJU code. The 11 SMBHs which form binaries and hierarchical triplets eventually merge after hardening through dynamical friction, stellar scattering, and gravitational wave (GW) emission. The binaries form at eccentricities of $e \sim 0.3$-$0.9$, with one system evolving to a very high eccentricity of $e = 0.998$, and merge on timescales of a few tens to several hundred megayears. During the simulation the merger-induced GW recoil kicks eject one SMBH remnant from the central host galaxy. This temporarily drives the galaxy off the $M_{\bullet}$-$σ_{\star}$ relation, however the galaxy returns to the relation due to subsequent galaxy mergers, which bring in new SMBHs. This showcases a possible mechanism contributing to the observed scatter of the $M_{\bullet}$-$σ_{\star}$ relation. Finally, we show that Pulsar Timing Arrays and LISA would be able to detect parts of the GW signals from the SMBH mergers that occur during the $\sim 4\,\mathrm{Gyr}$ time span simulated with KETJU.
△ Less
Submitted 11 April, 2022; v1 submitted 7 December, 2021;
originally announced December 2021.
-
The SAMI Galaxy Survey: rules of behaviour for spin-ellipticity radial tracks in galaxies
Authors:
Alexander Rawlings,
Caroline Foster,
Jesse van de Sande,
Dan S. Taranu,
Scott M. Croom,
Joss Bland-Hawthorn,
Sarah Brough,
Julia J. Bryant,
Matthew Colless,
Claudia del P. Lagos,
Iraklis S. Konstantopoulos,
Jon S. Lawrence,
Ángel R. López-Sánchez,
Nuria P. F. Lorente,
Anne M. Medling,
Sree Oh,
Matt S. Owers,
Samuel N. Richards,
Nicholas Scott,
Sarah M. Sweet,
Sukyoung K. Yi
Abstract:
We study the behaviour of the spin-ellipticity radial tracks for 507 galaxies from the Sydney AAO Multi-object Integral Field (SAMI) Galaxy Survey with stellar kinematics out to $\geq1.5R_\text{e}$. We advocate for a morpho-dynamical classification of galaxies, relying on spatially-resolved photometric and kinematic data. We find the use of spin-ellipticity radial tracks is valuable in identifying…
▽ More
We study the behaviour of the spin-ellipticity radial tracks for 507 galaxies from the Sydney AAO Multi-object Integral Field (SAMI) Galaxy Survey with stellar kinematics out to $\geq1.5R_\text{e}$. We advocate for a morpho-dynamical classification of galaxies, relying on spatially-resolved photometric and kinematic data. We find the use of spin-ellipticity radial tracks is valuable in identifying substructures within a galaxy, including embedded and counter-rotating discs, that are easily missed in unilateral studies of the photometry alone. Conversely, bars are rarely apparent in the stellar kinematics but are readily identified on images. Consequently, we distinguish the spin-ellipticity radial tracks of seven morpho-dynamical types: elliptical, lenticular, early spiral, late spiral, barred spiral, embedded disc, and 2-sigma galaxies. The importance of probing beyond the inner radii of galaxies is highlighted by the characteristics of galactic features in the spin-ellipticity radial tracks present at larger radii. The density of information presented through spin-ellipticity radial tracks emphasises a clear advantage to representing galaxies as a track, rather than a single point, in spin-ellipticity parameter space.
△ Less
Submitted 4 October, 2019; v1 submitted 3 October, 2019;
originally announced October 2019.