The formation and early evolution of low mass young stellar objects (YSOs) are investigated using three-dimensional non-ideal magneto-hydrodynamics simulations. We investigate the evolution of YSOs up to ~ 10^4 yr after protostar formation, at which protostellar mass reaches ~ 0.1 M_\odot . We particularly focus on the impact of the dust model on the evolution. We found that a circumstellar disk is formed in all simulations regardless of the dust model. Disk size is approximately 10 AU at the protostar formation epoch, and it increases to several tens of AU at ~ 10^4 yr after protostar formation. Disk mass is comparable to central protostellar mass and gravitational instability develops. In the simulations with small dust size, the warp of the pseudodisk develops ~ 10^4 yr after protostar formation. The warp strengthens magnetic braking in the disk and decreases disk size. Ion-neutral drift can occur in the infalling envelope under the conditions that the typical dust size is a \gtrsim 0.2{\mu}m and the protostar (plus disk) mass is M \gtrsim 0.1 M_\odot. The outflow activity is anti-correlated to the dust size and the strong outflow appears with small dust grains.

The formation process of Population III (PopIII) stars in the mass accretion phase is investigated by numerical experiments. The barotropic relation of primordial gas and artificial stiffening of the equation of state in very dense regions ($> 10^{15}~{\rm cm}^{-3}$) enables us to follow the fragmentation of PopIII circumstellar disks and the merging processes of the fragments. The disk becomes gravitationally unstable to fragmentation , followed by a rapid merger process typically within 100 yrs, which roughly corresponds to one orbital time of the circumstellar disk. We also find that the fragmentation of the gas disk around a multiple system, a circumbinary disk, is rare; however, it is frequent in the disk around an individual protostar. We also perform a simulation with standard sink particles, where the number and total mass of sink particles are in rough agreement with those of the stiff equation of state runs. Based on the results of these numerical results, we model the evolution of the number of fragments with a simple phenomenological equation. We find that the average number of fragments is roughly proportional to $t^{0.3}$, where $t$ is the elapsed time since the formation of the first protostar. Next, we compare this trend with a number of published numerical studies by scaling the elapsed time according to the scale-free nature of the system. As a result, we find most of the results in the literature agree well with the relation. Present results combined with the previous studies in the literature imply that the PopIII stars tend to be born not as single stars, but in multiple systems.

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.

The evolution of collapsing clouds embedded in different star-forming environments is investigated using three-dimensional non-ideal magnetohydrodynamics simulations considering different cloud metallicities ($Z/\thinspace Z_\odot$ = 0, 10$^{-5}$, 10$^{-4}$, 10$^{-3}$, 10$^{-2}$, 10$^{-1}$ and 1) and ionisation strengths ($C_\zeta$=0, 0.01, 1 and 10, where $C_\zeta$ is a coefficient controlling the ionisation intensity and $C_\zeta=1$ corresponds to the ionisation strength of nearby star-forming regions). With all combinations of these considered values of $Z/\thinspace Z_\odot$ and $C_\zeta$, 28 different star-forming environments are prepared and simulated. The cloud evolution in each environment is calculated until the central density reaches $n\approx10^{16}\,{\rm cm}^{-3}$ just before protostar formation, and the outflow driving conditions are derived. An outflow appears when the (first) adiabatic core forms in a magnetically active region where the magnetic field is well coupled with the neutral gas. In cases where outflows are driven, their momentum fluxes are always comparable to the observations of nearby star-forming regions. Thus, these outflows should control the mass growth of the protostars as in the local universe. Roughly, an outflow appears when $Z/\thinspace Z_\odot>10^{-4}$ and $C_\zeta \ge 0.01$. It is expected that the transition of the star formation mode from massive stars to normal solar-type stars occurs when the cloud metallicity is enhanced to the range of $Z/\thinspace Z_\odot\approx 10^{-4}$--$10^{-3}$, above which relatively low-mass stars would preferentially appear as a result of strong mass ejection.

In nearby star-forming clouds, amplification and dissipation of the magnetic field are known to play crucial roles in the star-formation process. The star-forming environment varies from place to place and era to era in galaxies. In the present study, amplification and dissipation of magnetic fields in star-forming clouds are investigated under different environments using magnetohydrodynamics (MHD) simulations. We consider various star-forming environments in combination with the metallicity and the ionization strength, and prepare prestellar clouds having two different mass-to-flux ratios. We calculate the cloud collapse until protostar formation using ideal and non-ideal (inclusion and exclusion of Ohmic dissipation and ambipolar diffusion) MHD calculations to investigate the evolution of the magnetic field. We perform 288 runs in total and show the diversity of the density range within which the magnetic field effectively dissipates, depending on the environment. In addition, the dominant dissipation process (Ohmic dissipation or ambipolar diffusion) is shown to strongly depend on the star-forming environment. Especially, for the primordial case, magnetic field rarely dissipates without ionization source, while it efficiently dissipates when very weak ionization sources exist in the surrounding environment. The results of the present study help to clarify star formation in various environments.

Metal enrichment by the first-generation (Pop III) stars is the very first step of the matter cycle in the structure formation and it is followed by the formation of extremely metal-poor (EMP) stars. To investigate the enrichment process by the Pop III stars, we carry out a series of numerical simulations including the feedback effects of photoionization and supernovae (SNe) of Pop III stars with a range of masses of minihaloes (MHs), M_halo , and Pop III stars, M_PopIII . We find that the metal-rich ejecta reaches neighbouring haloes and external enrichment (EE) occurs when the halo binding energy is sufficiently below the SN explosion energy, E_SN . The neighbouring haloes are only superficially enriched, and the metallicity of the clouds is [Fe/H] < -5. Otherwise, the SN ejecta falls back and recollapses to form enriched cloud, i.e. internal enrichment (IE) process takes place. In case that a Pop III star explodes as a core-collapse SNe (CCSNe), MHs undergo IE, and the metallicity in the recollapsing region is -5 < [Fe/H] < -3 in most cases. We conclude that IE from a single CCSN can explain the formation of EMP stars. For pair-instability SNe (PISNe), EE takes place for all relevant mass range of MHs, consistent with no observational sign of PISNe among EMP stars.

Low-mass population III (PopIII) stars of $\lesssim 0.8 M_{\odot}$ could survive up until the present. Non-detection of low-mass PopIII stars in our Galaxy has already put a stringent constraint on the initial mass function (IMF) of PopIII stars, suggesting that PopIII stars have a top-heavy IMF. On the other hand, some claims that the lack of such stars stems from metal enrichment of their surface by accretion of heavy elements from interstellar medium (ISM). We investigate effects of the stellar wind on the metal accretion onto low-mass PopIII stars because accretion of the local ISM onto the Sun is prevented by the solar wind even for neutrals. The stellar wind and radiation of low-mass PopIII stars are modeled based on knowledge of nearby low-mass stellar systems including our Sun. We find that low-mass PopIII stars traveling across the Galaxy forms the stellar magnetosphere in most of their life. Once the magnetosphere is formed, most of neutral interstellar particles are photoionized before reaching to the stellar surface and are blown away by the wind. Especially, the accretion abundance of iron will be reduced by a factor of $< 10^{-12}$ compared with Bondi-Hoyle-Lyttleton accretion. The metal accretion can enhance iron abundance [Fe/H] only up to $\sim -14$. This demonstrates that low-mass PopIII stars remain pristine and will be found as metal free stars and that further searches for them are valuable to constrain the IMF of PopIII stars.

We study the number and the distribution of low mass Pop III stars in the Milky Way. In our numerical model, hierarchical formation of dark matter minihalos and Milky Way sized halos are followed by a high resolution cosmological simulation. We model the Pop III formation in H2 cooling minihalos without metal under UV radiation of the Lyman-Werner bands. Assuming a Kroupa IMF from 0.15 to 1.0 Msun for low mass Pop III stars, as a working hypothesis, we try to constrain the theoretical models in reverse by current and future observations. We find that the survivors tend to concentrate on the center of halo and subhalos. We also evaluate the observability of Pop III survivors in the Milky Way and dwarf galaxies, and constraints on the number of Pop III survivors per minihalo. The higher latitude fields require lower sample sizes because of the high number density of stars in the galactic disk, the required sample sizes are comparable in the high and middle latitude fields by photometrically selecting low metallicity stars with optimized narrow band filters, and the required number of dwarf galaxies to find one Pop III survivor is less than ten at <100 kpc for the tip of redgiant stars. Provided that available observations have not detected any survivors, the formation models of low mass Pop III stars with more than ten stars per minihalo are already excluded. Furthermore, we discuss the way to constrain the IMF of Pop III star at a high mass range of > 10 Msun.

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

We perform cosmological hydrodynamics simulations with non-equilibrium primordial chemistry to obtain 59 minihalos that host first stars. The obtained minihalos are used as initial conditions of local three dimensional radiation hydrodynamics simulations to investigate the formation of the first stars. We find two-thirds of the minihalos host multiple stars, while the rest of them have single stars. The mass of the stars found in our simulations are in the range of 1 Msun \la M \la 300 Msun, peaking at several x 10 Msun. Most of the very massive stars of \ga 140 Msun are born as single stars, although not all of the single stars are very massive. We also find a few stars of \la 1 Msun that are kicked by the gravitational three body interactions to the position distant from the center of mass. The frequency that a star forming minihalo contains a binary system is \sim 50%. We also investigate the abundance pattern of the stellar remnants by summing up the contributions from the first stars in the simulations. Consequently, the pattern is compatible with that of the low metallicity Damped Lyman-alpha systems or the Extremely Metal Poor (EMP) stars, if the mass spectrum obtained in our experiment is shifted to the low mass side by 0.2 dex. If we consider the case that an EMP star is born in the remnant of the individual minihalo without mixing with others, the chemical signature of the pair instability supernova is more prominent, because most of them are born as single stars.

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.

We perform a three dimensional radiation hydrodynamics simulation to investigate the formation of first stars from initial collapse of a primordial gas cloud to formation and growth of protostars. The simulation is integrated until 0.1 Myrs after the formation of the primary protostar by which the protostars have already settled onto main sequence stars. This is the first attempt of simulating first star formation to take into account the ultraviolet radiative feedback effect by the multiple protostars as well as the three dimensional effects such as fragmentation of the accretion disk. We find that the mass accretions onto the population III protostars are significantly suppressed by the radiative feedback from themselves. As a result, we find five stars formed in this particular simulation, and that the final mass of the stars are < 60Msun, including a star of 4.4Msun. Formation of such a star hints at the existence of even lower-mass stars that would live today.

The disk wind, which is powered by the radiation force due to spectral lines (line force), is studied for broad absorption line (BAL) quasars. We investigate the structure of the disk wind based on the non-hydrodynamic method and compare with wind properties inferred from X-ray observations of BAL quasars. In this paper, we apply the stellar wind theory to the initial condition (the mass outflow rate at the base of the wind). We found the funnel-shaped winds with a half opening angle of 50^{circ} for the case of epsilon=0.3-0.9 and M_{BH}=10^{7-8.5}M_odot, where epsilon is the Eddington ratio and M_{BH} is the black hole mass. Thus, the absorption features are observed for an observer of which a viewing angle is around 50^{circ}. A probability of BAL quasars is 7-11%, which is roughly consistent the abundance ratio of BAL quasars, 10-15%. Here, the probability is estimated by the solid angle, that the absorbing features would be detected, divided by 4pi. In contrast, if the Eddington ratio is smaller than 0.01 or if the black hole is very massive, M_{BH} < 10^9M_{odot}, the disk wind is not launched due to the less effective line force. Then, the quasars are identified as non-BAL quasars independently of the observer's viewing angle.

First, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter haloes of ~10^4M_sun can produce first luminous objects with the aid of dark matter cusps. Therefore, the mass of first objects is smaller by roughly two orders of magnitude than in the previous prediction. This implies that the number of Pop III stars formed in the early universe could be significantly larger than hitherto thought. Secondly, the feedback by photo-ionization and photo-dissociation photons in the first objects is explored with radiation hydrodynamic simulations, and it is demonstrated that multiple stars can form in a 10^5M_sun halo. Thirdly, the fragmentation of an accretion disk around a primordial protostar is explored with photo-dissociation feedback. As a result, it is found that the photo-dissociation can reduce the mass accretion rate onto protostars. Also, protostars as small as 0.8M_sun may be ejected and evolve with keeping their mass, which might be detected as "real first stars" in the Galactic halo. Finally, state-of-the-art radiation hydrodynamic simulations are performed to investigate the internal ionization of first galaxies and the escape of ionizing photons. We find that UV feedback by forming massive stars enhances the escape fraction even in a halo as massive as > 6* 10^9M_sun, while it reduces the star formation rate significantly. This may have a momentous impact on the cosmic reionization.

We investigate generation processes of magnetic fields around first stars. Since the first stars are expected to form anisotropic ionization fronts in the surrounding clumpy media, magnetic fields are generated by effects of radiation force as well as the Biermann battery effect. We have calculated the amplitude of magnetic field generated by the effects of radiation force around the first stars in the preceding paper, in which the Biermann battery effects are not taken into account.In this paper, we calculate the generation of magnetic fields by the Biermann battery effect as well as the effects of radiation force, utilizing the radiation hydrodynamics simulations. As a result, we find that the generated magnetic field strengths are ~ 10^{-19}G-10^{-17}G at ~ 100pc-1kpc scale mainly by the Biermann battery, which is an order of magnitude larger than the results of our previous study. We also find that this result is insensitive to various physical parameters including the mass of the source star, distance between the source and the dense clump, unless we take unlikely values of these parameters.

We investigate seed magnetic field generation in the early universe by radiation force of first stars. In the previous study with steady assumption, large amplitudes(10^{-15} G for first stars, 10^{-11}G for QSOs) are predicted. In this study, we formulate this issue in unsteady framework. Then, we consider a specific model of magnetic field generation around a very massive first star. Consequently, 1) we find steady assumption is not valid in realistic situation, and 2) obtain much smaller magnetic field strength than predicted by Langer et al. (2003). In addition, we find momentum transfer process during photoionization is more important than Thomson scattering. The resultant magnetic flux density around the first star is < 10^{-19}G. These seed magnetic field will not affect subsequent star formation in the neighbor of first stars.

We explore the formation of secondary Population III (Pop III) stars under radiation hydrodynamic (RHD) feedback by a preformed massive star. To properly treat RHD feedback, we perform three-dimensional RHD simulations incorporating the radiative transfer of ionizing photons as well as H_2 dissociating photons from a preformed star. A collapsing gas cloud is settled at a given distance from a 120Msun Pop III star, and the evolution of the cloud is pursued including RHD feedback. We derive the threshold density depending on the distance, above which the cloud can keep collapsing owing to the shielding of H_2 dissociating radiation. We find that an H_2 shell formed ahead of an ionizing front works effectively to shield the H_2 dissociating radiation, leading to the positive feedback for the secondary Pop III star formation. Also, near the threshold density, the envelope of gas cloud is stripped significantly by a shock associated with an ionizing front. By comparing the mass accretion timescale with the Kelvin-Helmholtz timescale, we estimate the mass of secondary Pop III stars. It turns out that the stripping by a shock can reduce the mass of secondary Pop III stars down to \approx 20Msum.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

The first-generation stars in the $\Lambda$CDM universe are considered to have formed in dark halos with total masses in the range $\sim 10^{5}-10^{7}M_sun$ at $z \sim 20-50$. These stars expected to be very massive and in some cases they end their lives as the first supernovae (SNe). We explore the problem of whether star formation in low mass dark halos (< 10^{7} M_sun) was triggered or suppressed by the SN feedback from neighboring star-forming halos. We take into consideration mainly two effects by the SN shock: one is the evacuation of gas components from the halos and the other is the promotion of H_2 formation because of the enhanced ionization degree by shock heating. Combining above effects, we find that the star formation activities in the neighboring dark matter halos (M < 10^{7} M_sun) are basically suppressed in case they are located close to the SN center, because of the gas evacuation effect. The critical distance within which the gas is blown away falls within the range $\sim 0.3-1.5$kpc depending on the SN energy and the halo mass. In addition, we find there is very little window in the parameter space where star formation activities in dark halos are induced or promoted by neighboring SN.

We explore the impact of ultraviolet (UV) radiation from massive Population III (Pop III) stars of 25, 40, 80, and 120 M_sun on the subsequent Pop III star formation. In this paper, particular attention is paid to the dependence of radiative feedback on the mass of source Pop III star. UV radiation from the source star can work to impede the secondary star formation through the photoheating and photodissociation processes. Recently, Susa & Umemura (2006) have shown that the ionizing radiation alleviates the negative effect by H_2-dissociating radiation from 120$M_sun PopIII star, since an H_2 shell formed ahead of an ionizing front can effectively shield H_2-dissociating radiation. On the other hand, it is expected that the negative feedback by H_2-dissociating radiation can be predominant if a source star is less massive, since a ratio of the H_2-dissociating photon number to the ionizing photon number becomes higher. In order to investigate the radiative feedback effects from such less massive stars, we perform three-dimensional radiation hydrodynamic simulations, incorporating the radiative transfer effect of ionizing and H_2-dissociating radiation. As a result, we find that if a source star is less massive than ~25M_sun, the ionizing radiation cannot suppress the negative feedback of H_2-dissociating radiation. Therefore, the fate of the neighboring clouds around such less massive stars is determined solely by the flux of H_2-dissociating radiation from source stars. With making analytic estimates of H_2 shell formation and its shielding effect, we derive the criteria for radiation hydrodynamic feedback depending on the source star mass.

We perform radiation hydrodynamics simulations on the evolution of galactic gas disks irradiated by ultraviolet radiation background. We find gas disks with N_H > 10^21 cm^-2 exposed to ultraviolet radiation at a level of I_21=1 can be self-shielded from photoheating, whereas the disk with N_H < 10^21 cm^-2 cannot. We also find that the unshielded disks keep smooth density distribution without any sign of fragmentation, while the self-shielded disks easily fragment into small pieces by self-gravity, possibly followed by star formation. The suppression of star formation in unshielded disks is different from photoevaporation effect, since the assumed dark halo potential is deep enough to keep the photoheated gas. Presence of such critical threshold column density would be one of the reason for the so-called down-sizing feature of present-day galaxies.

We investigate the stability of black hole accretion disks in a primordial environment (POP III disks for short), by solving the vertical structure of optically thick disks, including convective energy transport, and by employing a one-zone model for optically thin isothermal disks. Because of the absence of metals in POP III disks, we find significant differences in stability associated with ionization between POP III disks and the disks of solar metallicity. An unstable branch in S-shaped equilibrium curves on the Mdot-Sigma (mass accretion rate - surface density) plane extends to a larger surface density compared with the case of disks of solar metallicity. The resulting equilibrium loci indicate that quasi-periodic oscillations in luminosity can also be driven in POP III disks, and their maximal luminosity is typically by an order of magnitude larger than that of the disks of solar metallicity. Such a strong outburst of POP III disks can be observed by future huge telescopes, in case that the mass is supplied onto the disks at the Bondi accretion rates in typical virialized small dark halos.

We investigate the dissipation of magnetic flux in primordial star-forming clouds throughout their collapse including the run-away collapse phase as well as the accretion phase. We solve the energy equation and the non-equilibrium chemical reactions in the collapsing gas, in order to obtain the radial distribution of the ionized fraction during the collapse. As a result, we find the ionized fraction is high enough for the magnetic field to couple with the gas throughout the evolution of the cloud. This result suggests that the jet formation from protostars as well as the activation of magneto-rotational instability in the accretion disk are enabled in the presence of the cosmological seed magnetic flux proposed by Langer et al.(2003).

We investigate the star formation process in primordial environment in the presence of radiative feedback by other population III stars formed earlier. In this paper, we focus our attention on the effects by photodissociative radiation toward the full understanding of the radiative feedback effects. We perform three dimensional radiation hydrodynamics simulations on this issue as well as analytic estimates, paying special attention on the self-shielding effect and dynamics of the star-forming cloud. As a result, we find that the ignition timing of the source star is crucial. If the ignition is later than the epoch when the central density of the collapsing cloud exceeds $\sim 10^3-10^4{\rm cm^{-3}}$, the collapse cannot be reverted, even if the source star is located at $\la$ 100pc. The uncertainty of the critical density comes from the variety of initial conditions of the collapsing cloud. We also find the analytic criteria for a cloud to collapse with given central density, temperature and the Lyman-Werner(LW) band flux which irradiates the cloud. Although we focus on the radiation from neighbor stars, this result can also be applied to the effects of diffuse LW radiation field, that is expected to be built up prior to the reionization of the universe. We find that self-gravitating clouds can easily self-shield from diffuse LW radiation and continue their collapse for densities larger than $\sim 10^3 {\rm cm^{-3}}$.

We explore the possibility of subsequent star formation after a first star forms in a Pop III object, by focusing on the radiation hydrodynamic (RHD) feedback brought by ionizing photons as well as H2 dissociating photons. For the purpose, we perform three-dimensional RHD simulations, where the radiative transfer of ionizing photons and H2 dissociating photons from a first star is self-consistently coupled with hydrodynamics based on a smoothed particle hydrodynamics method. As a result, it is shown that density peaks above a threshold density can keep collapsing owing to the shielding of H2 dissociating radiation by an H2 shell formed ahead of a D-type ionization front. But, below the threshold density, an M-type ionization front accompanied by a shock propagates, and density peaks are radiation hydrodynamically evaporated by the shock. The threshold density is dependent on the distance from a source star, which is $\approx 10^2 cm^{-3}$ for the source distance of 30pc. Taking into consideration that the extent of a Pop III object is $\approx 100$pc and density peaks within it have the density of $10^{2-4}$cm$^{-3}$, it is concluded that the secondary star formation is allowed in the broad regions in a Pop III object.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

We have been constructed a brand-new radiation hydrodynamics solver based upon Smoothed Particle Hydrodynamics (SPH), which works on parallel computer system. The code is designed to investigate the formation and evolution of the first generation objects at $z \gtrsim 10$, where the radiative feedback from various sources play important roles. The code can compute the fraction of chemical species e, H$^+$, H, H$^-$, H$_2$, and H$_2^+$ by fully implicit time integration. It also can deal with multiple sources of ionizing radiation, as well as the radiation at Lyman-Werner band. We compare the results for a few test calculations with the results of one dimensional simulations, in which we find good agreements with each other. We also evaluate the speedup by parallelization, that is found to be almost ideal, as far as the number of sources is comparable to the number of processors.

Recent observations on the cosmic microwave background by Wilkinson Microwave Anisotropy Probe(WMAP) strongly suggest that the reionization of the universe took place quite early (z \sim 17). On the other hand, it has been pointed out that the cold dark matter cosmology suffers from the substructure problem that subgalactic halos are overproduced than the observed dwarf galaxies in the Local Group. In this paper, as a potential mechanism to solve this problem, we consider the feedback effects of the early reionization on the formation of small-scale structures. For this purpose, we perform 3D radiation hydrodynamic simulations with incorporating the radiative transfer for ionizing photons. As a result, it is found that the early reionization is so devastating for low-mass systems with M_{vir}\la 10^8 M_\odot or v_{circ}\la 20km/s, and almost all gas is photo-evaporated in more than 95% of low-mass systems. Such a strong negative feedback on the formation of low-mass galaxies may solve the substructure problem and support the picture that Local Group dwarf galaxies are descendants of the more massive halos that experienced and survived tidal stripping.

We study the formation and evolution of HII regions around the first stars formed at z=10-30. We use a one-dimensional Lagrangian hydrodynamics code which self-consistently incorporates radiative transfer and non-equilibrium primordial gas chemistry. The star-forming region is defined as a spherical gas cloud with a Population III star embedded at the center. We explore a large parameter space of host halo mass, gas density profile, and stellar luminosity. The formation of the HII region is characterized by initial slow expansion of a weak D-type ionization front near the center, followed by rapid propagation of an R-type front throughout the outer gas envelope. We find that the transition between the two front types is indeed a critical condition for the complete ionization of halos of cosmological interest. In small mass (< 10^6 M_sun) halos, the transition takes place within a few 10^5 yr, yielding high escape fractions (>80%) of both ionizing and photodissociating photons. In larger mass (> 10^7 M_sun) halos, the ionization front remains to be of D-type over the lifetime of the massive star, the HII region is confined well inside the virial radius, and the escape fractions are essentially zero. We derive an analytic formula that reproduces well the results of our simulations. We discuss immediate implications of the present results for the star formation history and early reionization of the Universe.

We report the strength of seed magnetic flux of accretion disk surrounding the PopIII stars. The magnetic field in accretion disk might play an important role in the transport of angular momentum because of the turbulence induced by Magneto-Rotational Instability (MRI). On the other hand, since the primordial star-forming clouds contain no heavy elements and grains, they experience much different thermal history and different dissipation history of the magnetic field in the course of their gravitational contraction, from those in the present-day star-forming molecular clouds. In order to assess the magnetic field strength in the accretion disk of PopIII stars, we calculate the thermal history of the primordial collapsing clouds, and investigate the coupling of magnetic field with primordial gas. As a result, we find that the magnetic field strongly couple with primordial gas cloud throughout the collapse, i.e. the magnetic field are frozen to the gas as far as initial field strength satisfies $ B\lesssim 10^{-5}(n_{\rm H}/10^3{\rm cm^{-3}})^{0.55}{\rm G}$.

Recently, ``CDM crisis'' is under discussion. The main point of this crisis is that number of substructures presented by cosmological N-body simulations based on CDM scenario for structure formation is much larger than observed substructures. Therefore, it is crucial for this crisis to discriminate whether expected number of CDM substructures really exist but non-luminous or do not exist. In this paper, we present a new idea to detect such invisible substructures by utilizing a gravitational lensing. Here, we consider quasars that are gravitationally lensed by a foreground galaxy. A substructure around the lensing galaxy may superposed on one of the lensed images of such quasars. In this situation, additional image splitting should occur in the image behind the substructure, and further multiple images are created. This is ``quasar mesolensing''. We estimate separation and time delay between further multiple images due to quasar mesolensing. The expected value is $1 \sim 30$ milli-arcsecond for the separation and future fine resolution imaging enable us to find invisible substructures, and is $1 \sim 10^3$ second for the time delay and high-speed monitoring of such quasar will be able to find ``echo''-like variation due to quasar mesolensing in intrinsic variability of the quasar. Furthermore, we evaluate that the optical depth for the quasar mesolensing is $\sim 0.1$. Consequently, if we monitor a few multiple quasars, we can find ``echo''-like variation in one of the images after intrinsic flux variations of quasars.

We reanalyze the photoevaporation problem of subgalactic objects irradiated by UVB radiation in a reionized universe. For the purpose, we perform three-dimensional radiation smoothed-particle-hydrodynamics calculations, where the radiative transfer is solved by a direct method and also the non-equilibrium chemistry of primordial gas including H$_2$ molecules is incorporated. Attention is concentrated on radiative transfer effects against UVB for the formation of subgalactic objects with $T_{\rm vir}\la 10^4$K. We consider the reionization model with $z_{\rm reion} \approx 7$ and also the earlier reionization model ($z_{\rm reion} \approx 17$) inferred by the WMAP. We find that the star formation is suppressed appreciably by UVB, but baryons at high-density peaks are self-shielded even during the reionization, forming some amount of stars eventually. In that sense, the photoevaporation for subgalactic systems is not so perfect as argued by one-dimensional spherical calculations. The final stellar fraction depends on the collapse epoch and the mass of system, but almost regardless of the reionization epoch. For instance, a few tenths of formed stars are born after the cosmic reionization in $z_{\rm reion} \approx 7$ model, while more than 90% stars are born after the reionization in the WMAP reionization model. Thus, effects of UVB feedback on the substructure problem with CDM scenario should be evaluated with careful treatment of the radiative transfer. The star clusters formed at high-density peaks coalesce with each other in a dissipationless fashion in a dark matter potential, resultantly forming a spheroidal system. As a result, these low-mass galaxies have large mass-to-light ratios such as observed in dwarf spheroidals in the Local Group.

Full radiative transfer in 3D space has been solved to pursue the reionization history in an inhomogeneous universe. It has been shown that the reionization of an inhomogeneous universe is not a prompt event, but a fairly slow process. Also, QSO absorption line systems are simulated with using the results of radiative transfer calculations. Ly$\alpha$ continuum depression implies that the metagalactic UV intensity decreases rapidly with $z$ at $z>5$, and the reionization must have taken place between $z=6$ and 10. Finally, it is stressed that H$\alpha$ forest is a more powerful tool to probe the reionization history and the density fluctuations in the universe at $z>5$.

Using 3D radiative transfer calculations on the reionization of an inhomogeneous universe, QSO absorption line systems are simulated and they are compared with observations of Ly$\alpha$ continuum depression at high redshifts. By this comparison, it is found that the metagalactic UV intensity decreases rapidly with $z$ at $z>4$ as $I_{21}=0.5\exp [3(4-z)]$, and the reionization must have taken place between $z=6$ and 10. Based on this time-dependence of UV background intensity, we explore the collapse of pregalactic clouds in the UV background, and find that the self-shielding is prominent above a mass scale as $ M_{BIF}=3.0\times 10^{11} M_\odot [(1+z_c)/5]^{-4.2}(I_{21}/0.5)^{0.6}. $ This mass scale predicts the bifurcation of galactic morphology, and by confrontation with observations it turns out that the bifurcation mass successfully discriminates between elliptical and spiral galaxies.

We explore possibilities of collapse and star formation in Population III objects exposed to the external ultraviolet background (UVB) radiation. Assuming spherical symmetry, we solve self-consistently radiative transfer of photons, non-equilibrium H2 chemistry, and gas hydrodynamics. Although the UVB does suppress the formation of low mass objects, the negative feedback turns out to be weaker than previously suggested. In particular, the cut-off scale of collapse drops significantly below the virial temperature 10^4 K at weak UV intensities, due to both self-shielding of the gas and H2 cooling. Clouds above this cut-off tend to contract highly dynamically, further promoting self-shielding and H2 formation. For plausible radiation intensities and spectra, the collapsing gas can cool efficiently to temperatures well below 10^4 K before rotationally supported and the final H2 fraction reaches 10^{-3}. Our results imply that star formation can take place in low mass objects collapsing in the UVB. The threshold baryon mass for star formation is \sim 10^9 solar mass for clouds collapsing at redshifts z \simlt 3, but drops significantly at higher redshifts. In a conventional cold dark matter universe, the latter coincides roughly with that of the 1 \sigma density fluctuations. Objects near and above this threshold can thus constitute `building blocks' of luminous structures, and we discuss their links to dwarf spheroidal/elliptical galaxies and faint blue objects. These results suggest that the UVB can play a key role in regulating the star formation history of the Universe.

Based upon a novel paradigm of the galaxy formation under UV background, the evolutionary bifurcation of pregalactic clouds is confronted with observations on elliptical and spiral galaxies. The theory predicts that the dichotomy between the dissipational and dissipationless galaxy formation stems from the degree of self-shielding from the UV background radiation. This is demonstrated on a bifurcation diagram of collapse epochs versus masses of pregalactic clouds. Using the observed properties, the collapse epochs are assessed for each type of galaxies with attentive mass estimation. By the direct comparison of the theory with the observations, it turns out that the theoretical bifurcation branch successfully discriminates between elliptical and spiral galaxies. This suggests that the UV background radiation could play a profound role for the differentiation of the galactic morphology into the Hubble sequence.

The collapse of marginally Jeans unstable primordial gas clouds in the presence of UV radiation field is discussed. Assuming that the dynamical collapse proceeds approximately in an isothermal self-similar fashion, we investigate the thermal evolution of collapsing central core until H$_2$ cooling dominates photoheating and the temperature drops to below $10^4$K. Consequently, the mass of the cooled core is evaluated as $M_{cool}=3.6\times 10^6 M_\odot(I_{21}/1)^{-0.32}$. This scale depends only on the incident UV intensity, and provides a lower limit to the mass of collapsed objects in the UV radiation field.

To elucidate the effects of the UV background radiation on the collapse of pregalactic clouds, we implement a radiation-hydrodynamical calculation, combining one-dimensional spherical hydrodynamics with an accurate treatment of the radiative transfer of ionizing photons. Both absorption and scattering of UV photons are explicitly taken into account. It turns out that a gas cloud contracting within the dark matter potential does not settle into hydrostatic equilibrium, but undergoes run-away collapse even under the presence of the external UV field. The cloud center is shown to become self-shielded against ionizing photons by radiative transfer effects before shrinking to the rotation barrier. Based on our simulation results, we further discuss the possibility of H2 cooling and subsequent star formation in a run-away collapsing core. The present results are closely relevant to the survival of subgalactic Population III objects as well as to metal injection into intergalactic space.

The pancake collapse of pregalactic clouds under UV background radiation is explored with a one-dimensional sheet model. Here, attention is concentrated on elucidating the basic physics on the thermal evolution of pregalactic clouds exposed to diffuse UV radiation. So, we treat accurately the radiation transfer for the ionizing photons, with solving chemical reactions regarding hydrogen molecules as well as atoms. The self-shielding against UV radiation by H$_2$ Lyman-Werner bands, which regulates the photo-dissociation of hydrogen molecules, is also taken into account. As a result, it is found that when the UV background radiation is at a level of $10^{-22} (\nu/\nu_L)^{-1} {\rm erg~ s^{-1} cm^{-2} Hz^{-1} str^{-1}}$, the cloud evolution bifurcates with a critical mass as $M_{\rm SB} = 2.2\times 10^{11} M_\odot [(1+z_c)/5]^{-4.2}$, where $z_c$ is the final collapse epoch. A cloud more massive than $M_{\rm SB}$ cools below $5\times 10^3$K due to H$_2$ line emission at the pancake collapse and would undergo the initial starburst. The pancake possibly evolves into a virialized system in a dissipationless fashion. Consequently, this leads to the dissipationless galaxy formation at $3\lsim z_c \lsim 10$. A cloud less massive than $M_{\rm SB}$ cannot cool by H$_2$ emission shortly after the pancake collapse, but could cool in the course of shrinking to the rotation barrier. This is likely to lead to the dissipational galaxy formation at relatively low redshifts as $0\lsim z_c \lsim 4$. The present results provide a solid physical mechanism which controls the star formation efficiency in the pregalactic clouds. In the context of a standard CDM cosmology, $M_{\rm SB}$ lies between 1$\sigma$ and 2$\sigma$ density fluctuations.

We investigate the evolution of cosmological low mass (low virial temperature) objects and the formation of the first luminous objects. First, the `cooling diagram' for low mass objects is shown. We assess the cooling rate taking into account the contribution of H_2, which is not in chemical equilibrium generally, with a simple argument of time scales. The reaction rates and the cooling rate of H_2 are taken from the recent results by Galli & Palla (1998). Using this cooling diagram, we also estimate the formation condition of luminous objects taking into account the supernova (SN) disruption of virialized clouds. We find that the mass of the first luminous object is several times 10^7 solar mass, because smaller objects may be disrupted by the SNe before they become luminous. Metal pollution of low mass (Ly-alpha) clouds also discussed. The resultant metallicity of the clouds is about 1/1000 of the solar metallicity.

We investigate the thermal and dynamical evolution of primordial gas clouds in the universe after decoupling. Comparing the time-scale of dynamical evolution with that of fragmentation, we can estimate the typical fragmentation scale. We propose the following scenario of the formation process of first luminous objects consisting of large number stars. First, by pancake collapse of the overdensity regions in the expanding universe or collision between clouds in potential wells, quasi-plane shocks form. If the shock-heated temperature is higher than about $10^4$ K, the postshock gas cools down to several hundred K by \Hbun line cooling, and the shock-compressed layer fragments into filamentary clouds. The filamentary cloud collapses dynamically once more and fragments into cloud cores. Finally, a primordial star forms in a cloud core. We show that the minimum mass of the first star is essentially determined by the Chandrasekhar mass. Also, we investigate the dynamical collapse of cloud cores by numerical simulation and show that the evolution paths of the central regions of the cores depend only very weakly on the total core mass. After mass accretion, a massive star may be formed in a core, since the estimated mass accretion rate is very large. In such a case, it may be possible for many massive stars form almost simultaneously. Then the clouds can be luminous objects. On the other hand, if the shock-heated temperature is lower, effective star formation is delayed significantly.

The contraction processes of metal-free molecular clouds of starlike mass (or cloud cores) are investigated. We calculate radiative transfer of the H_2 lines and examine quasi-static contraction with radiative cooling. Comparing two time-scales, the free-fall time t_ff and the time-scale of quasi-static contraction t_qsc (nearly equal to t_cool, the cooling time) of these cores, we find that the ratio of the two time-scales t_ff/t_qsc, i.e., the efficiency of cooling, becomes larger with contraction even under the existence of cold and opaque envelopes. In particular, for fragments of primordial filamentary clouds, for which t_ff is nearly equal to t_qsc at the fragmentation epoch, they collapse dynamically in the free-fall time-scale. This efficiency of cooling is unique to line cooling.

We re-examine the thermal evolution of the postshock layer in primordial gas clouds. Comparing the time scales, we find that the evolutionary paths of postshock regions in primordial gas clouds can be basically understood in terms of the diagram drawn in the ionization degree vs temperature plane. The results obtained from the diagram are independent of the density in the case that we do not include photodissociation and photoionization. We also argue that the diagram is not only relevant to the case of the steady postshock flow, but also to the isochorically cooling gas.

We discuss the fragmentation of primordial gas clouds in the universe after decoupling. Comparing the time scale of collapse with that of fragmentation, we obtain the typical mass of a fragment both numerically and analytically. It is shown that the estimated mass gives the minimum mass of a fragment which is formed from the primordial gas cloud and is essentially determined by the Chandrasekhar mass.

We investigated the thermal evolution of non-spherical primordial gas clouds of mass $ M=10^6M_\odot $. We studied two cases: 1) prolate and oblate clouds without angular momentum, 2) initially spherical, oblate and prolate clouds with angular momentum. In the spherical case, the magnitude of the angular momentum is the key quantity which determines the fragment mass. The fragment mass is found to be $M \gsim 100 M_\odot$ for $10^{-2} \lsim \lambda \lsim 1$, where $\lambda$ is the cosmological spin parameter. For the oblate shape of initial gas clouds, the angular momentum is almost always not important, In the case of prolate clouds with some angular momentum, collapse proceeds in the same way as spherical case, and the bounce occurs to form a disk due to the angular momentum. In any case, the estimated ``fragment mass'', must be above $100 M_\odot$ for the primordial gas clouds with typical angular momentum and typical oblateness or prolateness.

We consider the joint probability distribution function for the mass contrast and angular momentum of over-density regions on the proto- galactic scale and investigate the formation of massive black holes at redshift $z\gsim10$. We estimate the growth rate of the angular momentum by the linear perturbation theory and the decay rate by the Compton drag and apply the Press-Schechter theory to obtain the formation rate of massive black holes, assuming the full reionization of the universe at $z=z_{ion}\gg 10$. We find the correlation between the mass contrast and angular momentum vanishes in the linear theory. However, application of the Press-Schechter theory introduces a correlation between the mass contrast and angular momentum of bound objects. Using thus obtained probability distribution, we calculate the mass fraction of black holes with $M\sim10^6-10^8M_{\odot}$ in the universe. We find that it crucially depends on the reionization epoch $z_{ion}$. Specifically, for the standard CDM power spectrum with the COBE normalization, the condition $z_{ion}\gsim 500$ must be satisfied to reproduce the observed number density of QSOs.