We explore the induction of circular polarization (CP) of near-infrared light in star-forming regions using three-dimensional radiative transfer calculations. The simulations trace the change of Stokes parameters at each scattering/absorption process in a dusty gas slab composed of aligned grains. We find that the CP degree enlarges significantly according as the size of dust grains increases and exceeds $\sim 20$ percent for micron-size grains. Therefore, if micron-size grains are dominant in a dusty gas slab, the high CP observed around luminous young stellar objects can be accounted for. The distributions of CP show the asymmetric quadrupole patters regardless of the grain sizes. Also, we find that the CP degree depends on the relative position of a dusty gas slab. If a dusty gas slab is located behind a star-forming region, the CP reaches $\sim 60$ percent in the case of $1.0~{\rm \mu m}$ dust grains. Hence, we suggest that the observed variety of CP maps can be explained by different size distributions of dust grains and the configuration of aligned grains.

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.

ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$\alpha$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$\alpha$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.

Recent observations have successfully detected [O III] $88.3\,{\rm \mu m}$ and [C II] $157.6\,{\rm \mu m}$ lines from galaxies in the early Universe with the Atacama Large Millimeter Array (ALMA). Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present relations between the metal line emission and galaxy evolution at $z=6-15$. We find that galaxies during their starburst phases have high [O III] luminosity of $\sim 10^{42}~\rm erg~s^{-1}$. Once supernova feedback quenches star formation, [O III] luminosities rapidly decrease and continue to be zero for $\sim 100\,{\rm Myr}$. The slope of the relation between $\log{(\rm SFR/M_{\odot}~yr^{-1})}$ and $\log{(L_{\rm [O_{III}]}/{\rm L_{\odot}})}$ at $z=6-9$ is 1.03, and 1.43 for $\log{(L_{\rm [C_{II}]}/{\rm L_{\odot}})}$. As gas metallicity increases from sub-solar to solar metallicity by metal enrichment from star formation and feedback, the line luminosity ratio $L_{\rm [O_{III}]} / L_{\rm [C_{II}]}$ decreases from $\sim 10$ to $\sim 1$ because the O/C abundance ratio decreases due to carbon-rich winds from AGB stars. Therefore, we suggest that the combination of [O III] and [C II] lines is a good probe to investigate the relative distribution of ionized and neutral gas in high-$z$ galaxies. In addition, we show that deep [C II] observations with a sensitivity of $\sim 10^{-2}~{\rm mJy~arcsec^{-2}}$ can probe the extended neutral gas disks of high-$z$ galaxies.

We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$\alpha$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$\alpha$ emitters with a bright ($\log L_{\rm Ly\alpha}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$\alpha$ emission, and present the photometric and spectroscopic properties of a total of seven LABs; the two new LABs and five previously-known LABs at $z=5.7-6.6$. The z70-1 LAB shows the extended Ly$\alpha$ emission with a scale length of $1.4\pm 0.2$ kpc, about three times larger than the UV continuum emission, making z70-1 the most distant LAB identified to date. All of the 7 LABs, except z49-1, exhibit no AGN signatures such as X-ray emission, {\sc Nv}$\lambda$1240 emission, or Ly$\alpha$ line broadening, while z49-1 has a strong {\sc Civ}$\lambda$1548 emission line indicating an AGN on the basis of the UV-line ratio diagnostics. We carefully model the point-spread functions of the HSC images, and conduct two-component exponential profile fitting to the extended Ly$\alpha$ emission of the LABs. The Ly$\alpha$ scale lengths of the core (star-forming region) and the halo components are $r_{\rm c}=0.6-1.2$ kpc and $r_{\rm h}=2.0-13.8$ kpc, respectively. The average $r_{\rm h}$ of the LABs falls on the extrapolation of the $r_{\rm h}$-Ly$\alpha$ luminosity relation of the Ly$\alpha$ halos around VLT/MUSE star-forming galaxies at the similar redshifts, suggesting that typical LABs at $z\gtrsim5$ are not special objects, but star-forming galaxies at the bright end.

We study various implementations of supernova feedback model and present the results of our `Osaka feedback model' using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {\small GADGET-3}. Our model is a modified version of Stinson et al.'s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical properties of SN remnants from the results of Chevalier and McKee \& Ostriker, such as the effective radius of SN bubble and the remnant life-time, in the form of Sedov-Taylor (ST)-like solutions with the effect of radiative cooling. Our model utilizes the local, physical parameters such as density and temperature of the ISM rather than galactic or halo properties to determine the galactic wind velocity or mass-loading factor. The Osaka model succeeds in self-regulating star formation, and naturally produces galactic outflow with variable velocities depending on the local environment and available SN energy as a function of time.An important addition to our previous work by Aoyama et al. is the implementation of the {\small CELib} chemistry library which allows us to deal with the time-dependent input of energy and metal yields for type Ia \& II supernovae (SNe) and asymptotic giant branch (AGB) stars. As initial tests of our model, we apply it to isolated galaxy simulations, and examine various galactic properties and compare with observational data including metal abundances.

Recent observations have successfully detected UV-bright and infrared-bright galaxies in the epoch of reionization. However, the origin of their radiative properties has not been understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present predictions of multi-wavelength radiative properties of the first galaxies at $z\sim 6-15$. Using zoom-in initial conditions, we investigate three massive galaxies and their satellites in different environment and halo masses: $M_{\rm h} = 2.4\times 10^{10}{\rm M_{\odot}}$ (Halo-10), $1.6\times 10^{11}{\rm M_{\odot}}$ (Halo-11) and $0.7\times 10^{12}{\rm M_{\odot}}$ (Halo-12) at $z=6$. We find that most of gas and dust are ejected from star-forming regions by supernova feedback, which allows UV photons to escape. We show that the peak of the spectral energy distribution (SED) rapidly changes between UV and infrared wavelengths on a time-scale of $\sim$\,100\,Myrs due to intermittent star formation and feedback, and the escape fraction of UV photons fluctuates in the range of $0.2-0.8$ at $z<10$ with a time-averaged value of 0.3. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. We predict the detectability of high-$z$ galaxies with the Atacama Large Millimeter Array (ALMA). For a sensitivity limit of $0.1\,{\rm mJy}$ at $850\,{\rm \mu m}$, the detection probability of galaxies in halos $M_{\rm h} \gtrsim 10^{11}\,{\rm M_{\odot}}$ at $z\lesssim 7$ exceeds fifty per cent. We argue that supernova feedback can produce the observed diversity of SEDs for high-$z$ galaxies.

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

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.

Detection of the redshifted 21cm-line signal from neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) is complicated by intense foregrounds such as galactic synchrotron and extragalactic radio galaxies. The 21cm-Lyman-$\alpha$ emitter(LAE) cross-correlation is one of the tools available to reduce the foreground effects because the foreground emission from such radio sources is statistically independent of LAE distribution. LAE surveys during the EoR at redshifts $z=6.6$ and $7.3$ are ongoing by the Subaru Hyper Suprime-Cam (HSC). Additionally, Prime Focus Spectrograph (PFS) will provide precise redshift information of the LAEs discovered by the HSC survey. In this paper, we investigate the detectability of the 21cm signal with the 21cm-LAE cross-correlation by using our improved reionization simulations. We also focus on the error budget and evaluate it quantitatively in order to consider a strategy to improve the signal-to-noise ratio. In addition, we explore an expansion of the LAE survey to suggest optimal survey parameters and show a potential to measure a characteristic size of ionized bubbles via the turnover scale of the cross-power spectrum. As a result, we find that the Murchison Widefield Array (MWA) has ability to detect the cross-power spectrum signal on large scales by combining LAE Deep field survey of HSC. We also show that the sensitivity is improved dramatically at small scales by adding redshift information from the PFS measurements. The Square Kilometre Array (SKA) has a potential to measure the turnover scale with an accuracy of $6\times10^{-3}~{\rm Mpc^{-1}}$.

The 21-cm signal in the vicinity of the first stars is expected to reflect properties of the first stars. In this study we pay special attention to tracing the time evolution of the ionizing photons' escape fraction, which affects the distribution of neutral hydrogen, by performing radiation hydrodynamics (RHD) simulations resolving dense gas in a halo. We find that the radial profile of 21-cm differential brightness temperature is quite sensitive to the stellar and halo masses, which reflects the time evolution of the escape fraction. In the case of a less massive star, ionizing photons hardly escape from its host halo due to the absorption by dense halo gas, thus an deep 21-cm absorption feature at just outside the halo lasts a long time. Whereas photons from a massive star well working to heat the ambient intergalactic medium turn out to cause a spatially extended 21-cm emission signature. Although individual signals are found to be undetectable with the Square Kilometre Array, our analysis using the results from the RHD simulations indicates that the properties of the first stars are imprinted on the 21-cm global signal: its amplitude depends not only on the cosmic star formation rate density, but also on the typical mass of the first stars due to the stellar-mass-dependent heating rate. Thus, we suggest that the initial mass function of the first stars is an essential factor in understanding the global signal.

The survey of Lyman $\alpha$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined $2^3$ models depending on the presence or absence of dispersion of halo Ly$\alpha$ emissivity, dispersion of the halo Ly$\alpha$ optical depth, $\tau_\alpha$, and halo mass dependence of $\tau_\alpha$. The unique free parameter in our model, a pivot value of $\tau_\alpha$, is calibrated so as to reproduce the $z=5.7$ Ly$\alpha$ luminosity function (LF). We compare our model predictions with Ly$\alpha$ LFs at $z=6.6$ and $7.3$, LAE angular auto-correlation functions (ACFs) at $z=5.7$ and $6.6$, and LAE fractions in Lyman break galaxies at $5<z<7$. The Ly$\alpha$ LFs and ACFs are reproduced by multiple models, but the LAE fraction turns out to be the most critical test. The dispersion of $\tau_\alpha$ and the halo mass dependence of $\tau_\alpha$ are essential to explain all observations reasonably. Therefore, a simple model of one-to-one correspondence between halo mass and Ly$\alpha$ luminosity with a constant Ly$\alpha$ escape fraction has been ruled out. Based on our best model, we present a formula to estimate the intergalactic neutral hydrogen fraction, $x_{\rm HI}$, from the observed Ly$\alpha$ luminosity density at $z\gtrsim6$. We finally obtain $x_{\rm HI}=0.5_{-0.3}^{+0.1}$ as a volume-average at $z=7.3$.

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

We study the impact of Ly$\alpha$ radiation feedback on globular cluster (GC) formation. In this Letter, we analytically derive the relation between star formation efficiency (SFE) and metallicity in spherical clouds with the Ly$\alpha$ radiation feedback. Our models show that the SFE becomes small as the metallicity decreases. In metal-poor gas clouds, Ly$\alpha$ photons are trapped for a long time and exert strong radiation force to the gas, resulting in the suppression of star formation. We find that bound star-clusters (SFE $\gtrsim$ 0.5) form only for the metallicity higher than $\sim 10^{-2.5}~Z_{\odot}$ in the case with the initial cloud mass $10^5~{\rm M}_\odot$ and the radius 5 pc. Our models successfully reproduce the lower bound of observed metallicity of GCs. Thus, we suggest that the Ly$\alpha$ radiation feedback can be essential in understanding the formation of GCs.

We present a novel Lyman alpha (Ly$\alpha$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$\alpha$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$\alpha$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$\alpha$ photons in the interstellar medium of young star-forming galaxies like Ly$\alpha$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$\alpha$ emission, which can be compared to the observations of LAEs.

Understanding Ly$\alpha$ emitting galaxies (LAEs) can be a key to reveal cosmic reionization and galaxy formation in the early Universe. Based on halo merger trees and Ly$\alpha$ radiation transfer calculations, we model redshift evolution of LAEs and their observational properties at $z \ge 6$. We consider ionized bubbles associated with individual LAEs and IGM transmission of Ly$\alpha$ photons. We find that Ly$\alpha$ luminosity tightly correlates with halo mass and stellar mass, while the relation with star formation rate has a large dispersion. Comparing our models with the observed luminosity function by Konno et al., we suggest that LAEs at $z \sim 7$ have galactic wind of $V_{\rm out} \gtrsim 100~\rm km\, s^{-1}$ and HI column density of $N_{\rm HI} \gtrsim 10^{20}~\rm cm^{-2}$. Number density of bright LAEs rapidly decreases as redshift increases, due to both lower star formation rate and smaller HII bubbles. Our model predicts future wide deep surveys with next generation telescopes, such as JWST, E-ELT and TMT, can detect LAEs at $z \sim 10$ with a number density of $n_{\rm LAE} \sim {\rm a~few~} \times10^{-6} ~\rm Mpc^{-3}$ for the flux sensitivity of $10^{-18} ~\rm erg\, cm^{-2}\, s^{-1}$. By combining these surveys with future 21-cm observations, it could be possible to detect both LAEs with $L_{Ly\alpha} \gtrsim 10^{42}~\rm erg~s^{-1}$ and their associated giant HII bubbles with the size $\gtrsim 250 ~\rm kpc$ at $z \sim 10$.

Star formation in high-redshift dwarf galaxies is a key to understand early galaxy evolution in the early Universe. Using the three-dimensional hydrodynamics code GIZMO, we study the formation mechanism of cold, high-density gas clouds in interacting dwarf galaxies with halo masses of $\sim 3 \times 10^{7}~M_{\odot}$, which are likely to be the formation sites of early star clusters. Our simulations can resolve both the structure of interstellar medium on small scales of $\lesssim 0.1$ pc and the galactic disk simultaneously. We find that the cold gas clouds form in the post-shock region via thermal instability due to metal-line cooling, when the cooling time is shorter than the galactic dynamical time. The mass function of cold clouds shows almost a power-law initially with an upper limit of thermally unstable scale. We find that some clouds merge into more massive ones with $\gtrsim 10^{4}~M_{\odot}$ within $\sim 2~{\rm Myr}$. Only the massive cold clouds with $\gtrsim 10^{3}~M_{\odot}$ can keep collapsing due to gravitational instability, resulting in the formation of star clusters. In addition, we investigate the dependence of cloud mass function on metallicity and ${\rm H_{2}}$ abundance, and show that the cases with low metallicities ($\lesssim 10^{-2}~Z_{\odot}$) or high ${\rm H_{2}}$ abundance ($\gtrsim 10^{-3}$) cannot form massive cold clouds with $\gtrsim 10^{3}~M_{\odot}$.

We present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution towards observable galaxies at $z \sim 6$. We focus on three different galaxies which end up in halos with masses $M_{h} = 2.4 \times10^{10}~h^{-1}\; M_{\odot}$ (Halo-10), $1.6 \times10^{11}~h^{-1}\; M_{\odot}$ (Halo-11) and $0.7 \times10^{12}~h^{-1} M_{\odot}$ (Halo-12) at z=6. Our simulations also probe impacts of different sub-grid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that star formation occurs intermittently due to supernova (SN) feedback at z > 10, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no-feedback or lower SF efficiency models can sustain galactic disk for long periods > 10 Myr. The expulsion of gas at the galactic center also affects the inner dark matter density profile. However, SN feedback does not seem to keep the shallow profile of dark matter for a long period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the star formation rates (SFR) and stellar masses of observed Lyman-$\alpha$ emitters (LAEs) at z = 7-8 fairly well given observational uncertainties. In addition, we investigate the effect of UV background radiation on star formation as an external feedback source, and find that earlier reionization extends the quenching time of star formation due to photo-ionization heating, but does not affect the stellar mass at z=6.

Evidence for dust around supermassive black holes (SMBHs) in the early Universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black-holes (IMBHs) of $\sim 10^{5}~M_{\odot}$ in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the accretion of dusty gas onto IMBHs proceeds gently with small fluctuations of the accretion rate, whereas that of pristine gas causes more violent periodic bursts. At dust-to-gas mass ratios similar to the solar neighborhood, the time averaged luminosity becomes smaller than that for primordial gas by one order of magnitude and the time-averaged Eddington ratio ranges from $\sim 10^{-4}$ to $\sim 10^{-2}$ in clouds with initial gas densities of $n_{\rm H} = 10 - 1000~\rm cm^{-3}$. Our calculations show that the effect of dust opacity alone is secondary compared to the radiation pressure on dust in regulating the BH growth. We also derive spectral energy distributions at IR bands by calculating dust thermal emission and show that the flux ratio between $\lambda \lesssim 20~\rm \mu m$ and $\gtrsim 100~\rm \mu m$ is closely related to the Eddington ratio. Thermal emission from hot dust near the BH dominates only during the high accretion phase, producing higher flux density at $\lesssim 20~\rm \mu m$. Therefore, we suggest that the combinations of MIR observations by JWST and FIR observation by ALMA or Spitzer can be used to estimate the Eddington ratio of massive BHs.

We here derive upper limits on the mass and luminosity of Population III (POPIII) dominated proto-galaxies based on the collapse of primordial gas under the effect of angular momentum loss via Ly$\alpha$ radiation drag and the gas accretion onto a galactic centre. Our model predicts that POPIII-dominated galaxies at z ~ 7 are hosted by haloes with $M_{\rm halo} \sim 1.5 \times 10^{8} - 1.1 \times 10^{9} \rm ~M_{\odot}$, that they have Ly$\alpha$ luminosities of $L_{\rm Ly\alpha} \sim 3.0 \times 10^{42} - 2.1 \times 10^{43}$ erg/s, stellar mass of $M_{\rm star} \sim 0.8 \times 10^{5} - 2.5 \times 10^{6} \rm ~M_{\odot}$, and outflowing gas with velocities $V_{\rm out} \sim 40$ km/s due to Ly$\alpha$ radiation pressure. We show that the POPIII galaxy candidate CR7 violates the derived limits on stellar mass and Ly$\alpha$ luminosity and thus is unlikely to be a POPIII galaxy. POPIII-dominated galaxies at z ~ 7 have HeII line emission that is ~1- 3 orders of magnitude lower then that of Ly$\alpha$, they have high Ly$\alpha$ equivalent width of > 300 $\AA$ and should be found close to bright star forming galaxies. The HeII 1640 $\AA$ line is in comfortable reach of next generation telescopes, like the JWST or TMT.

Cosmic reionization is known to be a major phase transition of the gas in the Universe. Since astronomical objects formed in the early Universe, such as the first stars, galaxies and black holes, are expected to have caused cosmic reionization, the formation history and properties of such objects are closely related to the reionization process. In spite of the importance of exploring reionization, our understandings regarding reionization is not sufficient yet. Square Kilometre Array (SKA) is a next-generation large telescope that will be operated in the next decade. Although several programs of next-generation telescopes are currently scheduled, the SKA will be the unique telescope with a potential to directly observe neutral hydrogen up to z~30, and provide us with valuable information on the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The early science with the SKA will start in a few years; it is thus the time for us to elaborate a strategy for CD/EoR Science with the SKA. The purpose of this document is to introduce Japanese scientific interests in the SKA project and to report results of our investigation.

We present calculations on the formation of massive black holes with 10^5 Msun at z > 6 that can be the seeds of supermassive black holes at z > 6. Under the assumption of compact star cluster formation in merging galaxies, star clusters in haloes of 10^8 ~ 10^9 Msun can undergo rapid core-collapse leading to the formation of very massive stars (VMSs) with ~1000 Msun which directly collapse into black holes with similar masses. Star clusters in halos of > 10^9 Msun experience type-II supernovae before the formation of VMSs due to long core-collapse time scales. We also model the subsequent growth of black holes via accretion of residual stars in clusters. 2-body relaxation efficiently re-fills the loss cones of stellar orbits at larger radii and resonant relaxation at small radii is the main driver for accretion of stars onto black holes. As a result, more than ninety percent of stars in the initial cluster are swallowed by the central black holes before z=6. Using dark matter merger trees we derive black hole mass functions at z=6-20. The mass function ranges from 10^3 to 10^5 Msun at z <~ 15. Major merging of galaxies of >~ 4*10^8 Msun at z ~ 20 successfully leads to the formation of >~ 10^5 Msun BHs by z >~ 10 which can be the potential seeds of supermassive black holes seen today.

We use high-resolution zoom-in cosmological simulations of galaxies of Romano-Diaz et al., post-processing them with a panchromatic three-dimensional radiation transfer code to obtain the galaxy UV luminosity function (LF) at z ~ 6-12. The galaxies are followed in a rare, heavily overdense region within a ~ 5-sigma density peak, which can host high-z quasars, and in an average density region, down to the stellar mass of M_star ~ 4* 10^7 Msun. We find that the overdense regions evolve at a substantially accelerated pace --- the most massive galaxy has grown to M_star ~ 8.4*10^10 Msun by z = 6.3, contains dust of M_dust~ 4.1*10^8 Msun, and is associated with a very high star formation rate, SFR ~ 745 Msun/yr.The attained SFR-M_star correlation results in the specific SFR slowly increasing with M_star. Most of the UV radiation in massive galaxies is absorbed by the dust, its escape fraction f_esc is low, increasing slowly with time. Galaxies in the average region have less dust, and agree with the observed UV LF. The LF of the overdense region is substantially higher, and contains much brighter galaxies. The massive galaxies are bright in the infrared (IR) due to the dust thermal emission, with L_IR~ 3.7*10^12 Lsun at z = 6.3, while L_IR < 10^11 Lsun for the low-mass galaxies. Therefore, ALMA can probe massive galaxies in the overdense region up to z ~ 10 with a reasonable integration time. The UV spectral properties of disky galaxies depend significantly upon the viewing angle.The stellar and dust masses of the most massive galaxy in the overdense region are comparable to those of the sub-millimetre galaxy (SMG) found by Riechers et al. at z = 6.3, while the modelled SFR and the sub-millimetre flux fall slightly below the observed one. Statistical significance of these similarities and differences will only become clear with the upcoming ALMA observations.

The Lyman-alpha (Lya) emission has played an important role in detecting high-redshift galaxies, including recently distant ones at redshift z > 7. It may also contain important information on the origin of these galaxies. Here, we investigate the formation of a typical L* galaxy and its observational signatures at the earliest stage, by combining a cosmological hydrodynamic simulation with three-dimensional radiative transfer calculations using the newly improved ART^2 code. Our cosmological simulation uses the Aquila initial condition which zooms in onto a Milky Way-like halo with high resolutions, and our radiative transfer couples multi-wavelength continuum, Lya line, and ionization of hydrogen. We find that the modeled galaxy starts to form at redshift z ~ 24 through efficient accretion of cold gas, which produces a strong Lya line with a luminosity of L(Lya) ~ 10^42 erg/s as early as z ~ 14. The Lya emission appears to trace the cold, dense gas. The lines exhibit asymmetric, single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, the contribution to the total Lya luminosity by excitation cooling increases with redshift, and it becomes dominant at z >~ 6. We predict that L* galaxies such as the modeled one may be detected at z <~ 8 by JWST and ALMA with a reasonable integration time. Beyond redshift 12, however, only Lya line may be observable by spectroscopic surveys. Our results suggest that Lya line is one of the most powerful tools to detect the first generation of galaxies, and to decipher their formation mechanism.

Using varying models for the star formation rate (SFR) of Population (Pop) III and II stars at z>6 we derive the expected redshift history of the global 21 cm signal from the inter-galactic medium (IGM). To recover the observed Thomson scattering optical depth of the cosmic microwave background (CMB) requires SFRs at the level of ~ 10^-3 Msun/yr/Mpc^3 at z ~ 15 from Pop III stars, or ~ 10^-1Msun/yr/Mpc^3 at z ~ 7 from Pop II stars. In the case the SFR is dominated by Pop III stars, the IGM quickly heats above the CMB at z > 12 due to heating from supernovae. In addition, Lyman-alpha photons from haloes hosting Pop III stars couple the spin temperature to that of the gas, resulting in a deep absorption signal. If the SFR is dominated by Pop II stars, the IGM slowly heats and exceeds the CMB temperature at z ~ 10. However, the larger and varying fraction of Pop III stars is able to break this degeneracy. We find that the impact of the initial mass function (IMF) of Pop III stars on the 21 cm signal results in an earlier change to a positive signal if the IMF slope is ~ -1.2. Measuring the 21 cm signal at z > 10 with next generation radio telescopes such as the Square Kilometre Array will be able to investigate the contribution from Pop III and Pop II stars to the global star formation rate.

Observations of the redshifted 21 cm line with upcoming radio telescopes promise to transform our understanding of the cosmic reionization. To unravel the underlying physical process, we investigate the 21 cm structures of three different ionizing sources, Population (Pop) III stars, the first galaxies, and the first quasars, by using radiative transfer simulations that include both ionization of neutral hydrogen and resonant scattering of Lya photons. We find that Pop III stars and quasars produce a smooth transition from an ionized and hot state to a neutral and cold one, owing to their hard spectral energy distribution with abundant ionizing photons, in contrast to the sharp transition in galaxies. Furthermore, Lya scattering plays a dominant role in producing the 21 cm signal as it determines the relation between hydrogen spin temperature and gas kinetic temperature. This effect, also called Wouthuysen-Field coupling, depends strongly on the ionizing source. It is the strongest around galaxies, where the spin temperature is highly coupled to that of the gas, resulting in extended absorption troughs in the 21 cm brightness temperature. On the other hand, in the case of Pop III stars, the 21 cm signal shows both emission and absorption regions around a small HII bubble. For quasars, a large emission region in the 21 cm signal is produced, and the absorption region decreases as the size of the HII bubble becomes large due to the limited traveling time of photons. We predict that future surveys from large radio arrays such as MWA, LOFAR and SKA may be able to detect the 21 cm signals of primordial galaxies and quasars, but not likely Pop III stars due to its small angular diameter.

We present results on the radiation drag exerted by an isotropic and homogenous background of Lya photons on neutral gas clouds orbiting within HII regions around Population III stars of different masses. The Doppler shift causes a frequency difference between photons moving in the direction of the cloud and opposite to it resulting in a net momentum loss of the cloud in the direction of motion. We find that half of the angular momentum of gas with v_theta <~ 20 km/s near (r <~ 3 kpc) a Population III star of 120 Msun at z=20 is lost within ~ 10^6 yr. The radiation drag is a strong function of cloud velocity that peaks at v ~ 20 km/s reflecting the frequency dependence of the photon cross section. Clouds moving with velocities larger than ~ 100 km/s loose their angular momentum on time scales of ~ 10^8 yr. At lower redshifts radiation drag becomes inefficient as the Lya photon density in HII regions decreases by a factor (1+z)^3 and angular momentum is lost on time scales > 10^8 yr even for low velocity clouds. Our results suggest that a sweet spot exists for the loss of angular momentum by radiation drag for gas clouds at z > 10 and with v ~ 20 km/s. Comparison to dynamical friction forces acting on typical gas clouds suggest that radiation drag is the dominant effect impacting the orbit. We propose that this effect can suppress the formation of extended gas discs in the first galaxies and help gas accretion near galactic centres and central black holes.

A large number of high-redshift galaxies have been discovered via their narrow-band Lya line or broad-band continuum colors in recent years. The nature of the escaping process of photons from these early galaxies is crucial to understanding galaxy evolution and the cosmic reionization. Here, we investigate the escape of Lya, non-ionizing UV-continuum (l = 1300 - 1600 angstrom in rest frame), and ionizing photons (l < 912 angstrom) from galaxies by combining a cosmological hydrodynamic simulation with three-dimensional multi-wavelength radiative transfer calculations. The galaxies are simulated in a box of 5^3 h^-3 Mpc^3 with high resolutions using the Aquila initial condition which reproduces a Milky Way-like galaxy at redshift z=0. We find that the escape fraction (fesc) of these different photons shows a complex dependence on redshift and galaxy properties: fesc(Lya) and fesc(UV) appear to evolve with redshift, and they show similar, weak correlations with galaxy properties such as mass, star formation, metallicity, and dust content, while fesc(Ion) remains roughly constant at ~ 0.2 from z ~ 0 - 10, and it does not show clear dependence on galaxy properties. fesc(Lya) correlates more strongly with fesc(UV) than with fesc(Ion). In addition, we find a relation between the emergent Lya luminosity and the ionizing photon emissivity of Lyman Alpha Emitters (LAEs). By combining this relation with the observed luminosity functions of LAEs at different redshift, we estimate the contribution from LAEs to the reionization of intergalactic medium (IGM). Our result suggests that ionizing photons from LAEs alone are not sufficient to ionize IGM at z > 6, but they can maintain the ionization of IGM at z ~ 0 - 5.

Recent observations have indicated the existence of dust in high-redshift galaxies, however, the dust properties in them are still unknown. Here we present theoretical constraints on dust properties in Lyman break galaxies (LBGs) at z=3 by post-processing a cosmological smoothed particle hydrodynamics simulation with radiative transfer calculations. We calculate the dust extinction in 2800 dark matter halos using the metallicity information of individual gas particles in our simulation. We use only bright galaxies with rest-frame UV magnitude M_1700 < -20 mag, and study the dust size, dust-to-metal mass ratio, and dust composition. From the comparison of calculated color excess between B and V-band (i.e., E(B-V)) and the observations, we constrain the typical dust size, and show that the best-fitting dust grain size is ~ 0.05 micron, which is consistent with the results of theoretical dust models for Type-II supernova. Our simulation with the dust extinction effect can naturally reproduce the observed rest-frame UV luminosity function of LBGs at z=3 without assuming an ad hoc constant extinction value. In addition, in order to reproduce the observed mean E(B-V), we find that the dust-to-metal mass ratio needs to be similar to that of the local galaxies, and that the graphite dust is dominant or at least occupy half of dust mass.

Recent observations show that a large number of Lyman-alpha emitters (LAEs) at high redshift z >= 3 have unusually high Lya equivalent widths (EW > 400 angstrom). However, the origin of these high EWs is an open question. Here, we investigate the impacts of photon trapping on the Lya EW and other properties by tracking the Lya radiative transfer in spherical galactic clouds. We find that the delayed escape of the Lya photons can change the Lya properties significantly. During the transition phase from optically thick to optically thin where the Lya photons can escape simultaneously, the EW can be boosted to ~ 1000 angstrom, the Lya luminosity can be increased by a factor of a few, and the line profile can be significantly broadened. The boost factor appears to depend on the galaxy properties such as mass and star formation rate and timescale, therefore future investigation combing 3D Lya RT calculations with cosmological simulations of galaxy formation and evolution is needed to fully understand the Lya properties of early star-forming galaxies.

Recent observations have discovered a population of extended Lya sources, dubbed Lyman-alpha blobs (LABs), at high redshift z ~ 3 - 6.6. These LABs typically have a luminosity of L ~ 10^42-10^44 erg/s, and a size of tens of kiloparsecs, with some giant ones reaching up to D ~ 100 kpc. However, the origin of these LABs is not well understood. In this paper, we investigate a merger model for the formation of LABs by studying Lya emission from interacting galaxies at high redshifts by means of a combination of hydrodynamics simulations with three dimensional radiative transfer calculations. Our galaxy simulations focus on a set of binary major mergers of galaxies with a mass range of 3-7 *10^12 Msun in the redshift range of z ~ 3 -7, and we use the newly improved ART^2 code to perform the radiative transfer calculations which couple multi-wavelength continuum, ionization of hydrogen, and Lya line emission. We find that intense star formation and enhanced cooling induced by gravitational interaction produce strong Lya emission from these merging galaxies. The Lya emission appears to be extended due to the extended distribution of sources and gas. During the close encounter of galaxy progenitors when the star formation rate peaks at ~ 10^3 Msun/yr, our model produces Lya blobs with luminosity of L ~ 10^42-10^44 erg/s, and size of D ~ 10-20 kpc at z>6 and D ~ 20-50 kpc at z ~ 3, in broad agreement with observations in the same redshift range. Our results suggest that merging galaxies may produce some typical LABs as observed, but the giant ones may be produced by mergers more massive than those in our model, or a combination of mergers and cold accretion from filaments on a large scale.

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.

We study the effect of local stellar radiation and UVB on the physical properties of DLAs and LLSs at z=3 using cosmological SPH simulations. We post-process our simulations with the ART code for radiative transfer of local stellar radiation and UVB. We find that the DLA and LLS cross sections are significantly reduced by the UVB, whereas the local stellar radiation does not affect them very much except in the low-mass halos. This is because clumpy high-density clouds near young star clusters effectively absorb most of the ionizing photons from young stars. We also find that the UVB model with a simple density threshold for self-shielding effect can reproduce the observed column density distribution function of DLAs and LLSs very well, and we validate this model by direct radiative transfer calculations of stellar radiation and UVB with high angular resolution. We show that, with a self-shielding treatment, the DLAs have an extended distribution around star-forming regions typically on ~ 10-30 kpc scales, and LLSs are surrounding DLAs on ~ 30-60 kpc scales. Our simulations suggest that the median properties of DLA host haloes are: Mh = 2.4*10^10 Msun, SFR = 0.3 Msun/yr, M* = 2.4*10^8 Msun, and Z/Zsun = 0.1. About 30 per cent of DLAs are hosted by haloes having SFR = 1 - 20 Msun/yr, which is the typical SFR range for LBGs. More than half of DLAs are hosted by the LBGs that are fainter than the current observational limit. Our results suggest that fractional contribution to LLSs from lower mass haloes is greater than for DLAs. Therefore the median values of LLS host haloes are somewhat lower with Mh = 9.6*10^9 Msun, SFR = 0.06 Msun/yr, M* = 6.5*10^7 Msun and Z/Zsun = 0.08. About 80 per cent of total LLS cross section are hosted by haloes with SFR < 1 Msun/yr, hence most LLSs are associated with low-mass halos with faint LBGs below the current detection limit.

Narrow-band Lya line and broad-band continuum have played important roles in the discovery of high-redshift galaxies in recent years. Hence, it is crucial to study the radiative transfer of both Lya and continuum photons in the context of galaxy formation and evolution in order to understand the nature of distant galaxies. Here, we present a three-dimensional Monte Carlo radiative transfer code, All-wavelength Radiative Transfer with Adaptive Refinement Tree (ART^2), which couples Lya line and multi-wavelength continuum, for the study of panchromatic properties of galaxies and interstellar medium. This code is based on the original version of Li et al., and features three essential modules: continuum emission from X-ray to radio, Lya emission from both recombination and collisional excitation, and ionization of neutral hydrogen. The coupling of these three modules, together with an adaptive refinement grid, enables a self-consistent and accurate calculation of the Lya properties. As an example, we apply ART^2 to a cosmological simulation that includes both star formation and black hole growth, and study in detail a sample of massive galaxies at redshifts z=3.1 - 10.2. We find that these galaxies are Lya emitters (LAEs), whose Lya emission traces the dense gas region, and that their Lya lines show a shape characteristic of gas inflow. Furthermore, the Lya properties, including photon escape fraction, emergent luminosity, and equivalent width, change with time and environment. Our results suggest that LAEs evolve with redshift, and that early LAEs such as the most distant one detected at z ~ 8.6 may be dwarf galaxies with a high star formation rate fueled by infall of cold gas, and a low Lya escape fraction.

The Lya emission has been observed from galaxies over a redshift span z ~ 0 - 8.6. However, the evolution of high-redshift Lya emitters (LAEs), and the link between these populations and local galaxies, remain poorly understood. Here, we investigate the Lya properties of progenitors of a local L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative transfer calculations using the new ART^2 code. We find that the main progenitor (the most massive one) of a Milky Way-like galaxy has a number of Lya properties close to those of observed LAEs at z ~ 2 - 6, but most of the fainter ones appear to fall below the detection limits of current surveys. The Lya photon escape fraction depends sensitively on a number of physical properties of the galaxy, such as mass, star formation rate, and metallicity, as well as galaxy morphology and orientation. Moreover, we find that high-redshift LAEs show blue-shifted Lya line profiles characteristic of gas inflow, and that the Lya emission by excitation cooling increases with redshift, and becomes dominant at z > 6. Our results suggest that some observed LAEs at z ~ 2-6 with luminosity of L_Lya ~ 10^{42-43} ergs/s may be similar to the main progenitor of the Milky Way at high redshift, and that they may evolve into present-day L* galaxies.

Based on a three-dimensional model of an early star-forming galaxy, we explore the evolution of the sub-millimeter brightness. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, where the SFR is ~ 10 Msun/yr at t<0.3 Gyr and several Msun/yr at t>0.3 Gyr. The former phase well reproduces the observed properties of LAEs and the latter does LBGs. We solve the three-dimensional radiative transfer in the clumpy interstellar media in this model galaxy, taking the size distributions of dust grains into account, and calculate the dust temperature as a function of galactic evolutionary time. We find that the clumpiness of interstellar media plays an important role for the sub-millimeter brightness. In the LAE phase, dust grains are concentrated on clumpy star-forming regions that are distributed all over the galaxy, and the grains can effectively absorb UV radiation from stars. As a result, the dust is heated up to T>35 K. In the LBG phase, the continuous supernovae drive dust grains far away from star-forming regions. Then, the grains cannot absorb much radiation from stars, and becomes into a cold state close to the CMB temperature. Consequently, the dust temperature decreases with the evolutionary time, where the mass-weighted mean temperature is T=26 K at t=0.1 Gyr and T=21 K at t=1.0 Gyr. By this analysis, it turns out that the sub-millimeter brightness is higher in the LAE phase than that in the LBG phase, although the dust-to-gas ratio increases monotonically as a function of time. We derive the spectral energy distributions by placing the model galaxy at a given redshift. The peak flux at 850 micron is found to be S_850 ~ 0.2 - 0.9 mJy if the model galaxy is placed at 6>z>2. This means that ALMA can detect an early star-forming galaxy with SFR of ~ 10 Msun/yr by less than one hour integration with 16 antennas.

We study the impact of ultraviolet background (UVB) radiation field and the local stellar radiation on the H_I column density distribution f(N_HI) of damped Ly-alpha systems (DLAs) and sub-DLAs at z=3 using cosmological smoothed particle hydrodynamics simulations. We find that, in the previous simulations with an optically thin approximation, the UVB was sinking into the H_I cloud too deeply, and therefore we underestimated the f(N_HI) at 19 < log(N_HI) < 21.2 compared to the observations. If the UVB is shut off in the high-density regions with n_gas > 6 x 10^{-3} cm^{-3}, then we reproduce the observed f(N_HI) at z=3 very well. We also investigate the effect of local stellar radiation by post-processing our simulation with a radiative transfer code, and find that the local stellar radiation does not change the f(N_HI) very much. Our results show that the shape of f(N_HI) is determined primarily by the UVB with a much weaker effect by the local stellar radiation and that the optically thin approximation often used in cosmological simulation is inadequate to properly treat the ionization structure of neutral gas in and out of DLAs. Our result also indicates that the DLA gas is closely related to the transition region from optically-thick neutral gas to optically-thin ionized gas within dark matter halos.

Combing the three-dimensional radiative transfer (RT) calculation and cosmological SPH simulations, we study the escape fraction of ionizing photons (f_esc) of high-redshift galaxies at z=3-6. Our simulations cover the halo mass range of M_h = 10^9 - 10^12 M_sun. We postprocess several hundred simulated galaxies with the Authentic Radiative Transfer (ART) code to study the halo mass dependence of f_esc. In this paper, we restrict ourselves to the transfer of stellar radiation from local stellar population in each dark matter halo. We find that the average f_esc steeply decreases as the halo mass increases, with a large scatter for the lower mass haloes. The low mass haloes with M_h ~ 10^9 M_sun have large values of f_esc (with an average of ~ 0.4), whereas the massive haloes with M_h ~ 10^11 M_sun show small values of f_esc (with an average of ~ 0.07). This is because in our simulations, the massive haloes show more clumpy structure in gas distribution, and star-forming regions are embedded inside these clumps, making it more difficult for the ionizing photons to escape. On the other hand, in low mass haloes, there are often conical regions of highly ionized gas due to the shifted location of young star clusters from the center of dark matter halo, which allows the ionizing photons to escape more easily than in the high-mass haloes. By counting the number of escaped ionizing photons, we show that the star-forming galaxies can ionize the intergalactic medium at z=3-6. The main contributor to the ionizing photons is the haloes with M_h < 10^10 M_sun owing to their high f_esc. The large dispersion in f_esc suggests that there may be various sizes of H{\sc ii} bubbles around the haloes even with the same mass in the early stages of reionization. We also examine the effect of UV background radiation field on f_esc using simple, four different treatment of UV background.

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) at redshifts 3<z<7 to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura (2006),which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is 10^{11}M_sun. We solve not only photo-ionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor 1.5-8.5 in the LAE phase and by a factor 2.5-11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor ~ 2. The resultant escape fractions are 0.07-0.47 in the LAE phase and 0.06-0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux ratio at 1500 angstrom to 900 angstrom of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high-z LAEs and LBGs can ionize the IGM at z=3-5. However, ionizing radiation from LAEs as well as LBGs falls short to ionize the IGM at z>6. That implies that additional ionization sources may required at z>6.

We study massive spin-1/2 fermion production in nonsingular superstring cosmology, taking into account one-loop quantum corrections to a superstring effective action with dilaton and modulus fields. While no creation occurs in the massless limit, massive fermions can be produced by the existence of a time-dependent frequency. Due to the increase of the Hubble expansion rate during the modulus-driven phase, the occupation of number of fermions continues to grow until the point of the graceful exit, after which fermion creation ceases with the decrease of the Hubble rate.

We study the dynamics and perturbations during inflation and reheating in a multi-field model where a second scalar field $\chi$ is nonminimally coupled to the scalar curvature $(\frac12 \xi R\chi^2$). When $\xi$ is positive, the usual inflationary prediction for large-scale anisotropies is hardly altered while the $\chi$ fluctuation in sub-Hubble modes can be amplified during preheating for large $\xi$. For negative values of $\xi$, however, long-wave modes of the $\chi$ fluctuation exhibit exponential increase during inflation, leading to the strong enhancement of super-Hubble metric perturbations even when $|\xi|$ is less than unity. This is because the effective $\chi$ mass becomes negative during inflation. We constrain the strength of $\xi$ and the initial $\chi$ by the amplitude of produced density perturbations. One way to avoid nonadiabatic growth of super-Hubble curvature perturbations is to stabilize the $\chi$ mass through a coupling to the inflaton. Preheating may thus be necessary in these models to protect the stability of the inflationary phase.