No Arabic abstract
The presence of young massive stars orbiting on eccentric rings within a few tenths of a parsec of the supermassive black hole in the Galactic centre is challenging for theories of star formation. The high tidal shear from the black hole should tear apart the molecular clouds that form stars elsewhere in the Galaxy, while transporting the stars to the Galactic centre also appears unlikely during their stellar lifetimes. We present numerical simulations of the infall of a giant molecular cloud that interacts with the black hole. The transfer of energy during closest approach allows part of the cloud to become bound to the black hole, forming an eccentric disc that quickly fragments to form stars. Compressional heating due to the black hole raises the temperature of the gas to 100-1000K, ensuring that the fragmentation produces relatively high stellar masses. These stars retain the eccentricity of the disc and, for a sufficiently massive initial cloud, produce an extremely top-heavy distribution of stellar masses. This potentially repetitive process can therefore explain the presence of multiple eccentric rings of young stars in the presence of a supermassive black hole.
We show that the rapid formation of super-massive black holes in quasars can indeed be understood in terms of major galaxy mergers followed by disk accretion. The necessary short disk evolution time can be achieved provided the disk viscosity is sufficiently large, which, for instance, is the case for hydrodynamic turbulence, unlimited by shock dissipation. We present numerical calculations for a representative case. This general picture can account for (a) the presence of highly luminous quasars at redshifts z > 6; (b) for the peak in quasar activity at z ~ 2; and (c) for a subsequent rapid disappearance of quasars at later epochs.
We addressed the so far unexplored issue of outflows induced by exponentially growing power sources, focusing on early supermassive black holes (BHs). We assumed that these objects grow to $10^9;M_{odot}$ by z=6 by Eddington-limited accretion and convert 5% of their bolometric output into a wind. We first considered the case of energy-driven and momentum-driven outflows expanding in a region where the gas and total mass densities are uniform and equal to the average values in the Universe at $z>6$. We derived analytic solutions for the evolution of the outflow, finding that, for an exponentially growing power with e-folding time $t_{Sal}$, the late time expansion of the outflow radius is also exponential, with e-folding time of $5t_{Sal}$ and $4t_{Sal}$ in the energy-driven and momentum-driven limit, respectively. We then considered energy-driven outflows produced by QSOs at the center of early dark matter halos of different masses and powered by BHs growing from different seeds. We followed the evolution of the source power and of the gas and dark matter density profiles in the halos from the beginning of the accretion until $z=6$. The final bubble radius and velocity do not depend on the seed BH mass but are instead smaller for larger halo masses. At z=6, bubble radii in the range 50-180 kpc and velocities in the range 400-1000 km s$^{-1}$ are expected for QSOs hosted by halos in the mass range $3times10^{11}-10^{13};M_{odot}$. By the time the QSO is observed, we found that the total thermal energy injected within the bubble in the case of an energy-driven outflow is $E_{th}sim5 times 10^{60}$ erg. This is in excellent agreement with the value of $E_{th}=(6.2pm 1.7)times 10^{60}$ erg measured through the detection of the thermal Sunyaev-Zeldovich effect around a large population of luminous QSOs at lower redshift. [abridged]
We study the angular resolution of the gravitational wave detector LISA and show that numerical relativity can drastically improve the accuracy of position location for coalescing Super Massive Black Hole (SMBH) binaries. For systems with total redshifted mass above $10^7 M_{odot}$, LISA will mainly see the merger and ring-down of the gravitational wave (GW) signal, which can now be computed numerically using the full Einstein equations. Using numerical waveforms that also include about ten GW cycles of inspiral, we improve inspiral-only position estimates by an order of magnitude. We show that LISA localizes half of all such systems at $z=1$ to better than 3 arcminutes and the best 20% to within one arcminute. This will give excellent prospects for identifying the host galaxy.
We report results of a project investigating the growth of super-massive black holes (BHs) by disk accretion. We find that the BH mass growth is quick enough to account for the inferred masses in the highest-redshift quasars, and the growth time is an inverse function of the final BH mass as seems to be required by recent X-ray surveys.
We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5times 10^8M_{odot}$ black hole with accretion rates varying from $sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $sim 0.1-0.4c$. When the accretion rate is smaller than $sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $sim 5%-7%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1%$. We always find the kinetic energy luminosity associated with the outflow is only $sim 15%-30%$ of the radiative luminosity. The mass flux lost in the outflow is $sim 15%-50%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks.