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.

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.

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 $\Sigma$. At $Z=1~Z_{\odot}$, for instance, the SFE increases from 0.02 to 0.3 with $\Sigma$ from $10$ to $300~{\rm 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 HII 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 is reduced by a factor of $\sim 0.3$ for the cases with $Z=10^{-2}~Z_{\odot}$ regardless of $\Sigma$. 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.

The star formation in molecular clouds is inefficient. The ionizing EUV radiation ($h \nu \geq 13.6$ eV) from young clusters has been considered as a primary feedback effect to limit the star formation efficiency (SFE). We here focus on effects of the stellar FUV radiation (6 eV $\leq h \nu \leq$ 13.6 eV) during the cloud disruption stage. The FUV radiation may further reduce the SFE via photoelectric heating, and it also affects the chemical states of the gas that is not converted to stars ("cloud remnants") via photodissociation of molecules. We have developed a one-dimensional semi-analytic model which follows the evolution of both the thermal and chemical structure of a photodissociation region (PDR) during the dynamical expansion of an HII region. We investigate how the FUV feedback limits the SFE, supposing that the star formation is quenched in the PDR where the temperature is above a threshold value (e.g., 100K). Our model predicts that the FUV feedback contributes to reduce the SFEs for the massive ($M_{\rm cl} \gtrsim 10^5 M_{\odot}$) clouds with the low surface densities ($\Sigma_{\rm cl} \lesssim 100$ M$_{\odot}$pc$^{-2}$). Moreover, we show that a large part of the H$_2$ molecular gas contained in the cloud remnants should be "CO-dark" under the FUV feedback for a wide range of cloud properties. Therefore, the dispersed molecular clouds are potential factories of the CO-dark gas, which returns into the cycle of the interstellar medium.

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.

We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.

Coalescence of intermediate-mass black holes (IMBHs) as a result of the migration toward galactic centers via dynamical friction may contribute to the formation of supermassive BHs. Here we reinvestigate the gaseous dynamical friction, which was claimed to be inefficient with radiative feedback from BHs in literature, by performing 3D radiation-hydrodynamics simulations that solve the flow structure in the vicinity of BHs. We consider a $10^4~M_\odot$ BH moving at the velocity $V_{\rm flow}$ through the homogeneous medium with metallicity $Z$ in the range of $0-0.1~Z_\odot$ and density $n_{\infty}$. We show that, if $n_{\infty} \lesssim 10^{6}~{\rm cm^{-3}}$ and $V_{\rm flow} \lesssim 60~{\rm km~s^{-1}}$, the BH is accelerated forward because of the gravitational pull from a dense shell ahead of an ionized bubble around the BH, regardless of the value of $Z$. If $n_{\infty} \gtrsim 10^{6}~{\rm cm^{-3}}$, however, our simulation shows the opposite result. The ionized bubble and associating shell temporarily appear, but immediately go downstream with significant ram pressure of the flow. They eventually converge into a massive downstream wake, which gravitationally drags the BH backward. The BH decelerates over the timescale of $\sim 0.01$~Myr, much shorter than the dynamical timescale in galactic disks. Our results suggest that IMBHs that encounter the dense clouds rapidly migrate toward galactic centers, where they possibly coalescence with others.

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

The super-Eddington accretion onto intermediate seed BHs is a potential formation mode of supermassive black holes exceeding $10^9~M_\odot$ in the early universe. We here investigate how such rapid accretion may occur with finite amounts of heavy elements contained in the gas and dust. In our 1D radiation-hydrodynamics simulations, the radiative transfer is solved for both the direct UV lights emitted by an accretion disk and the diffuse IR lights thermally emitted by dust grains. Our results show that the radiative force by the IR lights causes a strong feedback to regulate the mass accretion. The resulting mean accretion rate is lower with the higher metallicity, and there is the critical metallicity $Z \sim 10^{-2}~Z_\odot$, above which the super-Eddington accretion is prevented by the radiation pressure of the IR lights. With this taken into account, we examine if the dusty super-Eddington accretion occurs in young galaxies using a simple model. We show that a sufficient number of galaxies at $z \gtrsim 10$ can be such potential sites if BHs accrete the cold dense gas with $T \sim 10^2$ K, approximately the thermal equilibrium value at $Z = 10^{-2}~Z_\odot$. We argue that the efficiency of the BH growth via the rapid accretion depends on the metallicity, and that the metallicity slightly lower than $10^{-2}~Z_\odot$ provides a chance for the most efficient growth.

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.

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

Following our recent work devoted to the effect of accretion on the pre-main-sequence evolution of low-mass stars, we perform a detailed analysis of episodic excursions of low-mass protostars in the Hertzsprung-Russell (H-R) diagram triggered by strong mass accretion bursts typical of FU Orionis-type objects (FUors). These excursions reveal themselves as sharp increases in the stellar total luminosity and/or effective temperature of the protostar and can last from hundreds to a few thousands of years, depending on the burst strength and characteristics of the protostar. During the excursions, low-mass protostars occupy the same part of the H-R diagram as young intermediate-mass protostars in the quiescent phase of accretion. Moreover, the time spent by low-mass protostars in these regions is on average a factor of several longer than that spent by the intermediate-mass stars in quiescence. During the excursions, low-mass protostars pass close to the position of most known FUors in the H-R diagram, but owing to intrinsic ambiguity the model stellar evolutionary tracks are unreliable in determining the FUor properties. We find that the photospheric luminosity in the outburst state may dominate the accretion luminosity already after a few years after the onset of the outburst, meaning that the mass accretion rates of known FUors inferred from the bolometric luminosity may be systematically overestimated, especially in the fading phase.

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.

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.

We perform a suite of radiation hydrodynamics simulations of photoevaporating disks with varying the metallicity in a wide range of $10^{-3} \, Z_\odot \leq Z \leq 10^{0.5} \, Z_\odot $. We follow the disk evolution for over $\sim 5000$ years by solving hydrodynamics, radiative transfer, and non-equilibrium chemistry. Our chemistry model is updated from the first paper of this series by adding X-ray ionization and heating. We study the metallicity dependence of the disk photoevaporation rate and examine the importance of X-ray radiation. In the fiducial case with solar metallicity, including the X-ray effects does not significantly increase the photoevaporation rate when compared to the case with ultra-violet (UV) radiation only. At sub-solar metallicities in the range of $Z \gtrsim 10^{-1.5} \, Z_\odot $, the photoevaporation rate increases as metallicity decreases owing to the reduced opacity of the disk medium. The result is consistent with the observational trend that disk lifetimes are shorter in low metallicity environments. Contrastingly, the photoevaporation rate decreases at even lower metallicities of $Z \lesssim 10^{-1.5} \, Z_\odot $, because dust-gas collisional cooling remains efficient compared to far UV photoelectric heating whose efficiency depends on metallicity. The net cooling in the interior of the disk suppresses the photoevaporation. However, adding X-ray radiation significantly increases the photoevaporation rate, especially at $Z \sim 10^{-2}\, Z_\odot$. Although the X-ray radiation itself does not drive strong photoevaporative flows, X-rays penetrate deep into the neutral region in the disk, increase the ionization degree there, and reduce positive charges of grains. Consequently, the effect of photoelectric heating by far UV radiation is strengthened by the X-rays and enhances the disk photoevaporation.

We report ALMA Cycle 3 observations in CO isotopes toward a dense core, MC27/L1521F in Taurus, which is considered to be at an early stage of multiple star formation in a turbulent environment. Although most of the high-density parts of this core are considered to be as cold as $\sim$10 K, high-angular resolution ($\sim$20 au) observations in $^{12}$CO ($J$ = 3--2) revealed complex warm ($>$15--60 K) filamentary/clumpy structures with the sizes from a few tens of au to $\sim$1,000 au. The interferometric observations of $^{13}$CO and C$^{18}$O show that the densest part with arc-like morphologies associated with the previously identified protostar and condensations are slightly redshifted from the systemic velocity of the core. We suggest that the warm CO clouds may be consequences of shock heating induced by interactions among the different density/velocity components that originated from the turbulent motions in the core. However, such a small-scale and fast turbulent motion does not correspond to a simple extension of the line-width-size relation (i.e., Larson'{}s law), and thus the actual origin remains to be studied. The high-angular resolution CO observations are expected to be essential in detecting small-scale turbulent motions in dense cores and to investigate protostar formation therein.

The formation of supermassive stars (SMSs) is a possible pathway to seed supermassive black holes in the early universe. This chapter summarizes recent theoretical efforts to understand their evolution, highlighting effects of very rapid accretion at the rates of $\dot{M}_* \gtrsim 0.1~$M$_\odot$ yr$^{-1}$. Stellar evolution calculations predict that such an accreting SMS has a characteristic feature, i.e., a very large radius that monotonically increases as the stellar mass increases. The radius exceeds $7000$ R$_\odot$ (or $30$ AU) after the star accretes the gas of more than $10^3$ M$_\odot$. We show that the emergence of the "supergiant protostar" stage is a key for the formation of SMSs, because resulting radiative feedback against the accretion flow is substantially weakened during this stage. We also show that these SMSs end their lives while the accretion continues, and that their final fates vary with different accretion rates.

We conduct a theoretical study of the formation of massive stars over a wide range of metallicities from 1e-5 to 1Zsun and evaluate the star formation efficiencies (SFEs) from prestellar cloud cores taking into account multiple feedback processes. Unlike for simple spherical accretion, in the case of disk accretion feedback processes do not set upper limits on stellar masses. At solar metallicity, launching of magneto-centrifugally-driven outflows is the dominant feedback process to set SFEs, while radiation pressure, which has been regarded to be pivotal, has only minor contribution even in the formation of over-100Msun stars. Photoevaporation becomes significant in over-20Msun star formation at low metallicities of <1e-2Zsun, where dust absorption of ionizing photons is inefficient. We conclude that if initial prestellar core properties are similar, then massive stars are rarer in extremely metal-poor environments of 1e-5 - 1e-3Zsun. Our results give new insight into the high-mass end of the initial mass function and its potential variation with galactic and cosmological environments.

We present the first simulations of the formation and feedback of massive stars which account for radiation forces as well as photoionization feedback (along with protostellar outflows). In two different accretion scenarios modeled, we determine the relative strength of these feedback components and derive the size of the reservoir from which the forming stars gained their masses. We find that protostellar outflows alone limit the stellar mass growth only in an accretion scenario with a finite mass reservoir; when including accretion and ram pressure from large scales (> 0.1 pc), protostellar outflows do not limit stellar mass growth at all. Photoionization and HII regions dominate the feedback ladder only at later times and only on large scales. Specifically, photoionization yields a broadening of the bipolar outflow cavities and a reduction of the gravitational infall momentum by about 50 percent, but does not limit the stellar mass accretion. On the other hand, we find radiation forces restrain the gravitational infall toward the circumstellar disk, impact the gravito-centrifugal equilibrium at the outer edge of the disk, and eventually shut down stellar accretion completely. The most massive star formed in the simulations accreted 95 solar masses before disk destruction; this mass was drawn-in from an accretion reservoir of approximately 240 solar masses and 0.24 pc in radius. In the regime of very massive stars, the final mass of these stars is controlled by their own radiation force feedback.

Protoplanetary disks are thought to have lifetimes of several million years in the solar neighborhood, but recent observations suggest that the disk lifetimes are shorter in a low metallicity environment. We perform a suite of radiation hydrodynamics simulations of photoevaporation of protoplanetary disks to study the disk structure and its long-term evolution of $\sim 10000$ years, and the metallicity dependence of mass-loss rate. Our simulations follow hydrodynamics, extreme and far ultra-violet radiative transfer, and non-equilibrium chemistry in a self-consistent manner. Dust grain temperatures are also calculated consistently by solving the radiative transfer of the stellar irradiation and grain (re-)emission. We vary the disk gas metallicity over a wide range of $10^{-4}~ Z_\odot \leq Z \leq 10 ~Z_\odot$. The photoevaporation rate is lower with higher metallicity in the range of $10^{-1} \,Z_\odot \lesssim Z \lesssim 10 \,Z_\odot$, because dust shielding effectively prevents far-ultra violet (FUV) photons from penetrating into and heating the dense regions of the disk. The photoevaporation rate sharply declines at even lower metallicities in $10^{-2} \,Z_\odot \lesssim Z \lesssim 10^{-1}\,Z_\odot$, because FUV photoelectric heating becomes less effective than dust-gas collisional cooling. The temperature in the neutral region decreases, and photoevaporative flows are excited only in an outer region of the disk. At $10^{-4}\,Z_\odot \leq Z \lesssim 10^{-2}\,Z_\odot$, HI photoionization heating acts as a dominant gas heating process and drives photoevaporative flows with roughly a constant rate. The typical disk lifetime is shorter at $Z=0.3~Z_\odot$ than at $Z = Z_\odot$, being consistent with recent observations of the extreme outer galaxy.

The driving mechanism of massive outflows observed in high-mass star-forming regions is investigated using three-dimensional magnetohydrodynamics (MHD) and protostellar evolution calculations. In our previous paper, we showed that the mass outflow rate depends strongly on the mass accretion rate onto the circumstellar disk around a high-mass protostar, and massive outflows may be driven by the magnetic effect in high-mass star-forming cores. In the present study, in order to verify that the MHD disk wind is the primary driving mechanism of massive outflows, we quantitatively compare outflow properties obtained through simulations and observations. Since the outflows obtained through simulations are slightly younger than those obtained through observations, the time-integrated quantities of outflow mass, momentum, and kinetic energy are slightly smaller than those obtained through observations. On the other hand, time-derivative quantities of mass ejection rate, outflow momentum flux, and kinetic luminosity obtained through simulations are in very good agreement with those obtained through observations. This indicates that the MHD disk wind greatly contributes to the massive outflow driving from high-mass protostars, and the magnetic field might significantly control the high-mass star formation process.

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

We report ALMA observations in 0.87 mm continuum and $^{12}$CO ($J$ = 3--2) toward a very low-luminosity ($<$0.1 $L_{\odot}$) protostar, which is deeply embedded in one of the densest core MC27/L1521F, in Taurus with an indication of multiple star formation in a highly dynamical environment. The beam size corresponds to $\sim$20 AU, and we have clearly detected blueshifted/redshifted gas in $^{12}$CO associated with the protostar. The spatial/velocity distributions of the gas show there is a rotating disk with a size scale of $\sim$10 AU, a disk mass of $\sim$10$^{-4}$ $M_{\odot}$ and a central stellar mass of $\sim$0.2 $M_{\odot}$. The observed disk seems to be detached from the surrounding dense gas, although it is still embedded at the center of the core whose density is $\sim$10$^{6}$ cm$^{-3}$. The current low-outflow activity and the very low luminosity indicate that the mass accretion rate onto the protostar is extremely low in spite of a very early stage of star formation. We may be witnessing the final stage of the formation of $\sim$0.2 $M_{\odot}$ protostar. However, we cannot explain the observed low luminosity with the standard pre-main-sequence evolutionary track unless we assume cold accretion with an extremely small initial radius of the protostar ($\sim$0.65 $R_\odot$). These facts may challenge our current understanding of the low mass star formation, in particular the mass accretion process onto the protostar and the circumstellar disk.

The origin of super-massive black holes in the early universe remains poorly understood.Gravitational collapse of a massive primordial gas cloud is a promising initial process,but theoretical studies have difficulty growing the black hole fast enough.We report numerical simulations of early black hole formation starting from realistic cosmological conditions.Supersonic gas motions left over from the Big Bang prevent early gas cloud formation until rapid gas condensation is triggered in a proto-galactic halo. A protostar is formed in the dense, turbulent gas cloud, and it grows by sporadic mass accretion until it acquires 34,000 solar masses.The massive star ends its life with a catastrophic collapse to leave a black hole -- a promising seed for the formation of a monstrous black hole.

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 \mu {\rm m}$ at $1Z_{\odot}$, shifts to $< 3 \mu {\rm m}$ at $Z \lesssim 0.1 Z_\odot$. Second, a characteristic feature at $10 \mu \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.

The pre-main-sequence evolution of low-mass stars and brown dwarfs is studied numerically starting from the formation of a protostellar/proto-brown dwarf seed and taking into account the mass accretion onto the central object during the initial several Myr of evolution. The stellar evolution was computed using the STELLAR evolution code developed by Yorke & Bodenheimer with recent modifications by Hosokawa et al. The mass accretion rates were taken from numerical hydrodynamics models of Vorobyov & Basu computing the circumstellar disk evolution starting from the gravitational collapse of pre-stellar cloud cores of various mass and angular momentum. The resulting stellar evolution tracks were compared with the isochrones and isomasses calculated using non-accreting models. We find that mass accretion in the initial several Myr of protostellar evolution can have a strong effect on the subsequent evolution of young stars and brown dwarfs. The disagreement between accreting and non-accreting models in terms of the total stellar luminosity L_st, stellar radius R_st and effective temperature T_eff depends on the thermal efficiency of accretion, i.e., on the fraction of accretion energy absorbed by the central object. The largest mismatch is found for the cold accretion case, in which essentially all accretion energy is radiated away. The relative deviations in L_st and R_st in this case can reach 50% for 1.0-Myr-old objects and remain notable even for 10-Myr-old objects. In the hot and hybrid accretion cases, in which a constant fraction of accretion energy is absorbed, the disagreement between accreting and non-accreting models becomes less pronounced, but still remains notable for 1.0-Myr-old objects. These disagreements may lead to the wrong age estimate for objects of (sub-)solar mass when using the isochrones based on non-accreting models (abridged).

The relation between the mass accretion rate onto the circumstellar disc and the rate of mass ejection by magnetically driven winds is investigated using three-dimensional magnetohydrodynamics simulations. Using a spherical cloud core with a varying ratio of thermal to gravitational energy, which determines the mass accretion rate onto the disc, to define the initial conditions, the outflow propagation for approximately 10^4 yr after protostar formation is then calculated for several cloud cores. The mass ejection rate and accretion rate are comparable only when the magnetic energy of the initial cloud core is comparable to the gravitational energy. Consequently, in strongly magnetised clouds a higher mass accretion rate naturally produces both massive protostars and massive outflows. The simulated outflow mass, momentum, kinetic energy and momentum flux agree well with observations, indicating that massive stars form through the same mechanism as low-mass stars but require a significantly strong magnetic field to launch massive outflows.

We study formation and long-term evolution of a circumstellar disk in a collapsing molecular cloud core using a resistive magnetohydrodynamic simulation. While the formed circumstellar disk is initially small, it grows as accretion continues and its radius becomes as large as 200 AUs toward the end of the Class-I phase. A pair of grand-design spiral arms form due to gravitational instability in the disk, and they transfer angular momentum in the highly resistive disk. Although the spiral arms disappear in a few rotations as expected in a classical theory, new spiral arms form recurrently as the disk soon becomes unstable again by gas accretion. Such recurrent spiral arms persist throughout the Class-0 and I phase. We then perform synthetic observations and compare our model with a recent high-resolution observation of a young stellar object Elias 2-27, whose circumstellar disk has grand design spiral arms. We find good agreement between our theoretical model and the observation. Our model suggests that the grand design spiral arms around Elias 2-27 are consistent with material arms formed by gravitational instability. If such spiral arms commonly exist in young circumstellar disks, it implies that young circumstellar disks are considerably massive and gravitational instability is the key process of angular momentum transport.

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.

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

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.

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.

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

We describe an overall picture of galactic-scale star formation. Recent high-resolution magneto-hydrodynamical simulations of two-fluid dynamics with cooling/heating and thermal conduction have shown that the formation of molecular clouds requires multiple episodes of supersonic compression. This finding enables us to create a scenario in which molecular clouds form in interacting shells or bubbles on a galactic scale. First we estimate the ensemble-averaged growth rate of molecular clouds over a timescale larger than a million years. Next we perform radiation hydrodynamics simulations to evaluate the destruction rate of magnetized molecular clouds by the stellar FUV radiation. We also investigate the resultant star formation efficiency within a cloud which amounts to a low value (a few percent) if we adopt the power-law exponent -2.5 for the mass distribution of stars in the cloud. We finally describe the time evolution of the mass function of molecular clouds over a long timescale (>1Myr) and discuss the steady state exponent of the power-law slope in various environments.

A promising formation channel of SMBHs at redshift 6 is the so-called DC model, which posits that a massive seed BH forms through gravitational collapse of a $\sim 10^5~M_\odot$ SMS. We study the evolution of such a SMS growing by rapid mass accretion. In particular, we examine the impact of time-dependent mass accretion of repeating burst and quiescent phases that are expected to occur with a self-gravitating circumstellar disk. We show that the stellar evolution with such episodic accretion differs qualitatively from that expected with a constant accretion rate, even if the mean accretion rate is the same. Unlike the case of constant mass accretion, whereby the star expands roughly following $R_* \simeq 2.6 \times 10^3 R_\odot (M_*/100~M_\odot)^{1/2}$, the protostar can substantially contract during the quiescent phases between accretion bursts. The stellar effective temperature and ionizing photon emissivity increase accordingly as the star contracts, which can cause strong ionizing feedback and halt the mass accretion onto the star. With a fixed duration of the quiescent phase $\Delta t_{\rm q}$, such contraction occurs in early evolutionary phases, i.e. for $M_* \lesssim 10^3~M_\odot$ with $\Delta t_{\rm q} \simeq 10^3$ yr. For later epochs and larger masses but the same $\Delta t_{\rm q}$, contraction is negligible even during quiescent phases. With larger quiescent times $\Delta t_{\rm q}$, however, the star continues to contract during quiescent phases even for the higher stellar masses. We show that such behavior is well understood by comparing the interval time and the thermal relaxation time for a bloated surface layer. We conclude that the UV radiative feedback becomes effective if the quiescent phase associated by the burst accretion is longer than $\sim 10^3$ yr, which is possible in an accretion disk forming in the direct collapse model.

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.

We present radiation transfer (RT) simulations of evolutionary sequences of massive protostars forming from massive dense cores in environments of high surface densities. The protostellar evolution is calculated with a detailed multi-zone model, with the accretion rate regulated by feedback from an evolving disk-wind outflow cavity. Disk and envelope evolutions are calculated self-consistently. In this framework, an evolutionary track is determined by three environmental initial conditions: the initial core mass M_c, the mean surface density of the ambient star-forming clump Sigma_cl, and the rotational-to-gravitational energy ratio of the initial core, beta_c. Evolutionary sequences with various M_c, Sigma_cl, beta_c are constructed. We find that in a fiducial model with M_c=60Msun, Sigma_cl=1 g/cm^2 and beta_c=0.02, the final star formation efficiency >~0.43. For each evolutionary track, RT simulations are performed at selected stages, with temperature profiles, SEDs, and images produced. At a given stage the envelope temperature is highly dependent on Sigma_cl, but only weakly dependent on M_c. The SED and MIR images depend sensitively on the evolving outflow cavity, which gradually wides as the protostar grows. The fluxes at <~100 microns increase dramatically, and the far-IR peaks move to shorter wavelengths. We find that, despite scatter caused by different M_c, Sigma_cl, beta, and inclinations, sources at a given evolutionary stage appear in similar regions on color-color diagrams, especially when using colors at >~ 70 microns, where the scatter due to the inclination is minimized, implying that such diagrams can be useful diagnostic tools of evolutionary stages of massive protostars. We discuss how intensity profiles along or perpendicular to the outflow axis are affected by environmental conditions and source evolution.

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.

The generation process of magnetic field around a proto-first-star is studied. Utilizing the recent numerical result of proto-first-star formation based upon the radiation hydrodynamics simulations, we assess the magnetic field strength generated by the radiative force and the Biermann battery effect. We find that magnetic field of \sim 10^{-9} G is generated on the surface of the accretion disk around the proto-first-star. The field strength on the accretion disk is smaller by two orders of magnitude than the critical value, above which the gravitational fragmentation of the disk is suppressed. Thus, the generated seed magnetic field hardly affect the dynamics of on-site first star formation directly, unless efficient amplification process is taken into consideration. We also find that the generated magnetic field is continuously blown out from the disk on the outflows to the poles, that are driven by the thermal pressure of photoheated gas. The strength of the diffused magnetic field in low density regions is \sim 10^{-14}-10^{-13}G at n_H=10^3 cm^{-3} which could play important roles on the next generation star formation, as well as the seeds of magnetic field exist in present-day universe.

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.

The 6.7 GHz methanol maser emission, a tracer of forming massive stars, sometimes shows enigmatic periodic flux variations over several 10-100 days. In this Letter, we propose that this periodic variations could be explained by the pulsation of massive protostars growing under rapid mass accretion with rates of Mdot > 10^-3 Msun/yr. Our stellar evolution calculations predict that the massive protostars have very large radius exceeding 100 Rsun at maximum, and we here study the pulsational stability of such the bloated protostars by way of the linear stability analysis. We show that the protostar becomes pulsationally unstable with various periods of several 10-100 days, depending on different accretion rates. With the fact that the stellar luminosity when the star is pulsationally unstable also depends on the accretion rate, we derive the period-luminosity relation log (L/Lsun) = 4.62 + 0.98log(P/100 day), which is testable with future observations. Our models further show that the radius and mass of the pulsating massive protostar should also depend on the period. It would be possible to infer such protostellar properties and the accretion rate with the observed period. Measuring the maser periods enables a direct diagnosis of the structure of accreting massive protostars, which are deeply embedded in dense gas and inaccessible with other observations.

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.

The evolution of protostellar outflow is investigated with resistive magneto-hydrodynamic nested-grid simulations that cover a wide range of spatial scales (\sim 1AU - 1pc). We follow cloud evolution from the pre-stellar core stage until the infalling envelope dissipates long after the protostar formation. We also calculate protostellar evolution to derive protostellar luminosity with time-dependent mass accretion through a circumstellar disk. The protostellar outflow is driven by the first core prior to protostar formation and is directly driven by the circumstellar disk after protostar formation. The opening angle of the outflow is large in the Class 0 stage. A large fraction of the cloud mass is ejected in this stage, which reduces the star formation efficiency to \sim 50%. After the outflow breaks out from the natal cloud, the outflow collimation is gradually improved in the Class I stage. The head of the outflow travels more than \sim 10^5AU in \sim 10^5 yr. The outflow momentum, energy and mass derived in our calculations agree well with observations. In addition, our simulations show the same correlations among outflow momentum flux, protostellar luminosity and envelope mass as those in observations. These correlations differ between Class 0 and I stages, which is explained by different evolutionary stages of the outflow; in the Class 0 stage, the outflow is powered by the accreting mass and acquires its momentum from the infalling envelope; in the Class I stage, the outflow enters the momentum-driven snow-plough phase. Our results suggest that protostellar outflow should determine the final stellar mass and significantly affect the early evolution of low-mass protostars.

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.

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.

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.

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.

We revisit the problem of low-mass pre-main-sequence (PMS) stellar evolution and its observational consequences for where stars fall on the Hertzsprung-Russell diagram (HRD). In contrast to most previous work, our models follow stars as they grow from small masses via accretion, and we perform a systematic study of how the stars' HRD evolution is influenced by their initial radius, by the radiative properties of the accretion flow, and by the accretion history, using both simple idealized accretion histories and histories taken from numerical simulations of star cluster formation. We compare our numerical results to both non-accreting isochrones and to the positions of observed stars in the HRD, with a goal of determining whether both the absolute ages and the age dispersions inferred from non-accreting isochrones are reliable. We show that non-accreting isochrones can sometimes overestimate stellar ages for more massive stars (those with effective temperatures above \sim 3500 K), thereby explaining why non-accreting isochrones often suggest a systematic age difference between more and less massive stars in the same cluster. However, we also find the only way to produce a similar overestimate for the ages of cooler stars is by selecting parameters that are strongly inconsistent with both theoretical models and observational constraints. We therefore conclude that inferred ages and age spreads for cool stars are reliable, at least to the extent that the observed bolometric luminosities and temperatures are accurate. Finally, we note that the time-dependence of the mass accretion rate has remarkably little effect on low-mass stars' evolution on the HRD, and that such time-dependence may be neglected for all stars except those with effective temperatures above \sim 4000 K.

There is currently no accepted theoretical framework for the formation of the most massive stars, and the manner in which protostars continue to accrete and grow in mass beyond \sim10Msun is still a controversial topic. In this study we use several prescriptions of stellar accretion and a description of the Galactic gas distribution to simulate the luminosities and spatial distribution of massive protostellar population of the Galaxy. We then compare the observables of each simulation to the results of the Red MSX Source (RMS) survey, a recently compiled database of massive young stellar objects. We find that the observations are best matched by accretion rates which increase as the protostar grows in mass, such as those predicted by the turbulent core and competitive accretion (i.e. Bondi-Hoyle) models. These 'accelerating accretion' models provide very good qualitative and quantitative fits to the data, though we are unable to distinguish between these two models on our simulations alone. We rule out models with accretion rates which are constant with time, and those which are initially very high and which fall away with time, as these produce results which are quantitatively and/or qualitatively incompatible with the observations. To simultaneously match the low- and high-luminosity YSO distribution we require the inclusion of a 'swollen-star' pre-main-sequence phase, the length of which is well-described by the Kelvin-Helmholz timescale. Our results suggest that the lifetime of the YSO phase is \sim 10^5yrs, whereas the compact Hii-region phase lasts between \sim 2 - 4 \times 10^5yrs depending on the final mass of the star. Finally, the absolute numbers of YSOs are best matched by a globally averaged star-formation rate for the Galaxy of 1.5-2Msun/yr.

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.

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.

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.

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.

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)

To investigate the process of star formation triggered by the expansion of an HII region, we present a multi-wavelength analysis of the Galactic HII region RCW 120 and its surroundings. The collect and collapse model predicts that the layer of gas and dust accumulated between the ionization and shock fronts during the expansion of the HII region collapses and forms dense fragments, giving rise to potential sites of massive-star formation. The aim of our study is to look for such massive fragments and massive young stars on the borders of RCW 120. We mapped the RCW 120 region in the cold dust continuum emission at 1.2 mm to search for these fragments. We supplemented this study with the available near- (2MASS) and mid-IR (GLIMPSE) data to locate the IR sources observed towards this region and to analyse their properties. We then compared the observational results with the predictions of Hosokawa & Inutsuka's model (2005, 2006). At 1.2 mm we detected eight fragments towards this region, five located on its borders. The largest fragment has a mass of about 370 Msun. Class I and Class II young stellar objects are detected all over the region, with some observed far from the ionization front. This result emphasises the possible importance distant interactions between the radiation, escaping from the ionized region, and the surrounding medium.

We examine the triggering process of molecular cloud formation around diffuse HII regions. We calculate the time evolution of the shell as well as of the HII region in a two-phase neutral medium, solving the UV and FUV radiative transfer, the thermal and chemical processes in the time-dependent hydrodynamics code. In the cold neutral medium, the ambient gas is swept up in the cold (T \sim 30-40K) and dense (n \sim 10^3 cm^-3) shell around the HII region. In the shell, H_2 molecules are formed from the swept-up HI gas, but CO molecules are hardly formed. The reformation of H_2 molecules is more efficient with a higher-mass central star. The physical and chemical properties of gas in the shell are just intermediate between those of the neutral medium and molecular clouds observed by the CO emission. The dense shell with cold HI/H_2 gas easily becomes gravitationally unstable, and breaks up into small clouds. The cooling layer just behind the shock front also suffers from thermal instability, and will fragment into cloudlets with some translational motions. We suggest that the predicted cold ``dark'' HI/H_2 gas should be detected as the HI self-absorption (HISA) feature. We have sought such features in recent observational data, and found shell-like HISA features around the giant HII regions, W4 and W5. The shell-like HISA feature shows good spatial correlation with dust emission, but poor correlation with CO emission. Our quantitative analysis shows that the HISA cloud can be as cold as T \sim a few x 10K. (abridged)

We report the Chandra detection of hard X-ray emission from the Welch ring in W49A, an organized structure of ultra-compact (UC) HII regions containing a dozen nascent early-type stars. Two UC HII regions are associated with hard X-ray emission in a deep Advanced CCD Imaging Spectrometer image exposed for 96.7 ks. One of the two X-ray sources has no near-infrared counterpart and is extended by ~5 arcsec, or ~0.3 pc, at a distance of ~11.4 kpc, which is spatially aligned with the cometary radio continuum emission associated with the UC HII region. The X-ray spectrum of the emission, when fit with a thermal model, indicates a heavily absorbed plasma with extinction of \~5x10^{23}/cm^{2}, temperature of ~7 keV, and X-ray luminosity in the 3.0-8.0 keV band of ~3x10^{33} ergs/s. Both the luminosity and the size of the emission resemble the extended hard emission found in UC HII regions in Sagittarius B2, yet they are smaller by an order of magnitude than the emission found in massive star clusters such as NGC 3603. Three possibilities are discussed for the cause of the hard extended emission in the Welch ring: an ensemble of unresolved point sources, shocked interacting winds of the young O stars, and a wind-blown bubble interacting with ambient cold matter.

We study the expansion of the ionization and dissociation fronts (DFs) in a radially stratified molecular cloud, whose density distribution is represented as n(r) \propto r^-w. We focus on cases with w \leq 1.5, when the ionization front is ``trapped'' in the cloud and expands with the preceding shock front. The simultaneous evolution of the outer photodissociation region (PDR) is examined in detail. First, we analytically probe the time evolution of the column densities of the shell and envelope outside the HII region, which are key physical quantities for the shielding of dissociating photons. Next, we perform numerical calculations, and study how the thermal/chemical structure of the outer PDR changes with different density gradients. We apply our numerical model to the Galactic HII region, Sharpless 219 (Sh219). The time evolution of the column densities of the shell and outer envelope depends on w, and qualitatively changes across w = 1. In the cloud with w < 1, the shell column density increases as the HII region expands. The DFs are finally trapped in the shell, and the molecular gas gradually accumulates in the shell. The molecular shell and envelope surround the HII region. With w > 1, on the other hand, the shell column density initially increases, but finally decreases. The column density of the outer envelope also quickly decreases as the HII region swells up. It becomes easier and easier for the dissociating photons to penetrate the shell and envelope. The PDR broadly extends around the trapped HII region. A model with w = 1.5 successfully explains the observational properties of Sh219. Our model suggests that a density-bounded PDR surrounds the photon-bounded HII region in Sh219.

We quantitatively examine the significance of star formation triggered in the swept-up shell around an expanding HII region. If the swept-up molecular gas is sufficiently massive, new OB stars massive enough to repeat the triggering process will form in the shell. We determine the lower limit (M_thr) for the mass of the star that sweeps up the molecular gas, where at least one new star with mass M_* > M_thr forms after the shell fragmentation. To calculate the threshold stellar mass, M_thr, we examine how massive molecular shells can form around various central stars, by performing detailed numerical radiation hydrodynamics calculations. The mass of the photodissociated gas is generally larger than the mass of the photoionized gas. However, the swept-up molecular mass exceeds the photodissociated mass with a higher-mass star of M_* > 20 Msun. The accumulated molecular mass generally increases with the stellar mass, and amounts to 10^{4-5} Msun for M_* > 20 Msun with an ambient density of n \sim 100/cc. The threshold stellar mass is M_thr \sim 18 Msun with the star-formation efficiency of \epsilon \sim 0.1 and n \sim 100/cc. We examine the generality of this mode of run-away triggering for different sets of parameters, and found that M_thr \sim 15-20 Msun in various situations. If the ambient density is too high or the star-formation efficiency is too low, the triggering is not run-away, but a single event.

We analyze the dynamical expansion of the HII region, photodissociation region, and the swept-up shell, solving the UV- and FUV-radiative transfer, the thermal and chemical processes in the time-dependent hydrodynamics code. Following our previous paper, we investigate the time evolutions with various ambient number densities and central stars. Our calculations show that basic evolution is qualitatively similar among our models with different parameters. The molecular gas is finally accumulated in the shell, and the gravitational fragmentation of the shell is generally expected. The quantitative differences among models are well understood with analytic scaling relations. The detailed physical and chemical structure of the shell is mainly determined by the incident FUV flux and the column density of the shell, which also follow the scaling relations. The time of shell-fragmentation, and the mass of the gathered molecular gas are sensitive tothe ambient number density. In the case of the lower number density, the shell-fragmentation occurs over a longer timescale, and the accumulated molecular gas is more massive. The variations with different central stars are more moderate. The time of the shell-fragmentation differs by a factor of several with the various stars of M_* = 12-101 M_sun. According to our numerical results, we conclude that the expanding HII region should be an efficient trigger for star formation in molecular clouds if the mass of the ambient molecular material is large enough.

We analyze the dynamical expansion of the HII region and outer photodissociation region (PDR) around a massive star by solving the UV and FUV radiation transfer and the thermal and chemical processes in a time-dependent hydrodynamics code. We focus on the physical structure of the shell swept up by the shock front (SF) preceding the ionization front (IF). After the IF reaches the initial Stromgren radius, the SF emerges in front of the IF and the geometrically thin shell bounded with the IF and the SF is formed. The gas density inside the shell is about 10-100 times as high as the ambient gas density. Initially the dissociation fronts expands faster than IF and the PDR is formed outside the HII region. Thereafter the IF and SF gradually overtakes the proceeding dissociation fronts (DFs), and eventually DFs are taken in the shell. The chemical composition within the shell is initially atomic, but hydrogen and carbon monoxide molecules are gradually formed. This is partly because the IF and SF overtake DFs and SF enters the molecular region, and partly because the reformation timescales of the molecules become shorter than the dynamical timescale. The gas shell becomes dominated by the molecular gas by the time of gravitational fragmentation, which agrees with some recent observations. A simple estimation of star formation rate in the shell can provide a significant star formation rate in our galaxy.

We study the demographics of supermassive black holes (SMBHs) in the local and high-z universe with hard X-ray luminosity functions (HXLFs) of AGNs. First, we compare the mass accretion history at z>0 with optical luminosity functions (OLFs) and that with HXLFs. We consider the constraints on parameters of mass accretion (energy conversion efficiency and Eddington ratio) and conclude that the constraints based on HXLFs are more adequate rather than that based on OLFs. The sub-Eddington case is allowed only when we use HXLFs. Next, we estimate the upper limit of the cumulative mass density of SMBHs at any redshifts. For an application, we examine if Lyman-Break galaxies (LBGs) at z=3 already have SMBHs in their centers which is suggested by recent observations. If most of LBGs already has massive SMBHs at z=3, the resultant mass density of SMBHs at z=0 should exceed the observational estimate because such SMBHs should further grow by accretion. We show that the special conditions should be met in order that a large part of LBGs have SMBHs (for example, large energy conversion efficiency and frequent mergers and/or direct formations at z>3). The possibility that nearly all LBGs have SMBHs with large mass ratio, such as M_BH/M_stellar > 0.005, is reliably ruled out.(abridged)

The detectability of gravitational wave (GW) radiation from accretion disks is discussed based on various astrophysical contexts. In order to emit GW radiation, the disk shape should lose axial symmetry. We point out that a significant deformation is plausible in non-radiative hot accretion disks because of enhanced magnetic activity, whereas it is unlikely for standard-type cool disks. We have analyzed the 3D magnetohydrodynamical (MHD) simulation data of magnetized accretion flow, finding non-axisymmetric density patterns. The corresponding ellipticity is $\epsilon \sim 0.01$. The expected time variations of GW radiation are overall chaotic, but there is a hint of quasi-periodicity. GW radiation has no interesting consequence, however, in the case of close binaries, because of very tiny disk masses. GW radiation is not significant, either, for AGN because of very slow rotation velocities. The most promising case can be found in gamma-ray bursts or supernovae, in which a massive torus (or disk) with a solar mass or so may be formed around a stellar-mass compact object as the result of a merger of compact objects, or by the fallback of exploded material towards the center in a supernova. Although much more intense GW radiation is expected before the formation of the torus, the detection of GW radiation in the subsequent accretion phase is of great importance, since it will provide a good probe to investigating their central engines.

Using the theoretical models of the QSO formation, we can reproduce optical QSO luminosity functions (LFs) at high redshifts (z > 2.2). Two different models can reproduce LFs successfully, though the lifetime of QSOs, t_Q, and the relation between the black hole mass and the host halo mass are different each other; t_Q = 10^6yr, in one model, t_Q > 10^7yr, in other models. Here, we propose a method to break this degeneracy. We calculate the mass function of supermassive black holes (SMBHs) at z=2.5, and compare the result with the current mass function obtained by Salucci et al.(1999). In the shorter lifetime model, the mass function at z=2.5 exceeds that of z=0.0 by one order of magnitude, then it should be ruled out. We conclude that the lifetime is at least t_Q > 10^7yr. Next, we examine the difference of the formation epoch of SMBHs existing at z=3.0 for each model under the model assumptions. We simply discuss the difference of formation epoch as another possible model-discriminator.

We use a large, heterogeneous sample of local active galactic nuclei (AGN) that includes Seyfert 1s, Seyfert 2s and PG quasars to investigate for the first time the relation between black hole mass (M_BH) and mid-infrared nuclear emission. We find a clear relation between M_BH and the 10 micron nuclear luminosity for these local AGNs. There are no significant differences between type 1 and type 2 objects, implying that the reprocessing of the 10 micron nuclear emission is not severely affected by geometric and optical depth effects. We also confirm that M_BH is related to the 2-10keV X-ray luminosity, but only for the Compton thin galaxies. We present a theoretical basis for these empirical relations and discuss possible reasons for the observed scatter. Our results show that rest-frame 10 micron and hard X-ray luminosities (especially the former, which is applicable to all AGN types) can be powerful tools for conducting a census of black hole masses at high redshift and for probing their cosmological evolution.

Cosmological evolution of the QSO luminosity functions (LFs) at NIR/optical/X-ray bands for 1.3 < z < 3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that although QSO luminosities only depend on mass-accretion rate, \Mdot, QSO spectra have a dependence on black-hole mass, M_{BH}, as well. The smaller M_{BH} is and/or the larger \Mdot is, the harder becomes the QSO NIR/optical/UV spectrum. We model disk spectra which can reproduce these features and calculated LFs for redshift z ~ 3 with the assumption of new-born QSOs being shining at the Eddington luminosity. The main results are: (i) the observed LFs at optical and X-rays can be simultaneously reproduced. (ii) LFs at optical and X-ray bands are not sensitive to M_{BH}, while LFs at NIR bands are; about one order of magnitude difference is expected in volume number densities at L_{I, J} ~ 10^{46} erg s^{-1} between the case that all QSOs would have the same spectral shape as that of M_{BH} = 10^{9} M_sun and the case with M_{BH} = 10^{11} M_sun. (iii) The resultant LFs at NIR are dominated by 10^{7} M_sun black-holes at L_{I, J} ~ 10^{44} erg s^{-1}, and by 10^{11} M_sun black-holes at L_{I, J} \~ 10^{46} erg s^{-1}. Future infrared observations from space(e.g.NGST) will probe cosmological evolution of black hole masses. For redshift z < 3, on the other hand, the observed optical/X-ray LFs can be fitted, if the initial QSO luminosity L_0 is below the Eddington luminosity. Interestingly, the best fitting values of l = L_0/L_{Edd} are different in B- and X-ray bands; l_B ~ 2.5 l_X. The reason for this discrepancy is briefly discussed.