Search

Supermassive stars (SMSs) and heavy seed black holes, as their remnants, are promising candidates for Supermassive Black Hole (SMBH) progenitors, especially for ones observed in the early universe $ z\simeq 8.5-10$ by recent JWST observations. Expected cradles of SMSs are the atomic cooling halos ($M_{\rm halo}\simeq 10^7~{\rm M_\odot}$), where "cold accretion" emerges and possibly forms SMSs. We perform a suit of cosmological radiation hydrodynamics simulations and investigate star formation after the emergence of cold accretion, solving radiative feedback from stars inside the halo. We follow the mass growth of the protostars for $\sim 3~{\rm Myr}$, resolving the gas inflow down to $\sim 0.1~{\rm pc}$ scales. We discover that, after cold accretion emerges, multiple SMSs of $m_\star \gtrsim 10^5{\rm M_\odot}$ form at the halo centre with the accretion rates maintained at $\dot{m}_\star \simeq 0.04~{\rm M_\odot}~{\rm yr^{-1}}$ for $\lesssim 3~{\rm Myr}$. Cold accretion supplies gas at a rate of $\dot{M}_{\rm gas}\gtrsim 0.01-0.1~{\rm M_\odot}~{\rm yr^{-1}}$ from outside the halo virial radius to the central gas disc. Gravitational torques from spiral arms transport gas further inward, which feeds the SMSs. Radiative feedback from stars suppresses H$_2$ cooling and disc fragmentation, while photoevaporation is prevented by a dense envelope, which attenuates ionising radiation. Our results suggest that cold accretion can bring efficient BH mass growth after seed formation in the later universe. Moreover, cold accretion and gas migration inside the central disc increase the mass concentration and provide a promising formation site for the extremely compact stellar clusters observed by JWST., Comment: 19 pages, 15 figures, accepted by MNRAS

The formation of massive objects via gravitational collapse is relevant both for explaining the origin of the first supermassive black holes and in the context of massive star formation. Here, we analyze simulations of the formation of massive objects pursued by different groups and in various environments, concerning the formation of supermassive black holes, primordial stars, as well as present-day massive stars. We focus particularly on the regime of small virial parameters, i.e., low ratios of the initial kinetic to gravitational energy, low to moderate Mach numbers, and the phase before feedback is very efficient. We compare the outcomes of collapse under different conditions using dimensionless parameters, particularly the star formation efficiency \epsilon_*, the fraction f_* of mass in the most massive object relative to the total stellar mass, and the fraction f_{\rm tot} of mass of the most massive object as a function of the total mass. We find that in all simulations analyzed here, f_{\rm tot} increases as a function of \epsilon_*, although the steepness of the increase depends on the environment. The relation between f_* and \epsilon_* is found to be more complex and also strongly depends on the number of protostars present at the beginning of the simulations. We show that a collision parameter, estimated as the ratio of the system size divided by the typical collision length, allows us to approximately characterize whether collisions will be important. We analyze the statistical correlation between the dimensionless quantities using the Spearman coefficient and confirm via a machine learning analysis that good predictions of f_* can be obtained from \epsilon_* together with a rough estimate of the collision parameter. This suggests that a good estimate of the mass of the most massive object can be obtained once the maximum efficiency for a given environment is known., Comment: 13 pages, 8 figures, 2 tables. Submitted to A&A

The stellar initial mass function (IMF) in the early universe is essential to understand the formation of ancient galaxies. To this end, we conduct a series of long-term radiation hydrodynamic simulations following star cluster formation, varying the metallicity from $Z/Z_\odot = 10^{-4}$ to $1$. We particularly consider the effects of protostellar radiative feedback, which modify the exact shape of the IMF and determine the star formation efficiency (SFE), i.e. the ratio between the mass in stars and the initial gas mass in the parental cloud. Our results show that the IMF changes from a Salpeter-type to a top-heavy function as the metallicity decreases. When $Z/Z_\odot \lesssim 10^{-2}$, the IMF becomes log-flat and distinct from a Salpeter-like IMF. Stellar feedback is effective in shaping both the low- and high-mass ends of the IMF. Heating of dust grains by stellar radiation suppresses small-scale fragmentation and reduces the number of low-mass stars with $M_* \lesssim 1~M_\odot$ at all metallicities. The ionizing radiation hinders the growth of massive stars, steepening the slope of the IMF at the high-mass end. The resulting feedback is more effective at lower metallicity, and star formation is regulated by stellar radiative feedback, with the SFE decreasing with decreasing metallicity. We suggest that the unexpectedly large number of UV-bright galaxies at $z>10$ reported by JWST observations can be explained by considering star cluster formation at $Z/Z_\odot \sim 10^{-2}$ or $10^{-3}$, where the IMF is top-heavy, but the SFE is not too low due to stellar feedback., Comment: 22 pages, 19 figures, submitted to MNRAS

We use the Cosmic Archaeology Tool (CAT) semi-analytical model to explore the contribution of Population (Pop) III/II stars and active galactic nuclei (AGNs) to the galaxy UV luminosity function (LF) evolution at $4 \leq z \leq 20$. We compare in particular with recent JWST data in order to explore the apparent tension between observations and theoretical models in the number density of bright galaxies at $z \gtrsim 10$. The model predicts a star formation history dominated by UV faint ($M_{\rm UV} > - 18$) galaxies, with a Pop III contribution of $\lesssim 10\%$ ($\lesssim 0.5\%$) at $z \simeq 20$ ($z \simeq 10$). Stars are the primary sources of cosmic reionization, with $5 - 10 \%$ of ionizing photons escaping into the intergalatic medium at $5 \leq z \leq 10$, while the contribution of unobscured AGNs becomes dominant only at $z \lesssim 5$. The predicted stellar and AGN UV LFs reproduce the observational data at $5 \lesssim z \lesssim 9 - 10$. At higher redshift, CAT predicts a steeper evolution in the faint-end slope ($M_{\rm UV} > - 18$), and a number density of bright galaxies ($M_{\rm UV} \simeq -20$) consistent with data at $z \sim 10 - 11$, but smaller by 0.8 dex at $z \sim 12 - 13$, and 1.2 dex at $z \sim 14 - 16$, when compared to the values estimated by recent studies. Including the AGN emission does not affect the above findings, as AGNs contribute at most to $\lesssim 10 \%$ of the total UV luminosity at $M_{\rm UV} < - 19$ and $z \gtrsim 10$. Interestingly, considering a gradual transition in the stellar IMF, modulated by metallicity and redshift as suggested by recent simulations, the model agrees with JWST data at $z \sim 12 - 13$, and the disagreement at $z \sim 14 - 16$ is reduced to 0.5 dex., Comment: 20 pages, 12 figures. Submitted to MNRAS

We estimate the fraction of halos that host supermassive black holes (SMBHs) forming through the direct collapse (DC) scenario by using cosmological N -body simulations combined with a semi-analytic model for galaxy evolution. While in most of earlier studies the occurrence of the DC is limited only in chemically pristine halos, we here suppose that the DC can occur also in halos with metallicity below a threshold value $Z_{\rm th} = 0$--$10^{-3}~{\rm Z}_{\bigodot}$, considering the super-competitive accretion pathway for DC black hole (DCBH) formation. In addition, we consider for the first time the effect of Lyman-Werner (LW) radiation from stars within host halos, i.e., internal radiation. We find that, with low threshold metallicities of $Z_{\rm th} \leq 10^{-4}~{\rm Z}_{\bigodot}$, the inclusion of internal radiation rather reduces the number density of DCBHs from $0.2$--$0.3$ to $0.03$--$0.06~{\rm Mpc}^{-3}$. This is because star formation is suppressed due to self-regulation, and the LW flux emitted by neighboring halos is reduced. Only when $Z_{\rm th}$ is as high as $10^{-3}~{\rm Z}_{\bigodot}$, internal radiation enhances the number density of DCBHs from $0.4$ to $1~{\rm Mpc}^{-3}$, thereby decreasing the threshold halo mass above which at least one DCBH forms from $2\times 10^{9}$ to $9\times 10^{8}~{\rm M}_{\bigodot}$. We also find that halos with $M_{\rm halo} \gtrsim 10^{11}$--$10^{12}~{\rm M}_{\bigodot}$ can host more than one DCBH at $z = 0$. This indicates that the DC scenario alone can explain the observed number of SMBH-hosting galaxies., Comment: 17 pages, 12 figures, 2 tables, accepted by MNRAS

We investigate the first emergence of the so-called cold accretion, the accretion flows deeply penetrating a halo, in the early universe with cosmological N-body/SPH simulations. We study the structure of the accretion flow and its evolution within small halos with $\lesssim 10^8~{\rm M}_\odot$ with sufficiently high spatial resolutions down to $\sim 1 \ {\rm pc}$ scale. While previous studies only follow the evolution for a short period after the primordial cloud collapse, we follow the long-term evolution until the cold accretion first appears, employing the sink particle method. We show that the cold accretion emerges when the halo mass exceeds $\sim 2.2\times 10^7 \ {\rm M}_\odot\left\{\left(1+z\right)/15 \right\}^{-3/2}$, the ${\it minimum}$ halo masses above which the accretion flow penetrates halos. We further continue simulations to study whether the cold accretion provides the dense shock waves, which have been proposed to give birth to supermassive stars (SMSs). We find that the accretion flow eventually hits a compact disc near the halo centre, creating dense shocks over a wide area of the disc surface. The resulting post-shock gas becomes dense and hot enough with its mass comparable to the Jeans mass $M_{\rm J}\sim 10^{4-5} \ {\rm M}_\odot$, a sufficient amount to induce the gravitational collapse, leading to the SMS formation., Comment: 17 pages, 17 figures, accepted for publication in MNRAS

H1821+643 is the most X-ray luminous non-beamed AGN of $L_\mathrm{14-150 keV}= 5.2\times 10^{45}$ erg s$^{-1}$ in the Swift/BAT ultra-hard X-ray survey and it is also a hyper-luminous infrared (IR) galaxy $L_\mathrm{IR} = 10^{13.2} L_\odot$ residing in the center of a massive galaxy cluster, which is a unique environment achieving the rapid mass assembly of black holes (BH) and host galaxies in the local universe. We decompose the X-ray to IR spectral energy distribution (SED) into the AGN and starburst component using the SED fitting tool CIGALE-2022.0 and show that H1821+643 consumes a large amount of cold gas ($\dot{M}_\mathrm{con}$) with star-formation rate of $\log ( \mathrm{SFR}/M_{\odot}~\mathrm{yr}^{-1}) = 3.01 \pm 0.04$ and BH accretion rate of $\log (\dot{M}_\mathrm{BH}/M_{\odot}~\mathrm{yr}^{-1}) = 1.20 \pm 0.05$. This high $\dot{M}_\mathrm{con}$ is larger than the cooling rate ($\dot{M}_\mathrm{cool}$) of the intra-cluster medium (ICM), $\dot{M}_\mathrm{con}/\dot{M}_\mathrm{cool} \gtrsim 1$, which is one to two order magnitude higher than the typical value of other systems, indicating that H1821 provides the unique and extreme environment of rapid gas consumption. We also show that H1821+643 has an efficient cooling path achieving from $10^7$ K to $10^2$ K thanks to [OIII] 63 $\mu \mathrm{m}$, which is a main coolant in low temperature range ($10^4$ K to $10^2$ K) with a cooling rate of $\dot{M}_{\mathrm{cool}}=3.2\times 10^5\ M_{\odot}\mathrm{~yr^{-1}}$, and the star-forming region extends over 40 kpc scale., Comment: 20 pages, 8 figures, accepted for publication in ApJ

We study star cluster formation at low metallicities of $Z/Z_\odot=10^{-4}$--$10^{-1}$ using three-dimensional hydrodynamics simulations. Particular emphasis is put on how the stellar mass distribution is affected by the cosmic microwave background radiation (CMB), which sets the temperature floor to the gas. Starting from the collapse of a turbulent cloud, we follow the formation of a protostellar system resolving $\sim$au scale. In relatively metal-enriched cases of $Z/Z_\odot \gtrsim 10^{-2}$, where the mass function resembles the present-day one in the absence of the CMB, high temperature CMB suppresses cloud fragmentation and reduces the number of low-mass stars, making the mass function more top-heavy than in the cases without CMB heating at $z\gtrsim10$. In lower-metallicity cases with $Z/Z_\odot \lesssim 10^{-3}$, where the gas temperature is higher than the CMB value due to inefficient cooling, the CMB has only a minor impact on the mass distribution, which is top-heavy regardless of the redshift. In cases either with a low metallicity of $Z/Z_\odot \lesssim 10^{-2}$ or at a high redshift $z\gtrsim10$, the mass spectrum consists of a low-mass Salpeter-like component, peaking at $0.1~M_\odot$, and a top-heavy component with $10$--$50~M_\odot$, with the fraction in the latter increasing with increasing redshift. In galaxies forming at $z\gtrsim10$, the major targets of the future instruments including JWST, CMB heating makes the stellar mass function significantly top-heavy, enhancing the number of supernova explosions by a factor of $1.4$ ($2.8$) at $z=10$ ($20$, respectively) compared to the prediction by Chabrier initial mass function when $Z/Z_\odot=0.1$., Comment: 18 pages, 12 figures, MNRAS submitted

We study star cluster formation in a low-metallicity environment using three dimensional hydrodynamic simulations. Starting from a turbulent cloud core, we follow the formation and growth of protostellar systems with different metallicities ranging from $10^{-6}$ to $0.1~Z_{\odot}$. The cooling induced by dust grains promotes fragmentation at small scales and the formation of low-mass stars with $M_{*} \sim 0.01$--$0.1~M_{\odot}$ when $Z/Z_{\odot} \gtrsim 10^{-5}$. While the number of low-mass stars increases with metallicity, the stellar mass distribution is still top-heavy for $Z/Z_{\odot} \lesssim 10^{-2}$ compared to the Chabrier initial mass function (IMF). In these cases, star formation begins after the turbulent motion decays and a single massive cloud core monolithically collapses to form a central massive stellar system. The circumstellar disk preferentially feeds the mass to the central massive stars, making the mass distribution top-heavy. When $Z/Z_{\odot}=0.1$, collisions of the turbulent flows promote the onset of the star formation and a highly filamentary structure develops owing to efficient fine-structure line cooling. In this case, the mass supply to the massive stars is limited by the local gas reservoir and the mass is shared among the stars, leading to a Chabrier-like IMF. We conclude that cooling at the scales of the turbulent motion promotes the development of the filamentary structure and works as an important factor leading to the present-day IMF., Comment: 20 pages, 14 figures, accepted for publication in MNRAS

Supermassive stars (SMSs) with $\sim10^{4-5}~\mathrm{M}_{\odot}$ are candidate objects for the origin of supermassive black holes observed at redshift $z$>6. They are supposed to form in primordial-gas clouds that provide the central stars with gas at a high accretion rate, but their growth may be terminated in the middle due to the stellar ionizing radiation if the accretion is intermittent and its quiescent periods are longer than the Kelvin-Helmholtz (KH) timescales at the stellar surfaces. In this paper, we examine the role of the ionizing radiation feedback based on the accretion history in two possible SMS-forming clouds extracted from cosmological simulations, following their evolution with vertically-integrated two-dimensional hydrodynamic simulations with detailed thermal and chemical models. The consistent treatment of the gas thermal evolution is crucial for obtaining the realistic accretion history, as we demonstrate by performing an additional run with a barotropic equation of state, in which the fluctuation of the accretion rate is artificially suppressed. We find that although the accretion becomes intermittent due to the formation of spiral arms and clumps in gravitationally unstable disks, the quiescent periods are always shorter than the KH timescales, implying that SMSs can form without affected by the ionizing radiation., Comment: 14 pages, 11 figures, accepted for publication in MNRAS

The direct collapse (DC) is a promising mechanism that provides massive seed black holes (BHs) with $\sim 10^{5}~M_{\odot}$ in the early universe. To study a long-term accretion growth of a DCBH thus formed, we perform cosmological radiation-hydrodynamics simulations, extending our previous work where we investigated its formation stage. With a high spatial resolution down below the Bondi radius, we show that the accretion rate onto the BH is far below the Eddington value. Such slow mass growth is partly because of the strong radiative feedback from the accreting BH. Moreover, we find that the BH has a large velocity of $\sim 100~{\rm km~s^{-1}}$ relative to the gas after it falls into the first galaxy, which substantially reduces the accretion rate. The latter effect stems from the fact that the DCBHs form in metal-free environments typically at $\sim 1~$kpc from the galaxy. The BH accelerates as it approaches the galactic center due to the gravity. The relative velocity never damps after that, and the BH does not settle down to the galactic center but continues to wander around it. An analytic estimate predicts that the DCBH formation within $\sim 100$~pc around the galactic center is necessary to decelerate the BH with dynamical friction before $z=7$. Since metal enrichment with $Z \sim 10^{-5} - 10^{-3}~Z_\odot$ is expected in such a case, the formation of DCBHs in the metal-poor environments is preferable for the subsequent rapid growth., Comment: 13 pages, 8 figures, submitted to MNRAS

Direct collapse black hole (DCBH) formation with mass $\gtrsim 10^{5}~M_{\odot}$ is a promising scenario for the origin of high-redshift supermassive black holes. It has usually been supposed that the DCBH can only form in the primordial gas since the metal enrichment enhances the cooling ability and causes the fragmentation into smaller pieces. What actually happens in such an environment, however, has not been explored in detail. Here, we study the impact of the metal enrichment on the clouds, conducting hydrodynamical simulations to follow the cloud evolution in cases with different degree of the metal enrichment $Z/Z_{\odot}=10^{-6}-10^{-3}$. Below $Z/Z_{\odot}=10^{-6}$, metallicity has no effect and supermassive stars form along with a small number of low-mass stars. With more metallicity $Z/Z_{\odot} \gtrsim 5 \times 10^{-6}$, although the dust cooling indeed promotes fragmentation of the cloud core and produces about a few thousand low-mass stars, the accreting flow preferentially feeds the gas to the central massive stars, which grows supermassive as in the primordial case. We term this formation mode as the {\it super competitive accretion}, where only the central few stars grow supermassive while a large number of other stars are competing for the gas reservoir. Once the metallicity exceeds $10^{-3}~Z_{\odot}$ and metal-line cooling becomes operative, the central star cannot grow supermassive due to lowered accretion rate. Supermassive star formation by the super competitive accretion opens up a new window for seed BHs, which relaxes the condition on metallicity and enhances the seed BH abundance., Comment: 10 pages, 6figures, submitted to MNRAS

The disk fragmentation is a possible process leading to the formation of Population III stellar binary systems. However, numerical simulations show diverse fates of the fragments; some evolve into stable binaries and others merge away with a central star. To clarify the physics behind such diversity, we perform a series of three dimensional hydrodynamics simulations in a controlled manner. We insert a point particle mimicking a fragment in a self-gravitating disk, where the initial mass and position are free parameters, and follow the orbital evolution for several tens of orbits. The results show great diversity even with such simple experiments. Some particles shortly merge away after migrating inward, but others survive as the migration stalls with the gap-opening in the disk. We find that our results are well interpreted postulating that the orbital angular momentum is extracted by (i) the gravitational torque from the disk spiral structure, and (ii) tidal disruption of a gravitationally-bound envelope around the particle. Our analytic evaluations show the processes (i) and (ii) are effective in an outer and inner part of the disk respectively. There is a window of the gap-opening in the middle, if the envelope mass is sufficiently large. These all agree with our numerical results. We further show that the binaries, which appear for the "survival" cases, gradually expand while accreting the disk gas. Our theoretical framework is freely scalable to be applied for the present-day star and planet formation., Comment: 16 pages, 14 figures, published in MNRAS

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes.", Comment: Solicited review article (accepted) for Publications of the Astronomical Society of Australia. 38 pages, 15 figures, 1 table

Formation of supermassive stars (SMSs) with mass ~10^4 Msun is a promising pathway to seed the formation of supermassive black holes in the early universe. The so-called direct-collapse (DC) model postulates that such an SMS forms in a hot gas cloud irradiated by a nearby star-forming galaxy. We study the DC SMS formation in a fully cosmological context using three-dimensional radiation hydrodynamics simulations. We initialize our simulations using the outputs of the cosmological simulation of Chon et al. (2016), where two DC gas clouds are identified. The long-term evolution over a hundred thousand years is followed from the formation of embryo protostars through their growth to SMSs. We show that the strength of the tidal force by a nearby galaxy determines the multiplicity of the formed stars and affects the protostellar growth. In one case, where a collapsing cloud is significantly stretched by strong tidal force, multiple star-disk systems are formed via filament fragmentation. Small-scale fragmentation occurs in each circumstellar disk, and more than 10 stars with masses of a few times 10^3 Msun are finally formed. Interestingly, about a half of them are found as massive binary stars. In the other case, the gas cloud collapses nearly spherically under a relatively weak tidal field, and a single star-disk system is formed. Only a few SMSs with masses ~ 10^4 Msun are found already after evolution of a hundred thousand years, and the SMSs are expected to grow further by gas accretion and to leave massive blackholes at the end of their lives., Comment: 19 pages, 20 figures, Submitted to MNRAS

A massive primordial halo near an intensely star-forming galaxy may collapse into a supermassive star (SMS) and leave a massive black hole seed of about $10^5~M_{\odot}$. To investigate the impact of ionizing radiation on the formation of an SMS from a nearby galaxy, we perform three-dimensional radiation hydrodynamical simulations by selecting a pair of massive dark matter halos forming at $z >10$. We find that rich structures such as clumps and filaments around the source galaxy shield the cloud from ionizing radiation. In fact, in some cases, cloud collapse is accelerated under ionizing radiation. This fact suggests that the ionization of the cloud's surroundings helps its collapse. Only strong radiation at the early stage of structure formation can halt the cloud collapse, but this is much stronger than observationally allowed value. We also explored the effect of ionizing radiation on a sample of 68 halos by employing an analytical model and found that increase in the mean density of the gas between the SMS forming cloud and the source galaxy protects the gas cloud from ionizing radiation as they approach each other. Thus, we conclude that ionizing radiation does not prevent the formation of an SMS in most of the cases., Comment: 11 pages, 9 figures, to be submitted to MNRAS

Gravitational collapse of a massive primordial gas cloud is thought to be a promising path for the formation of supermassive blackholes in the early universe. We study conditions for the so-called direct collapse (DC) blackhole formation in a fully cosmological context. We combine a semianalytic model of early galaxy formation with halo merger trees constructed from dark matter $N$-body simulations. We locate a total of 68 possible DC sites in a volume of $20\;h^{-1}\;\mathrm{Mpc}$ on a side. We then perform hydrodynamics simulations for 42 selected halos to study in detail the evolution of the massive clouds within them. We find only two successful cases where the gas clouds rapidly collapse to form stars. In the other cases, gravitational collapse is prevented by the tidal force exerted by a nearby massive halo, which otherwise should serve as a radiation source necessary for DC. Ram pressure stripping disturbs the cloud approaching the source. In many cases, a DC halo and its nearby light source halo merge before the onset of cloud collapse. Only when the DC halo is assembled through major mergers, the gas density increases rapidly to trigger gravitational instability. Based on our cosmological simulations, we conclude that the event rate of DC is an order of magnitude smaller than reported in previous studies, although the absolute rate is still poorly constrained. It is necessary to follow the dynamical evolution of a DC cloud and its nearby halo(s) in order to determine the critical radiation flux for DC., Comment: 22 pages, 22 figures, accepted for publication in ApJ

Starting from the geometrical concept of a 4-dimensional de-Sitter configuration of spheres in Euclidean 3-space and modelling voids in the Universe as spheres, we show that a uniform distribution over this configuration space implies a power-law for the void number density which is consistent with results from the excursion set formalism and with data, for an intermediate range of void volumes. The scaling dimension of the large scale structure can be estimated as well. We also discuss the effect of restricting the survey geometry on the void statistics. This work is a new application of de-Sitter geometry to cosmology and also provides a new geometrical perspective on self-similarity in cosmology., Comment: 8 pages, 4 figures, accepted by MNRAS. Minor changes, appendix added