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

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

We construct helium (He) star models with optically thick winds and compare them with the properties of Galactic Wolf-Rayet (WR) stars. Hydrostatic He-core solutions are connected smoothly to trans-sonic wind solutions that satisfy the regularity conditions at the sonic point. Velocity structures in the supersonic parts are assumed by a simple beta-type law. By constructing a center-to-surface structure, a mass-loss rate can be obtained as an eigenvalue of the equations. Sonic points appear at temperatures ~ 1.8e5 - 2.8e5 K below the Fe-group opacity peak, where the radiation force becomes comparable to the local gravity. Photospheres are located at radii 3-10 times larger than sonic points. The obtained mass-loss rates are comparable to those of WR stars. Our mass-loss rate - luminosity relation agrees well with the relation recently obtained by Graefener et al. (2017). Photospheric temperatures of WR stars tend to be cooler than our predictions. We discuss the effects of stellar evolution, detailed radiation transfer, and wind clumping, which are ignored in this paper.

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.

Supermassive black holes observed at high redshift $z\gtrsim6$ could grow from direct collapse black holes (DCBHs) with mass $\sim10^5\,M_{\odot}$, which result from the collapse of supermassive stars (SMSs). If a relativistic jet is launched from a DCBH, it can break out of the collapsing SMS and produce a gamma-ray burst (GRB). Although most of the GRB jets are off-axis from our line of sight, we show that the energy injected from the jet into a cocoon is huge $\sim10^{55-56}\,{\rm{erg}}$, so that the cocoon fireball is observed as ultra-luminous supernovae of $\sim10^{45-46}\rm{\,erg\,s^{-1}}$ for $\sim5000 [(1+z)/16] \rm{\,days}$. They are detectable by the future telescopes with near infrared bands, such as, $Euclid$, $WFIRST$, $WISH$, and $JWST$ up to $z\sim20$ and $\sim 100$ events per year, providing a direct evidence of the DCBH scenario.

The existence of black holes (BHs) of mass ~ 10^{9} M_sun at z > 6 is a big puzzle in astrophysics because even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 M_sun to grow into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive stars with mass M ~ 10^{5} M_sun into massive seed BHs. We find that if a jet is launched from the accretion disk around the central BH, the jet can break out the star because of the structure of the radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^{4} - 10^{6} s and flux of 10^{-11} - 10^{-8} erg s^{-1} cm^{-2} could be detectable by Swift. We estimate an event rate of < 1yr^{-1}. The total explosion energy is > 10^{55} - 10^{56} erg. The resulting negative feedback delays the growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely.

We examine the bright radio synchrotron counterparts of low-luminosity gamma-ray bursts (llGRBs) and relativistic supernovae (SNe) and find that they can be powered by spherical hypernova (HN) explosions. Our results imply that radio-bright HNe are driven by relativistic jets that are choked deep inside the progenitor stars or quasi-spherical magnetized winds from fast-rotating magnetars. We also consider the optical synchrotron counterparts of radio-bright HNe and show that they can be observed as precursors several days before the SN peak with an r-band absolute magnitude of M_r ~ -14 mag. While previous studies suggested that additional trans-relativistic components are required to power the bright radio emission, we find that they overestimated the energy budget of the trans-relativistic component by overlooking some factors related to the minimum energy of non-thermal electrons. If an additional trans-relativistic component exists, then a much brighter optical precursor with M_r ~ -20 mag can be expected. Thus, the scenarios of radio-bright HNe can be distinguished by using optical precursors, which can be detectable from < 100 Mpc by current SN surveys like the Kiso SN Survey, Palomar Transient Factory, and Panoramic Survey Telescope & Rapid Response System.

We investigate the possible origin of extended emissions (EEs) of short gamma-ray bursts with an isotropic energy of ~ 10^(50-51) erg and a duration of a few 10 s to ~ 100 s, based on a compact binary (neutron star (NS)-NS or NS-black hole (BH)) merger scenario. We analyze the evolution of magnetized neutrino-dominated accretion disks of mass ~ 0.1 M_sun around BHs formed after the mergers, and estimate the power of relativistic outflows via the Blandford-Znajek (BZ) process. We show that a rotation energy of the BH up to > 10^52 erg can be extracted with an observed time scale of > 30 (1+z) s with a relatively small disk viscosity parameter of alpha < 0.01. Such a BZ power dissipates by clashing with non-relativistic pre-ejected matter of mass M ~ 10^-(2-4) M_sun, and forms a mildly relativistic fireball. We show that the dissipative photospheric emissions from such fireballs are likely in the soft X-ray band (1-10 keV) for M ~ 10^-2 M_sun possibly in NS-NS mergers, and in the BAT band (15-150 keV) for M ~ 10^-4 M_sun possibly in NS-BH mergers. In the former case, such soft EEs can provide a good chance of ~ 6 yr^-1 for simultaneous detections of the gravitational waves with a ~ 0.1 deg angular resolution by soft X-ray survey facilities like Wide-Field MAXI.

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

We perform population synthesis simulations for Population III (Pop III) coalescing binary neutron stars (NS-NSs), neutron star - black hole binaries (NS-BHs), and binary black holes (BH-BHs) which merge within the age of the universe. We found that the typical mass of Pop III BH-BHs is $\sim 30\ \rm{M}_{\odot}$ so that the inspiral chirp signal of gravitational waves can be detected up to z=0.28 by KAGRA, Adv. LIGO, Adv. Virgo and GEO network. Our simulations suggest that the detection rate of the coalescing Pop III BH-BHs is $140 (68) \cdot ({\rm SFR}_{\rm p}/10^{-2.5} \rm{M}_{\odot} {\rm yr}^{-1} {\rm Mpc}^{-3}) \cdot {\rm Err}_{\rm sys} ~{\rm events} ~{\rm yr}^{-1}$ for the flat (Salpeter) initial mass function (IMF), respectively, where $\rm SFR_p$ and $\rm Err_{sys}$ are the peak value of the Pop III star formation rate and the possible systematic errors due to the assumptions in Pop III population synthesis, respectively. $\rm Err_{sys}=1$ correspond to conventional parameters for Pop I stars. From the observation of the chirp signal of the coalescing Pop III BH-BHs, we can determine both the mass and the redshift of the binary for the cosmological parameters determined by Planck satellite. Our simulations suggest that the cumulative redshift distribution of the coalescing Pop III BH-BHs depends almost only on the cosmological parameters. We might be able to confirm the existence of Pop III massive stars of mass $\sim 30~\rm M_{\odot}$ by the detections of gravitational waves if the merger rate of the Pop III massive BH-BHs dominates that of Pop I BH-BHs.

Long GRBs (LGRBs) have typical duration of ~ 30 s and some of them are associated with hypernovae, like Type Ic SN 1998bw. Wolf-Rayet stars are the most plausible LGRB progenitors, since the free-fall time of the envelope is consistent with the duration, and the natural outcome of the progenitor is a Type Ic SN. While a new population of ultra-long GRBs (ULGRBs), GRB 111209A, GRB 101225A, and GRB 121027A, has a duration of ~ 10^4 s, two of them are accompanied by superluminous-supernova (SLSN) like bumps, which are <~ 10 times brighter than typical hypernovae. Wolf-Rayet progenitors cannot explain ULGRBs because of too long duration and too bright SN-like bump. A blue supergiant (BSG) progenitor model, however, can explain the duration of ULGRBs. Moreover, SLSN-like bump can be attributed to the so-called cocoon-fireball photospheric emissions (CFPEs). Since a large cocoon is inevitably produced during the relativistic jet piercing though the BSG envelope, this component can be a smoking-gun evidence of BSG model for ULGRBs. In this paper, we examine u, g, r, i, and J-band light curves of three ULGRBs and demonstrate that they can be fitted quite well by our BSG model with the appropriate choices of the jet opening angle and the number density of the ambient gas. In addition, we predict that for 121027A, SLSN-like bump could have been observed for ~ 20 - 80 days after the burst. We also propose that some SLSNe might be CFPEs of off-axis ULGRBs without visible prompt emission.

Metal-poor massive stars may typically end up their lives as blue supergiants (BSGs). Gamma-ray bursts (GRBs) from such progenitors could have ultra-long duration of relativistic jets. For example Population III (Pop III) GRBs at z ~ 10-20 might be observable as X-ray rich events with a typical duration of T_90 ~ 10^4(1+z) sec. Recent GRB111209A at z = 0.677 has an ultra long duration of T_90 ~ 2.5*10^4 sec so that it have been suggested that the progenitor might be a metal-poor BSGs in the local universe. Here, we suggest luminous UV/optical/infrared emissions associated with such a new class of GRB from metal poor BSGs. Before the jet head breaks out the progenitor envelope, the energy injected by the jet is stored in a hot-plasma cocoon, which finally emerges and expands as a baryon-loaded fireball. We show that the photospheric emissions from the cocoon fireball could be intrinsically very bright (L_peak ~ 10^(42-44) erg/sec) in UV/optical bands (E_peak ~ 10 eV) with a typical duration of ~ 100 days in the rest frame. Such cocoon emissions from Pop III GRB might be detectable in infrared bands at ~ years after Pop III GRBs at up to z ~ 15 by up-coming facilities like JWST. We also suggest that GRB111209A might have been rebrightening in UV/optical bands up to an AB magnitude of < 26. The cocoon emissions from local metal-poor BSGs might have been already observed as luminous supernovae without GRB since they can be seen from the off-axis direction of the jet.

Recent numerical simulations suggest that Population III (Pop III) stars were born with masses not larger than $\sim 100 M_{\odot}$ but typically $\sim 40M_{\odot}$. By self-consistently considering the jet generation and propagation in the envelope of these low mass Pop III stars, we find that a Pop III blue super giant star has the possibility to raise a gamma-ray burst (GRB) even though it keeps a massive hydrogen envelope. We evaluate observational characters of Pop III GRBs and predict that Pop III GRBs have the duration of $\sim 10^5$ sec in the observer frame and the peak luminosity of $\sim 5 \times 10^{50} {\rm erg} {\rm sec}^{-1}$. Assuming that the $E_p-L_p$ (or $E_p-E_{\gamma, \rm iso}$) correlation holds for Pop III GRBs, we find that the spectrum peak energy falls $\sim$ a few keV (or $\sim 100$ keV) in the observer frame. We discuss the detectability of Pop III GRBs by future satellite missions such as EXIST and Lobster. If the $E_p-E_{\gamma, \rm iso}$ correlation holds, we have the possibility to detect Pop III GRBs at $z \sim 9$ as long duration X-ray rich GRBs by EXIST. On the other hand, if the $E_p-L_p$ correlation holds, we have the possibility to detect Pop III GRBs up to $z \sim 19$ as long duration X-ray flashes by Lobster.