No Arabic abstract
Primordial supermassive stars (SMSs) formed in atomic-cooling halos at z ~ 15 - 20 are leading candidates for the seeds of the first quasars. Past numerical studies of the evolution of SMSs have typically assumed constant accretion rates rather than the highly variable flows in which they form. We model the evolution of SMSs in the cosmological flows that create them using the Kepler stellar evolution and implicit hydrodynamics code. We find that they reach masses of 1 - 2 x $10^5 M_{odot}$ before undergoing direct-collapse to black holes (DCBHs) during or at the end of their main-sequence hydrogen burning, at 1 - 1.5 Myr, regardless of halo mass, spin, or merger history. We also find that realistic, highly-variable accretion histories allow for a much greater diversity of supermassive stellar structures, including in some cases largely thermally relaxed objects, which may provide a significant source of radiative feedback. Our models indicate that the accretion histories predicted for purely atomic-cooling halos may impose a narrow spectrum of masses on the seeds of the first massive quasars, however further studies incorporating realistic feedback will be essential in order to confirm whether or not this holds true in all cases. Our results also indicate that multiple SMSs at disparate stages of evolution can form in these halos, raising the possibility of SMS binaries and supermassive X-ray binaries (SMXBs), as well as DCBH mergers which could be detected by LISA.
The collapse of supermassive primordial stars in hot, atomically-cooled halos may have given birth to the first quasars at $z sim$ 15 - 20. Recent numerical simulations of these rapidly accreting stars reveal that they are cool, red hypergiants shrouded by dense envelopes of pristine atomically-cooled gas at 6,000 - 8,000 K, with luminosities $L$ $gtrsim$ 10$^{10}$ L$_{odot}$. Could such luminous but cool objects be detected as the first stage of quasar formation in future near infrared (NIR) surveys? We have now calculated the spectra of supermassive primordial stars in their birth envelopes with the Cloudy code. We find that some of these stars will be visible to JWST at $z lesssim$ 20 and that with modest gravitational lensing Euclid and WFIRST could detect them out to $z sim$ 10 - 12. Rather than obscuring the star, its accretion envelope enhances its visibility in the NIR today by reprocessing its short-wavelength flux into photons that are just redward of the Lyman limit in the rest frame of the star.
Supermassive primordial stars are suspected to be the progenitors of the most massive quasars at z~6. Previous studies of such stars were either unable to resolve hydrodynamical timescales or considered stars in isolation, not in the extreme accretion flows in which they actually form. Therefore, they could not self-consistently predict their final masses at collapse, or those of the resulting supermassive black hole seeds, but rather invoked comparison to simple polytropic models. Here, we systematically examine the birth, evolution and collapse of accreting non-rotating supermassive stars under accretion rates of 0.01-10 solar masses per year, using the stellar evolution code KEPLER. Our approach includes post-Newtonian corrections to the stellar structure and an adaptive nuclear network, and can transition to following the hydrodynamic evolution of supermassive stars after they encounter the general relativistic instability. We find that this instability triggers the collapse of the star at masses of 150,000-330,000 solar masses for accretion rates of 0.1-10 solar masses per year, and that the final mass of the star scales roughly logarithmically with the rate. The structure of the star, and thus its stability against collapse, is sensitive to the treatment of convection, and the heat content of the outer accreted envelope. Comparison with other codes suggests differences here may lead to small deviations in the evolutionary state of the star as a function of time, that worsen with accretion rate. Since the general relativistic instability leads to the immediate death of these stars, our models place an upper limit on the masses of the first quasars at birth.
Supermassive black holes (SMBHs) are ubiquitous in galaxies with a sizable mass. It is expected that a pair of SMBHs originally in the nuclei of two merging galaxies would form a binary and eventually coalesce via a burst of gravitational waves. So far theoretical models and simulations have been unable to predict directly the SMBH merger timescale from ab-initio galaxy formation theory, focusing only on limited phases of the orbital decay of SMBHs under idealized conditions of the galaxy hosts. The predicted SMBH merger timescales are long, of order Gyrs, which could be problematic for future gravitational wave searches. Here we present the first multi-scale $Lambda$CDM cosmological simulation that follows the orbital decay of a pair of SMBHs in a merger of two typical massive galaxies at $zsim3$, all the way to the final coalescence driven by gravitational wave emission. The two SMBHs, with masses $sim10^{8}$ M$_{odot}$, settle quickly in the nucleus of the merger remnant. The remnant is triaxial and extremely dense due to the dissipative nature of the merger and the intrinsic compactness of galaxies at high redshift. Such properties naturally allow a very efficient hardening of the SMBH binary. The SMBH merger occurs in only $sim10$ Myr after the galactic cores have merged, which is two orders of magnitude smaller than the Hubble time.
The next generation of electromagnetic and gravitational wave observatories will open unprecedented windows to the birth of the first supermassive black holes. This has the potential to reveal their origin and growth in the first billion years, as well as the signatures of their formation history in the local Universe. With this in mind, we outline three key focus areas which will shape research in the next decade and beyond: (1) What were the seeds of the first quasars; how did some reach a billion solar masses before z$sim7$? (2) How does black hole growth change over cosmic time, and how did the early growth of black holes shape their host galaxies? What can we learn from intermediate mass black holes (IMBHs) and dwarf galaxies today? (3) Can we unravel the physics of black hole accretion, understanding both inflows and outflows (jets and winds) in the context of the theory of general relativity? Is it valid to use these insights to scale between stellar and supermassive BHs, i.e., is black hole accretion really scale invariant? In the following, we identify opportunities for the Canadian astronomical community to play a leading role in addressing these issues, in particular by leveraging our strong involvement in the Event Horizon Telescope, the {it James Webb Space Telescope} (JWST), Euclid, the Maunakea Spectroscopic Explorer (MSE), the Thirty Meter Telescope (TMT), the Square Kilometer Array (SKA), the Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR), and more. We also discuss synergies with future space-based gravitational wave (LISA) and X-ray (e.g., Athena, Lynx) observatories, as well as the necessity for collaboration with the stellar and galactic evolution communities to build a complete picture of the birth of supermassive black holes, and their growth and their influence over the history of the Universe.
The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have less concentrated structure than fully-convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit $gtrsim 10^5~M_odot$ derived for the fully-convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with $dot M lesssim 0.1~M_odot~{rm yr}^{-1}$. With $dot{M} simeq 0.3 - 1~M_odot~{rm yr}^{-1}$, the star becomes GR-unstable during the helium-burning stage at $M simeq 2 - 3.5~times 10^5~M_odot$. In an extreme case with $10~M_odot~{rm yr}^{-1}$, the star does not collapse until the mass reaches $simeq 8.0times 10^5~M_odot$, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.