Source author record

Kazuyuki Omukai

Kazuyuki Omukai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

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

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$.

preprint2022arXiv

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

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.

preprint2022arXiv

Protostellar-disc fragmentation across all metallicities

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}$.

preprint2021arXiv

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

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.

preprint2020arXiv

Cosmological DCBH formation sites hostile for their growth

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.

preprint2020arXiv

Eccentric Black Hole Mergers in Active Galactic Nuclei

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.

preprint2020arXiv

Mass-gap Mergers in Active Galactic Nuclei

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

The Birth of a Massive First-Star Binary

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Formation of Massive Primordial Stars: Intermittent UV Feedback with Episodic Mass Accretion

We present coupled stellar evolution (SE) and 3D radiation-hydrodynamic (RHD) simulations of the evolution of primordial protostars, their immediate environment, and the dynamic accretion history under the influence of stellar ionizing and dissociating UV feedback. Our coupled SE-RHD calculations result in a wide diversity of final stellar masses covering 10 Msun $\lesssim M_* \lesssim$ 1000 Msun. The formation of very massive ($\gtrsim$ 250 Msun) stars is possible under weak UV feedback, whereas ordinary massive (a few x 10 Msun) stars form when UV feedback can efficiently halt the accretion. This may explain the peculiar abundance pattern of a Galactic metal-poor star recently reported by Aoki et al. (2014), possibly the observational signature of very massive precursor primordial stars. Weak UV feedback occurs in cases of variable accretion, in particular when repeated short accretion bursts temporarily exceed 0.01 Msun/yr, causing the protostar to inflate. In the bloated state, the protostar has low surface temperature and UV feedback is suppressed until the star eventually contracts, on a thermal adjustment timescale, to create an HII region. If the delay time between successive accretion bursts is sufficiently short, the protostar remains bloated for extended periods, initiating at most only short periods of UV feedback. Disk fragmentation does not necessarily reduce the final stellar mass. Quite the contrary, we find that disk fragmentation enhances episodic accretion as many fragments migrate inward and are accreted onto the star, thus allowing continued stellar mass growth under conditions of intermittent UV feedback. This trend becomes more prominent as we improve the resolution of our simulations. We argue that simulations with significantly higher resolution than reported previously are needed to derive accurate gas mass accretion rates onto primordial protostars.

preprint2016arXiv

From the first stars to the first black holes

The growth of the first super massive black holes (SMBHs) at z > 6 is still a major challenge for theoretical models. If it starts from black hole (BH) remnants of Population III stars (light seeds with mass ~ 100 Msun) it requires super-Eddington accretion. An alternative route is to start from heavy seeds formed by the direct collapse of gas onto a ~ 10^5 Msun BH. Here we investigate the relative role of light and heavy seeds as BH progenitors of the first SMBHs. We use the cosmological, data constrained semi-analytic model GAMETE/QSOdust to simulate several independent merger histories of z > 6 quasars. Using physically motivated prescriptions to form light and heavy seeds in the progenitor galaxies, we find that the formation of a few heavy seeds (between 3 and 30 in our reference model) enables the Eddington-limited growth of SMBHs at z > 6. This conclusion depends sensitively on the interplay between chemical, radiative and mechanical feedback effects, which easily erase the conditions that allow the suppression of gas cooling in the low metallicity gas (Z < Zcr and JLW > Jcr). We find that heavy seeds can not form if dust cooling triggers gas fragmentation above a critical dust-to-gas mass ratio (D > Dcr). In addition, the relative importance of light and heavy seeds depends on the adopted mass range for light seeds, as this dramatically affects the history of cold gas along the merger tree, by both SN and AGN-driven winds.

preprint2016arXiv

The Final Fates of Accreting Supermassive Stars

The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have less concentrated structure than fully-convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit $\gtrsim 10^5~M_\odot$ derived for the fully-convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with $\dot M \lesssim 0.1~M_\odot~{\rm yr}^{-1}$. With $\dot{M} \simeq 0.3 - 1~M_\odot~{\rm yr}^{-1}$, the star becomes GR-unstable during the helium-burning stage at $M \simeq 2 - 3.5~\times 10^5~M_\odot$. In an extreme case with $10~M_\odot~{\rm yr}^{-1}$, the star does not collapse until the mass reaches $\simeq 8.0\times 10^5~M_\odot$, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.

preprint2015arXiv

Primordial Star Formation under the Influence of Far Ultraviolet Radiation: 1540 Cosmological Halos and the Stellar Mass Distribution

We perform a large set of cosmological simulations of early structure formation and follow the formation and evolution of 1540 star-forming gas clouds to derive the mass distribution of primordial stars. The star formation in our cosmological simulations is characterized by two distinct populations, the so-called Population III.1 stars and primordial stars formed under the influence of far ultraviolet (FUV) radiation (Population III.2D stars). In this work, we determine the stellar masses by using the dependences on the physical properties of star-forming cloud and/or the external photodissociating intensity from nearby primordial stars, which are derived from the results of two-dimensional radiation hydrodynamic simulations of protostellar feedback. The characteristic mass of the Pop III stars is found to be a few hundred solar masses at z ~ 25, and it gradually shifts to lower masses with decreasing redshift. At high redshifts z > 20, about half of the star-forming gas clouds are exposed to intense FUV radiation and thus give birth to massive Pop III.2D stars. However, the local FUV radiation by nearby Pop III stars becomes weaker at lower redshifts, when typical Pop III stars have smaller masses and the mean physical separation between the stars becomes large owing to cosmic expansion. Therefore, at z < 20, a large fraction of the primordial gas clouds host Pop III.1 stars. At z =< 15, the Pop III.1 stars are formed in relatively cool gas clouds due to efficient radiative cooling by H_2 and HD molecules; such stars have masses of a few x 10 Msun. Since the stellar evolution and the final fate are determined by the stellar mass, Pop III stars formed at different epochs play different roles in the early universe.

preprint2015arXiv

Supermassive star formation via episodic accretion: protostellar disc instability and radiative feedback efficiency

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

preprint2014arXiv

Conditions for HD Cooling in the First Galaxies Revisited: Interplay between Far-Ultraviolet and Cosmic Ray Feedback in Population III Star Formation

HD dominates the cooling of primordial clouds with enhanced ionization, e.g. shock-heated clouds in structure formation or supernova remnants, relic HII regions of Pop III stars, and clouds with cosmic-ray (CR) irradiation. There, the temperature decreases to several 10 K and the characteristic stellar mass decreases to $\sim 10\ {\rm M}_{\odot}$, in contrast with first stars formed from undisturbed pristine clouds ($\sim 100\ {\rm M}_{\odot}$). However, without CR irradiation, even weak far ultra-violet (FUV) irradiation suppresses HD formation/cooling. Here, we examine conditions for HD cooling in primordial clouds including both FUV and CR feedback. At the beginning of collapse, the shock-compressed gas cools with its density increasing, while the relic HII region gas cools at a constant density. Moreover, shocks tend to occur in denser environments than HII regions. Owing to the higher column density and the more effective shielding, the critical FUV intensity for HD cooling in a shock-compressed gas becomes $\sim 10$ times higher than in relic HII regions. Consequently, in the shock-compressed gas, the critical FUV intensity exceeds the background level for most of the redshift we consider ($6 \lesssim z \lesssim 15$), while in relic HII regions, HD cooling becomes effective after the CR intensity increases enough at $z \lesssim 10$. Our result suggests that less massive ($\sim 10\ {\rm M}_{\odot}$) Pop III stars may be more common than previously considered and could be the dominant population of Pop III stars.

preprint2014arXiv

Dissipation of magnetic fields in star-forming clouds with different metallicities

We study dissipation process of magnetic fields in the metallicity range $0-1 Z_{\odot}$ for contracting prestellar cloud cores. By solving non-equilibrium chemistry for important charged species including charged grains, we evaluate the drift velocity of the magnetic-field lines with respect to the gas. We find that the magnetic flux dissipates in the density range $10^{12}{\rm cm^{-3}} \lesssim n_{\rm H} \lesssim 10^{17}{\rm cm^{-3}}$ for the solar-metallicity case at the scale of the core, which is assumed to be the Jeans scale. The dissipation density range becomes narrower for lower metallicity. The magnetic field is always frozen to the gas below metallicity $\lesssim 10^{-7}-10^{-6}Z_\odot$, depending on the ionization rate by cosmic rays and/or radioactivity. With the same metallicity, the dissipation density range becomes wider for lower ionization rate. The presence of such a dissipative regime is expected to cause various dynamical phenomena in protostellar evolution such as the suppression of jet/outflow launching and fragmentation of the circumstellar disks depending on the metallicity.

preprint2014arXiv

Dust grain growth and the formation of the extremely primitive star SDSS J102915+172927

Dust grains in low-metallicity star-forming regions may be responsible for the formation of the first low-mass stars. The minimal conditions to activate dust-induced fragmentation require the gas to be pre-enriched above a critical dust-to-gas mass ratio Dcr=[2.6--6.3]x10^-9 with the spread reflecting the dependence on the grain properties. The recently discovered Galactic halo star SDSS J102915+172927 has a stellar mass of 0.8 Msun and a metallicity of Z=4.5x10^-5 Zsun and represents an optimal candidate for the dust-induced low-mass star formation. Indeed, for the two most plausible Population III supernova progenitors, with 20 Msun and 35 Msun, the critical dust-to-gas mass ratio can be overcome provided that at least 0.4 Msun of dust condenses in the ejecta, allowing for moderate destruction by the reverse shock. Here we show that even if dust formation in the first supernovae is less efficient or strong dust destruction does occur, grain growth during the collapse of the parent gas cloud is sufficiently rapid to activate dust cooling and likely fragmentation into low-mass and long-lived stars. Silicates and magnetite grains can experience significant grain growth in the density range 10^9 /cc < nH<10^12 /cc by accreting gas-phase species (SiO, SiO2, and Fe) until their gas-phase abundance drops to zero, reaching condensation efficiencies =1. The corresponding increase in the dust-to-gas mass ratio allows dust-induced cooling and fragmentation to be activated at 10^12 /cc < nH < 10^14 /cc, before the collapsing cloud becomes optically thick to continuum radiation. We show that for all the initial conditions that apply to the parent cloud of SDSS J102915+172927, dust-driven fragmentation is able to account for the formation of the star.

preprint2014arXiv

Formation of an embryonic supermassive star in the first galaxy

We studied the gravitational collapse of a warm (~8000 K) primordial-gas cloud as a candidate progenitor for a supermassive star (SMS; >10^5 Msun) using a three-dimensional hydrodynamical simulation, including all the relevant cooling processes of both H_2 and H, which can potentially induce cloud fragmentation. This is the first simulation of this kind to resolve protostar formation. We find that the cloud undergoes runaway collapse without a major episode of fragmentation. Although the H_2 fraction jumps by a large factor via the three-body reaction at ~10^-13 g/cm^3, its cooling remains inefficient due to the optical thickness, and the temperature remains >3000 K. When the central core of the cloud becomes opaque to continuum radiation at ~10^-8 g/cm^3, a hydrostatic protostar with ~0.2 Msun is formed. The protostar grows to the mass ~1 Msun and the radius ~2 AU within ~1 yr via rapid accretion of dense filamentary flows. With high accretion rate ~2 Msun/yr, the protostar is expected to turn into a SMS within its lifetime, eventually collapsing to a seed for the supermassive black hole observed in the early Universe at z~7.

preprint2014arXiv

Gravitational instability in protostellar disks at low metallicities

Fragmentation of protostellar disks controls the growth of protostars and plays a key role in determining the final mass of newborn stars. In this paper, we investigate the structure and gravitational stability of the protostellar disks in the full metallicity range between zero and the solar value. Using the mass-accretion rates evaluated from the thermal evolution in the preceding collapse phase of the pre-stellar cores, we calculate disk structures and their evolution in the framework of the standard steady disks. Overall, with higher metallicity, more efficient cooling results in the lower accretion rate and lower temperature inside the disk: at zero metallicity, the accretion rate is ~ 1e-3Msun/yr and the disk temperature is ~ 1000 K, while at solar metallicity, ~ 1e-6Msun/yr and 10 K. Despite the large difference in these values, the zero- and solar-metallicity disks have similar stability properties: the Toomre parameter for the gravitational stability, which can be written using the ratio of temperatures in the disk and in the envelope as Q ~ (T_disk/T_env)^3/2, is > 1, i.e., marginally stable. At intermediate metallicities of 1e-5--1e-3Zsun, however, the disks are found to be strongly unstable with Q ~ 0.1--1 since dust cooling, which is effective only in the disks due to their high density (> 1e10 cm^-3), makes the temperature in the disks lower than that in the envelopes. This indicates that masses of the individual stars formed as a result of the protostellar disk fragmentation can be significantly smaller than their parent core in this metallicity range. The typical stellar mass in this case would be a few Msun, which is consistent with the observationally suggested mass-scale of extremely metal-poor stars.

preprint2014arXiv

One Hundred First Stars : Protostellar Evolution and the Final Masses

We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in each gas cloud and the subsequent protostellar evolution until the gas mass accretion onto the protostar is halted by stellar radiative feedback. The accretion rates differ significantly among the primordial gas clouds which largely determine the final stellar masses. For low accretion rates the growth of a protostar is self-regulated by radiative feedback effects and the final mass is limited to several tens of solar masses. At high accretion rates the protostar's outer envelope continues to expand and the effective surface temperature remains low; such protostars do not exert strong radiative feedback and can grow in excess to one hundred solar masses. The obtained wide mass range suggests that the first stars play a variety of roles in the early universe, by triggering both core-collapse supernovae and pair-instability supernovae as well as by leaving stellar mass black holes. We find certain correlations between the final stellar mass and the physical properties of the star-forming cloud. These correlations can be used to estimate the mass of the first star from the properties of the parent cloud or of the host halo, without following the detailed protostellar evolution.

preprint2014arXiv

Supernova dust formation and the grain growth in the early universe: The critical metallicity for low-mass star formation

We investigate the condition for the formation of low-mass second-generation stars in the early universe. It has been proposed that gas cooling by dust thermal emission can trigger fragmentation of a low-metallicity star-forming gas cloud. In order to determine the critical condition in which dust cooling induces the formation of low-mass stars, we follow the thermal evolution of a collapsing cloud by a one-zone semi-analytic collapse model. Earlier studies assume the dust amount in the local universe, where all refractory elements are depleted onto grains, and/or assume the constant dust amount during gas collapse. In this paper, we employ the models of dust formation and destruction in early supernovae to derive the realistic dust compositions and size distributions for multiple species as the initial conditions of our collapse calculations. We also follow accretion of heavy elements in the gas phase onto dust grains, i.e., grain growth, during gas contraction. We find that grain growth well alters the fragmentation property of the clouds, and that this still does not approach to the value in the local universe. The critical conditions can be written by the gas metallicity Zcr and the initial depletion efficiency fdep,0 of gas-phase metal onto grains, or dust-to-metal mass ratio, as (Zcr/10^{-5.5} Zsun) = (fdep,0/0.18)^{-0.44} with small scatters in the range of Zcr = [0.06--3.2]x10^{-5} Zsun. We also show that the initial dust composition and size distribution are important to determine Zcr.

preprint2014arXiv

The critical radiation intensity for direct collapse black hole formation: dependence on the radiation spectral shape

It has been proposed that supermassive black holes (SMBHs) are originated from direct-collapse black holes (DCBHs) that are formed at z gtrsim 10 in the primordial gas in the case that H2 cooling is suppressed by strong external radiation. In this work, we study the critical specific intensity J^crit required for DCBH formation for various radiation spectral shapes by a series of one-zone calculations of a collapsing primordial- gas cloud. We calculate the critical specific intensity at the Lyman-Werner (LW) bands J^crit_LW,21 (in units of 10^-21 erg s^-1 Hz^-1 sr^-1 cm^-2) for realistic spectra of metal-poor galaxies. We find J^crit is not sensitive to the age or metallicity for the constant star formation galaxies with J^crit_LW,21 = 1300-1400, while J^crit decreases as galaxies become older or more metal-enriched for the instantaneous starburst galaxies. However, such dependence for the instantaneous starburst galaxies is weak for the young or extremely metal-poor galaxies: J^crit_LW,21 = 1000-1400 for the young galaxies and J^crit_LW,21 approx 1400 for the extremely metal-poor galaxies. The typical value of J^crit for the realistic spectra is higher than those expected in the literature, which affects the estimated DCBH number density n_DCBH. By extrapolating the result of Dijkstra, Ferrara and Mesinger, we obtain n_DCBH sim 10^-10 cMpc^-3 at z = 10, although there is still large uncertainty in this estimation.

preprint2014arXiv

The origin of the most iron-poor star

We investigate the origin of carbon-enhanced metal-poor (CEMP) stars starting from the recently discovered $\rm [Fe/H]<-7.1$ star SMSS J031300 (Keller et al. 2014). We show that the elemental abundances observed on the surface of SMSS J031300 can be well fit by the yields of faint, metal free, supernovae. Using properly calibrated faint supernova explosion models, we study, for the first time, the formation of dust grains in such carbon-rich, iron-poor supernova ejecta. Calculations are performed assuming both unmixed and uniformly mixed ejecta and taking into account the partial destruction by the supernova reverse shock. We find that, due to the paucity of refractory elements beside carbon, amorphous carbon is the only grain species to form, with carbon condensation efficiencies that range between (0.15-0.84), resulting in dust yields in the range (0.025-2.25)M$_{\odot}$. We follow the collapse and fragmentation of a star forming cloud enriched by the products of these faint supernova explosions and we explore the role played by fine structure line cooling and dust cooling. We show that even if grain growth during the collapse has a minor effect of the dust-to-gas ratio, due to C depletion into CO molecules at an early stage of the collapse, the formation of CEMP low-mass stars, such as SMSS J031300, could be triggered by dust cooling and fragmentation. A comparison between model predictions and observations of a sample of C-normal and C-rich metal-poor stars supports the idea that a single common pathway may be responsible for the formation of the first low-mass stars.

preprint2014arXiv

Thermal instability and multi-phase interstellar medium in the first galaxies

We examine the linear stability and nonlinear growth of the thermal instability in isobarically contracting gas with various metallicities and FUV field strengths. When the H2 cooling is suppressed by FUV fields (G_0>10^-3) or the metallicity is high enough (Z/Zs>10^-3), the interstellar medium is thermally unstable in the temperature range 100-7000 K owing to the cooling by CII and OI fine-structure lines. In this case, a bi-phasic medium with a bimodal density probability distribution is formed as a consequence of the thermal instability. The characteristic scales of the thermal instability become smaller with increasing metallicity. Comparisons of the nonlinear simulations with different resolution indicates that the maximum scale of the thermal instability should be resolved with more than 60 cells to follow runaway cooling driven by the thermal instability. Under sufficiently weak FUV fields and with low metallcity, the density range of the thermal instability shrinks owing to dominance of H2 cooling. As the FUV intensity is reduced, bi-phasic structure becomes less remarkable and disappears eventually. Our basic results suggest that, in early galaxies, i) the thermal instability has little effect for the medium with Z/Zs<10^-4, ii) fragmentation by the thermal instability could determine mass spectrum of star clusters for 10^-4<Z/Zs<0.04, and iii) thermally bistable turbulent interstellar medium like our galaxy becomes ubiquitous for Z/Zs>0.04, although the threshold metallicity depends on the conditions such as thermal pressure, FUV strength and redshift.

preprint2013arXiv

Formation of Primordial Supermassive Stars by Rapid Mass Accretion

Supermassive stars (SMSs) forming via very rapid mass accretion (Mdot >~ 0.1 Msun/yr) could be precursors of supermassive black holes observed beyond redshift of about 6. Extending our previous work, we here study the evolution of primordial stars growing under such rapid mass accretion until the stellar mass reaches 10^{4 - 5} Msun. Our stellar evolution calculations show that a star becomes supermassive while passing through the "supergiant protostar" stage, whereby the star has a very bloated envelope and a contracting inner core. The stellar radius increases monotonically with the stellar mass, until =~ 100 AU for M_* >~ 10^4 Msun, after which the star begins to slowly contract. Because of the large radius the effective temperature is always less than 10^4 K during rapid accretion. The accreting material is thus almost completely transparent to the stellar radiation. Only for M_* >~ 10^5 Msun can stellar UV feedback operate and disturb the mass accretion flow. We also examine the pulsation stability of accreting SMSs, showing that the pulsation-driven mass loss does not prevent stellar mass growth. Observational signatures of bloated SMSs should be detectable with future observational facilities such as the James Webb Space Telescope. Our results predict that an inner core of the accreting SMS should suffer from the general relativistic instability soon after the stellar mass exceeds 10^5 Msun. An extremely massive black hole should form after the collapse of the inner core.

preprint2013arXiv

Photoevaporation of Circumstellar Disks Revisited: The Dust-Free Case

Photoevaporation by stellar ionizing radiation is believed to play an important role in the dispersal of disks around young stars. The mass loss model for dust-free disks developed by Hollenbach et al. is currently regarded as a conventional one and has been used in a wide variety of studies. However, the rate in this model was derived by the crude so-called 1+1D approximation of ionizing radiation transfer, which assumes that diffuse radiation propagates in a direction vertical to the disk. In this study, we revisit the photoevaporation of dust-free disks by solving the 2D axisymmetric radiative transfer for steady-state disks. Unlike that solved by the conventional model, we determine that direct stellar radiation is more important than the diffuse field at the disk surface. The radial density distribution at the ionization boundary is represented by the single power-law with an index -3/2 in contrast to the conventional double power-law. For this distribution, the photoevaporation rate from the entire disk can be written as a function of the ionizing photon emissivity, Phi_EUV, from the central star and the disk outer radius, r_d, as follows: Mdot_PE = 5.4 x 10^-5 x (Phi_EUV/10^49 sec^-1)^1/2 x (r_d/1000 AU)^1/2 Msun/yr. This new rate depends on the outer disk radius rather than on the gravitational radius as in the conventional model, caused by the enhanced contribution to the mass loss from the outer disk annuli. In addition, we discuss its applications to present-day as well as primordial star formation.

preprint2013arXiv

Pulsational instability of supergiant protostars: Do they grow supermassive by accretion?

Supermassive stars (SMSs; M>10^5 Msun) and their remnant black holes are promising progenitors for supermassive black holes (SMBHs) observed in the early universe at z>7. It has been postulated that SMSs forms through very rapid mass accretion onto a protostar at a high rate exceeding 0.01 Msun/yr. According to recent studies, such rapidly accreting protostars evolve into "supergiant protostars", i.e. protostars consisting of a bloated envelope and a contracting core, similar to giant star. However, like massive stars as well as giant stars, both of which are known to be pulsationally unstable, supergiant protostars may also be also unstable to launch strong pulsation-driven outflows. If this is the case, the stellar growth via accretion will be hindered by the mass loss. We here study the pulsational stability of the supergiant protostars in the mass range M<10^3 Msun through the method of the linear perturbation analysis. We find that the supergiant protostars with M>600 Msun and very high accretion rate Mdot>1.0 Msun/yr are unstable due to the kappa mechanism. The pulsation is excited in the He^+ ionization layer in the envelope. Even under a conservative assumption that all the pulsation energy is converted into the kinetic energy of the outflows, the mass-loss rate is ~10^-3 Msun/yr, which is lower than the accretion rate by more than two orders of magnitude. We thus conclude that the supergiant protostars should grow stably via rapid accretion at least in the mass range we studied. As long as the rapid accretion is maintained in the later stage, protostars will become SMSs, which eventually produce seeds for the high-z SMBHs.

preprint2012arXiv

Do the environmental conditions affect the dust-induced fragmentation in low-metallicity clouds ?: Effect of pre-ionization and far-ultraviolet/cosmic-ray fields

We study effects of the fully ionized initial state, or pre-ionization, on the subsequent thermal evolution of low-metallicity clouds under various intensities of the external far-ultraviolet(FUV) and cosmic-ray(CR) fields. The pre-ionization significantly affects the thermal and dynamical evolution of metal-free clouds without FUV/CRs by way of efficient HD formation. On the other hand, the pre-ionization effect on the thermal evolution is limited in very low-density regime for more metal-enriched clouds ([Z/H] >~ -4) or those under modest FUV (>10^{-3}) or CR field (>0.1 of the present-day Galactic disk levels). In any case, for >10^8 cm^{-3}, neither the initial ionization state nor the irradiating FUV strength affect the thermal evolution. The dust cooling is an important mechanism for making sub-solar mass fragments in low-metallicity gas. Since this fragmentation occurs at the temperature minimum by the dust cooling at >10^{10} cm^{-3}, this process is not vulnerable either to initial ionization state or external radiation.

preprint2012arXiv

Protostellar Feedback and Final Mass of the Second-Generation Primordial Stars

The first stars in the universe ionized the ambient primordial gas through various feedback processes. "Second-generation" primordial stars potentially form from this disturbed gas after its recombination. In this Letter, we study the late formation stage of such second-generation stars, where a large amount of gas accretes onto the protostar and the final stellar mass is determined when the accretion terminates. We directly compute the complex interplay between the accretion flow and stellar ultraviolet (UV) radiation, performing radiation-hydrodynamic simulations coupled with stellar evolution calculations. Because of more efficient H2 and HD cooling in the pre-stellar stage, the accretion rates onto the star are ten times lower than in the case of the formation of the first stars. The lower accretion rates and envelope density result in the occurrence of an expanding bipolar HII region at a lower protostellar mass M_* \simeq 10Msun, which blows out the circumstellar material, thereby quenching the mass supply from the envelope to the accretion disk. At the same time the disk loses mass due to photoevaporation by the growing star. In our fiducial case the stellar UV feedback terminates mass accretion onto the star at M_* \simeq 17Msun. Although the derived masses of the second-generation primordial stars are systematically lower than those of the first generation, the difference is within a factor of only a few. Our results suggest a new scenario, whereby the majority of the primordial stars are born as massive stars with tens of solar masses, regardless of their generations.

preprint2012arXiv

Rapidly Accreting Supergiant Protostars: Embryos of Supermassive Black Holes?

Direct collapse of supermassive stars (SMSs) is a possible pathway for generating supermassive black holes in the early universe. It is expected that an SMS could form via very rapid mass accretion with Mdot ~ 0.1 - 1 Msun/yr during the gravitational collapse of an atomic-cooling primordial gas cloud. In this paper we study how stars would evolve under such extreme rapid mass accretion, focusing on the early evolution until the stellar mass reaches 1000 Msun. To this end we numerically calculate the detailed interior structure of accreting stars with primordial element abundances. Our results show that for accretion rates higher than 0.01 Msun/yr, stellar evolution is qualitatively different from that expected at lower rates. While accreting at these high rates the star always has a radius exceeding 100 Rsun, which increases monotonically with the stellar mass. The mass-radius relation for stellar masses exceeding ~ 100 Msun follows the same track with R_* \propto M_*^0.5 in all cases with accretion rates > 0.01 Msun/yr; at a stellar mass of 1000 Msun the radius is about 7000 Rsun (~= 30 AU). With higher accretion rates the onset of hydrogen burning is shifted towards higher stellar masses. In particular, for accretion rates exceeding Mdot > 0.1 Msun/yr, there is no significant hydrogen burning even after 1000 Msun have accreted onto the protostar. Such "supergiant" protostars have effective temperatures as low as Teff ~= 5000 K throughout their evolution and because they hardly emit ionizing photons, they do not create an HII region or significantly heat their immediate surroundings. Thus, radiative feedback is unable to hinder the growth of rapidly accreting stars to masses in excess of 1000 Msun, as long as material is accreted at rates Mdot > 0.01 Msun/yr.

preprint2012arXiv

Supermassive black hole formation by the cold accretion shocks in the first galaxies

We propose a new scenario for supermassive star (SMS;>10^5Msun) formation in shocked regions of colliding cold accretion flows near the centers of first galaxies. Recent numerical simulations indicate that assembly of a typical first galaxy with virial temperature (~10^4K) proceeds via cold and dense flows penetrating deep to the center, where the supersonic streams collide each other to develop a hot and dense (~10^4K, ~10^3/cc) shocked gas. The post-shock layer first cools by efficient Ly alpha emission and contracts isobarically until 8000K. Whether the layer continues the isobaric contraction depends on the density at this moment: if the density is high enough for collisionally exciting H2 rovibrational levels (>10^4/cc), enhanced H2 collisional dissociation suppresses the gas to cool further. In this case, the layer fragments into massive (>10^5Msun) clouds, which collapse isothermally (~8000K) by the Ly alpha cooling without subsequent fragmentation. As an outcome, SMSs are expected to form and evolve eventually to seeds of supermassive black holes (SMBH). By calculating thermal evolution of the post-shock gas, we delimit the range of post-shock conditions for the SMS formation, which can be expressed as: T>6000K/(n/10^4/cc) for n<10^4/cc and T>5000-6000K for n>10^4/cc, depending somewhat on initial ionization degree. We found that metal enrichment does not affect the above condition for metallicity below 10^-3Zsun if metals are in the gas phase, while condensation of several percent of metals into dust decreases this critical value of metallicity by an order of magnitude. Unlike the previously proposed scenario for SMS formation, which postulates extremely strong ultraviolet radiation to quench H2 cooling, our scenario naturally explains the SMBH seed formation in the assembly process of the first galaxies, even without such a strong radiation.

preprint2012arXiv

The formation of the extremely primitive star SDSS J102915+172927 relies on dust

The relative importance of metals and dust grains in the formation of the first low-mass stars has been a subject of debate. The recently discovered Galactic halo star SDSS J102915+172927 (Caffau et al. 2011) has a mass less than 0.8 Msun and a metallicity of Z = 4.5 10^{-5} Zsun. We investigate the origin and properties of this star by reconstructing the physical conditions in its birth cloud. We show that the observed elemental abundance trend of SDSS J102915+172927 can be well fitted by the yields of core-collapse supernovae with metal-free progenitors of 20 Msun and 35 Msun. Using these selected supernova explosion models, we compute the corresponding dust yields and the resulting dust depletion factor taking into account the partial destruction by the supernova reverse shock. We then follow the collapse and fragmentation of a star forming cloud enriched by the products of these SN explosions at the observed metallicity of SDSS J102915+172927. We find that [0.05 - 0.1] Msun mass fragments, which then lead to the formation of low-mass stars, can occur provided that the mass fraction of dust grains in the birth cloud exceeds 0.01 of the total mass of metals and dust. This, in turn, requires that at least 0.4 Msun of dust condense in the first supernovae, allowing for moderate destruction by the reverse shock. If dust formation in the first supernovae is less efficient or strong dust destruction does occur, the thermal evolution of the SDSS J102915+172927 birth cloud is dominated by molecular cooling, and only > 8 Msun fragments can form. We conclude that the observed properties of SDSS J102915+172927 support the suggestion that dust must have condensed in the ejecta of the first supernovae and played a fundamental role in the formation of the first low-mass stars.

preprint2012arXiv

Variable Accretion Rates and Fluffy First Stars

We combine the output of hydrodynamical simulations of Population III star cluster formation with stellar evolution models, and calculate the evolution of protostars experiencing variable mass accretion rates due to interactions within a massive disk. We find that the primordial protostars are extended 'fluffy' objects for the bulk of their pre-main-sequence lifetimes. Accretion luminosity feedback from such objects is high, but as shown in previous work, has a minimal effect on the star cluster. The extended radii of the protostars, combined with the observation of close encounters in the simulations, suggests that mergers will occur in such systems. Furthermore, mass transfer between close protostellar binaries with extended radii could lead to massive tight binaries, which are a possible progenitor of gamma ray bursts.

preprint2011arXiv

Effect of cosmic ray/X-ray ionization on supermassive black hole formation

We study effects of external ionization by cosmic rays (CRs) and X-rays on the thermal evolution of primordial clouds under strong far-ultraviolet (FUV) radiation. A strong FUV radiation dissociates H2 and quenches its cooling. Even in such an environment, a massive cloud with Tvir>10^4 K can contract isothermally at 8000 K by Lyman alpha cooling. This cloud collapses monolithically without fragmentation, and a supermassive star (>10^5 Msun) is believed to form at the center, which eventually evolves to a supermassive black hole (SMBH). However, candidates of FUV sources, including star-forming galaxies, are probably sources of strong CRs and X-rays, as well. We find that the external ionization promotes H2 production and elevates the threshold FUV intensity Jcr needed for the SMBH formation for CR energy density U_CR>10^-14 erg/cm^3 or X-ray intensity J_X>10^-24 erg/s/cm^2/sr/Hz at 1 keV. The critical FUV flux increases in proportion to U_CR^{1/2} (J_X^{1/2}) in the high CR (X-ray, respectively) limit. With the same value of FUV intensity at the Lyman limit (13.6 eV), the H^- photodissociation rate, with threshold of 0.755 eV, increases and the H2 molecules decrease with decreasing effective temperature of the FUV sources T*. The lower value of T* thus results in the lower critical FUV flux at 13.6 eV. Using an empirical relation between intensities of FUV and CR/X-ray from nearby star-forming galaxies, we find that external ionization effect remarkably enhances the critical FUV flux for sources with T* as high as 10^5 K and composed of stars with <100 Msun to a level that is not realized in any halo. This indicates that, to induce SMBH formation, the FUV sources must be either Pop II/I galaxies with low brightness temperature, Pop III galaxies with a very top-heavy IMF, or Pop III galaxies too young to harbor sources of CR/X-ray, e.g., supernova remnants or high-mass X-ray binaries.

preprint2011arXiv

Protostellar Feedback Halts the Growth of the First Stars in the Universe

The first stars fundamentally transformed the early universe by emitting the first light and by producing the first heavy elements. These effects were predetermined by the mass distribution of the first stars, which is thought to have been fixed by a complex interplay of gas accretion and protostellar radiation. We performed radiation-hydrodynamics simulations that followed the growth of a primordial protostar through to the early stages as a star with thermo-nuclear burning. The circumstellar accretion disk was evaporated by ultraviolet radiation from the star when its mass was 43 times that of the Sun. Such massive primordial stars, in contrast to the often postulated extremely massive stars, may help explain the fact that there are no signatures of the pair-instability supernovae in abundance patterns of metal-poor stars in our galaxy.

preprint2011arXiv

The first low-mass stars: critical metallicity or dust-to-gas ratio?

We explore the minimal conditions which enable the formation of metal-enriched solar and sub-solar mass stars. We find that in the absence of dust grains, gas fragmentation occurs at densities nH ~ [10^4-10^5]cm^{-3} when the metallicity exceeds Z ~ 10^{-4} Zsun. The resulting fragmentation masses are > 10 Msun. The inclusion of Fe and Si cooling does not affect the thermal evolution as this is dominated by molecular cooling even for metallicities as large as Z = 10^{-2} Zsun. The presence of dust is the key driver for the formation of low-mass stars. We focus on three representative core-collapse supernova (SN) progenitors, and consider the effects of reverse shocks of increasing strength: these reduce the depletion factors, fdep = Mdust/(Mdust+Mmet), alter the shape of the grain size distribution function and modify the relative abundances of grain species and of metal species in the gas phase. We find that the lowest metallicity at which fragmentation occurs is Z=10^{-6} Zsun for gas pre-enriched by the explosion of a 20 Msun primordial SN (fdep > 0.22) and/or by a 35 Msun, Z=10^{-4} Zsun SN (fdep > 0.26); it is ~ 1 dex larger, when the gas is pre-enriched by a Z = 10^{-4} Zsun, 20 Msun SN (fdep > 0.04). Cloud fragmentation depends on the depletion factor and it is suppressed when the reverse shock leads to a too large destruction of dust grains. These features are all consistent with the existence of a minimum dust-to-gas ratio, Dcr, above which fragmentation is activated. We derive a simple analytic expression for Dcr which, for grain composition and properties explored in the present study, reads Dcr = [2.6 - 6.3] x 10^{-9}. When the dust-to-gas ratio of star forming clouds exceeds this value, the fragmentation masses range between 0.01 Msun and 1 Msun, thus enabling the formation of the first low-mass stars.

preprint2010arXiv

Evolution of Massive Protostars via Disk Accretion

Mass accretion onto (proto-)stars at high accretion rates > 10^-4 M_sun/yr is expected in massive star formation. We study the evolution of massive protostars at such high rates by numerically solving the stellar structure equations. In this paper we examine the evolution via disk accretion. We consider a limiting case of "cold" disk accretion, whereby most of the stellar photosphere can radiate freely with negligible backwarming from the accretion flow, and the accreting material settles onto the star with the same specific entropy as the photosphere. We compare our results to the calculated evolution via spherically symmetric accretion, the opposite limit, whereby the material accreting onto the star contains the entropy produced in the accretion shock front. We examine how different accretion geometries affect the evolution of massive protostars. For cold disk accretion at 10^-3 M_sun/yr the radius of a protostar is initially small, about a few R_sun. After several solar masses have accreted, the protostar begins to bloat up and for M \simeq 10 M_sun the stellar radius attains its maximum of 30 - 400 R_sun. The large radius about 100 R_sun is also a feature of spherically symmetric accretion at the same accreted mass and accretion rate. Hence, expansion to a large radius is a robust feature of accreting massive protostars. At later times the protostar eventually begins to contract and reaches the Zero-Age Main-Sequence (ZAMS) for M \simeq 30 M_sun, independent of the accretion geometry. For accretion rates exceeding several 10^-3 M_sun/yr the protostar never contracts to the ZAMS. The very large radius of several 100s R_sun results in a low effective temperature and low UV luminosity of the protostar. Such bloated protostars could well explain the existence of bright high-mass protostellar objects, which lack detectable HII regions.

preprint2010arXiv

Low-Metallicity Star Formation : Prestellar Collapse and Protostellar Accretion in the Spherical Symmetry

The collapse of dense cores with different metallicities is studied by hydrodynamical calculations coupled with detailed chemical and radiative processes. For this purpose, we construct a simple chemical network with non-equilibrium reactions among 15 chemical species, which reproduces the abundance of important molecular coolants by more detailed network very well. The evolution is followed until the formation of a hydrostatic protostar at the center. In a lower-metallicity gas cloud, the temperature during the collapse remains high owing to less efficient cooling. Using the temperature evolution at the center as a function the density, we discuss the possibility of fragmentation during the dust-cooling phase. The critical metallicity for the fragmentation is 10^{-5}Z_sun assuming moderate elongation of the cloud cores at the onset of this phase. From the density and velocity distributions at the time of protostar formation, we evaluate the mass accretion rate in the subsequent accretion phase. Using these accretion rates, we also calculate the evolution of the protostars under the assumption of stationary accretion flow. Finally, we discuss possible suppression of fragmentation by heating of the ambient gas by protostellar radiation, which is considered important in the contemporary star formation. We argue that it is negligible for <10^{-2}Zsun.

preprint2009arXiv

Binary Formation with Different Metallicities: Dependence on Initial Conditions

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower-metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for < 10^-4 Z_sun. With initial mass < 0.1 M_sun, if fragments are ejected in embryo from the host clouds by multi-body interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.

preprint2009arXiv

Dust coagulation in star formation with different metallicities

Dust grains coagulate into larger aggregates in dense gas. This changes their size distribution and possibly affects the thermal evolution of star-forming clouds. We here investigate dust coagulation in collapsing pre-stellar cores with different metallicities by considering the thermal motions of grains. We show that coagulation does occur even at low metallicity $\sim 10^{-6} \Zsun$. However, we also find (i) that the H$_2$ formation rate on dust grains is reduced only after the majority of H$_2$ is formed; and (ii) that the dust opacity is modified only after the core becomes optically thick. Therefore, we conclude that the effects of dust coagulation can safely be neglected in discussing the temperature evolution of the pre-stellar cores for any metallicity as long as the grain motions are thermal.

preprint2009arXiv

Low-Metallicity Protostars and the Maximum Stellar Mass Resulting from Radiative Feedback

The final mass of a newborn star is set at the epoch when the mass accretion onto the star is terminated. We study the evolution of accreting protostars and the limits of accretion in low metallicity environments. Accretion rates onto protostars are estimated via the temperature evolution of prestellar cores with different metallicities. The derived rates increase with decreasing metallicity, from 10^-6 M_sun/yr at Z = Z_sun to 10^-3 M_sun/yr at Z = 0. With the derived accretion rates, the protostellar evolution is numerically calculated. We find that, at lower metallicity, the protostar has a larger radius and reaches the zero-age main-sequence (ZAMS) at higher stellar mass. Using this protostellar evolution, we evaluate the upper stellar mass limit where the mass accretion is hindered by radiative feedback. We consider the effects of radiation pressure exerted on the accreting envelope, and expansion of the HII region. The mass accretion is finally terminated by radiation pressure on dust grains in the envelope for Z > 10^-3 Z_sun and by the expanding HII region for lower metallicty. The mass limit from these effects increases with decreasing metallicity from 10 M_sun at Z = Z_sun to about 300 M_sun at Z = 10^-6 Z_sun. The termination of accretion occurs after the central star arrives at the ZAMS at all metallicities, which allows us to neglect protostellar evolution effects in discussing the upper mass limit by stellar feedback. The fragmentation induced by line cooling in low-metallicity clouds yields prestellar cores with masses large enough that the final stellar mass is set by the feedback effects.

preprint2009arXiv

Metals, dust and the cosmic microwave background: fragmentation of high-redshift star-forming clouds

We investigate the effects of the Cosmic Microwave Background (CMB) radiation field on the collapse of prestellar clouds. Using a semi-analytic model to follow the thermal evolution of clouds with varying initial metallicities and dust contents at different redshifts, we study self-consistently the response of the mean Jeans mass at cloud fragmentation to metal line-cooling, dust-cooling and the CMB. In the absence of dust grains, at redshifts z < 10 moderate characteristic masses (of 10s of Msun) are formed when the metallicity is 10^{-4} Zsun < Z < 10^{-2.5} Zsun; at higher metallicities, the CMB inhibits fragmentation and only very large masses (of ~ 100s of Msun) are formed. These effects become even more dramatic at z > 10 and the fragmentation mass scales are always > 100s of Msun, independent of the initial metallicity. When dust grains are present, sub-solar mass fragments are formed at any redshift for metallicities Z > 10^{-6} Zsun because dust-cooling remains relatively insensitive to the presence of the CMB. When Z > 10^{-3} Zsun, heating of dust grains by the CMB at z > 5 favors the formation of larger masses, which become super-solar when Z > 10^{-2} Zsun and z > 10. Finally, we discuss the implications of our result for the interpretation of the observed abundance patterns of very metal-poor stars in the galactic halo.

preprint2009arXiv

Star Formation in Relic HII Regions of the First Stars: Binarity and Outflow Driving

Star formation in relic HII regions of the first stars is investigated using magneto-hydrodynamical simulations with a nested grid method that covers 10 orders of magnitude in spatial scale and 20 orders of magnitude in density contrast. Due to larger fraction of H_2 and HD molecules, its prestellar thermal evolution is considerably different from that in the first star formation. Reflecting the difference, two hydrostatic cores appear in a nested manner: a protostar is enclosed by a transient hydrostatic core, which appears during the prestellar collapse. If the initial natal core rotates fast at a rate with rotational to gravitational energy ratio β_0 > 0.01-0.1, the transient hydrostatic core fragments to \sim 10 M_\odot sub-cores at density \sim 10^9 cm^-3. With smaller rotation energy, fragmentation occurs at higher density while a single protostar forms without fragmentation if rotation is extremely slow with β_0 < 10^-6 -10^-5. If magnetic field is present, these threshold values of β_0 is boosted owing to angular momentum transport by the magnetic breaking. Magnetic field also drives the protostellar outflows. With strong magnetic field, two distinct outflows are observed: The slower one emanates from the transient hydrostatic core, while the faster one from the protostar. These flows may affect the final stellar mass by ejecting some of masses in the initial core, and also may play some role in driving and maintenance of interstellar turbulence in young galaxies.

preprint2009arXiv

Star Formation Triggered by Supernova Explosions in Young Galaxies

We study the evolution of supernova remnants in a low-metallicity medium $Z/Z_{\odot} = 10^{-4}$ -- $10^{-2}$ in the early universe, using one-dimensional hydrodynamics with non-equilibrium chemistry. Once a post-shock layer is able to cool radiatively, a dense shell forms behind the shock. If this shell becomes gravitationally unstable and fragments into pieces, next-generation stars are expected to form from these fragments. To explore the possibility of this triggered star formation, we apply a linear perturbation analysis of an expanding shell to our results and constrain the parameter range of ambient density, explosion energy, and metallicity where fragmentation of the shell occurs. For the explosion energy of $10^{51}{\rm ergs} (10^{52}{\rm ergs})$, the shell fragmentation occurs for ambient densities higher than $\gtrsim 10^{2} {\rm cm^{-3}}$ (10 ${\rm cm^{-3}}$, respectively). This condition depends little on the metallicity in the ranges we examined. We find that the mode of star formation triggered occurs only in massive ($\gtrsim 10^{8}M_{\odot}$) haloes.

preprint2008arXiv

Evolution of Massive Protostars with High Accretion Rates

Formation of massive stars by accretion requires a high accretion rate of > 10^-4 M_sun/yr to overcome the radiation pressure barrier of the forming stars. Here, we study evolution of protostars accreting at such high rates, by solving the structure of the central star and the inner accreting envelope simultaneously. The protostellar evolution is followed starting from small initial cores until their arrival at the stage of the Zero-Age Main Sequence (ZAMS) stars. An emphasis is put on evolutionary features different from those with a low accretion rate of 10^-5 M_sun/yr, which is presumed in the standard scenario for low-mass star formation. With the high accretion rate of 10^-3 M_sun/yr, the protostellar radius becomes very large and exceeds 100 R_sun. It is not until the stellar mass reaches 40 M_sun that hydrogen burning begins and the protostar reaches the ZAMS phase, and this ZAMS arrival mass increases with the accretion rate. At a very high accretion rate of > 3 x 10^-3 M_sun/yr, the total luminosity of the protostar becomes so high that the resultant radiation pressure inhibits the growth of the protostars under steady accretion before reaching the ZAMS stage. Therefore, the evolution under the critical accretion rate 3 x 10^-3 M_sun/yr gives the upper mass limit of possible pre-main-sequence stars at 60 M_sun. The upper mass limit of MS stars is also set by the radiation pressure onto the dusty envelope under the same accretion rate at 250 M_sun. We also propose that the central source enshrouded in the Orion KL/BN nebula has effective temperature and luminosity consistent with our model, and is a possible candidate for such protostars growing under the high accretion rate. (abridged)

preprint2006arXiv

Formation of Primordial Stars in a LCDM Universe

We study the formation of the first generation of stars in the standard cold dark matter model, using a very high-resolution hydordynamic simulations. Our simulation achieves a dynamic range of 10^{10} in length scale. With accurate treatment of atomic and molecular physics, it allows us to study the chemo-thermal evolution of primordial gas clouds to densities up to n = 10^{16}/cc without assuming any a priori equation of state; a six orders of magnitudes improvement over previous three-dimensional calculations. All the relevant atomic and molecular cooling and heating processes, including cooling by collision-induced continuum emission, are implemented. For calculating optically thick H2 cooling at high densities, we use the Sobolev method. To examine possible gas fragmentation owing to thermal instability, we compute explicitly the growth rate of isobaric perturbations. We show that the cloud core does not fragment in either the low-density or high-density regimes. We also show that the core remains stable against gravitational deformation and fragmentation. We obtain an accurate gas mass accretion rate within a 10 Msun innermost region around the protostar. The protostar is accreting the surrounding hot gas at a rate of 0.001-0.01 Msun/yr. From these findings we conclude that primordial stars formed in early minihalos are massive. We carry out proto-stellar evolution calculations using the obtained accretion rate. The resulting mass of the first star is M_ZAMS = 60-100 Msun, with the exact mass dependent on the actual accretion rate.

preprint2005arXiv

Formation of Population III Stars in Fossil HII Regions: Significance of HD

We study the evolution of gas in HII regions around the first stars after the death of the exciting stars. If the first star in a small halo dies without supernova (SN), subsequent star formation is possible in the same halo. We thus investigate the effect of ionization of the gas on subsequent star formation within small halos in the early universe using one-dimensional hydrodynamics with spherical symmetry along with non-equilibrium primordial gas chemistry. We find that the enhanced electron fraction facilitates the formation of molecular hydrogen at the cores of these halos. The low temperature circumstances produced by the H_2 cooling is suitable for HD formation and the resultant cooling further drops the temperature below 100 K. Consequently, low-mass stars with primordial abundances can form even in a small halo. After accreting the interstellar metals, these stars might resemble low-mass ultra metal-poor stars discovered in the present Galactic halo.

preprint2005arXiv

Primordial Molecular Emission in Population III Galaxies

We study formation of molecules in primordial prestellar clumps and evaluate the line luminosities to assess detectability by next-generation facilities. If the initial H_2 fraction is sufficiently high, HD becomes an important coolant in the clumps. The luminosity from such HD cooling clumps is lower than that from H_2 cooling ones because of lower temperature (<100K). As for Li reactions, we include the three-body LiH formation approximately. The Li molecular fraction remains very low (<10^{-3}) throughout the evolution owing to the high dissociative reaction rate of LiH +H -> Li + H_2. LiH does not become an important coolant in any density range. The luminous emission lines from the prestellar cores include H_2 rovibrational lines: 1-0 Q(1), 1-0 O(3), 1-0 O(5), and pure rotational lines: 0-0 S(3), 0-0 S(4), 0-0 S(5). The next-generation facilities SPICA and JWST are able to detect H_2 emission in a large pre-galactic cloud that forms metal-free stars at a high rate of \sim 10^3 M_s/yr at redshift z<10. We also derive an analytical expression for the luminosity that reproduces the numerical results.

preprint2004arXiv

H_2 Line Emission Associated with the Formation of the First Stars

Molecular hydrogen line radiation emitted in formation events of first-generation stars are evaluated in a discussion of its detectability by future observational facilities. H_2 luminosity evolution from the onset of prestellar collapse until the formation of a \sim 100 M_{\odot} protostar is followed. Calculations are extended not only to the early phase of the runaway collapse but also to the later phase of accretion, whose observational features have not been studied before. Contrary to the runaway collapse phase, where the pure-rotational lines are always dominant, in the accretion phase rovibrational line emission becomes prominent. The maximum luminosity is also attained in the accretion phase for strong emission lines. The peak intensity of the strongest rovibrational line reaches \sim 10^{-29} (W/m^2), corresponding to the flux density of 10^{-5} (μJy), for a source at the typical redshift of first-generation star formation, 1+z=20. Although the redshifted rovibrational H_2 emission from such an epoch falls in the wavelength range of the next-generation infrared satellite, Space Infrared Telescope for Cosmology and Astrophysics, for exceeding the detection threshold 10^7 such protostars are required to reach the maximum luminosity simultaneously in a pregalactic cloud. It is improbable that this condition is satisfied in a realistic scenario of early structure formation.

preprint1998arXiv

Formation of Primordial Protostars

The evolution of collapsing metal free protostellar clouds is investigated for various masses and initial conditions. We perform hydrodynamical calculations for spherically symmetric clouds taking account of radiative transfer of the molecular hydrogen lines and the continuum, as well as of chemistry of the molecular hydrogen. The collapse is found to proceed almost self-similarly like Larson-Penston similarity solution. In the course of the collapse, efficient three-body processes transform atomic hydrogen in an inner region of $\sim 1 M_{\sun}$ entirely into molecular form. However, hydrogen in the outer part remains totally atomic although there is an intervening transitional layer of several solar masses, where hydrogen is in partially molecular form. No opaque transient core is formed although clouds become optically thick to H$_{2}$ collision-induced absorption continuum, since H$_{2}$ dissociation follows successively. When the central part of the cloud reaches stellar densities ($\sim 10^{-2} {\rm g cm^{-3}}$), a very small hydrostatic core ($\sim 5 \times 10^{-3} M_{\sun}$) is formed and subsequently grows in mass as the ambient gas accretes onto it. The mass accretion rate is estimated to be $3.7 \times 10^{-2} M_{\sun} {\rm yr^{-1}} (M_{\ast}/M_{\sun})^{-0.37}$, where $M_{\ast}$ is instantaneous mass of the central core, by using a similarity solution which reproduces the evolution of the cloud before the core formation.