{"id":850027,"date":"2026-09-16T17:54:12","date_gmt":"2026-09-16T17:54:12","guid":{"rendered":"https:\/\/www.newsbeep.com\/us\/850027\/"},"modified":"2026-09-16T17:54:12","modified_gmt":"2026-09-16T17:54:12","slug":"overmassive-black-holes-and-little-red-dots-naturally-form-in-simulations","status":"publish","type":"post","link":"https:\/\/www.newsbeep.com\/us\/850027\/","title":{"rendered":"Overmassive black holes and little red dots naturally form in simulations"},"content":{"rendered":"

Our cosmological initial conditions are based on the simulation Phi-4096 of ref.\u200951<\/a>. In this work, a dark-matter-only N-body simulation was performed in a comoving box of side length 16\u2009h\u22121\u2009Mpc. The initial condition is generated by MUSIC52<\/a> at a redshift of 127. The cosmological parameters used follow the latest measurement by Planck53<\/a>: \u03a9m\u2009=\u20090.31, \u03a9\u039b\u2009=\u20090.69, \u03a9b\u2009=\u20090.048, h\u2009=\u20090.68, ns\u2009=\u20090.96 and \u03c38\u2009=\u20090.83. The simulation used 4,0963 dark-matter particles, corresponding to a particle mass of 5.13\u2009\u00d7\u2009103\u2009h\u22121\u2009M\u2299. This high resolution allows the identification of mini-halos with masses down to 105\u2009h\u22121\u2009M\u2299, sufficient to resolve the sites of Population\u2009III (hereafter Pop\u2009III) star formation and to track chemical enrichment across cosmic time. Halo merger trees were constructed from simulation snapshots between z=\u200935 and z\u2009=\u20097.5 using the ROCKSTAR phase space halo\/subhalo finder54<\/a> and the CONSISTENT-TREES merger tree code55<\/a>. On top of these, a semi-analytic galaxy-formation model was implemented to follow gas cooling, Pop\u2009II and Pop\u2009III star formation, supernova feedback, metal enrichment and local Lyman\u2013Werner (LW) radiation fields17<\/a>. We use the result of their model with no baryon streaming motion (0\u03c3vbc).<\/p>\n

The details of the semi-analytic model largely follow those described in ref.\u200956<\/a> with recent updates of the treatment of early star formation57<\/a>. Star formation in Pop\u2009II halos follows the standard prescription used in galaxy-formation models, in which the cold-gas component is converted into stars on a timescale tSF\u2009=\u2009tdyn\/\u03b1* with efficiency \u03b1*\u2009=\u20090.03 (ref.\u200958<\/a>). Along with star formation, the model assumes that the rate at which cold gas is reheated into the hot-gas phase by associated supernova feedback is proportional to the star-formation rate, \\(\\gamma {\\dot{M}}_{* }\\), in which \u03b3\u2009=\u2009(Vhalo\/110\u2009km\u2009s\u22121)\u22121.74 (refs.\u200956<\/a>,59<\/a>). When the progenitor halos never formed any stars, the stellar population follows a Pop\u2009III initial stellar mass function (IMF). The characteristic Pop\u2009III stellar mass is determined by the halo mass growth rate, which correlates with the final stellar mass60<\/a>. Stellar populations emit LW radiation that photodissociates molecular hydrogen and delays star formation in nearby halos. The LW radiation intensity is calculated separately for Pop\u2009II and Pop\u2009III stars61<\/a>: <\/p>\n

$${J}_{21,{\\rm{III}}}=\\sum _{i}15{\\left(\\frac{{r}_{i}}{1{\\rm{kpc}}}\\right)}^{-2}\\left(\\frac{{M}_{{\\rm{PopIII}},i}}{1000\\,{M}_{\\odot }}\\right),$$<\/p>\n

\n (1)\n <\/p>\n

$${J}_{21,{\\rm{II}}}=\\sum _{i}3{\\left(\\frac{{r}_{i}}{1{\\rm{kpc}}}\\right)}^{-2}\\left(\\frac{{M}_{{\\rm{PopII}},i}}{1000\\,{M}_{\\odot }}\\right),$$<\/p>\n

\n (2)\n <\/p>\n

in which ri is the distance to halo i and MPopIII,i and MPopII,i are the masses of Pop\u2009III and Pop\u2009II stars formed within the past 5\u2009Myr in halo i, respectively. The sums run over all halos that host Pop\u2009III or Pop\u2009II star formation. The coefficient in the Pop\u2009II expression is derived assuming a Scalo IMF (ref.\u200962<\/a>) and a stellar metallicity of Z\u2009=\u20090.001 (refs.\u200963<\/a>,64<\/a>).<\/p>\n

The LW background is primarily contributed by Pop\u2009II stars and its amplitude depends on assumptions about the IMF and the stellar models used. We note that uncertainties in the LW intensity modelling have only a small impact on our results, as we later test explicitly. A halo begins to form Pop\u2009III stars once its mass exceeds the critical threshold determined by the local LW intensity65<\/a>. Under strong LW irradiation, star formation is delayed until the halo virial temperature reaches about 8,000\u2009K, at which point Ly\u03b1 cooling triggers rapid collapse. Such halos typically experience high mass accretion rates and host more massive Pop\u2009III stars, extending to supermassive stars with M*\u2009\u2248\u2009105\u2009M\u2299, which subsequently collapse into heavy BH seeds. Indeed, Ishiyama and Hirano17<\/a> identified more than about 104 supermassive stars with masses larger than 105\u2009M\u2299 as heavy-seed candidates in their simulation volume.<\/p>\n

From this dataset, we select one representative halo predicted to form a massive seed, which serves as the target of our follow-up radiation-hydrodynamic simulations described in this work. The halo is chosen according to two criteria: (1) the local LW intensity at the time of seed formation exceeds the critical value of J21,crit\u2009=\u20091,000 and (2) the halo is not tidally disrupted by nearby massive galaxies. Here J21 denotes the FUV intensity normalized to 10\u221221\u2009erg\u2009s\u22121\u2009Hz\u22121\u2009cm\u22122\u2009sr\u22121. The second condition is necessary because candidate halos are often located near luminous galaxies that provide strong LW irradiation, but their gravitational collapse may be suppressed by external tidal fields66<\/a>. To exclude such cases, we select halos whose distances from the nearest massive galaxy satisfy rdist\u2009>\u200910dtidal, in which dtidal is defined as the separation at which the tidal radius imposed by the neighbouring galaxy equals the virial radius of the candidate halo. Within the simulated volume, we identify 60 halos that meet both criteria. Among them, we focus on the one that exhibits the most rapid mass growth, reaching a virial temperature of Tvir\u2009\u2243\u20098\u2009\u00d7\u2009103\u2009K, which marks the onset of atomic cooling and subsequent collapse.<\/p>\n

We trace back the initial position of the selected candidate halo to the cosmological initial conditions at redshift z\u2009=\u2009127. Extended Data Fig.\u20091<\/a> shows the time evolution of the halo mass, in which MBH1 forms. The white star marks the redshift at which the virial temperature of the halo reaches 8,000\u2009K, when it is identified as a DCBH halo. A zoom-in region is defined to be 40 times larger than the Lagrangian radius of the target halo, corresponding to a comoving size of about 400\u2009kpc. We confirm that this region is sufficiently large to encompass all material relevant to the formation of the heavy-seed BH in the selected halo. We measured the overdensity of the specified zoom-in region by computing the mean density within its bounding box of about 400 comoving kpc and found it to correspond to a 3.8\u03c3 fluctuation.<\/p>\n

Radiation-hydrodynamic calculation<\/p>\n

The radiation-hydrodynamic calculation is performed within this zoom-in region using the moving-mesh code AREPO (ref.\u200916<\/a>), extended to include prescriptions for star formation and BH accretion physics. The resolution of the zoom-in region at the initial snapshot is identical to that used in ref.\u200951<\/a>, in which the dark-matter particle and baryonic cell masses are 4.33\u2009\u00d7\u2009103 and 7.94\u2009\u00d7\u2009102\u2009h\u22121\u2009M\u2299, respectively. At this resolution, star formation within mini-halos can be reliably resolved. Adaptive mesh refinement is applied whenever the local cell size falls below 16 times the local Jeans length, ensuring that gravitational collapse is properly captured and that artificial fragmentation is avoided67<\/a>. Unlike the semi-analytic model in ref.\u200917<\/a>, the formation and properties of primordial stars are directly followed by the high-resolution radiation-hydrodynamic simulations described here. Uncertainties in the stellar mass and the resulting FUV background intensity in the semi-analytic model may affect whether a heavy-seed BH forms or not, but we later show that these uncertainties do not substantially alter our main conclusion about the rapid emergence of overmassive BHs.<\/p>\n

We solve a non-equilibrium primordial chemical network consisting of eight species: e\u2212, H, H+, H2, H\u2212, D, D+ and HD. The network and associated cooling processes follow the implementation in ref.\u200968<\/a> and include H2 rovibrational line cooling, atomic hydrogen line cooling (including Ly\u03b1) and free\u2013free and free\u2013bound emission of H and H\u2212. The ionization of H, the dissociation of H2 by LW radiation through the Solomon process and the photodetachment of H\u2212 are also taken into account. The chemical network and reaction rates used are similar to those in the minimal chemistry model presented in ref.\u200969<\/a>, except that we neglect \\({{\\rm{H}}}_{2}^{+}\\) and He. \\({{\\rm{H}}}_{2}^{+}\\) can catalyse H2 formation in primordial gas but it is less important in environments exposed to the background radiation field assumed in this study, namely a 105\u2009K blackbody spectrum69<\/a>. Heating and chemical reactions induced by FUV, extreme ultraviolet and X-ray radiation from nearby stars and BHs are also included, as described later.<\/p>\n

Helium chemistry becomes important at gas temperatures of several 104\u2013105\u2009K, at which it provides further cooling channels and increases the free-electron fraction. Although we do not solve the full non-equilibrium helium chemistry, we estimate the ionization state of helium by assuming chemical equilibrium and include the line-cooling rates of He\u2009I and He\u2009II. This effect is relevant for gas irradiated by both Pop\u2009III stars and accreting BHs. Radiation from the accreting BH raises the gas temperature well above 104\u2009K. The elevated temperature already increases the electron fraction and thereby promotes H2 formation; including non-equilibrium helium chemistry would further enhance the H2 abundance to some extent. However, this does not affect our main conclusion. Even if the enhanced H2 abundance lowers the cloud temperature, the large-scale inflow can be sustained. This is demonstrated by the formation of MBH2, in which a heavy-seed BH forms under a massive inflow despite efficient H2 cooling that lowers the gas temperature below 1,000\u2009K (Extended Data Fig.\u20097<\/a>).<\/p>\n

To avoid numerical contamination from the coarse outer region, gas cells outside the zoom-in volume are kept at fixed resolution and both radiative cooling and refinement are disabled. In this calculation, we treat the FUV radiation from the source halos separately from the radiation emitted by stars and BHs formed in the target region. Here we define source halos as halos that have already been labelled as star-forming when the target halo satisfies the DCBH formation criteria in the semi-analytic model. These halos host star-forming galaxies and provide the strong FUV radiation required for DCBH formation. Instead of explicitly following the star formation within the galaxies, we use a time-dependent LW intensity derived from the semi-analytic model in ref.\u200917<\/a>, which self-consistently evolves the star formation and feedback in the same cosmological volume. The resulting LW flux is applied uniformly to all gas cells within the zoom-in region as a uniform background field. We separately solve the ionizing radiation and X-rays emitted by stars and BHs formed inside the zoom-in region, but outside the source halos, using the ray-tracing scheme implemented in regularized smoothed particle hydrodynamic70<\/a>,71<\/a>. We include photoionization of H and the associated photoheating, photodissociation of H2 and HD and photodetachment of H\u2212. The photodetachment rate of H\u2212 is calculated under the optically thin approximation, whereas the other photoreaction rates are calculated using regularized smoothed particle hydrodynamics. We do not include the secondary ionization and heating. The luminosities and spectra of these sources are described later in the section \u2018Modelling Pop\u2009III stars and BHs\u2019.<\/p>\n

Extended Data Fig.\u20092<\/a> shows the time evolution of the LW intensity contributed by nearby source galaxies. The intensity increases rapidly and exceeds J21\u2009=\u20091,000 around redshift z\u2009\u2243\u200918, the critical level required to suppress H2 cooling and enable heavy-seed BH formation (for example, refs.\u200972<\/a>,73<\/a>,74<\/a>,75<\/a>). The spectrum of the background radiation is modelled as a blackbody with an effective temperature of TBB\u2009=\u2009105\u2009K, which reproduces the relative rates of H2 photodissociation and H\u2212 photodetachment expected from young star-forming galaxies73<\/a>. We consider the self-shielding of the background LW radiation using the prescription described in ref.\u200976<\/a>. We assume that the cloud extends over the local Jeans length, lJeans, and estimate the H2 column density as \\({N}_{{{\\rm{H}}}_{2}}={n}_{{{\\rm{H}}}_{2}}{l}_{{\\rm{Jeans}}}\\), in which \\({n}_{{{\\rm{H}}}_{2}}\\) is the number density of H2 in the cell. We neglect the contribution of ionizing photons from neighbouring galaxies because their mean free path is short and they are strongly attenuated by the surrounding intergalactic medium71<\/a>.<\/p>\n

During the radiation-hydrodynamic simulation, we do not assume any specific IMF but instead resolve individual star-formation events directly using a sink-particle method. A sink particle is introduced once the local gas density exceeds 2\u2009\u00d7\u2009106\u2009cm\u22123, with an accretion radius set to ten times the local mesh size. The sink accretes surrounding gas following the local gravitational potential and its properties are updated accordingly at each time step. We allow mergers of the sink particles once the distance of two sink particles becomes smaller than the sum of the sink radii.<\/p>\n

Modelling Pop\u2009III stars and BHs<\/p>\n

The luminosity and effective temperature of each protostar are determined from its instantaneous mass and accretion rate by interpolating results from one-dimensional stellar evolution models assuming constant accretion rates77<\/a>. If the stellar accretion rate exceeds the critical accretion rate of \\({\\dot{M}}_{{\\rm{c}}{\\rm{r}}{\\rm{i}}{\\rm{t}}}\\equiv 0.02\\,{M}_{\\odot }\\,{{\\rm{y}}{\\rm{r}}}^{-1}\\), we assume that the stellar radius expands owing to the injection of the large amount of entropy into the stellar envelope28<\/a>,78<\/a>,79<\/a>,80<\/a>. During this phase, we assume the surface temperature of the star to be 6,000\u2009K and the luminosity to be the Eddington value following the results of detailed stellar evolution calculations81<\/a>. After the accretion rate becomes smaller than \\({\\dot{M}}_{{\\rm{crit}}}\\), we gradually shrink the stellar radius to the radius expected for main-sequence stars on the timescale of the surface Kelvin\u2013Helmholtz time, which is given by ten times the stellar Kelvin\u2013Helmholtz time82<\/a>. Once the stellar age exceeds the lifetime predicted by these models, the star is converted into a BH particle if its final mass exceeds 260\u2009M\u2299 (ref.\u200983<\/a>). We note here that all of the stars formed during the simulation (the progenitor stars of LBH, MBH1 and MBH2) have masses greater than 260\u2009M\u2299 and collapse into the BHs without any supernova feedback. The mass accretion rate is measured by the sink method, in which the gas inside the sink radius is assimilated and added to the mass of the star particle. The sink radius is set to be 100 times the cell size, which will be converted into the sink particle, which ranges from 1 to 20\u2009pc. Because the sink radius is much larger than the physical stellar radius, we estimate the final stellar mass from the sink accretion rate using the empirical relation derived from previous high-resolution simulations of Pop\u2009III star formation60<\/a>,84<\/a>,<\/p>\n

$${M}_{\\ast }=250\\,{M}_{\\odot }{\\left(\\frac{{\\dot{M}}_{\\ast }}{2.8\\times 1{0}^{-3}{M}_{\\odot }{{\\rm{y}}{\\rm{r}}}^{-1}}\\right)}^{0.7},$$<\/p>\n

\n (3)\n <\/p>\n

in which \\({\\dot{M}}_{\\ast }\\) denotes the time-averaged accretion rate onto the growing protostar. When the stellar age reaches its lifetime, any excess gas mass accreted by the sink that exceeds the estimated stellar mass is returned to the surrounding gas cells to conserve mass. Throughout the simulation, we assign the sink-particle mass as the instantaneous stellar mass when evaluating radiative feedback. This approach overestimates the stellar luminosity and gives stronger feedback effects compared with the actual stellar mass, so the final stellar mass represents a lower limit to the actual value.<\/p>\n

We assume a stellar lifetime of 2\u2009Myr, typical for very massive stars with masses greater than 100\u2009M\u2299 (ref.\u200928<\/a>). After the stellar age exceeds 2\u2009Myr, we convert the stars into BHs. However, the lifetime of such stars\u2014particularly supermassive stars\u2014remains uncertain. Stars with masses greater than several 105\u2009M\u2299 become unstable owing to general-relativistic effects and collapse into BHs during the hydrogen-burning phase28<\/a>,85<\/a>,86<\/a>,87<\/a>,88<\/a>,89<\/a>. This implies that collapse may occur earlier than the canonical lifetime assumed for massive stars. Recent calculations90<\/a> show that the threshold mass for general relativistic instability depends on the mass accretion rate, spanning 2\u20138\u2009\u00d7\u2009105\u2009M\u2299. If the stars collapse into BHs earlier, the resulting BHs would remain embedded in dense gas for a longer period, thereby extending the phase of super-Eddington growth.<\/p>\n

After the stars collapse into BHs, we convert the corresponding sink particles into BH particles (LBH, MBH1 and MBH2). We keep their original accretion radii if the accretion radius is larger than the Bondi radius of the BHs for the ionized gas, otherwise the accretion radii are set to be the Bondi radii. Because the typical accretion radius is around a parsec, the simulations resolve the Bondi radius for most of the massive-seed BHs, allowing us to directly follow the gravitational capture of gas onto the BHs. The spectral energy distribution of accreting BHs is modelled as a double power-law continuum, F\u03bd\u2009\u221d\u2009\u03bd\u22120.6, extending from 1\u2009eV to 10\u2009eV, and F\u03bd\u2009\u221d\u2009\u03bd\u22121.5, extending from 10\u2009eV to 1\u2009keV (ref.\u200991<\/a>). This represents the integrated emission from the accretion disk and its associated corona. We allow the accretion rate to exceed the classical Eddington limit, consistent with theoretical models of super-Eddington accretion flows (for example, refs.\u200923<\/a>,25<\/a>,92<\/a>,93<\/a>). We assume the slim disk solution to give a bolometric luminosity when the accretion rate is higher than the Eddington rate, <\/p>\n

$${L}_{\\mathrm{bol}}=\\{\\begin{array}{cc}{{\\epsilon }}_{{\\rm{r}}}{\\dot{M}}_{\\mathrm{acc}}{c}^{2} & ({\\dot{M}}_{\\mathrm{acc}}\\le 2{\\dot{M}}_{\\mathrm{Edd}}),\\\\ 2\\left[1+\\log \\left(\\frac{{\\dot{M}}_{\\mathrm{acc}}}{2{\\dot{M}}_{\\mathrm{Edd}}}\\right)\\right]{L}_{\\mathrm{Edd}} & ({\\dot{M}}_{\\mathrm{acc}} > 2{\\dot{M}}_{\\mathrm{Edd}}),\\end{array}$$<\/p>\n

\n (4)\n <\/p>\n

in which \u03f5r\u2009=\u20090.1 is the radiative efficiency, \\({\\dot{M}}_{{\\rm{acc}}}\\) is the instantaneous gas accretion rate onto the BH, \\({\\dot{M}}_{{\\rm{Edd}}}\\) is the Eddington accretion rate and \\({L}_{{\\rm{Edd}}}\\equiv {\\dot{M}}_{{\\rm{Edd}}}\/{{\\epsilon }}_{{\\rm{r}}}{c}^{2}\\) is the Eddington luminosity94<\/a>. The emitted radiation is coupled to the surrounding gas through photoionization heating, which regulates the accretion flow and drives intermittent outflows around the BHs. We also allow mergers of BHs once the distance between two BHs becomes smaller than ten physical parsec.<\/p>\n

We do not include kinetic feedback from accreting BHs, such as jets or winds. Recent general relativistic magnetohydrodynamic simulations have shown that BHs accreting at super-Eddington rates can launch magnetically arrested jets and disk winds. If such jets break out and transport mass and momentum to larger scales in the surrounding medium, the accretion efficiency may decrease to less than the Eddington value95<\/a>,96<\/a>. The efficiency of jet launching should depend on the BH spin and magnetic field strength, both of which are highly uncertain. This uncertainty arises because modelling them would require following the stellar rotation at the birth of DCBHs and the detailed magnetic dynamo processes during accretion. The jet opening angle and its interaction with the accreting gas are also uncertain. We therefore consider our simulation as a fiducial model that provides an upper limit on BH growth under the assumption that kinetic feedback is absent.<\/p>\n

Zoom-in calculation of seed BH formation in the isolated cloud<\/p>\n

To assess the feasibility of heavy-seed BH formation, we perform a higher-resolution follow-up simulation that resolves the individual protostars expected to be the progenitors of the seed BHs. We focus on the formation site of the first heavy-seed BH, MBH1, identified in the cosmological radiation-hydrodynamic run. This seed, with a final mass of 6\u2009\u00d7\u2009105\u2009M\u2299, forms at redshift z\u2009\u2243\u200914. We extract a cubic region of one comoving kpc on a side centred on the host halo, shortly before the onset of collapse when the central gas density reaches 103\u2009cm\u22123. At this point, the central region of the cloud, of size a few tens of pc, becomes Jeans-unstable (Extended Data Fig.\u20096<\/a>). The re-simulation is performed using the same AREPO framework but with enhanced spatial and mass resolution to capture the internal fragmentation of the collapsing gas cloud. Sink particles are introduced once the local gas density exceeds 2\u2009\u00d7\u2009108\u2009cm\u22123, representing the formation of individual protostars. The sink radius is set to be four times the cell size, which will be converted into the sink particle, which ranges from 3,000 to 8,000\u2009au. We allow mergers of the sink particles once the distance of two sink particles becomes smaller than the sum of the sink radii.<\/p>\n

The numerical resolution used here is insufficient to resolve the opacity limit at n\u2009\u2248\u20091016\u2009cm\u22123, above which the gas becomes optically thick and protostars form97<\/a>. Indeed, the physical sizes of growing protostars are expected to be about 0.05\u2013100\u2009au, depending on the stellar mass and accretion rate81<\/a>, which is below our numerical resolution. We can, nevertheless, compare our results with higher-resolution studies. Becerra et al.98<\/a> studied the protostellar evolution changing nadib, the density above which the gas becomes adiabatic, from 108\u2009cm\u22123 to 1012\u2009cm\u22123. They have found that the mass evolution of supermassive protostars is not greatly affected by the numerical resolution, because the growth rate is mainly regulated by the large-scale accretion flow. They followed the early stellar evolution for 104\u2013105\u2009years and found that the accretion rate remains at about 1\u2009M\u2299\u2009yr\u22121, comparable to our results. Similarly, Chon and Omukai99<\/a> followed protostellar evolution in a DCBH halo for 2\u2009Myr with nadib\u2009=\u20091011\u2009cm\u22123. They found that the initial mass accretion rate remains around 1\u2009M\u2299\u2009yr\u22121. In their simulation, however, the accretion rate drops below 0.01\u2009M\u2299\u2009yr\u22121 at 104\u2013105\u2009years after protostar formation because of the limited gas supply. This contrasts our results, in which efficient accretion continues until around 2\u2009Myr after protostar formation.<\/p>\n

Zoom-in calculation of the BH accretion phase<\/p>\n

To reproduce the density structure around the massive-seed BH, we performed a high-resolution simulation that resolves gas densities up to 1012\u2009cm\u22123 and spatial scales down to 500\u2009au around the central BH. We take the snapshot at the moment when MBH2 forms at z\u2009=\u200913.1 and follow the evolution for 30\u2009kyr after its emergence. To capture the later evolution up to 0.5\u2009Myr after MBH2 formation, we also run a complementary simulation with a larger sink radius of 5,000\u2009au. Starting from the snapshot at 0.48\u2009Myr, we then reduce the sink radius back to 500\u2009au and continue the calculation to follow the detailed structure of the circum-BH disk for a further 0.55\u2009kyr. This procedure allows us to track both the long-term disk evolution and the small-scale accretion flow onto the BH.<\/p>\n

Host halo evolution of massive BHs<\/p>\n

Extended Data Fig.\u20093<\/a> shows the redshift evolution of the distance between the source galaxy and the host halos of MBH1, MBH2 and LBH. All BHs in our simulation form at distances of about ten physical kpc from the source galaxy, outside the virial radius of the source halo (dashed line). After their formation, the BHs migrate towards the source galaxy and enter the virial radius of the source halo at z\u2009\u2248\u200912. Initially, the BHs follow eccentric orbits, as indicated by the oscillatory behaviour of their distances. Over time, however, their orbits gradually decay and the BHs settle towards the galaxy centre by around z\u2009\u2248\u20098. We also find that, around this epoch, their accretion rates increase to about 0.1\u20131 times the Eddington value (Fig.\u20091f<\/a>).<\/p>\n

The host halos of the BHs may avoid metal enrichment because they are well separated, by about ten physical kpc, from the source galaxy that drives winds and pollutes the surrounding intergalactic medium with metals. Analytic estimates of galactic-wind expansion in ref.\u2009100<\/a> show that such winds expand only to around 1\u20132 physical kpc within 300\u2009Myr. Similarly, Ventura et al.101<\/a> statistically studied the distribution of superbubbles around star-forming halos in a 10\u2009h\u22121\u2009Mpc cosmological region and found that the bubble size is at most about 5\u2009kpc. These scales are much smaller than the separation found in our simulation, suggesting that the BH host halos can remain metal-poor despite their proximity to the source galaxy.<\/p>\n

We have shown that the key to forming an overmassive BH population is the formation of DCBHs in massive halos, with virial temperatures reaching about 40,000\u2009K. We attribute this to the delayed onset of protostar formation after Ly\u03b1 cooling becomes effective, caused by the finite collapse timescale of the halo and by dynamical heating. Extended Data Fig.\u20094<\/a> shows the time evolution of the halo growth rate and the peak gas density in the host halo of MBH1. The gas density increases steadily but protostar formation (at 300\u2009Myr) occurs only around 90\u2009Myr after the halo virial temperature first reaches 8,000\u2009K (at 210\u2009Myr), which is often assumed to mark the onset of DCBH formation. This delay arises for two reasons. First, the collapse timescale is on the order of the free-fall time, <\/p>\n

$${t}_{{\\rm{f}}{\\rm{f}}}\\approx \\sqrt{\\frac{1}{G\\rho }}=27.4\\,{\\rm{M}}{\\rm{y}}{\\rm{r}}{\\left(\\frac{n}{10{{\\rm{c}}{\\rm{m}}}^{-3}}\\right)}^{-1\/2},$$<\/p>\n

\n (5)\n <\/p>\n

in which n is the gas density of the gravitationally unstable region. This implies that DCBH formation should occur only after a timescale on the order of approximately 10\u2009Myr once the halo virial temperature exceeds 8,000\u2009K. The halo mass can grow during the initial free-fall phase to exceed Tvir\u2009=\u20098,000\u2009K. Second, dynamical heating caused by halo mergers further delays the collapse. During the first roughly 30\u2009Myr after the virial temperature reaches 8,000\u2009K, the peak gas density increases only weakly. Around this epoch, the host halo merges with nearby halos or clumps, which deposits further gravitational and kinetic energy into the collapsing gas and suppresses runaway collapse102<\/a>,103<\/a>. After these clump mergers subside and the halo growth rate decreases, the peak gas density begins to rise rapidly, eventually leading to the formation of the protostar that later collapses into MBH1. During this delay, the host halo grows sufficiently massive to enable subsequent super-Eddington accretion.<\/p>\n

Formation of extremely massive-seed BHs<\/p>\n

Here we show how extremely massive-seed BHs form\u2014through the formation of heavy seeds followed by a brief phase of super-Eddington accretion. Extended Data Fig.\u20091<\/a> shows the time evolution of the mass of the halo that later hosts the heavy-seed BH (MBH1). At the time of protostar formation, the halo mass is about 2\u2009\u00d7\u2009108\u2009M\u2299, corresponding to a virial temperature of 4\u2009\u00d7\u2009104\u2009K. This is much higher than the canonical virial temperature of 8,000\u2009K at which halo collapse is typically expected. Indeed, the semi-analytic model predicts that this halo would have collapsed at z\u2009\u2248\u200918.1, when its mass was about 107\u2009M\u2299, nearly an order of magnitude smaller than the host halo of MBH1. During the delayed collapse, more gas accumulates within the halo, leading to the formation of more massive stars and enhanced accretion onto the resulting BHs.<\/p>\n

Once the gas within a halo becomes gravitationally unstable, it collapses to form stars. The final stellar and BH masses are determined by how much gas is accreted onto the central protostar during its growth phase. Extended Data Fig.\u20095<\/a> shows the radial profiles of gas density (Extended Data Fig.\u20095a<\/a>), temperature (Extended Data Fig.\u20095b<\/a>) and escape velocity (Extended Data Fig.\u20095c<\/a>) at the time of protostar formation. The grey line represents the light-seed case (LBH) and the purple and green lines correspond to the two massive seeds, MBH1 and MBH2, respectively. The density structure differs markedly between light and heavy seeds. In the massive-seed cases, the high-density core extends to radii an order of magnitude larger than in the light-seed case, a consequence of the higher gas temperature in the collapsing cloud. Extended Data Fig.\u20095b<\/a> shows that the temperature in heavy-seed formation remains near 104\u2009K\u2014an order of magnitude higher than in the light-seed case\u2014owing to strong ultraviolet irradiation that suppresses molecular cooling. The collapse of this larger, hotter core also raises the escape velocity, which exceeds 10\u2009km\u2009s\u22121. This high binding energy facilitates continued mass accretion, as the ionized gas remains gravitationally bound and can feed the central protostar efficiently.<\/p>\n

To visualize how gravitational collapse proceeds, we compare the enclosed mass with the critical Bonnor\u2013Ebert (BE) mass, the threshold above which a cloud becomes gravitationally unstable to collapse104<\/a>,105<\/a>. The top panels in Extended Data Fig.\u20096<\/a> show the BE mass (green) and the enclosed mass (purple) for LBH, MBH1 and MBH2, from left to right. In the light-seed case, only the innermost 0.1\u20131\u2009pc region is gravitationally unstable, whereas in the heavy-seed cases, the unstable region extends to 10\u2013100\u2009pc. The bottom panels show the ratio of enclosed to BE mass as a function of enclosed mass; regions in which this ratio exceeds unity are gravitationally unstable. Only gas within about 103\u2009M\u2299 becomes unstable in the light-seed case, whereas up to about 106\u2009M\u2299, it is unstable in the heavy-seed cases. These characteristic mass scales correspond roughly to the resulting BH masses: about 800\u2009M\u2299 for LBH and 3\u2009\u00d7\u2009105\u2009M\u2299 and 6\u2009\u00d7\u2009105\u2009M\u2299 for MBH1 and MBH2, respectively. The dynamical time at the edge of the gravitationally unstable region is <\/p>\n

$${t}_{{\\rm{dyn}}}\\approx \\sqrt{\\frac{{R}^{3}}{G{M}_{{\\rm{enc}}}}}=1.5\\times 1{0}^{7}\\,{\\rm{years}}{\\left(\\frac{R}{100{\\rm{pc}}}\\right)}^{3\/2}{\\left(\\frac{{M}_{{\\rm{enc}}}}{1{0}^{6}{M}_{\\odot }}\\right)}^{-1\/2},$$<\/p>\n

\n (6)\n <\/p>\n

which is longer than the lifetime of massive stars (about 2\u2009Myr). This implies that the dense envelope remains bound and should continue to fall onto the BHs after the central protostars collapse, unless feedback from the protostars clears it away.<\/p>\n

The temperature decrease at r\u2009\u2272\u20091\u2009pc in MBH2 is caused by the enhanced electron fraction and the subsequent increase in molecular hydrogen formation, similar to the Pop\u2009III star formation in a fossil H\u2009II region106<\/a>. Extended Data Fig.\u20097<\/a> shows two-dimensional histograms of temperature (left), molecular hydrogen fraction (middle) and electron fraction (right) as functions of gas density, for MBH1 (top) and MBH2 (bottom). Gas at temperature T\u2009\u2248\u2009100\u20131,000\u2009K in MBH2 shows an enhanced H2 fraction, indicating that molecular formation and the associated cooling are responsible for the temperature drop. The electron fraction in the outer envelope is almost fully ionized but decreases towards higher densities as recombination proceeds. The increased electron fraction during collapse promotes H2 formation through the H\u2212 channel, the dominant pathway in primordial gas, which is catalysed by free electrons. The ionization source is the strong radiation from MBH1, which fully ionizes and heats the surrounding gas. Indeed, the temperature increase at r\u2009\u2273\u2009103\u2009pc in MBH2 indicates that intense radiation from MBH1 photoionizes the intergalactic medium. Despite the enhanced H2 formation, a massive-seed BH still forms in MBH2, as large-scale gravitational instability drives a strong gas inflow that overwhelms the cooling. This behaviour is analogous to the heavy-seed formation in the super-competitive accretion scenario, in which fine-structure lines and dust cooling reduce the gas temperature, whereas large-scale collapse collects a substantial amount of gas to form massive-seed BHs99<\/a>,107<\/a>.<\/p>\n

Observability of growing massive BHs<\/p>\n

To estimate the optical depth and H\u03b1 luminosity, we restrict the analysis to gas cells located within 30\u00b0 of the disk plane. We estimate the electron column density, Ne, by summing the total number of electrons contained within each spherical shell and dividing by the surface area of the shell, yielding an average column density at radius r. The Thomson optical depth is then calculated as \u03c4\u2009=\u2009Ne\u03c3T, in which \u03c3T\u2009=\u20096.65\u2009\u00d7\u200910\u221225\u2009cm2 is the Thomson scattering cross-section. In estimating the H\u03b1 luminosity, we have assumed case B recombination. This may underestimate H\u03b1 luminosity because resonance scattering will increase the H\u03b1 emissivity at densities of n\u2009\u2248\u20091010\u20131011\u2009cm\u22123 (ref.\u2009108<\/a>). The detailed radiative-transfer calculations for similar conditions indicate that the emergent H\u03b1\/H\u03b2 is enhanced to values of about 6\u201310, consistent with the observed values for LRDs32<\/a>.<\/p>\n

Present Chandra constraints indicate that LRDs are generally X-ray weak. Only three LRDs have secure individual X-ray detections so far47<\/a>,109<\/a>, whereas the remaining LRD population is largely undetected on an individual basis. Stacking analyses have also yielded mostly tentative signals or non-detections, including stacks of X-ray-undetected sources110<\/a>,111<\/a>,112<\/a>. To assess the detectability of X-rays emitted by the central BH in our simulation, we evaluated the hydrogen column density, NH, at two different snapshots, t\u2009=\u200926 and 523\u2009kyr. Extended Data Fig.\u20098<\/a> shows the distribution of NH measured along 64 different lines of sight evenly spaced around the accreting MBH2. During the early super-Eddington accretion phase at t\u2009=\u200926\u2009kyr, the column density exceeds 1026\u2009cm\u22122 in all directions. This implies that X-rays emitted from the BH are strongly Compton scattered and therefore undetectable from any viewing angle, consistent with the non-detections obtained from stacking Chandra data112<\/a>. At the later epoch, nearly Eddington phase at t\u2009=\u2009523\u2009kyr, the BH remains embedded in dense gas but the column density decreases to NH\u2009=\u20091025\u20131026\u2009cm\u22122, lower than in the early super-Eddington phase. Even in this phase, however, the X-ray emission is still expected to be difficult to observe. This analysis indicates that the LRD-like system found in our study remains X-ray dark for at least the first approximately 0.5\u2009Myr, in agreement with present observational constraints.<\/p>\n

Inefficient growth of Pop III remnant BH<\/p>\n

We find that LBH undergoes almost negligible growth during the simulation. Its mass remains close to 800\u2009M\u2299 until z\u2009\u2248\u200912, followed only by modest growth thereafter. We mainly attribute this inefficient growth to the shallow gravitational potential of the host halo at the time of Pop\u2009III star formation. Extended Data Fig.\u20095c<\/a> shows that the escape velocity is only about 2\u20135\u2009km\u2009s\u22121, smaller than the sound speed of ionized gas, which is about 10\u2009km\u2009s\u22121. This indicates that the halo potential is too shallow to retain ionized gas, so the gas inside the halo is rapidly expelled by photoevaporation40<\/a>,113<\/a>,114<\/a>,115<\/a>. Indeed, the gas density around LBH remains \u22721\u2009cm\u22123, with a temperature of approximately 104\u2009K (Extended Data Fig.\u20095a,b<\/a>). Under these conditions, the Bondi accretion rate is only about 10\u22128\u201310\u22127\u2009M\u2299\u2009yr\u22121, much lower than the Eddington accretion rate of LBH.<\/p>\n

Around z\u2009\u2248\u200912, the halo potential becomes sufficiently deep for gas to begin accumulating at the halo centre. However, the inflowing gas is predominantly supplied to MBH1 rather than to LBH. MBH1 formed in a halo located close to the original host halo of LBH and the two host halos subsequently merged. After the merger, MBH1 becomes the primary accretor because of its larger mass, whereas LBH remains less efficiently fed. This is partly because LBH has a smaller mass and therefore a longer dynamical-friction timescale than the massive BHs, further suppressing its orbital decay and growth26<\/a>,40<\/a>. Although stochastic interactions with nearby gas increase the mass of LBH to several thousand solar masses, its accretion rate remains much lower than those of the massive BHs.<\/p>\n

DCBH formation in lower-z Universe<\/p>\n

To test how ubiquitous DCBH formation and LRD-like systems are, and to examine whether this phenomenon can occur in the later Universe, we performed radiation hydrodynamics calculations for another candidate region selected from ref.\u200917<\/a>. To identify a lower-redshift sample, we chose a candidate halo that satisfies the DCBH formation criteria at z\u2009\u2248\u200914, the latest formation redshift among our 62 samples.<\/p>\n

We find eight seed BHs forming in the region around the target halo. Extended Data Fig.\u20099a<\/a> shows the spatial distribution of the massive BHs formed in this region at z\u2009=\u20097.011. At this snapshot, the lowest-redshift DCBH, labelled MBH8, has already formed. Some of the BHs have migrated to the centres of massive halos, whereas other clumps around the target halo are still forming BHs. Extended Data Fig.\u20099b,c<\/a> shows the redshift evolution of the BH mass and the accretion rate normalized by the Eddington rate, respectively. All seed BHs grow beyond 106\u2009M\u2299 after experiencing brief episodes of super-Eddington accretion. Some BHs, such as MBH1 and MBH3, repeatedly increase their accretion rates and undergo several super-Eddington phases. These results demonstrate that DCBH formation followed by later super-Eddington accretion is not confined to the earliest cosmic epochs but can also occur at later times, overlapping with the redshift range in which LRDs have been detected.<\/p>\n

Robustness tests<\/p>\n

To assess numerical robustness and model sensitivity, we performed one higher-resolution run (lv13) and three simulations with reduced FUV background intensities (w2, w5 and w10). In the last three runs, we kept the parameters of the semi-analytic model fixed but reduced the LW intensity from its fiducial value obtained from the model. In the higher-resolution run, we use an effective resolution of 8,1923 in the zoom-in region, corresponding to dark-matter and baryonic particle masses of 542 and 99.2\u2009h\u22121\u2009M\u2299, respectively. This resolution is sufficient to resolve mini-halos in which the first stars form102<\/a>,116<\/a>. The run yields results that are nearly identical to the fiducial case: one LBH forms at z\u2009\u2248\u200922, followed by the formation of two massive seeds. The blue line in Extended Data Fig.\u200910<\/a> shows the evolution of the most massive BH in the higher-resolution run, demonstrating that the mass growth is largely insensitive to the numerical resolution achieved in our default simulation.<\/p>\n

To examine uncertainties in the semi-analytic predictions, we also ran simulations with lower FUV background intensities. In our semi-analytic model, the stellar component of the source galaxy is populated using standard galaxy-formation prescriptions, but the star-formation rate and efficiency are subject to substantial uncertainty. A lower star-formation efficiency would result in a smaller luminosity and therefore a reduced FUV intensity incident on the massive-BH-forming halo. The reduced-FUV runs also capture possible spatial variations in the radiation field, which are not explicitly modelled in our semi-analytic treatment. If heavy seeds still form under weaker FUV irradiation, this supports both the robustness of the semi-analytic predictions and the insensitivity to spatial variations.<\/p>\n

The green, yellow and red lines in Extended Data Fig.\u200910<\/a> show the BH mass evolution when the background intensity is reduced by factors of two, five and ten, respectively. The dependence on the seed-formation epoch and subsequent growth is weak. As the FUV intensity decreases, the most massive BH forms slightly earlier, but the difference in formation redshift between the fiducial run and the lowest-intensity case is less than \u0394z\u2009\u2243\u20090.1. This behaviour is consistent with the expectation that weaker FUV irradiation allows more efficient H2 cooling, leading to earlier collapse of the target halo. The later mass growth is also only mildly affected: lower FUV intensity results in a smaller final BH mass but the variation remains within a factor of three. Even in the weakest-FUV case, the final BH reaches 1.9\u2009\u00d7\u2009107\u2009M\u2299 by z\u2009=\u20097, highlighting the robustness of overmassive BH formation despite model uncertainties.<\/p>\n","protected":false},"excerpt":{"rendered":"Our cosmological initial conditions are based on the simulation Phi-4096 of ref.\u200951. In this work, a dark-matter-only N-body…\n","protected":false},"author":2,"featured_media":850028,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[49],"tags":[93312,44060,93313,1159,1160,199,79],"class_list":["post-850027","post","type-post","status-publish","format-standard","has-post-thumbnail","category-physics","tag-computational-astrophysics","tag-early-universe","tag-galaxies-and-clusters","tag-humanities-and-social-sciences","tag-multidisciplinary","tag-physics","tag-science"],"_links":{"self":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts\/850027","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/comments?post=850027"}],"version-history":[{"count":0,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/posts\/850027\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/media\/850028"}],"wp:attachment":[{"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/media?parent=850027"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/categories?post=850027"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.newsbeep.com\/us\/wp-json\/wp\/v2\/tags?post=850027"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}