Showing 1–50 of 117 results for author: Omukai, K

Impact of Ly$α$ radiation force on super-Eddington accretion onto a massive black hole

Authors: Takuya Mushano, Takumi Ogawa, Ken Ohsuga, Hidenobu Yajima, Kazuyuki Omukai

Abstract: The viability of super-Eddington accretion remains a topic of intense debate, crucial for understanding the formation of supermassive black holes in the early universe. However, the impact of Ly$α$ radiation force on this issue remains poorly understood. We investigate the propagation of the Ly$α$ photons and evaluate the Ly$α$ radiation force within a spherically symmetric accreting HI gas onto t… ▽ More The viability of super-Eddington accretion remains a topic of intense debate, crucial for understanding the formation of supermassive black holes in the early universe. However, the impact of Ly$α$ radiation force on this issue remains poorly understood. We investigate the propagation of the Ly$α$ photons and evaluate the Ly$α$ radiation force within a spherically symmetric accreting HI gas onto the central black hole. We solve the radiation transfer equation, incorporating the destruction processes of Ly$α$ photons through two-photon decay and collisional de-excitation. We find that the Ly$α$ photons, originating in the HII region around black holes, suffer from multiple resonance scattering before being destroyed via two-photon decay and collisional de-excitation. Hence, the Ly$α$ radiation force undergoes a significant amplification, surpassing gravity at the innermost section of the HI region. This amplification, quantified as the force multiplier, reaches approximately 130 and remains nearly constant, regardless of the optical depth at the line center, provided the optical thickness of the flow is within the range of $10^{10-14}$. The requisite lower limit of the product of gas density and black hole mass to realize the super-Eddington accretion is found to be in the range $(2-{40}) \times 10^9 M_\odot\,{\rm cm}^{-3}$, which is a few to tens of times larger than the minimum value obtained without accounting for the Ly$α$ radiation force. The pronounced amplification of the Ly$α$ radiation force poses a substantial challenge to the feasibility of super-Eddington accretion. △ Less

Submitted 6 October, 2024; originally announced October 2024.

Comments: Accepted for publication in Publications of the Astronomical Society of Japan (PASJ), 10 pages, 5 figures, 1 table

Gravitational collapse at low to moderate Mach numbers: The relationship between star formation efficiency and the fraction of mass in the massive object

Authors: Jorge Saavedra-Bastidas, Dominik R. G. Schleicher, Ralf S. Klessen, Sunmyon Chon, Kazuyuki Omukai, Thomas Peters, Lewis R. Prole, Bastián Reinoso, Rafeel Riaz, Paulo Solar

Abstract: 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-… ▽ More 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 ε_*, 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 ε_*, although the steepness of the increase depends on the environment. The relation between f_* and ε_* 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 ε_* 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. △ Less

Submitted 2 August, 2024; originally announced August 2024.

Comments: 13 pages, 8 figures, 2 tables. Submitted to A&A

Journal ref: A&A 690, A186 (2024)

The Impact of Stellar Radiative Feedback on Formation of Young Massive Clusters via Fast HI Gas Collisions

Authors: Ryunosuke Maeda, Tsuyoshi Inoue, Kazuyuki Omukai, Yasuo Fukui, Kisetsu Tsuge

Abstract: Young massive clusters (YMCs) are dense aggregates of young stars and are often speculated as potential precursors to globular clusters. However, the formation mechanism of massive and compact gas clumps that precede YMCs remains unknown. In this paper, we study the formation of such massive clumps via fast HI gas collisions (~100 km/s) as suggested by recent observations and their subsequent evol… ▽ More Young massive clusters (YMCs) are dense aggregates of young stars and are often speculated as potential precursors to globular clusters. However, the formation mechanism of massive and compact gas clumps that precede YMCs remains unknown. In this paper, we study the formation of such massive clumps via fast HI gas collisions (~100 km/s) as suggested by recent observations and their subsequent evolution into YMCs by using three-dimensional magnetohydrodynamics simulations involving self-gravity and detailed thermal/chemical processes. In particular, the impact of ionization feedback from stellar radiation is included in an approximate fashion where the temperature within the HII regions is elevated to 10,000 K, while supernova feedback is not included. We examine whether the resulting massive clumps can survive this ionization feedback and evolve into YMCs. Our simulations reveal the emergence of gas clumps that do not only possess substantial mass (~10^5 M_sun) but also sufficient compactness (~5 pc). Notably, these clumps exhibit significantly higher escape velocities compared to the sound speed of the HII region, indicating effective gravitational retention of gas against feedback-induced evaporation. Consequently, these conditions foster efficient star formation within the massive gas clumps, ultimately leading to their evolution into YMCs. We also perform simulations involving lower-velocity gas collisions, approximately 15 km/s, typical shock velocities induced by galactic superbubbles. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 15 pages, 13 figures, and 1 table; Accepted for publication in ApJ

Impact of Turbulent Magnetic Fields on Disk Formation and Fragmentation in First Star Formation

Authors: Kenji Eric Sadanari, Kazuyuki Omukai, Kazuyuki Sugimura, Tomoaki Matsumoto, Kengo Tomida

Abstract: Recent cosmological hydrodynamic simulations have suggested that the first stars in the universe often form as binary or multiple systems. However, previous studies typically overlooked the potential influence of magnetic fields during this process, assuming them to be weak and minimally impactful. Emerging theoretical investigations, however, propose an alternative perspective, suggesting that tu… ▽ More Recent cosmological hydrodynamic simulations have suggested that the first stars in the universe often form as binary or multiple systems. However, previous studies typically overlooked the potential influence of magnetic fields during this process, assuming them to be weak and minimally impactful. Emerging theoretical investigations, however, propose an alternative perspective, suggesting that turbulent dynamo effects within first-star forming clouds can generate strong magnetic fields. In this study, we perform three-dimensional ideal magnetohydrodynamics simulations, starting from the gravitational collapse of a turbulent cloud core to the early accretion phase, where disk fragmentation frequently occurs. Our findings reveal that turbulent magnetic fields, if they reach an equipartition level with turbulence energy across all scales during the collapse phase, can significantly affect the properties of the multiple systems. Specifically, both magnetic pressure and torques contribute to disk stabilization, leading to a reduction in the number of fragments, particularly for low-mass stars. Additionally, our observations indicate the launching of protostellar jets driven by magnetic pressure of toroidal fields, although their overall impact on star formation dynamics appears to be minor. Given the case with which seed magnetic fields amplify to the full equipartition level, our results suggest that magnetic fields likely play a significant role in shaping the initial mass function of the first stars, highlighting the importance of magnetic effects on star formation in the early universe. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 21 pages, 13 figures, 2 tables, Accepted for publication in PASJ

First star formation in extremely early epochs

Authors: Mana Ito, Kazuyuki Omukai

Abstract: First stars play crucial roles in development of the universe, influencing events like cosmic reionization and the chemical enrichment. While first stars are conventionally thought to form at around $z \sim 20-30$ in the standard $Λ$ Cold Dark Matter ($Λ$CDM) cosmology, observational constraints on small-scale density fluctuations remain limited, possibly differing significantly from the scale-inv… ▽ More First stars play crucial roles in development of the universe, influencing events like cosmic reionization and the chemical enrichment. While first stars are conventionally thought to form at around $z \sim 20-30$ in the standard $Λ$ Cold Dark Matter ($Λ$CDM) cosmology, observational constraints on small-scale density fluctuations remain limited, possibly differing significantly from the scale-invariant fluctuations assumed in the $Λ$CDM model. Should this be the case, the formation of first stars could occur much earlier than typically predicted. In this study, we investigate the formation process of first stars in the extremely early epochs of $z \gtrsim 100$ in the post-recombination universe. At such early times, the effects of the warm cosmic microwave background (CMB) become significant. We calculate the collapse of primordial star-forming clouds using a one-zone thermo-chemical model that accounts for CMB influences on radiative heating, Compton cooling, and photodissociation reactions. We found that the impact of the CMB on the evolution is limited at $z \lesssim 100$, with the temperature evolution closely resembling the conventional model. However, within the range $100 \lesssim z \lesssim 400$, the formation of H$_2$ via the H$^-$ channel is impeded by H$^-$ photodetachment induced by the CMB, leading to higher temperatures compared to standard one. Consequently, first stars with masses exceeding $1000 ~\mathrm{M}_\odot$ can emerge at $z \gtrsim 100$. Furthermore, at $z \gtrsim 500$, the temperature evolution becomes nearly isothermal solely due to atomic cooling, as H$_2$ formation is entirely suppressed. In such cases, supermassive stars with masses around $\sim 10^5 ~\mathrm{M}_\odot$ are expected to form solely via atomic cooling. These findings emphasize the significant variation in the typical mass of the first stars depending on the epoch of formation. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Impact of radiative feedback on the initial mass function of metal-poor stars

Authors: Sunmyon Chon, Takashi Hosokawa, Kazuyuki Omukai, Raffaella Schneider

Abstract: 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… ▽ More 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. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 22 pages, 19 figures, submitted to MNRAS

Formation of Massive and Wide First-star Binaries in Radiation Hydrodynamics Simulations

Authors: Kazuyuki Sugimura, Tomoaki Matsumoto, Takashi Hosokawa, Shingo Hirano, Kazuyuki Omukai

Abstract: We study the formation of Pop III stars by performing radiation hydrodynamics simulations for three different initial clouds extracted from cosmological hydrodynamics simulations. Starting from the cloud collapse stage, we follow the growth of protostars by accretion for $\sim 10^5$ yr until the radiative feedback from the protostars suppresses the accretion and the stellar properties are nearly f… ▽ More We study the formation of Pop III stars by performing radiation hydrodynamics simulations for three different initial clouds extracted from cosmological hydrodynamics simulations. Starting from the cloud collapse stage, we follow the growth of protostars by accretion for $\sim 10^5$ yr until the radiative feedback from the protostars suppresses the accretion and the stellar properties are nearly fixed. We find that the Pop III stars form in massive and wide binaries/small-multiple stellar systems, with masses $>30\,M_\odot$ and separations $>2000$ au. We also find that the properties of the final stellar system correlate with those of the initial clouds: the total mass increases with the cloud-scale accretion rate, and the angular momentum of the binary orbit matches that of the initial cloud. While the total mass of the system in our simulations is consistent with our previous single-star formation simulations, individual masses are lower due to mass sharing, suggesting potential modification in the extent of feedback from Pop III stars in the subsequent evolution of the Universe. We also identify such systems as mini-binaries embedded in a wider outer multiple-star system, which could evolve into progenitors for observed gravitational wave events. △ Less

Submitted 23 January, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Comments: 30 pages, 20 figures, published in ApJ

Journal ref: ApJ, 959, 17 (2023)

Metallicity Dependence of Molecular Cloud Hierarchical Structure at Early Evolutionary Stages

Authors: Masato I. N. Kobayashi, Kazunari Iwasaki, Kengo Tomida, Tsuyoshi Inoue, Kazuyuki Omukai, Kazuki Tokuda

Abstract: The formation of molecular clouds out of HI gas is the first step toward star formation. Its metallicity dependence plays a key role to determine star formation through the cosmic history. Previous theoretical studies with detailed chemical networks calculate thermal equilibrium states and/or thermal evolution under one-zone collapsing background. The molecular cloud formation in reality, however,… ▽ More The formation of molecular clouds out of HI gas is the first step toward star formation. Its metallicity dependence plays a key role to determine star formation through the cosmic history. Previous theoretical studies with detailed chemical networks calculate thermal equilibrium states and/or thermal evolution under one-zone collapsing background. The molecular cloud formation in reality, however, involves supersonic flows, and thus resolving the cloud internal turbulence/density structure in three dimension is still essential. We here perform magnetohydrodynamics simulations of 20 km s^-1 converging flows of Warm Neutral Medium (WNM) with 1 micro Gauss mean magnetic field in the metallicity range from the Solar (1.0 Zsun) to 0.2 Zsun environment. The Cold Neutral Medium (CNM) clumps form faster with higher metallicity due to more efficient cooling. Meanwhile, their mass functions commonly follow dn/dm proportional to m^-1.7 at three cooling times regardless of the metallicity. Their total turbulence power also commonly shows the Kolmogorov spectrum with its 80 percent in the solenoidal mode, while the CNM volume alone indicates the transition towards the Larson's law. These similarities measured at the same time in the unit of the cooling time suggest that the molecular cloud formation directly from the WNM alone requires a longer physical time in a lower metallicity environment in the 1.0--0.2 Zsun range. To explain the rapid formation of molecular clouds and subsequent massive star formation possibly within less than 10 Myr as observed in the Large/Small Magellanic Clouds (LMC/SMC), the HI gas already contains CNM volume instead of pure WNM. △ Less

Submitted 23 July, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 23 pages, 11 figures. Accepted for publication in ApJ

Exploring the nature of UV-bright $z \gtrsim 10$ galaxies detected by JWST: star formation, black hole accretion, or a non-universal IMF?

Authors: Alessandro Trinca, Raffaella Schneider, Rosa Valiante, Luca Graziani, Arianna Ferrotti, Kazuyuki Omukai, Sunmyon Chon

Abstract: 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 gal… ▽ More 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. △ Less

Submitted 12 February, 2024; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 20 pages, 12 figures. Submitted to MNRAS

Stability of relativistic tidal response against small potential modification

Authors: Takuya Katagiri, Hiroyuki Nakano, Kazuyuki Omukai

Abstract: The tidal response of compact objects in an inspiraling binary system is measured by a set of tidal Love and dissipation numbers imprinted in the gravitational waveforms. While a four-dimensional black hole in vacuum within General Relativity has vanishing Love numbers, a black hole in alternative theories of gravity can acquire non-vanishing Love numbers. The dissipation numbers may quantify Plan… ▽ More The tidal response of compact objects in an inspiraling binary system is measured by a set of tidal Love and dissipation numbers imprinted in the gravitational waveforms. While a four-dimensional black hole in vacuum within General Relativity has vanishing Love numbers, a black hole in alternative theories of gravity can acquire non-vanishing Love numbers. The dissipation numbers may quantify Planckian corrections at the horizon scale. These properties will allow a test of classical theories of gravity in the strong-field regime with gravitational-wave observation. Since black holes are not in the exact vacuum environment in astrophysical situations, the following question arises: can the environment affect the tidal response? In this paper, we investigate the stability of the tidal response of a Schwarzschild black hole for frequency-dependent tidal-field perturbations against a small modification of the background. Our analysis relies on the scattering theory, which overcomes difficulties in computing the relativistic tidal Love numbers. The tidal Love and dissipation numbers can be understood from the property of sufficiently low-frequency scattering waves. We show that the tidal Love numbers are sensitive to the property of the modification. Therefore, we need careful consideration of the environment around the black hole in assessing the deviation of the underlying theory of gravity from General Relativity with the Love numbers. The modification has less impact on the dissipation numbers, indicating that quantifying the existence of the event horizon with them is not spoiled. We also demonstrate that in a composite system, i.e., a compact object with environmental effects, the Love and dissipation numbers are approximately determined by the sum of the numbers of each component. △ Less

Submitted 12 October, 2023; v1 submitted 10 April, 2023; originally announced April 2023.

Comments: 44 pages, 19 figures, v2: minor corrections, v3: minor corrections, accepted for publication in PRD

Direct-collapse black hole formation induced by internal radiation of host halos

Authors: Gen Chiaki, Sunmyon Chon, Kazuyuki Omukai, Alessandro Trinca, Raffaella Schneider, Rosa Valiante

Abstract: 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… ▽ More 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. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 17 pages, 12 figures, 2 tables, accepted by MNRAS

A needle in a haystack? Catching Pop III stars in the Epoch of Reionization: I. Pop III star forming environments

Authors: Alessandra Venditti, Luca Graziani, Raffaella Schneider, Laura Pentericci, Claudia Di Cesare, Umberto Maio, Kazuyuki Omukai

Abstract: Despite extensive search efforts, direct observations of the first (Pop III) stars have not yet succeeded. Theoretical studies have suggested that late Pop III star formation is still possible in pristine clouds of high-mass galaxies, coexisting with Pop II stars, down to the Epoch of Reionization (EoR). Here we reassess this finding by exploring Pop III star formation in six $50h^{-1} ~$cMpc simu… ▽ More Despite extensive search efforts, direct observations of the first (Pop III) stars have not yet succeeded. Theoretical studies have suggested that late Pop III star formation is still possible in pristine clouds of high-mass galaxies, coexisting with Pop II stars, down to the Epoch of Reionization (EoR). Here we reassess this finding by exploring Pop III star formation in six $50h^{-1} ~$cMpc simulations performed with the hydrodynamical code dustyGadget. We find that Pop III star formation ($\sim 10^{-3.4} - 10^{-3.2} ~ \mathrm{M_\odot yr^{-1} cMpc^{-3}}$) is still occurring down to $z \sim 6 - 8$, i.e. well within the reach of deep JWST surveys. At these epochs, $\gtrsim 10 \%$ of the rare massive galaxies with $M_\star \gtrsim 3 \times 10^9 ~ \mathrm{M_\odot}$ are found to host Pop III stars, although with a Pop III/Pop II mass fraction $\lesssim 0.1 \%$. Regardless of their mass, Pop III hosting galaxies are mainly found on the main sequence, at high star formation rates, probably induced by accretion of pristine gas. This scenario is also supported by their increasing star formation histories and their preferential location in high-density regions of the cosmic web. Pop III stars are found both in the outskirts of metal-enriched regions and in isolated, pristine clouds. In the latter case, their signal may be less contaminated by Pop IIs, although its detectability will strongly depend on the specific line-of-sight to the source, due to the complex morphology of the host galaxy and its highly inhomogeneous dust distribution. △ Less

Submitted 20 April, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: 24 pages, 13 figures, 4 tables. Accepted for publication in MNRAS. Main changes: improved statistical sample with analysis of the spatial distribution of Pop III regions; corrections to the reported values of 12+Log(O/H); comparison observational data updated

Non-ideal magnetohydrodynamic simulations of the first star formation: the effect of ambipolar diffusion

Authors: Kenji Eric Sadanari, Kazuyuki Omukai, Kazuyuki Sugimura, Tomoaki Matsumoto, Kengo Tomida

Abstract: In the present-day universe, magnetic fields play such essential roles in star formation as angular momentum transport and outflow driving, which control circumstellar disc formation/fragmentation and also the star formation efficiency. While only a much weaker field has been believed to exist in the early universe, recent theoretical studies find that strong fields can be generated by turbulent d… ▽ More In the present-day universe, magnetic fields play such essential roles in star formation as angular momentum transport and outflow driving, which control circumstellar disc formation/fragmentation and also the star formation efficiency. While only a much weaker field has been believed to exist in the early universe, recent theoretical studies find that strong fields can be generated by turbulent dynamo during the gravitational collapse. Here, we investigate the gravitational collapse of a cloud core ($\sim 10^{3}\ \rm cm^{-3}$) up to protostar formation ($\sim 10^{20}\ \rm cm^{-3}$) by non-ideal magnetohydrodynamics (MHD) simulations considering ambipolar diffusion (AD), the dominant non-ideal effects in the primordial-gas. We systematically study rotating cloud cores either with or without turbulence and permeated with uniform fields of different strengths. We find that AD can slightly suppress the field growth by dynamo especially on scales smaller than the Jeans-scale at the density range $10^{10}-10^{14}\ \rm cm^{-3}$, while we could not see the AD effect on the temperature evolution, since the AD heating rate is always smaller than compression heating. The inefficiency of AD makes the field as strong as $10^{3}-10^{5} \rm\ G$ near the formed protostar, much stronger than in the present-day cases, even in cases with initially weak fields. The magnetic field affects the inflow motion when amplified to the equipartition level with turbulence on the Jeans-scale, although disturbed fields do not launch winds. This might suggest that dynamo amplified fields have smaller impact on the dynamics in the later accretion phase than other processes such as ionisation feedback. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 15 pages, 11 figures, 2 tables. Accepted for publication in MNRAS

Vanishing Love of Black Holes in General Relativity: From Spacetime Conformal Symmetry of a Two-dimensional Reduced Geometry

Authors: Takuya Katagiri, Masashi Kimura, Hiroyuki Nakano, Kazuyuki Omukai

Abstract: We study the underlying structure of the vanishing of the Love numbers of both Schwarzschild and Kerr black holes in terms of spacetime conformal symmetry in a unified manner for the static spin-$s$ fields. The perturbations can be reduced with the harmonic decomposition to a set of infinite static scalar fields in a two-dimensional anti-de Sitter spacetime~$({\rm AdS}_2)$. In the reduced system,… ▽ More We study the underlying structure of the vanishing of the Love numbers of both Schwarzschild and Kerr black holes in terms of spacetime conformal symmetry in a unified manner for the static spin-$s$ fields. The perturbations can be reduced with the harmonic decomposition to a set of infinite static scalar fields in a two-dimensional anti-de Sitter spacetime~$({\rm AdS}_2)$. In the reduced system, each scalar field is paired with another, implying that all multipole modes of the perturbation can be regarded as symmetric partners, which can be understood from the property of the supersymmetry algebra. The generator of the supersymmetric structure is constructed from a closed conformal Killing vector field of ${\rm AdS}_2$. The associated conserved quantity allows one to show no static response, i.e., vanishing of the Love and dissipation numbers. We also discuss the vanishing Love numbers of the Kerr black hole with the nonzero dissipation numbers for the non-axisymmetric perturbations in terms of a radial constant found in a parallel manner as the axisymmetric field case even though the interpretation for the structure is controversial. The symmetric structure corresponds to the ``ladder'' symmetry in Hui et al. [JCAP 01, no.01, 032 (2022)] although the geometrical origin is different. Our ladder operator includes the generators of hidden symmetries in previous works. △ Less

Submitted 6 June, 2023; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: 41 pages, v2: added discussion on the result of the Kerr black hole in terms of the Teukolsky equation, v3: minor corrections, accepted for publication in PRD

Protostellar-disc fragmentation across all metallicities

Authors: Ryoki Matsukoba, Kei E. I. Tanaka, Kazuyuki Omukai, Eduard I. Vorobyov, Takashi Hosokawa

Abstract: Cosmic metallicity evolution possibly creates the diversity of star formation modes at different epochs. Gravitational fragmentation of circumstellar discs provides an important formation channel of multiple star systems, including close binaries. We here study the nature of disc fragmentation, systematically performing a suite of two-dimensional radiation-hydrodynamic simulations, in a broad rang… ▽ More Cosmic metallicity evolution possibly creates the diversity of star formation modes at different epochs. Gravitational fragmentation of circumstellar discs provides an important formation channel of multiple star systems, including close binaries. We here study the nature of disc fragmentation, systematically performing a suite of two-dimensional radiation-hydrodynamic simulations, in a broad range of metallicities, from the primordial to the solar values. In particular, we follow relatively long-term disc evolution over 15 kyr after the disc formation, incorporating the effect of heating by the protostellar irradiation. Our results show that the disc fragmentation occurs at all metallicities $1$--$0$ $Z_{\odot}$, yielding self-gravitating clumps. Physical properties of the clumps, such as their number and mass distributions, change with the metallicity due to different gas thermal evolution. For instance, the number of clumps is the largest for the intermediate metallicity range of $10^{-2}$--$10^{-5}$ $Z_{\odot}$, where the dust cooling is effective exclusively in a dense part of the disc and causes the fragmentation of spiral arms. The disc fragmentation is more modest for $1$--$0.1$ $Z_{\odot}$ thanks to the disc stabilization by the stellar irradiation. Such metallicity dependence agrees with the observed trend that the close binary fraction increases with decreasing metallicity in the range of $1$--$10^{-3}$ $Z_{\odot}$. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 17 pages, 16 figures, submitted to MNRAS

Impact of the cosmic background radiation on the initial mass function of metal-poor stars

Authors: Sunmyon Chon, Haruka Ono, Kazuyuki Omukai, Raffaella Schneider

Abstract: 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 protostella… ▽ More 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$. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 18 pages, 12 figures, MNRAS submitted

Merger rate density of binary black holes through isolated Population I, II, III and extremely metal-poor binary star evolution

Authors: Ataru Tanikawa, Takashi Yoshida, Tomoya Kinugawa, Alessandro A. Trani, Takashi Hosokawa, Hajime Susa, Kazuyuki Omukai

Abstract: We investigate the formation of merging binary black holes (BHs) through isolated binary evolution, performing binary population synthesis calculations covering an unprecedentedly wide metallicity range of Population (Pop) I, II, III, and extremely metal-poor (EMP) binary stars. We find that the predicted merger rate density and primary BH mass ($m_1$) distribution are consistent with the gravitat… ▽ More We investigate the formation of merging binary black holes (BHs) through isolated binary evolution, performing binary population synthesis calculations covering an unprecedentedly wide metallicity range of Population (Pop) I, II, III, and extremely metal-poor (EMP) binary stars. We find that the predicted merger rate density and primary BH mass ($m_1$) distribution are consistent with the gravitational wave (GW) observations. Notably, Pop III and EMP ($< 10^{-2}$ $Z_\odot$) binary stars yield most of the pair instability (PI) mass gap events with $m_1 = 65$--$130$ $M_\odot$. Pop III binary stars contribute more to the PI mass gap events with increasing redshift, and all the PI mass gap events have the Pop III origin at redshifts $\gtrsim 8$. Our result can be assessed by future GW observations in the following two points. First, there are no binary BHs with $m_1=100$--$130$ $M_\odot$ in our result, and thus the $m_1$ distribution should suddenly drop in the range of $m_1=100$--$130$ $M_\odot$. Second, the PI mass gap event rate should increase toward higher redshift up to $\sim 11$, since those events mainly originate from the Pop III binary stars. We find that the following three assumptions are needed to reproduce the current GW observations: a top-heavy stellar initial mass function and the presence of close binary stars for Pop III and EMP binary stars, and inefficient convective overshoot in the main-sequence phase of stellar evolution. Without any of the above, the number of PI mass gap events becomes too low to reproduce current GW observations. △ Less

Submitted 2 March, 2022; v1 submitted 20 October, 2021; originally announced October 2021.

Comments: 26 pages, 13 figures, 2 tables, ApJ accepted. BSEEMP, a binary population synthesis code used in this paper, is published on github: https://github.com/atrtnkw/bseemp

Journal ref: The Astrophysical Journal, 926, 83, 2022 February

Radiation hydrodynamics simulations of line-driven AGN disc winds: metallicity dependence and black hole growth

Authors: Mariko Nomura, Kazuyuki Omukai, Ken Ohsuga

Abstract: Growth of the black holes (BHs) from the seeds to supermassive BHs (SMBHs, $\sim\!10^9\,M_\odot$) is not understood, but the mass accretion must have played an important role. We performed two-dimensional radiation hydrodynamics simulations of line-driven disc winds considering the metallicity dependence in a wide range of the BH mass, and investigated the reduction of the mass accretion rate due… ▽ More Growth of the black holes (BHs) from the seeds to supermassive BHs (SMBHs, $\sim\!10^9\,M_\odot$) is not understood, but the mass accretion must have played an important role. We performed two-dimensional radiation hydrodynamics simulations of line-driven disc winds considering the metallicity dependence in a wide range of the BH mass, and investigated the reduction of the mass accretion rate due to the wind mass loss. Our results show that denser and faster disc winds appear at higher metallicities and larger BH masses. The accretion rate is suppressed to $\sim\! 0.4$--$0.6$ times the mass supply rate to the disc for the BH mass of $M_{\rm BH}\gtrsim 10^5\,M_{\odot}$ in high-metallicity environments of $Z\gtrsim Z_\odot$, while the wind mass loss is negligible when the metallicity is sub-solar ($\sim 0.1Z_\odot$). By developing a semi-analytical model, we found that the metallicity dependence of the line force and the BH mass dependence of the surface area of the wind launch region are the cause of the metallicity dependence ($\propto\! Z^{2/3}$) and BH mass dependencies ($\propto\! M_{\rm BH}^{4/3}$ for $M_{\rm BH}\leq 10^6\,M_\odot$ and $\propto\! M_{\rm BH}$ for $M_{\rm BH}\geq 10^6\,M_\odot$) of the mass-loss rate. Our model suggests that the growth of BHs by the gas accretion effectively slows down in the regime $\gtrsim 10^{5}M_\odot$ in metal-enriched environments $\gtrsim Z_\odot$. This means that the line-driven disc winds may have an impact on late evolution of SMBHs. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Light, medium-weight or heavy? The nature of the first supermassive black hole seeds

Authors: F. Sassano, R. Schneider, R. Valiante, K. Inayoshi, S. Chon, K. Omukai, L. Mayer, P. R. Capelo

Abstract: Observations of hyper-luminous quasars at $z>6$ reveal the rapid growth of supermassive black holes (SMBHs $>10^9 \rm M_{\odot}$) whose origin is still difficult to explain. Their progenitors may have formed as remnants of massive, metal free stars (light seeds), via stellar collisions (medium-weight seeds) and/or massive gas clouds direct collapse (heavy seeds). In this work we investigate for th… ▽ More Observations of hyper-luminous quasars at $z>6$ reveal the rapid growth of supermassive black holes (SMBHs $>10^9 \rm M_{\odot}$) whose origin is still difficult to explain. Their progenitors may have formed as remnants of massive, metal free stars (light seeds), via stellar collisions (medium-weight seeds) and/or massive gas clouds direct collapse (heavy seeds). In this work we investigate for the first time the relative role of these three seed populations in the formation of $z>6$ SMBHs within an Eddington-limited gas accretion scenario. To this aim, we implement in our semi-analytical data-constrained model a statistical description of the spatial fluctuations of Lyman-Werner (LW) photo-dissociating radiation and of metal/dust enrichment. This allows us to set the physical conditions for BH seeds formation, exploring their relative birth rate in a highly biased region of the Universe at $z>6$. We find that the inclusion of medium-weight seeds does not qualitatively change the growth history of the first SMBHs: although less massive seeds ($<10^3 \rm M_\odot$) form at a higher rate, the mass growth of a $\sim 10^9 \rm M_\odot$ SMBH at $z<15$ is driven by efficient gas accretion (at a sub-Eddington rate) onto its heavy progenitors ($10^5 \rm M_\odot$). This conclusion holds independently of the critical level of LW radiation and even when medium-weight seeds are allowed to form in higher metallicity galaxies, via the so-called super-competitive accretion scenario. Our study suggests that the genealogy of $z \sim 6$ SMBHs is characterized by a rich variety of BH progenitors, which represent only a small fraction ($< 10 - 20\%$) of all the BHs that seed galaxies at $z > 15$. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: (21 pages, 18 figures, 2 tables. Accepted for publication in MNRAS)

Magnetohydrodynamic effect on first star formation: prestellar core collapse and protostar formation

Authors: Kenji Eric Sadanari, Kazuyuki Omukai, Kazuyuki Sugimura, Tomoaki Matsumoto, Kengo Tomida

Abstract: Recent theoretical studies have suggested that a magnetic field may play a crucial role in the first star formation in the universe. However, the influence of the magnetic field on the first star formation has yet to be understood well. In this study, we perform three-dimensional magnetohydrodynamic simulations taking into account all the relevant cooling processes and non-equilibrium chemical rea… ▽ More Recent theoretical studies have suggested that a magnetic field may play a crucial role in the first star formation in the universe. However, the influence of the magnetic field on the first star formation has yet to be understood well. In this study, we perform three-dimensional magnetohydrodynamic simulations taking into account all the relevant cooling processes and non-equilibrium chemical reactions up to the protostar density, in order to study the collapse of magnetized primordial gas cores with self-consistent thermal evolution. Our results show that the thermal evolution of the central core is hardly affected by a magnetic field, because magnetic forces do not prevent the contraction along the fields lines. We also find that the magnetic braking extracts the angular momentum from the core and suppresses fragmentation depending on the initial strength of the magnetic field. The angular momentum transport by the magnetic outflows is less effective than that by the magnetic braking because the outflows are launched only in a late phase of the collapse. Our results indicate that the magnetic effects become important for the field strength $B> 10^{-8}(n_{\rm H}/1\ \rm cm^{-3})^{2/3}\ \rm G$, where $n_{\rm H}$ is the number density, during the collapse phase. Finally, we compare our results with simulations using a barotropic approximation and confirm that this approximation is reasonable at least for the collapse phase. Nevertheless, self-consistent treatment of the thermal and chemical processes is essential for extending simulations to the accretion phase, in which radiative feedback by protostars plays a crucial role. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 18 pages, 15 figures, 2 tables. Accepted for publication in MNRAS

Formation of the first galaxies in the aftermath of the first supernovae

Authors: Makito Abe, Hidenobu Yajima, Sadegh Khochfar, Claudio Dalla Vecchia, Kazuyuki Omukai

Abstract: We perform high-resolution cosmological hydrodynamic simulations to study the formation of the first galaxies that reach the masses of $10^{8-9}~h^{-1}~M_\odot$ at $z=9$. The resolution of the simulations is high enough to resolve minihaloes and allow us to successfully pursue the formation of multiple Population (Pop) III stars, their supernova (SN) explosions, resultant metal-enrichment of the i… ▽ More We perform high-resolution cosmological hydrodynamic simulations to study the formation of the first galaxies that reach the masses of $10^{8-9}~h^{-1}~M_\odot$ at $z=9$. The resolution of the simulations is high enough to resolve minihaloes and allow us to successfully pursue the formation of multiple Population (Pop) III stars, their supernova (SN) explosions, resultant metal-enrichment of the inter-galactic medium (IGM) in the course of the build-up of the system. Metals are ejected into the IGM by multiple Pop III SNe, but some of the metal-enriched gas falls back onto the halo after $\gtrsim 100~\rm Myr$. The star formation history of the first galaxy depends sensitively on the initial mass function (IMF) of Pop III stars. The dominant stellar population transits from Pop III to Pop II at $z\sim 12-15$ in the case of power-law Pop III IMF, ${\rm d}n/{\rm d}M \propto M^{-2.35}$ with the mass range $10-500~M_\odot$. At $z\lesssim 12$, stars are stably formed in the first galaxies with a star formation rate of $\sim 10^{-3}$-$10^{-1}~M_\odot/{\rm yr}$. In contrast, for the case with a flat IMF, gas-deprived first galaxies form due to frequent Pop III pair-instability SNe, resulting in the suppression of subsequent Pop II star formation. In addition, we calculate UV continuum, Ly$α$- and H$α$-line fluxes from the first galaxies. We show that the James Webb Space Telescope will be able to detect both UV continuum, Ly$α$ and H$α$ line emission from first galaxies with halo mass $\gtrsim 10^{9}~M_\odot$ at $z \gtrsim 10$. △ Less

Submitted 7 September, 2021; v1 submitted 6 May, 2021; originally announced May 2021.

Comments: 14 pages, 15 figures, accepted for publication in MNRAS

Signatures of Hierarchical Mergers in Black Hole Spin and Mass distribution

Authors: Hiromichi Tagawa, Zoltán Haiman, Imre Bartos, Bence Kocsis, Kazuyuki Omukai

Abstract: Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($χ_\mathrm{eff}$), the precession spin parameter (… ▽ More Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($χ_\mathrm{eff}$), the precession spin parameter ($χ_\mathrm{p}$), and the chirp mass ($m_\mathrm{chirp}$) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with ${\bar χ}_\mathrm{typ}\equiv \overline{(χ_\mathrm{eff}^2+χ_\mathrm{p}^2)^{1/2}}$, up to roughly the maximum $m_\mathrm{chirp}$ among first-generation (1g) BHs, and (ii) a plateau at ${\bar χ}_\mathrm{typ}\sim 0.6$ at higher $m_\mathrm{chirp}$. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from ${\bar χ}_\mathrm{typ}$ as a function of $m_\mathrm{chirp}$. The GW data observed in LIGO/Virgo O1--O3a prefers an increase in ${\bar χ}_\mathrm{typ}$ at low $m_\mathrm{chirp}$, which is consistent with the growth of the BH spin magnitude by hierarchical mergers, at $\sim 2 σ$ confidence. A Bayesian analysis suggests that 1g BHs have the maximum mass of $\sim 15$--$30\,M_\odot$ if the majority of mergers are of high-generation BHs (not among 1g-1g BHs), which is consistent with mergers in active galactic nucleus disks and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favored to have non-zero spin magnitudes of $\sim 0.3$. We also forecast that signatures for hierarchical mergers in the ${\bar χ}_\mathrm{typ}$ distribution can be confidently recovered once the number of GW events increases to $\gtrsim O(100)$. △ Less

Submitted 10 August, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 20 pages, 10 figures, accepted in MNRAS

Gravitational Wave Physics and Astronomy in the nascent era

Authors: Makoto Arimoto, Hideki Asada, Michael L. Cherry, Michiko S. Fujii, Yasushi Fukazawa, Akira Harada, Kazuhiro Hayama, Takashi Hosokawa, Kunihito Ioka, Yoichi Itoh, Nobuyuki Kanda, Koji S. Kawabata, Kyohei Kawaguchi, Nobuyuki Kawai, Tsutomu Kobayashi, Kazunori Kohri, Yusuke Koshio, Kei Kotake, Jun Kumamoto, Masahiro N. Machida, Hideo Matsufuru, Tatehiro Mihara, Masaki Mori, Tomoki Morokuma, Shinji Mukohyama , et al. (28 additional authors not shown)

Abstract: The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion,… ▽ More The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion, we review the fundamental outcomes and prospects of gravitational wave physics and astronomy. We survey the development focusing on representative sources of gravitational waves: binary black holes, binary neutron stars, and supernovae. We also summarize the role of gravitational wave observations as a probe of new physics. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 89pages, Prog. Theor. Exp. Phys. accepted

Transition of the initial mass function in the metal-poor environments

Authors: Sunmyon Chon, Kazuyuki Omukai, Raffaella Schneider

Abstract: 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… ▽ More 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. △ Less

Submitted 25 August, 2021; v1 submitted 8 March, 2021; originally announced March 2021.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

Ionization degree and magnetic diffusivity in star-forming clouds with different metallicities

Authors: Daisuke Nakauchi, Kazuyuki Omukai, Hajime Susa

Abstract: Magnetic fields play such essential roles in star formation as transporting angular momentum and driving outflows from a star-forming cloud, thereby controlling the formation efficiency of a circumstellar disc and also multiple stellar systems. The coupling of magnetic fields to the gas depends on its ionization degree. We calculate the temperature evolution and ionization degree of a cloud for va… ▽ More Magnetic fields play such essential roles in star formation as transporting angular momentum and driving outflows from a star-forming cloud, thereby controlling the formation efficiency of a circumstellar disc and also multiple stellar systems. The coupling of magnetic fields to the gas depends on its ionization degree. We calculate the temperature evolution and ionization degree of a cloud for various metallicities of Z/Zsun = 1e-6, 1e-5, 1e-4, 1e-3, 1e-2, 1e-1, and 1. We update the chemical network by reversing all the gas-phase processes and by considering grain-surface chemistry, including grain evaporation, thermal ionization of alkali metals, and thermionic emission from grains. The ionization degree at nH ~ 1e15-1e19 /cm^3 becomes up to eight orders of magnitude higher than that obtained in the previous model, owing to the thermionic emission and thermal ionization of K and Na, which have been neglected so far. Although magnetic fields dissipate owing to ambipolar diffusion or Ohmic loss at nH < 1e15 /cm^3, the fields recover strong coupling to the gas at nH ~ 1e15 /cm^3, which is lower by a few orders of magnitude compared to the previous work. We develop a reduced chemical network by choosing processes relevant to major coolants and charged species. The reduced network consists of 104 (161) reactions among 28 (38) species in the absence (presence, respectively) of ionization sources. The reduced model includes H2 and HD formation on grain surfaces as well as the depletion of O, C, OH, CO, and H2O on grain surfaces. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 23 pages, 12 figures, 3 tables, accepted for publication in MNRAS

Mass-gap Mergers in Active Galactic Nuclei

Authors: Hiromichi Tagawa, Bence Kocsis, Zoltan Haiman, Imre Bartos, Kazuyuki Omukai, Johan Samsing

Abstract: The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins that could be outside of the range expected from isolated stellar evolution. These merging objects could have undergone previous mergers. Such hierarchical mergers are predicted to be frequent in active galactic nuclei (AGN) disks, where binaries form and evolve efficien… ▽ More The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins that could be outside of the range expected from isolated stellar evolution. These merging objects could have undergone previous mergers. Such hierarchical mergers are predicted to be frequent in active galactic nuclei (AGN) disks, where binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Here we compare the properties of these observed events to the theoretical models of mergers in AGN disks, which are obtained by performing one-dimensional $N$-body simulations combined with semi-analytical prescriptions. The high BH masses in GW190521 are consistent with mergers of high-generation (high-g) BHs where the initial progenitor stars had high metallicity, 2g BHs if the original progenitors were metal-poor, or 1g BHs that had gained mass via super-Eddington accretion. Other measured properties related to spin parameters in GW190521 are also consistent with mergers in AGN disks. Furthermore, mergers in the lower mass gap or those with low mass ratio as found in GW190814 and GW190412 are also reproduced by mergers of 2g-1g or 1g-1g objects with significant accretion in AGN disks. Finally, due to gas accretion, the massive neutron star merger reported in GW190425 can be produced in an AGN disk. △ Less

Submitted 20 December, 2020; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 12 pages, 4 figures, accepted in ApJ

Disk fragmentation and intermittent accretion onto supermassive stars

Authors: Ryoki Matsukoba, Eduard I. Vorobyov, Kazuyuki Sugimura, Sunmyon Chon, Takashi Hosokawa, Kazuyuki Omukai

Abstract: 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 i… ▽ More 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. △ Less

Submitted 4 November, 2020; originally announced November 2020.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

Eccentric Black Hole Mergers in Active Galactic Nuclei

Authors: Hiromichi Tagawa, Bence Kocsis, Zoltan Haiman, Imre Bartos, Kazuyuki Omukai, Johan Samsing

Abstract: The astrophysical origin of gravitational wave (GW) transients is a timely open question in the wake of discoveries by LIGO/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and t… ▽ More The astrophysical origin of gravitational wave (GW) transients is a timely open question in the wake of discoveries by LIGO/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and the number of suspected hierarchical mergers observed by LIGO/Virgo. The binary eccentricity distribution can provide further information to distinguish between astrophysical models. Here we derive the eccentricity distribution of BH mergers in AGN disks. We find that eccentricity is mainly due to binary-single (BS) interactions, which lead to most BH mergers in AGN disks having a significant eccentricity at $0.01\,\mathrm{Hz}$, detectable by LISA. If BS interactions occur in isotropic-3D directions, then $8$--$30\%$ of the mergers in AGN disks will have eccentricities at $10\,\mathrm{Hz}$ above $e_{10\,\rm Hz}\gtrsim 0.03$, detectable by LIGO/Virgo/KAGRA, while $5$--$17\%$ of mergers have $e_{10\,\rm Hz}\geq 0.3$. On the other hand, if BS interactions are confined to the AGN-disk plane due to torques from the disk, with 1-20 intermediate binary states during each interaction, or if BHs can migrate to $\lesssim10^{-3}\,\mathrm{pc}$ from the central supermassive black hole, then $10$--$70\%$ of the mergers will be highly eccentric ($e_{10\,\rm Hz} \geq 0.3$), consistent with the possible high eccentricity in GW190521. △ Less

Submitted 18 December, 2020; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 9 pages, 3 figures, accepted in ApJL

Pulsation-driven mass loss from massive stars behind stellar mergers in metal-poor dense clusters

Authors: Daisuke Nakauchi, Kohei Inayoshi, Kazuyuki Omukai

Abstract: The recent discovery of high-redshift (z > 6) supermassive black holes (SMBH) favors the formation of massive seed BHs in protogalaxies. One possible scenario is formation of massive stars ~ 1e3-1e4 Msun via runaway stellar collisions in a dense cluster, leaving behind massive BHs without significant mass loss. We study the pulsational instability of massive stars with the zero-age main-sequence (… ▽ More The recent discovery of high-redshift (z > 6) supermassive black holes (SMBH) favors the formation of massive seed BHs in protogalaxies. One possible scenario is formation of massive stars ~ 1e3-1e4 Msun via runaway stellar collisions in a dense cluster, leaving behind massive BHs without significant mass loss. We study the pulsational instability of massive stars with the zero-age main-sequence (ZAMS) mass Mzams/Msun = 300-3000 and metallicity Z/Zsun = 0-0.1, and discuss whether or not pulsation-driven mass loss prevents massive BH formation. In the MS phase, the pulsational instability excited by the epsilon-mechanism grows in ~ 1e3 yrs. As the stellar mass and metallicity increase, the mass-loss rate increases to < 1e-3 Msun/yr. In the red super-giant (RSG) phase, the instability is excited by the kappa-mechanism operating in the hydrogen ionization zone and grows more rapidly in ~ 10 yrs. The RSG mass-loss rate is almost independent of metallicity and distributes in the range of ~ 1e-3-1e-2 Msun/yr. Conducting the stellar structure calculations including feedback due to pulsation-driven winds, we find that the stellar models of Mzams/Msun = 300-3000 can leave behind remnant BHs more massive than ~ 200-1200 Msun. We conclude that massive merger products can seed monster SMBHs observed at z > 6. △ Less

Submitted 31 August, 2020; originally announced August 2020.

Comments: 15 pages, 10 figures, accepted for publication in the Astrophysical Journal

Cosmological DCBH formation sites hostile for their growth

Authors: Sunmyon Chon, Takashi Hosokawa, Kazuyuki Omukai

Abstract: 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,… ▽ More 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. △ Less

Submitted 10 September, 2020; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: 13 pages, 8 figures, submitted to MNRAS

Accretion bursts in low-metallicity protostellar disks

Authors: E. I. Vorobyov, V. G. Elbakyan, K. Omukai, T. Hosokawa, R. Matsukoba, M. Guedel

Abstract: The early evolution of protostellar disks with metallicities in the $Z=1.0-0.01~Z_\odot$ range was studied with a particular emphasis on the strength of gravitational instability and the nature of protostellar accretion in low-metallicity systems. Numerical hydrodynamics simulations in the thin-disk limit were employed that feature separate gas and dust temperatures, and disk mass-loading from the… ▽ More The early evolution of protostellar disks with metallicities in the $Z=1.0-0.01~Z_\odot$ range was studied with a particular emphasis on the strength of gravitational instability and the nature of protostellar accretion in low-metallicity systems. Numerical hydrodynamics simulations in the thin-disk limit were employed that feature separate gas and dust temperatures, and disk mass-loading from the infalling parental cloud cores. Models with cloud cores of similar initial mass and rotation pattern, but distinct metallicity were considered to distinguish the effect of metallicity from that of initial conditions. The early stages of disk evolution in low-metallicity models are characterized by vigorous gravitational instability and fragmentation. Disk instability is sustained by continual mass-loading from the collapsing core. The time period that is covered by this unstable stage is much shorter in the $Z=0.01~Z_\odot$ models as compared to their higher metallicity counterparts thanks to the higher mass infall rates caused by higher gas temperatures (that decouple from lower dust temperatures) in the inner parts of collapsing cores. Protostellar accretion rates are highly variable in the low-metallicity models reflecting a highly dynamical nature of the corresponding protostellar disks. The low-metallicity systems feature short, but energetic episodes of mass accretion caused by infall of inward-migrating gaseous clumps that form via gravitational fragmentation of protostellar disks. These bursts seem to be more numerous and last longer in the $Z=0.1~Z_\odot$ models in comparison to the $Z=0.01~Z_\odot$ case. Variable protostellar accretion with episodic bursts is not a particular feature of solar metallicity disks. It is also inherent to gravitationally unstable disks with metallicities up to 100 times lower than solar. △ Less

Submitted 30 June, 2020; originally announced June 2020.

Comments: Accepted for publication by Astronomy & Astrophysics

Star cluster formation and cloud dispersal by radiative feedback: dependence on metallicity and compactness

Authors: Hajime Fukushima, Hidenobu Yajima, Kazuyuki Sugimura, Takashi Hosokawa, Kazuyuki Omukai, Tomoaki Matsumoto

Abstract: We study star cluster formation in various environments with different metallicities and column densities by performing a suite of three-dimensional radiation hydrodynamics simulations. We find that the photoionization feedback from massive stars controls the star formation efficiency (SFE) in a star-forming cloud, and its impact sensitively depends on the gas metallicity $Z$ and initial cloud sur… ▽ More We study star cluster formation in various environments with different metallicities and column densities by performing a suite of three-dimensional radiation hydrodynamics simulations. We find that the photoionization feedback from massive stars controls the star formation efficiency (SFE) in a star-forming cloud, and its impact sensitively depends on the gas metallicity $Z$ and initial cloud surface density $Σ$. At $Z=1~Z_{\odot}$, SFE increases as a power law from 0.03 at $Σ= 10~M_{\odot}{\rm pc^{-2}}$ to 0.3 at $Σ= 300~M_{\odot}{\rm pc^{-2}}$. In low-metallicity cases $10^{-2}- 10^{-1} Z_{\odot}$, star clusters form from atomic warm gases because the molecule formation time is not short enough with respect to the cooling or dynamical time. In addition, the whole cloud is disrupted more easily by expanding H{\sc ii} bubbles which have higher temperature owing to less efficient cooling. With smaller dust attenuation, the ionizing radiation feedback from nearby massive stars is stronger and terminate star formation in dense clumps. These effects result in inefficient star formation in low-metallicity environments: the SFE drops by a factor of $\sim 3$ at $Z=10^{-2}~Z_{\odot}$ compared to the results for $Z=1~Z_{\odot}$, regardless of $Σ$. Newborn star clusters are also gravitationally less bound. We further develop a new semi-analytical model that can reproduce the simulation results well, particularly the observed dependencies of the SFEs on the cloud surface densities and metallicities. △ Less

Submitted 24 August, 2020; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS

Thermal evolution of protoplanetary disks: from $β$-cooling to decoupled gas and dust temperatures

Authors: Eduard I. Vorobyov, Ryoki Matsukoba, Kazuyuki Omukai, Manuel Guedel

Abstract: Aims: We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods: Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified $β$-cooling approach with and without irradia… ▽ More Aims: We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods: Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified $β$-cooling approach with and without irradiation, in which the rate of disk cooling is proportional to the local dynamical time, a fiducial model with equal dust and gas temperatures calculated taking viscous heating, irradiation, and radiative cooling into account, and also a more sophisticated approach allowing decoupled dust and gas temperatures. Results: We found that the gas temperature may significantly exceed that of dust in the outer regions of young disks thanks to additional compressional heating caused by the infalling envelope material in the early stages of disk evolution and slow collisional exchange of energy between gas and dust in low-density disk regions. The outer envelope however shows an inverse trend with the gas temperatures dropping below that of dust. The global disk evolution is only weakly sensitive to temperature decoupling. Nevertheless, separate dust and gas temperatures may affect the chemical composition, dust evolution, and disk mass estimates. Constant-$β$ models without stellar and background irradiation fail to reproduce the disk evolution with more sophisticated thermal schemes because of intrinsically variable nature of the $β$-parameter. Constant-$β$ models with irradiation can better match the dynamical and thermal evolution, but the agreement is still incomplete. Conclusions: Models allowing separate dust and gas temperatures are needed when emphasis is placed on the chemical or dust evolution in protoplanetary disks, particularly in sub-solar metallicity environments. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Comments: 19 pages, 21 Figures, accepeted for publication in Astronomy & Astrophysics

Formation of massive stars under protostellar radiation feedback: Very metal-poor stars

Authors: Hajime Fukushima, Takashi Hosokawa, Gen Chiaki, Kazuyuki Omukai, Naoki Yoshida, Rolf Kuiper

Abstract: We study the formation of very metal-poor stars under protostellar radiative feedback effect. We use cosmological simulations to identify low-mass dark matter halos and star-forming gas clouds within them. We then follow protostar formation and the subsequent long-term mass accretion phase of over one million years using two-dimensional radiation-hydrodynamics simulations. We show that the critica… ▽ More We study the formation of very metal-poor stars under protostellar radiative feedback effect. We use cosmological simulations to identify low-mass dark matter halos and star-forming gas clouds within them. We then follow protostar formation and the subsequent long-term mass accretion phase of over one million years using two-dimensional radiation-hydrodynamics simulations. We show that the critical physical process that sets the final mass is formation and expansion of a bipolar HII region. The process is similar to the formation of massive primordial stars, but radiation pressure exerted on dust grains also contributes to halting the accretion flow in the low-metallicity case. We find that the net feedback effect in the case with metallicity $Z = 10^{-2}~Z_{\odot}$ is stronger than in the case with $Z \sim 1~Z_{\odot}$. With decreasing metallicity, the radiation pressure effect becomes weaker, but photoionization heating of the circumstellar gas is more efficient owing to the reduced dust attenuation. In the case with $Z = 10^{-2}~Z_{\odot}$, the central star grows as massive as 200 solar-masses, similarly to the case of primordial star formation. We conclude that metal-poor stars with a few hundred solar masses can be formed by gas accretion despite the strong radiative feedback. △ Less

Submitted 22 April, 2021; v1 submitted 5 April, 2020; originally announced April 2020.

Comments: 18 pages, 13 figures, accepted for publication in MNRAS

The Birth of a Massive First-Star Binary

Authors: Kazuyuki Sugimura, Tomoaki Matsumoto, Takashi Hosokawa, Shingo Hirano, Kazuyuki Omukai

Abstract: We study the formation of massive Population III binary stars using a newly developed radiation hydrodynamics code with the adaptive mesh refinement and adaptive ray-tracing methods. We follow the evolution of a typical primordial star-forming cloud obtained from a cosmological hydrodynamics simulation. Several protostars form as a result of disk fragmentation and grow in mass by the gas accretion… ▽ More We study the formation of massive Population III binary stars using a newly developed radiation hydrodynamics code with the adaptive mesh refinement and adaptive ray-tracing methods. We follow the evolution of a typical primordial star-forming cloud obtained from a cosmological hydrodynamics simulation. Several protostars form as a result of disk fragmentation and grow in mass by the gas accretion, which is finally quenched by the radiation feedback from the protostars. Our code enables us, for the first time, to consider the feedback by both the ionizing and dissociating radiation from the multiple protostars, which is essential for self-consistently determining their final masses. At the final step of the simulation, we observe a very wide ($\gtrsim 10^4\,\mathrm{au}$) binary stellar system consisting of $60$ and $70\,M_\odot$ stars. One of the member stars also has two smaller mass ($10\,M_\odot$) companion stars orbiting at $200$ and $800\,\mathrm{au}$, making up a mini-triplet system. Our results suggest that massive binary or multiple systems are common among Population III stars. △ Less

Submitted 13 June, 2020; v1 submitted 31 January, 2020; originally announced February 2020.

Comments: 6 pages, 3 figures, published in ApJL, movie available at https://www.youtube.com/watch?v=794O0yGWGp0

Journal ref: ApJL 892, L14 (2020)

Supermassive Star Formation via Super Competitive Accretion in Slightly Metal-enriched Clouds

Authors: Sunmyon Chon, Kazuyuki Omukai

Abstract: 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, howeve… ▽ More 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. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 10 pages, 6figures, submitted to MNRAS

Ionization degree and magnetic diffusivity in the primordial star-forming clouds

Authors: Daisuke Nakauchi, Kazuyuki Omukai, Hajime Susa

Abstract: Magnetic fields play such roles in star formation as the angular momentum transport in star-forming clouds, thereby controlling circumstellar disc formation and even binary star formation efficiency. The coupling between the magnetic field and gas is determined by the ionization degree in the gas. Here, we calculate the thermal and chemical evolution of the primordial gas by solving chemical react… ▽ More Magnetic fields play such roles in star formation as the angular momentum transport in star-forming clouds, thereby controlling circumstellar disc formation and even binary star formation efficiency. The coupling between the magnetic field and gas is determined by the ionization degree in the gas. Here, we calculate the thermal and chemical evolution of the primordial gas by solving chemical reaction network where all the reactions are reversed. We find that at ~ 10^14-10^18 /cm^3, the ionization degree becomes 100-1000 times higher than the previous results due to the lithium ionization by thermal photons trapped in the cloud, which has been omitted so far. We construct the minimal chemical network which can reproduce correctly the ionization degree as well as the thermal evolution by extracting 36 reactions among 13 species. Using the obtained ionization degree, we evaluate the magnetic field diffusivity. We find that the field dissipation can be neglected for global fields coherent over > a tenth of the cloud size as long as the field is not so strong as to prohibit the collapse. With magnetic fields strong enough for ambipolar diffusion heating to be significant, the magnetic pressure effects to slow down the collapse and to reduce the compressional heating become more important, and the temperature actually becomes lower than in the no-field case. △ Less

Submitted 17 April, 2019; originally announced April 2019.

Comments: 17 pages, 13 figures, 2 tables, submitted to MNRAS

Spectral Energy Distribution of the First Galaxies: Contribution from Pre-Main-Sequence Stars

Authors: Hiroto Mitani, Naoki Yoshida, Kazuyuki Omukai, Takashi Hosokawa

Abstract: One of the major goals of next-generation space-borne and ground-based telescopes is to detect and characterize the first galaxies that were in place in the first few hundred million years after the big bang. We study the spectral energy distribution (SED) of the first galaxies and discuss the prospects for detection and identification. We consider very young star-forming galaxies at $z=15$ and in… ▽ More One of the major goals of next-generation space-borne and ground-based telescopes is to detect and characterize the first galaxies that were in place in the first few hundred million years after the big bang. We study the spectral energy distribution (SED) of the first galaxies and discuss the prospects for detection and identification. We consider very young star-forming galaxies at $z=15$ and incorporate the contribution from pre-main-sequence (PMS) stars. Unlike in the present-day galaxies, primordial protostars are not embedded in dusty gas clouds, and hence the light from them can be visible at a wide range of wavelengths. We use MESA code to follow the PMS evolution and use the BT-Settl model to calculate the SED of individual PMS stars. We show that PMS stars contribute to boost the flux in the mid-infrared, and that the galaxy SED at very early evolutionary phases is overall redder than at later phases. The infrared flux contribution is comparable to that caused by emission lines powered by massive stars. We argue that the contribution from PMS stars is important for characterizing young galaxies in the early Universe and also for target selection with future deep galaxy surveys. △ Less

Submitted 4 July, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Comments: 5 pages, 4 figures, published online in MNRAS letters

Gravitational stability and fragmentation condition for discs around accreting supermassive stars

Authors: Ryoki Matsukoba, Sanemichi Z. Takahashi, Kazuyuki Sugimura, Kazuyuki Omukai

Abstract: Supermassive stars (SMSs) with mass $\sim10^{5}~\rm{M}_{\odot}$ are promising candidates for the origin of supermassive black holes observed at redshift $\gtrsim6$. They are supposed to form as a result of rapid accretion of primordial gas, although it can be obstructed by the time variation caused by circum-stellar disc fragmentation due to gravitational instability. To assess the occurrence of f… ▽ More Supermassive stars (SMSs) with mass $\sim10^{5}~\rm{M}_{\odot}$ are promising candidates for the origin of supermassive black holes observed at redshift $\gtrsim6$. They are supposed to form as a result of rapid accretion of primordial gas, although it can be obstructed by the time variation caused by circum-stellar disc fragmentation due to gravitational instability. To assess the occurrence of fragmentation, we study the structure of marginally gravitationally unstable accretion discs, by using a steady one-dimensional thin disc model with detailed treatment of chemical and thermal processes. Motivated by two SMS formation scenarios, i.e., those with strong ultraviolet radiation background or with large velocity difference between the baryon and the dark matter, we consider two types of flows, i.e., atomic and molecular flows, respectively, for a wide range of the central stellar mass $10-10^5~\rm{M}_{\odot}$ and the accretion rate $10^{-3}-1~\rm{M}_{\odot}~\rm{yr}^{-1}$. In the case of a mostly atomic gas flowing to the disc outer boundary, the fragmentation condition is expressed as the accretion rate being higher than the critical value of $10^{-1}~\rm{M}_{\odot}~\rm{yr}^{-1}$ regardless of the central stellar mass. On the other hand, in the case of molecular flows, there is a critical disc radius outside of which the disc becomes unstable. Those conditions appears to be marginally satisfied according to numerical simulations, suggesting that disc fragmentation can be common during SMS formation. △ Less

Submitted 31 December, 2018; originally announced January 2019.

Comments: 17 pages, 10 figures, accepted for publication in MNRAS

Probing the High-Redshift Universe with SPICA: Toward the Epoch of Reionization and Beyond

Authors: E. Egami, S. Gallerani, R. Schneider, A. Pallottini, L. Vallini, E. Sobacchi, A. Ferrara, S. Bianchi, M. Bocchio, S. Marassi, L. Armus, L. Spinoglio, A. W. Blain, M. Bradford, D. L. Clements, H. Dannerbauer, J. A. Fernández-Ontiveros, E. González-Alfonso, M. J. Griffin, C. Gruppioni, H. Kaneda, K. Kohno, S. C. Madden, H. Matsuhara, P. Najarro , et al. (14 additional authors not shown)

Abstract: With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper… ▽ More With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $μ$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI. △ Less

Submitted 25 September, 2018; v1 submitted 20 September, 2018; originally announced September 2018.

Comments: 24 pages, 8 figures, accepted for publication in PASA (SPICA special issue)

Condition for low-mass star formation in shock-compressed metal-poor clouds

Authors: Daisuke Nakauchi, Kazuyuki Omukai, Raffaella Schneider

Abstract: Shocks may have been prevalent in the early Universe, associated with virialization and supernova explosions, etc. Here, we study thermal evolution and fragmentation of shock-compressed clouds, by using a one-zone model with detailed thermal and chemical processes. We explore a large range of initial density (1-1e5 /cm^3), metallicity (0-1e-2 Z_sun), UV strength (0-500 times Galactic value), and c… ▽ More Shocks may have been prevalent in the early Universe, associated with virialization and supernova explosions, etc. Here, we study thermal evolution and fragmentation of shock-compressed clouds, by using a one-zone model with detailed thermal and chemical processes. We explore a large range of initial density (1-1e5 /cm^3), metallicity (0-1e-2 Z_sun), UV strength (0-500 times Galactic value), and cosmic microwave background temperature (10 and 30 K). Shock-compressed clouds contract isobarically via atomic and molecular line cooling, until self-gravitating clumps are formed by fragmentation. If the metals are only in the gas-phase, the clump mass is higher than ~ 3 M_sun in any conditions we studied. Although in some cases with a metallicity higher than ~ 1e-3 Z_sun, re-fragmentation of a clump is caused by metal-line cooling, this fragment mass is higher than ~ 30 M_sun. On the other hand, if about half the mass of metals is condensed in dust grains, as in the Galactic interstellar medium, dust cooling triggers re-fragmentation of a clump into sub-solar mass pieces, for metallicities higher than ~ 1e-5 Z_sun. Therefore, the presence of dust is essential in low-mass (< M_sun) star formation from a shock-compressed cloud. △ Less

Submitted 14 July, 2018; originally announced July 2018.

Comments: 15 pages, 8 figures, accepted for publication in MNRAS

Condition for dust evacuation from the first galaxies

Authors: Hajime Fukushima, Hidenobu Yajima, Kazuyuki Omukai

Abstract: Dust enables low-mass stars to form from low-metallicity gas by inducing fragmentation of clouds via the cooling by its thermal emission. Dust may, however, be evacuated from star-forming clouds due to radiation force from massive stars. We here study the condition for the dust evacuation by comparing the dust evacuation time with the time of cloud destruction due to either expansion of HII region… ▽ More Dust enables low-mass stars to form from low-metallicity gas by inducing fragmentation of clouds via the cooling by its thermal emission. Dust may, however, be evacuated from star-forming clouds due to radiation force from massive stars. We here study the condition for the dust evacuation by comparing the dust evacuation time with the time of cloud destruction due to either expansion of HII regions or supernovae. The cloud destruction time has weak dependence on the cloud radius, while the dust evacuation time becomes shorter for a cloud with the smaller radius. The dust evacuation thus occurs in compact star-forming clouds whose column density is $N_{\rm H} \simeq 10^{24} - 10^{26} ~{\rm cm^{-2}}$. The critical halo mass above which the dust evacuation occurs becomes lower for higher formation redshift, e.g., $\sim 10^{9}~M_{\odot}$ at redshift $z \sim 3$ and $\sim 10^{7}~M_{\odot}$ at $z \sim 9$. In addition, metallicity of the gas should be less than $\sim 10^{-2} ~ Z_{\odot}$. Otherwise the dust attenuation reduces the radiation force significantly. From the dust-evacuated gas, massive stars are likely to form even with metallicity above $\sim 10^{-5}~Z_{\odot}$, the critical value for low-mass star formation due to the dust cooling. This can explain the dearth of ultra-metal poor stars with the metallicity lower than $\sim 10^{-4}~Z_{\odot}$. △ Less

Submitted 23 March, 2018; originally announced March 2018.

Comments: 16 pages, 18 figures, accepted for publication in MNRAS

Stunted accretion growth of black holes by combined effect of the flow angular momentum and radiation feedback

Authors: Kazuyuki Sugimura, Takashi Hosokawa, Hidenobu Yajima, Kohei Inayoshi, Kazuyuki Omukai

Abstract: Accretion on to seed black holes (BHs) is believed to play a crucial role in formation of supermassive BHs observed at high-redshift (z>6). Here, we investigate the combined effect of gas angular momentum and radiation feedback on the accretion flow, by performing 2D axially symmetric radiation hydrodynamics simulations that solve the flow structure across the Bondi radius and the outer part of th… ▽ More Accretion on to seed black holes (BHs) is believed to play a crucial role in formation of supermassive BHs observed at high-redshift (z>6). Here, we investigate the combined effect of gas angular momentum and radiation feedback on the accretion flow, by performing 2D axially symmetric radiation hydrodynamics simulations that solve the flow structure across the Bondi radius and the outer part of the accretion disc simultaneously. The accreting gas with finite angular momentum forms a rotationally-supported disc inside the Bondi radius, where the accretion proceeds by the angular momentum transport due to assumed alpha-type viscosity. We find that the interplay of radiation and angular momentum significantly suppresses accretion even if the radiative feedback is weakened in an equatorial shadowing region. The accretion rate is O(alpha)\sim O(0.01-0.1) times the Bondi value, where alpha is the viscosity parameter. By developing an analytical model, we show that such a great reduction of the accretion rate persists unless the angular momentum is so small that the corresponding centrifugal radius is \lesssim 0.04 times the Bondi radius. We argue that BHs are hard to grow quickly via rapid mass accretion considering the angular momentum barrier presented in this paper. △ Less

Submitted 21 September, 2018; v1 submitted 20 February, 2018; originally announced February 2018.

Comments: 17 pages, 12 figures, published in MNRAS

Journal ref: MNRAS, 478, 3961-3975 (2018)

Upper stellar mass limit by radiative feedback at low-metallicities: metallicity and accretion rate dependence

Authors: Hajime Fukushima, Kazuyuki Omukai, Takashi Hosokawa

Abstract: We investigate the upper stellar mass limit set by radiative feedback by the forming star with various accretion rates and metallicities. To this end, we numerically solve the structures of both a protostar and its surrounding accretion envelope assuming a spherical symmetric and steady flow. The optical depth of the dust cocoon, a dusty part of the accretion envelope, differs among the direct lig… ▽ More We investigate the upper stellar mass limit set by radiative feedback by the forming star with various accretion rates and metallicities. To this end, we numerically solve the structures of both a protostar and its surrounding accretion envelope assuming a spherical symmetric and steady flow. The optical depth of the dust cocoon, a dusty part of the accretion envelope, differs among the direct light from the stellar photosphere and the diffuse light re-emitted as dust thermal emission. As a result, varying the metallicity qualitatively changes the way that the radiative feedback suppresses the accretion flow. With a fixed accretion rate of $10^{-3}M_{\odot} {\rm yr^{-1}}$, the both direct and diffuse lights jointly operate to prevent the mass accretion at $Z \gtrsim 10^{-1}Z_{\odot}$. At $Z \lesssim 10^{-1}Z_{\odot}$, the diffuse light is no longer effective, and the direct light solely limits the mass accretion. At $Z \lesssim 10^{-3}Z_{\odot}$, the HII region formation plays an important role in terminating the accretion. The resultant upper mass limit increases with decreasing metallicity, from a few $\times~10~M_\odot$ to $\sim 10^3~M_\odot$ over $Z = 1Z_{\odot} - 10^{-4}~Z_\odot$. We also illustrate how the radiation spectrum of massive star-forming cores changes with decreasing metallicity. First, the peak wavelength of the spectrum, which is located around $30 μ{\rm m}$ at $1Z_{\odot}$, shifts to $< 3 μ{\rm m}$ at $Z \lesssim 0.1 Z_\odot$. Second, a characteristic feature at $10 μ\rm{m}$ due to the amorphous silicate band appears as a dip at $1Z_{\odot}$, but changes to a bump at $Z \lesssim 0.1 Z_{\odot}$. Using these spectral signatures, we can search massive accreting protostars in nearby low-metallicity environments with up-coming observations. △ Less

Submitted 1 October, 2017; originally announced October 2017.

Comments: 21 pages, 20 figures, accepted for publication in MNRAS

Primordial protostars accreting beyond the $ΩΓ$-limit: radiation effect around the star-disk boundary

Authors: Sanemichi Z. Takahashi, Kazuyuki Omukai

Abstract: We consider whether the maximum mass of first stars is imposed by the protostellar spin, i.e., by the so-called $ΩΓ$-limit, which requires the sum of the radiation and centrifugal forces at the stellar surface be smaller than the inward pull of the gravity. Once the accreting protostar reaches such a marginal state, the star cannot spin up more and is not allowed to accrete more gas with inward an… ▽ More We consider whether the maximum mass of first stars is imposed by the protostellar spin, i.e., by the so-called $ΩΓ$-limit, which requires the sum of the radiation and centrifugal forces at the stellar surface be smaller than the inward pull of the gravity. Once the accreting protostar reaches such a marginal state, the star cannot spin up more and is not allowed to accrete more gas with inward angular momentum flux. So far, however, the effect of stellar radiation on the structure of the accretion disk has not been properly taken into account in discussing the effect of $ΩΓ$-limit on the first star formation. Here, we obtain a series of the steady accretion-disk solutions considering such effect and find solutions without net angular momentum influx to the stars with arbitrary rotation rates, in addition to those with finite angular momentum flows. The accretion of positive angular momentum flows pushes the star beyond the $ΩΓ$-limit, which is allowed only with the external pressure provided by the circumstellar disk. On the other hand, the accretion with no net angular momentum influx does not result in the spin-up of the star. Thus, the existence of the solution with no net angular momentum influx indicates that the protostars can keep growing in mass by accretion even after they reach the $ΩΓ$-limit. △ Less

Submitted 1 August, 2017; originally announced August 2017.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Fates of the dense cores formed by fragmentation of filaments: do they fragment again or not?

Authors: Kazuyuki Sugimura, Yurina Mizuno, Tomoaki Matsumoto, Kazuyuki Omukai

Abstract: Fragmentation of filaments into dense cores is thought to be an important step in forming stars. The bar-mode instability of spherically collapsing cores found in previous linear analysis invokes a possibility of re-fragmentation of the cores due to their ellipsoidal (prolate or oblate) deformation. To investigate this possibility, here we perform three-dimensional self-gravitational hydrodynamics… ▽ More Fragmentation of filaments into dense cores is thought to be an important step in forming stars. The bar-mode instability of spherically collapsing cores found in previous linear analysis invokes a possibility of re-fragmentation of the cores due to their ellipsoidal (prolate or oblate) deformation. To investigate this possibility, here we perform three-dimensional self-gravitational hydrodynamics simulations that follow all the way from filament fragmentation to subsequent core collapse. We assume the gas is polytropic with index γ, which determines the stability of the bar-mode. For the case that the fragmentation of isolated hydrostatic filaments is triggered by the most unstable fragmentation mode, we find the bar mode grows as collapse proceeds if γ< 1.1, in agreement with the linear analysis. However, it takes more than ten orders-of-magnitude increase in the central density for the distortion to become non-linear. In addition to this fiducial case, we also study non-fiducial ones such as the fragmentation is triggered by a fragmentation mode with a longer wavelength and it occurs during radial collapse of filaments and find the distortion rapidly grows. In most of astrophysical applications, the effective polytropic index of collapsing gas exceeds 1.1 before ten orders-of-magnitude increase in the central density. Thus, supposing the fiducial case of filament fragmentation, re-fragmentation of dense cores would not be likely and their final mass would be determined when the filaments fragment. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 13 pages, 11 figures, 4 tables, submitted to MNRAS

Triggering the formation of direct collapse black holes by their congeners

Authors: Bin Yue, Andrea Ferrara, Fabio Pacucci, Kuzuyuki Omukai

Abstract: Direct collapse black holes (DCBHs) are excellent candidates as seeds of supermassive black holes (SMBHs) observed at $z \gsim 6$. The formation of a DCBH requires a strong external radiation field to suppress $\rm H_2$ formation and cooling in a collapsing gas cloud. Such strong field is not easily achieved by first stars or normal star-forming galaxies. Here we investigate a scenario in which th… ▽ More Direct collapse black holes (DCBHs) are excellent candidates as seeds of supermassive black holes (SMBHs) observed at $z \gsim 6$. The formation of a DCBH requires a strong external radiation field to suppress $\rm H_2$ formation and cooling in a collapsing gas cloud. Such strong field is not easily achieved by first stars or normal star-forming galaxies. Here we investigate a scenario in which the previously-formed DCBH can provide the necessary radiation field for the formation of additional ones. Using one-zone model and the simulated DCBH Spectral Energy Distributions (SEDs) filtered through absorbing gas initially having column density $N_{\rm H}$, we derive the critical field intensity, $J_{\rm LW}^{\rm crit}$, to suppress $\rm H_2$ formation and cooling. For the SED model with $N_{\rm H}=1.3\times10^{25}$ cm$^{-2}$, $8.0\times10^{24}$ cm$^{-2}$ and $5.0\times10^{24}$ cm$^{-2}$, we obtain $J_{\rm LW}^{\rm crit}\approx22$, 35 and 54, all much smaller than the critical field intensity for normal star-forming galaxies ($J_{\rm LW}^{\rm crit}\simgt 1000$). X-ray photons from previously-formed DCBHs build up a high-$z$ X-ray background (XRB) that may boost the $J_{\rm LW}^{\rm crit}$. However, we find that in the three SED models $J_{\rm LW}^{\rm crit}$ only increases to $\approx80$, 170 and 390 respectively even when $\dtρ_\bullet$ reaches the maximum value allowed by the present-day XRB level ($0.22, 0.034, 0.006~M_\odot$yr$^{-1}$Mpc$^{-3}$), still much smaller than the galactic value. Although considering the XRB from first galaxies may further increase $J_{\rm LW}^{\rm crit}$, we conclude that our investigation supports a scenario in which DCBH may be more abundant than predicted by models only including galaxies as external radiation sources. △ Less

Submitted 24 March, 2017; v1 submitted 23 December, 2016; originally announced December 2016.

Comments: 18 pages, 14 figures, 5 tables, ApJ in press

Do Stellar Winds Prevent the Formation of Supermassive Stars by Accretion?

Authors: Daisuke Nakauchi, Takashi Hosokawa, Kazuyuki Omukai, Hideyuki Saio, Ken'ichi Nomoto

Abstract: Supermassive stars (SMS; ~ 10^5 M_sun) formed from metal-free gas in the early Universe attract attention as progenitors of supermassive black holes observed at high redshifts. To form SMSs by accretion, central protostars must accrete at as high rates as ~ 0.1-1 M_sun/yr. Such protostars have very extended structures with bloated envelopes, like super-giant stars, and are called super-giant proto… ▽ More Supermassive stars (SMS; ~ 10^5 M_sun) formed from metal-free gas in the early Universe attract attention as progenitors of supermassive black holes observed at high redshifts. To form SMSs by accretion, central protostars must accrete at as high rates as ~ 0.1-1 M_sun/yr. Such protostars have very extended structures with bloated envelopes, like super-giant stars, and are called super-giant protostars (SGPSs). Under the assumption of hydrostatic equilibrium, SGPSs have density inverted layers, where the luminosity becomes locally super-Eddington, near the surface. If the envelope matter is allowed to flow out, however, a stellar wind could be launched and hinder the accretion growth of SGPSs before reaching the supermassive regime. We examine whether radiation-driven winds are launched from SGPSs by constructing steady and spherically symmetric wind solutions. We find that the wind velocity does not reach the escape velocity in any case considered. This is because once the temperature falls below ~ 10^4 K, the opacity plummet drastically owing to the recombination of hydrogen and the acceleration ceases suddenly. This indicates that, in realistic non-steady cases, even if outflows are launched from the surface of SGPSs, they would fall back again. Such a "wind" does not result in net mass loss and does not prevent the growth of SGPSs. In conclusion, SGPSs will grow to SMSs and eventually collapse to massive BHs of ~ 10^5 M_sun, as long as the rapid accretion is maintained. △ Less

Submitted 29 November, 2016; originally announced November 2016.

Comments: 11 pages, 13 figures, accepted for publication in MNRAS

Limits on Pop III star formation with the most iron-poor stars

Authors: M. de Bennassuti, S. Salvadori, R. Schneider, R. Valiante, K. Omukai

Abstract: We study the impact of star-forming mini-haloes, and the Initial Mass Function (IMF) of Population III (Pop III) stars, on the Galactic halo Metallicity Distribution Function (MDF) and on the properties of C-enhanced and C-normal stars at [Fe/H]<-3. For our investigation we use a data-constrained merger tree model for the Milky Way formation, which has been improved to self-consistently describe t… ▽ More We study the impact of star-forming mini-haloes, and the Initial Mass Function (IMF) of Population III (Pop III) stars, on the Galactic halo Metallicity Distribution Function (MDF) and on the properties of C-enhanced and C-normal stars at [Fe/H]<-3. For our investigation we use a data-constrained merger tree model for the Milky Way formation, which has been improved to self-consistently describe the physical processes regulating star-formation in mini-haloes, including the poor sampling of the Pop III IMF. We find that only when star-forming mini-haloes are included the low-Fe tail of the MDF is correctly reproduced, showing a plateau that is built up by C-enhanced metal-poor (CEMP) stars imprinted by primordial faint supernovae. The incomplete sampling of the Pop III IMF in inefficiently star-forming mini-haloes (< $10^{-3}$ $M_\odot$/yr) strongly limits the formation of Pair Instability Supernovae (PISNe), with progenitor masses $m_{\rm popIII}$=[140-260] $M_\odot$, even when a flat Pop III IMF is assumed. Second-generation stars formed in environments polluted at >50% level by PISNe are thus extremely rare, corresponding to $\approx$ 0.25% of the total stellar population at [Fe/H]<-2, which is consistent with recent observations. The low-Fe tail of the MDF strongly depends on the Pop III IMF shape and mass range. Given the current statistics, we find that a flat Pop III IMF model with $m_{\rm popIII}$=[10-300] $M_\odot$ is disfavoured by observations. We present testable predictions for Pop III stars extending down to lower masses, with $m_{\rm popIII}$=[0.1-300] $M_\odot$. △ Less

Submitted 20 October, 2016; v1 submitted 18 October, 2016; originally announced October 2016.

Comments: 15 pages, 11 figures. Accepted for publication in MNRAS. The only change is the correction of a mistake in the list of authors

Rapid Black Hole Growth under Anisotropic Radiation Feedback

Authors: Kazuyuki Sugimura, Takashi Hosokawa, Hidenobu Yajima, Kazuyuki Omukai

Abstract: Discovery of high-redshift (z > 6) supermassive black holes (BHs) may indicate that the rapid (or super-Eddington) gas accretion has aided their quick growth. Here, we study such rapid accretion of the primordial gas on to intermediate-mass (10^2 - 10^5 M_sun) BHs under anisotropic radiation feedback. We perform two-dimensional radiation hydrodynamics simulations that solve the flow structure acro… ▽ More Discovery of high-redshift (z > 6) supermassive black holes (BHs) may indicate that the rapid (or super-Eddington) gas accretion has aided their quick growth. Here, we study such rapid accretion of the primordial gas on to intermediate-mass (10^2 - 10^5 M_sun) BHs under anisotropic radiation feedback. We perform two-dimensional radiation hydrodynamics simulations that solve the flow structure across the Bondi radius, from far outside of the Bondi radius down to a central part which is larger than a circum-BH accretion disc. The radiation from the unresolved circum-BH disc is analytically modeled considering self-shadowing effect. We show that the flow settles into a steady state, where the flow structure consists of two distinct parts: (1) bipolar ionized outflowing regions, where the gas is pushed outward by thermal gas pressure and super-Eddington radiation pressure, and (2) an equatorial neutral inflowing region, where the gas falls toward the central BH without affected by radiation feedback. The resulting accretion rate is much higher than that in the case of isotropic radiation, far exceeding the Eddington-limited rate to reach a value slightly lower than the Bondi one. The opening angle of the equatorial inflowing region is determined by the luminosity and directional dependence of the central radiation. We find that photoevaporation from its surfaces set the critical opening angle of about ten degrees below which the accretion to the BH is quenched. We suggest that the shadowing effect allows even stellar-remnant BHs to grow rapidly enough to become high-redshift supermassive BHs. △ Less

Submitted 21 September, 2018; v1 submitted 11 October, 2016; originally announced October 2016.

Comments: 19 pages, 12 figures, 8 tables, corrected typos in the appendix

Journal ref: MNRAS, 469, 62-79 (2017)