The rapid assembly of the massive black holes that power the luminous quasars observed at z ∼ 6–7 remains a puzzle. Various direct collapse models have been proposed to head-start black hole growth ...from initial seeds with masses ∼105 M⊙, which can then reach a billion solar mass while accreting at the Eddington limit. Here, we propose an alternative scenario based on radiatively inefficient supercritical accretion of stellar-mass holes embedded in the gaseous circumnuclear discs (CNDs) expected to exist in the cores of high-redshift galaxies. Our sub-pc resolution hydrodynamical simulations show that stellar-mass holes orbiting within the central 100 pc of the CND bind to very high density gas clumps that arise from the fragmentation of the surrounding gas. Owing to the large reservoir of dense cold gas available, a stellar-mass black hole allowed to grow at super-Eddington rates according to the ‘slim-disc’ solution can increase its mass by three orders of magnitudes within a few million years. These findings are supported by simulations run with two different hydro codes, ramses based on the Adaptive Mesh Refinement technique and gizmo based on a new Lagrangian Godunov-type method, and with similar, but not identical, sub-grid recipes for star formation, supernova feedback, black hole accretion and feedback. The low radiative efficiency of supercritical accretion flows are instrumental to the rapid mass growth of our black holes, as they imply modest radiative heating of the surrounding nuclear environment.
Missing [C ii] emission from early galaxies Carniani, S; Ferrara, A; Maiolino, R ...
Monthly notices of the Royal Astronomical Society,
12/2020, Letnik:
499, Številka:
4
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
ALMA observations have revealed that C ii 158 μm line emission in high-z galaxies is ≈2–3 × more extended than the UV continuum emission. Here we explore whether surface brightness dimming ...(SBD) of the C ii line is responsible for the reported C ii deficit, and the large $L_{\rm O\, \small {III}}/L_{\rm C\, \small {II}}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z > 6 galaxies observed in both C ii and O iii. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect C ii emission in the whole sample, with an extent systematically larger than the O iii emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40 per cent of the extended C ii component might be missed at an angular resolution of 0.8 arcsec, implying that $L_{\rm C\, \small {II}}$ is underestimated by a factor ≈2 in data at low (<7) signal-to-noise ratio. By combining these results, we conclude that $L_{\rm C\, \small {II}}$ of z > 6 galaxies lies, on average, slightly below the local $L_{\rm C\, \small {II}}-\mathrm{ SFR}$ relation (Δz = 6–9 = −0.07 ± 0.3), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm O\, \small {III}}/L_{\rm C\, \small {II}}\lt 10$, i.e. more in line with current hydrodynamical simulations.
In this paper, we explore a possible route of black hole seed formation that appeals to a model by Davies, Miller & Bellovary who considered the case of the dynamical collapse of a dense cluster of ...stellar black holes subjected to an inflow of gas. Here, we explore this case in a broad cosmological context. The working hypotheses are that (i) nuclear star clusters form at high redshifts in pre-galactic discs hosted in dark matter haloes, providing a suitable environment for the formation of stellar black holes in their cores, (ii) major central inflows of gas occur on to these clusters due to instabilities seeded in the growing discs and/or to mergers with other gas-rich haloes and (iii) following the inflow, stellar black holes in the core avoid ejection due to the steepening to the potential well, leading to core collapse and the formation of a massive seed of ≲ 1000 M⊙. We simulate a cosmological box tracing the build-up of the dark matter haloes and their embedded baryons, and explore cluster evolution with a semi-analytical model. We show that this route is feasible, peaks at redshifts z ≲ 10 and occurs in concomitance with the formation of seeds from other channels. The channel is competitive relative to others, and is independent of the metal content of the parent cluster. This mechanism of gas-driven core collapse requires inflows with masses at least 10 times larger than the mass of the parent star cluster, occurring on time-scales shorter than the evaporation/ejection time of the stellar black holes from the core. In this respect, the results provide upper limit to the frequency of this process.
ABSTRACT
We study the effect of stellar feedback (photodissociation/ionization, radiation pressure, and winds) on the evolution of a Giant Molecular Cloud (GMC), by means of a 3D radiative transfer, ...hydrosimulation implementing a complex chemical network featuring H2 formation and destruction. We track the formation of individual stars with mass $M\gt 1\, {\rm M}_{\odot }$ with a stochastic recipe. Each star emits radiation according to its spectrum, sampled with 10 photon bins from near-infrared to extreme ultraviolet bands; winds are implemented by energy injection in the neighbouring cells. We run a simulation of a GMC with mass $M=10^5\, {\rm M}_{\odot }$, following the evolution of different gas phases. Thanks to the simultaneous inclusion of different stellar feedback mechanisms, we identify two stages in the cloud evolution: (1) radiation and winds carve ionized, low-density bubbles around massive stars, while FUV radiation dissociates most H2 in the cloud, apart from dense, self-shielded clumps; (2) rapid star formation (SFR$\simeq 0.1\, {\rm M}_{\odot }\, {\rm yr}^{-1}$) consumes molecular gas in the dense clumps, so that UV radiation escapes and ionizes the remaining $\mathrm{H\,{\small I}}$ gas in the GMC. H2 is exhausted in 1.6 Myr, yielding a final star formation efficiency of 36 per cent. The average intensity of FUV and ionizing fields increases almost steadily with time; by the end of the simulation (t = 2.5 Myr) we find 〈G0〉 ≃ 103 (in Habing units), and a ionization parameter 〈Uion〉 ≃ 102, respectively. The ionization field has also a more patchy distribution than the FUV one within the GMC. Throughout the evolution, the escape fraction of ionizing photons from the cloud is fion, esc ≲ 0.03.
Disentangling the different stages of the star formation process, in particular in the high-mass regime, is a challenge in astrophysics. Chemical clocks could help alleviate this problem, but their ...evolution strongly depends on many parameters, leading to degeneracy in the interpretation of the observational data. One of these uncertainties is the degree of CO depletion. We present here the first self-consistent magneto-hydrodynamic simulations of high-mass, star-forming regions at different scales, fully coupled with a nonequilibrium chemical network, which includes C-N-O bearing molecules. Depletion and desorption processes are treated time dependently. The results show that full CO depletion (i.e., all gas-phase CO frozen-out on the surface of dust grains) can be reached very quickly, in one-third or even smaller fractions of the freefall time, whether the collapse proceeds on slow or fast timescales. This leads to a high level of deuteration in a short time, both for typical tracers like N2H+, as well as for the main ion H 3 + , the latter being in general larger and more extended. N2 depletion is slightly less efficient, and no direct effects on N-bearing molecules and deuterium fractionation are observed. We show that CO depletion is not the only driver of deuteration, and that there is a strong impact on Dfrac when changing the grain size. We finally apply a two-dimensional Gaussian point-spread function to our results to mimic observations with single-dish and interferometers. Our findings suggest that the low-values observed in high-mass star-forming clumps are in reality masking a full-depletion stage in the inner 0.1 pc region.
Abstract
Black holes with masses of $\rm 10^6\text{-}10^9\,M_{\odot }$ dwell in the centres of most galaxies, but their formation mechanisms are not well known. A subdominant dissipative component of ...dark matter with similar properties to the ordinary baryons, known as mirror dark matter, may collapse to form massive black holes during the epoch of first galaxies formation. In this study, we explore the possibility of massive black hole formation via this alternative scenario. We perform three-dimensional cosmological simulations for four distinct haloes and compare their thermal, chemical, and dynamical evolution in both the ordinary and the mirror sectors. We find that the collapse of haloes is significantly delayed in the mirror sector due to the lack of $\rm H_2$ cooling and only haloes with masses above $\rm \ge\!10^7\, M_{\odot }$ are formed. Overall, the mass inflow rates are $\rm \ge\!10^{-2}\,M_{\odot }\,yr^{ -1}$ and there is less fragmentation. This suggests that the conditions for the formation of massive objects, including black holes, are more favourable in the mirror sector.
ABSTRACT
We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z = 6.14 ...thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified ‘delayed cooling’ feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n ∼ 103–105 cm−3). The sites where star formation ignites are characterized by maximum densities of the order of 105 cm−3 and gravitational pressures Pgrav/k >107 K cm−3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z = 6.14 is 3.4$\times 10^7~\rm M_{\odot }$, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of $\sim 10^6~\rm M_{\odot }$ and half-mass sizes of ∼100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one orders of magnitude lower than the observed one. In the size–mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.
ABSTRACT
Cosmic rays are a global source of ionization, and the ionization fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, and affect the ...chemistry and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic ray ionization rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from ${\mathrm{ H}_3}^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic ray ionization rate from observational quantities. The theoretical predictions are validated with high-resolution 3D numerical simulations and applied to the well-known core L1544; we obtained an estimate of ζ2 ∼ 2–3 × 10−17 s−1. Our results and the analytical formulae provided represent the first model-independent robust tool to probe the cosmic ray ionization rate in the densest part of star-forming regions (on spatial scales of R ≤ 0.05 pc). An error analysis is presented to give statistical relevance to our study.
Context . Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, the characterisation of their ionisation rate continues to pose a challenge from ...the observational point of view. Aims . In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionisation rate, ζ 2 , from molecular line observations. These have been derived from the chemical kinetics of the involved species, but they have not yet been validated using synthetic data processed with a standard observative pipeline. In this work, we aim to bridge this gap. Methods . We performed a radiative transfer on a set of three-dimensional magneto-hydrodynamical simulations of prestellar cores, exploring different initial ζ 2 , evolutionary stages, types of radiative transfer (for instance assuming local-thermodynamic-equilibrium conditions), and telescope responses. We then computed the column densities of the involved tracers to determine ζ 2 , employing a recently proposed method based on the detection of H 2 D + . We compared this approach with a previous method, based on more common tracers. Both approaches are commonly used. Results . Our results confirm that the equation based on the detection of H 2 D + accurately retrieves the actual ζ 2 within a factor of two to three in the physical conditions explored in our tests. Since we have also explored a non-local thermodynamic equilibrium (non-LTE) radiative transfer, this work indirectly offers insights into the excitation temperatures of common transitions at moderate volume densities (n ≈ 10 5 cm −3 ). We also performed a few tests using a previous methodology that is independent of H 2 D + , which overestimates the actual ζ 2 by at least two orders of magnitude. We considered a new derivation of this method, however, we found that it still leads to high over-estimations. Conclusions . The method based on H 2 D + is further validated in this work and demonstrates a reliable method for estimating ζ 2 in cold and dense gas. On the contrary, the former analytical equation, as already pointed out by its authors, has no global domain of application. Thus, we find that it ought to be employed with caution.
ABSTRACT
We use zoom-in, hydrodynamical, cosmological N-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged ...project, to study key structural properties of dark matter haloes when the Universe was only $0.92\, {\rm Gyr}$ old. The very high resolution (maximum physical resolution $0.3\, {h}^{-1}\, {\rm pc}$ at z = 6.14, smallest dark matter particle mass $164\, {\rm M}_{\odot }$) allows us to reach the very low mass end of the stellar-to-halo mass relation ($M_{\rm vir}=10^{7.5{\!-\!}9.5}\, {\rm M}_{\odot }$) to study the processes that mould dark matter haloes during the first stages of structure formation. We investigate the role of baryonic cooling and stellar feedback, modelled from individual stars, in shaping haloes, and of environmental effects as accretion of dark matter along cosmic filaments and mergers. We find that the onset of star formation (typically for $\log M_{\rm vir}/\, {\rm M}_{\odot }\simeq 7.6$) causes the inner cusp in the haloes’ density profile to flatten into a core with constant density and size proportionally to the halo virial mass. Even at these mass scales, we confirm that baryons make haloes that have formed stars rounder in the central regions than haloes that have not formed stars yet, with median minor-to-major 〈q〉 and intermediate-to-major 〈s〉 axes 0.66 and 0.84, respectively. Our morphological analysis shows that, at z = 6.14, haloes are largely prolate in the outer parts, with the major axis aligned along filaments of the cosmic web or towards smaller sub-haloes, with the degree of elongation having no significant dependence on the halo mass.