Supersonic Gas Streams Enhance the Formation of Massive Black Holes in the Early Universe

The origin of super-massive black holes in the early universe remains poorly understood.Gravitational collapse of a massive primordial gas cloud is a promising initial process,but theoretical studies have difficulty growing the black hole fast enough.We report numerical simulations of early black hole formation starting from realistic cosmological conditions.Supersonic gas motions left over from the Big Bang prevent early gas cloud formation until rapid gas condensation is triggered in a proto-galactic halo. A protostar is formed in the dense, turbulent gas cloud, and it grows by sporadic mass accretion until it acquires 34,000 solar masses.The massive star ends its life with a catastrophic collapse to leave a black hole -- a promising seed for the formation of a monstrous black hole.Comment: Published in Science, combined with updated SOM, additional images and movies are available at http://www-utap.phys.s.u-tokyo.ac.jp/naoki.yoshida/Blackhole/0929e.htm

Low density, radiatively inefficient rotating-accretion flow onto a black hole

We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $\alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $r\sim R_{\rm B}$, where the density follows $\rho \propto (1+R_{\rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $\rho \propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($\dot{M}\propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $\dot{M}_{\rm B}$ The net accretion rate depends on the viscous strength, following $\dot{M}/\dot{M}_{\rm B}\propto (\alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $\dot{M}_{\rm B}/\dot{M}_{\rm Edd}< 10^{-3}$ having minimal radiation cooling, where $\dot{M}_{\rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{\rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $\dot{M}/\dot{M}_{\rm B} \sim 10^{-3}-10^{-2}$. This rate is as low as $\dot{M}/\dot{M}_{\rm Edd} \sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.Comment: 16 pages, 17 figures, published in MNRA

Radiative magneto-hydrodynamics in massive star formation and accretion disks

We briefly overview our newly developed radiation transport module for MHD simulations and two actual applications. The method combines the advantage of the speed of the Flux-Limited Diffusion approximation and the high accuracy obtained in ray-tracing methods.Comment: 2 pages, 1 figure, Proceedings of the IAU Symposium 259, Cosmic Magnetic Fields: From Planets, to Stars and Galaxie

Hydrodynamics of embedded planets' first atmospheres - III. The role of radiation transport for super-Earth planets

The population of close-in super-Earths, with gas mass fractions of up to 10% represents a challenge for planet formation theory: how did they avoid runaway gas accretion and collapsing to hot Jupiters despite their core masses being in the critical range of $M_\mathrm{c} \simeq 10 M_\mathrm{\oplus}$? Previous three-dimensional (3D) hydrodynamical simulations indicate that atmospheres of low-mass planets cannot be considered isolated from the protoplanetary disc, contrary to what is assumed in 1D-evolutionary calculations. This finding is referred to as the recycling hypothesis. In this Paper we investigate the recycling hypothesis for super-Earth planets, accounting for realistic 3D radiation hydrodynamics. Also, we conduct a direct comparison in terms of the evolution of the entropy between 1D and 3D geometries. We clearly see that 3D atmospheres maintain higher entropy: although gas in the atmosphere loses entropy through radiative cooling, the advection of high entropy gas from the disc into the Bondi/Hill sphere slows down Kelvin-Helmholtz contraction, potentially arresting envelope growth at a sub-critical gas mass fraction. Recycling, therefore, operates vigorously, in line with results by previous studies. However, we also identify an "inner core" -- in size $\approx$ 25% of the Bondi radius -- where streamlines are more circular and entropies are much lower than in the outer atmosphere. Future studies at higher resolutions are needed to assess whether this region can become hydrodynamically-isolated on long time-scales.Comment: 16 pages, 12 figures, accepted for publication at MNRA