No Arabic abstract
(Abridged) We review the results of the first multi-scale, hydrodynamical simulations of mergers between galaxies with central supermassive black holes (SMBHs) to investigate the formation of SMBH binaries in galactic nuclei. We demonstrate that strong gas inflows produce nuclear disks at the centers of merger remnants whose properties depend sensitively on the details of gas thermodynamics. In numerical simulations with parsec-scale spatial resolution in the gas component and an effective equation of state appropriate for a starburst galaxy, we show that a SMBH binary forms very rapidly, less than a million years after the merger of the two galaxies. Binary formation is significantly suppressed in the presence of a strong heating source such as radiative feedback by the accreting SMBHs. We also present preliminary results of numerical simulations with ultra-high spatial resolution of 0.1 pc in the gas component. These simulations resolve the internal structure of the resulting nuclear disk down to parsec scales and demonstrate the formation of a central massive object (~ 10^8 Mo) by efficient angular momentum transport. This is the first time that a radial gas inflow is shown to extend to parsec scales as a result of the dynamics and hydrodynamics involved in a galaxy merger, and has important implications for the fueling of SMBHs. Due to the rapid formation of the central clump, the density of the nuclear disk decreases significantly in its outer region, reducing dramatically the effect of dynamical friction and leading to the stalling of the two SMBHs at a separation of ~1 pc. We discuss how the orbital decay of the black holes might continue in a more realistic model which incorporates star formation and the multi-phase nature of the ISM.
It has been recently suggested that supermassive black holes at z = 5-6 might form from super-fast (dot M > 10^4 Msun/yr) accretion occurring in unstable, massive nuclear gas disks produced by mergers of Milky-Way size galaxies. Interestingly, such mechanism is claimed to work also for gas enriched to solar metallicity. These results are based on an idealized polytropic equation of state assumption, essentially preventing the gas from cooling. We show that under more realistic conditions, the disk rapidly (< 1 yr) cools, the accretion rate drops, and the central core can grow only to approx 100 Msun. In addition, most of the disk becomes gravitationally unstable in about 100 yr, further quenching the accretion. We conclude that this scenario encounters a number of difficulties that possibly make it untenable.
Supermassive black hole (SMBH) binaries represent the main target for missions such as the Laser Interferometer Space Antenna and Pulsar Timing Arrays. The understanding of their dynamical evolution prior to coalescence is therefore crucial to improving detection strategies and for the astrophysical interpretation of the gravitational wave data. In this paper, we use high-resolution $N$-body simulations to model the merger of two equal-mass galaxies hosting a central SMBH. In our models, all binaries are initially prograde with respect to the galaxy sense of rotation. But, binaries that form with a high eccentricity, $egtrsim 0.7$, quickly reverse their sense of rotation and become almost perfectly retrograde at the moment of binary formation. The evolution of these binaries proceeds towards larger eccentricities, as expected for a binary hardening in a counter-rotating stellar distribution. Binaries that form with lower eccentricities remain prograde and at comparatively low eccentricities. We study the origin of the orbital flip by using an analytical model that describes the early stages of binary evolution. This model indicates that the orbital plane flip is due to the torque from the triaxial background mass distribution that naturally arises from the galactic merger process. Our results imply the existence of a population of SMBH binaries with a high eccentricity and could have significant implications for the detection of the gravitational wave signal emitted by these systems.
We study the collapse of rapidly rotating supermassive stars that may have formed in the early Universe. By self-consistently simulating the dynamics from the onset of collapse using three-dimensional general-relativistic hydrodynamics with fully dynamical spacetime evolution, we show that seed perturbations in the progenitor can lead to the formation of a system of two high-spin supermassive black holes, which inspiral and merge under the emission of powerful gravitational radiation that could be observed at redshifts z>10 with the DECIGO or Big Bang Observer gravitational-wave observatories, assuming supermassive stars in the mass range 10^4-10^6 Msol. The remnant is rapidly spinning with dimensionless spin a^*=0.9. The surrounding accretion disk contains ~10% of the initial mass.
When galaxies collide, dynamical friction drives their central supermassive black holes close enought to each other such that gravitational radiation becomes the leading dissipative effect. Gravitational radiation takes away energy, momentum and angular momentum from the compact binary, such that the black holes finally merge. In the process, the spin of the dominant black hole is reoriented. On observational level, the spins are directly related to the jets, which can be seen at radio frequencies. Images of the X-shaped radio galaxies together with evidence on the age of the jets illustrate that the jets are reoriented, a phenomenon known as spin-flip. Based on the galaxy luminosity statistics we argue here that the typical galaxy encounters involve mass ratios between 1:3 to 1:30 for the central black holes. Based on the spin-orbit precession and gravitational radiation we also argue that for this typical mass ratio in the inspiral phase of the merger the initially dominant orbital angular momentum will become smaller than the spin, which will be reoriented. We prove here that the spin-flip phenomenon typically occurs already in the inspiral phase, and as such is describable by post-Newtonian techniques.
The binding energy liberated by the coalescence of supermassive black hole (SMBH) binaries during galaxy mergers is thought to be responsible for the low density cores often found in bright elliptical galaxies. We use high-resolution $N$-body and Monte Carlo techniques to perform single and multi-stage galaxy merger simulations and systematically study the dependence of the central galaxy properties on the binary mass ratio, the slope of the initial density cusps, and the number of mergers experienced. We study both the amount of depleted stellar mass (or ``mass deficit), $M_{rm def}$, and the radial extent of the depleted region, $r_{rm b}$. We find that $r_{rm b}simeq r_{rm SOI}$ and that $M_{rm def}$ varies in the range $0.5$ to $4M_{bullet}$, with $r_{rm SOI}$ the influence radius of the remnant SMBH and $M_{bullet}$ its mass. The coefficients in these relations depend weakly on the binary mass ratio and remain remarkably constant through subsequent mergers. We conclude that the core size and mass deficit do not scale linearly with the number of mergers, making it hard to infer merger histories from observations. On the other hand, we show that both $M_{rm def}$ and $r_{rm b}$ are sensitive to the morphology of the galaxy merger remnant, and that adopting spherical initial conditions, as done in early work, leads to misleading results. Our models reproduce the range of values for $M_{rm def}$ found in most observational work, but span nearly an order of magnitude range around the true ejected stellar mass.