We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find.

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$\alpha$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$\alpha$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

Estimating the temperature and mass of dust in high-$z$ galaxies is essential for discussions of the origin of dust in the early Universe. However, this suffers from limited sampling of the infrared spectral-energy distribution. Here we present an algorithm for deriving the temperature and mass of dust in a galaxy, assuming dust to be in radiative equilibrium. We formulate the algorithm for three geometries: a thin spherical shell, a homogeneous sphere, and a clumpy sphere. We also discuss effects of the mass absorption coefficients of dust at ultraviolet and infrared wavelengths, $\kappa_{\rm UV}$ and $\kappa_{\rm IR}$, respectively. As an example, we apply the algorithm to a normal, dusty star-forming galaxy at $z=7.5$, A1689zD1, for which three data points in the dust continuum are available. Using $\kappa_{\rm UV}=5.0\times10^4$ cm$^2$ g$^{-1}$ and $\kappa_{\rm IR}=30(\lambda/100\mu m)^{-\beta}$ cm$^2$ g$^{-1}$ with $\beta=2.0$, we obtain dust temperatures of 38--70~K and masses of $10^{6.5-7.3}$ M$_\odot$ for the three geometries considered. We obtain similar temperatures and masses from just a single data point in the dust continuum, suggesting the usefulness of the algorithm for high-$z$ galaxies with limited infrared observations. In the clumpy-sphere case, the temperature becomes equal to that of the usual modified black-body fit, because an additional parameter describing the clumpiness works as an adjuster. The best-fit clumpiness parameter is $\xi_{\rm cl}=0.1$, corresponding to $\sim10$\% of the volume filling factor of the clumps in this high-$z$ galaxy if the clump size is $\sim10$ pc, similar to that of giant molecular clouds in the local Universe.

We present our analysis of the LyC emission and escape fraction of 111 spectroscopically verified galaxies with and without AGN from $2.26<z<4.3$. We extended our ERS sample from Smith et al. (2018; arXiv:1602.01555) with 64 galaxies in the GOODS North and South fields using WFC3/UVIS F225W, F275W, and F336W mosaics we independently drizzled using the HDUV, CANDELS, and UVUDF data. Among the 17 AGN from the 111 galaxies, one provided a LyC detection in F275W at $m_{AB}=23.19$ mag (S/N $\simeq$ 133) and $GALEX$ NUV at $m_{AB}=23.77$ mag (S/N $\simeq$ 13). We simultaneously fit $SDSS$ and $Chandra$ spectra of this AGN to an accretion disk and Comptonization model and find $f_{esc}$ values of $f_{esc}^{F275W}\simeq 28^{+20}_{-4}$% and $f_{esc}^{NUV}\simeq 30^{+22}_{-5}$%. For the remaining 110 galaxies, we stack image cutouts that capture their LyC emission using the F225W, F275W, and F336W data of the GOODS and ERS samples, and both combined, as well as subsamples of galaxies with and without AGN, and $all$ galaxies. We find the stack of 17 AGN dominate the LyC production from $\langle z\rangle\simeq 2.3-4.3$ by a factor of $\sim$10 compared to all 94 galaxies without AGN. While the IGM of the early universe may have been reionized mostly by massive stars, there is evidence that a significant portion of the ionizing energy came from AGN.

Detecting $\rm H_I$ 21cm line in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) suffers from foreground contamination such as Galactic synchrotron and extragalactic radio sources. Cross-correlation between the 21cm line and Lyman-$\alpha$ emitter (LAE) galaxies is a powerful tool to identify the 21cm signal since the 21cm line emission has correlation with LAEs while the LAEs are statistically independent of the foregrounds. So far, the detectability of 21cm-LAE cross-power spectrum has been investigated with simple LAE models where the observed Ly$\alpha$ luminosity is proportional to the dark matter halo mass. However, the previous models were inconsistent with the latest observational data of LAEs obtained with Subaru/Hyper Suprime-Cam (HSC). Here, we revisit the detectability of 21cm-LAE cross-power spectrum adopting a state-of-the-art LAE model consistent with all Subaru/HSC observations such as the Ly$\alpha$ luminosity function, LAE angular auto-correlation, and the LAE fractions in the continuum selected galaxies. We find that resultant cross-power spectrum with the updated LAE model is reduced at small scales ($k\sim 1\ \rm Mpc^{-1}$) compared to the simple models, while the amplitudes at large scales ($k \lesssim 0.2 \ \rm Mpc^{-1}$) are not affected so much. We conclude that the large-scale signal would be detectable with Square Kilometre Array (SKA) and HSC LAE cross-correlation but detecting the small scale signal would require an extended HSC LAE survey with an area of $\sim 75\ \rm deg^2$ or 3000 hrs observation time of 21cm line with SKA.

Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$\alpha)\leq50$\AA\ and we find no correlation or anti-correlation between EW (Ly$\alpha$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$.

We calculate the pulsed fraction (PF) of the super-critical column accretion flows onto magnetized neutron stars (NSs), of which the magnetic axis is misaligned with the rotation axis, based on the simulation results by Kawashima et al.(2016, PASJ, 68, 83). Here, we solve the geodesic equation for light in the Schwarzschild spacetime in order to take into account the light bending effect. The gravitational redshift and the relativistic doppler effect from gas motions of the accretion columns are also incorporated. The pulsed emission appears since the observed luminosity, which exceeds the Eddington luminosity for the stellar-mass black holes, periodically changes via precession of the column caused by the rotation of the NS. The PF tends to increase as $\theta_{\rm obs}$ approaching to $\theta_{\rm B}$, where $\theta_{\rm obs}$ and $\theta_{\rm B}$ are the observer's viewing angle and the polar angle of the magnetic axis measured from the rotation axis. The maximum PF is around 50 %. Also, we find that the PF becomes less than 5 % for $\theta_{\rm obs} \lesssim 5^\circ$ or for $\theta_{\rm B} \lesssim 5^\circ$. Our results are consistent with observations of ultraluminous X-ray pulsars (ULXPs) with few exceptions, since the ULXPs mostly exhibit the PF of $\lesssim$ 50 %. Our present study supports the hypothesis that the ULXPs are powered by the super-critical column accretion onto NSs.

Using a sample of galaxies at $z\approx 3$ with detected Lyman Continuum (LyC) leakage in the SSA22 field, we attempt to verify a proposed indirect method for identifying cases with high LyC escape fraction \fesc\ based on measurements of the H$\beta$ equivalent width (EW) and the $\beta$ slope of the UV continuum. To this end, we present Keck/MOSFIRE H$\beta$ emission line flux measurements of LyC galaxies (LCGs) at spectroscopic redshifts $z_{\mathrm{spec}} \sim 3.3$, Lyman break galaxies (LBGs) at photometric redshifts $z_{\mathrm{phot}} = 2.7$--$3.7$, and Ly$\alpha$ emitters at $z_{\mathrm{phot}} = 3.1$. We also reconfirm the spectroscopic redshifts and measure the H$\beta$ emission line fluxes from 2 LCGs and 6 LBGs. For the LCG in our sample with the most extreme \fesc, as revealed by the direct detection of LyC photons, we find that the EW(H$\beta$)--$\beta$ method gives a broadly consistent estimate for \fesc, although the error bars remain very large. We also discuss how a combination of $f_\mathrm{esc}$ measurements based on direct and indirect methods can shed light on the LyC escape mechanism and the anisotropy of the leakage.

We report on the detection of the [CII] 157.7 $\mu$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $\mu$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $\mu$m dust emission (90 $\mu$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $\mu$m rest-frame) down to 18 $\mu$Jy (3$\sigma$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $\mu$Jy continuum emission at 850 $\mu$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($\beta_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $\beta_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization.

The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity star-forming clouds which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses. We combine the observational information from galaxies at redshifts 5 < z < 10 to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle. Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we build a dust emission template for this population of galaxies in the epoch of reionization. Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media. We have found evidence for the first dust grains formed in the universe. But, above all, this letter underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the epoch of reionization.

We search for galaxies with a strong Balmer break (Balmer Break Galaxies; BBGs) at $z \sim 6$ over a 0.41 deg$^2$ effective area in the COSMOS field. Based on rich imaging data, including data obtained with the Atacama Large Millimeter/submillimeter Array (ALMA), three candidates are identified by their extremely red $K - [3.6]$ colors as well as by non-detection in X-ray, optical, far-infrared (FIR), and radio bands. The non-detection in the deep ALMA observations suggests that they are not dusty galaxies but BBGs at $z \sim 6$, although contamination from Active Galactic Nuclei (AGNs) at $z \sim 0$ cannot be completely ruled out for the moment. Our spectral energy distribution (SED) analyses reveal that the BBG candidates at $z \sim 6$ have stellar masses of $\approx 5 \times 10^{10} M_{\odot}$ dominated by old stellar populations with ages of $\gtrsim 700$ Myr. Assuming that all the three candidates are real BBGs at $z \sim 6$, we estimate the stellar mass density (SMD) to be $2.4^{+2.3}_{-1.3} \times 10^{4} M_{\odot}$ Mpc$^{-3}$. This is consistent with an extrapolation from the lower redshift measurements. The onset of star formation in the three BBG candidates is expected to be several hundred million years before the observed epoch of $z \sim 6$. We estimate the star-formation rate density (SFRD) contributed by progenitors of the BBGs to be 2.4 -- 12 $\times 10^{-5} M_{\odot}$ yr$^{-1} $Mpc$^{-3}$ at $z > 14$ (99.7\% confidence range). Our result suggests a smooth evolution of the SFRD beyond $z = 8$.

Photometric observations of the spectroscopically confirmed $z\approx 9.1$ galaxy MACS1149-JD1 have indicated the presence of a prominent Balmer break in its spectral energy distribution, which may be interpreted as due to very large fluctuations in its past star formation activity. In this paper, we investigate to what extent contemporary simulations of high-redshift galaxies produce star formation rate variations sufficiently large to reproduce the observed Balmer break of MACS1149-JD1. We find that several independent galaxy simulations are unable to account for Balmer breaks of the inferred size, suggesting that MACS1149-JD1 either must be a very rare type of object or that our simulations are missing some key ingredient. We present predictions of spectroscopic Balmer break strength distributions for $z\approx 7-9$ galaxies that may be tested through observations with the upcoming James Webb Space Telescope and also discuss the impact that various assumptions on dust reddening, Lyman continuum leakage and deviations from a standard stellar initial mass function would have on the results.

With the Atacama Large Millimeter/sub-millimeter Array (ALMA), we report detections of the far-infrared (FIR) [OIII] 88 $\mu$m line and the underlying dust continuum in the two quasars in the reionization epoch, J205406.48-000514.8 (hereafter J2054-0005) at $z=6.0391\pm0.0002$ and J231038.88+185519.7 (hereafter J2310+1855) at $z=6.0035\pm0.0007$. The [OIII] luminosity of J2054-0005 and J2310+1855 are $L_{\rm [OIII]} = (6.8\pm0.6) \times 10^{9}$ and $(2.4\pm0.6) \times 10^{9}$ $L_{\rm \odot}$, corresponding to $\approx$ 0.05\%\ and 0.01\%\ of the total infrared luminosity, $L_{\rm TIR}$, respectively. Combining these [OIII] luminosities with [CII] 158 $\mu$m luminosities in the literature, we find that J2054-0005 and J2310+1855 have the [OIII]-to-[CII] luminosity ratio of $2.1\pm0.4$ and $0.3\pm0.1$, respectively, the latter of which is the lowest among objects so far reported at $z>6$. Combining [OIII] observations at $z\approx6-9$ from our study and the literature, we identify the [OIII] line deficit: objects with larger $L_{\rm TIR}$ have lower $L_{\rm [OIII]}$-to-$L_{\rm TIR}$ ratios. Furthermore, we also find that the anti-correlation is shifted toward higher $L_{\rm TIR}$ value when compared to the local [OIII] line deficit.

We report results of a search for galaxies at z>3 with Lyman continuum (LyC) emission using a narrow-band filter NB359 with Subaru / Suprime-Cam in a ~800 arcmin$^2$ blank field around the GOODS-N. We use 103 star-forming galaxies (SFGs) and 8 AGNs with spectroscopic redshifts in a range between 3.06 and 3.5, and 157 photometrically selected z=3.1 Lyman $\alpha$ emitter (LAE) candidates as the targets. After removing galaxies spectroscopically confirmed to be contaminated by foreground sources, we found two SFGs and one AGN as candidate LyC emitting sources among the targets with spectroscopic redshifts. Among LAE candidates, five sources are detected in the NB359 image, and three among them may be contaminated by foreground sources. We compare the sample galaxies in the GOODS-N with those in the SSA22, where a prominent protocluster at z=3.1 is known and the LyC search using the same NB359 filter has been made. Frequency of galaxies with LyC leakage in the SSA22 field may be about two times higher than the galaxies in the GOODS-N with the sample UV magnitude range, although the numbers of LyC detections in these fields are too small to make a statistically significant conclusion. By combining the sample galaxies in these fields, we place the 3$\sigma$ upper limits of the observed LyC-to-UV flux density ratio and LyC escape fraction for galaxies at z=3.1 with absolute UV magnitude $M_{UV} < -18.8$ as ($f_{LyC}/f_{UV})_{obs} < 0.036$ and $f_{esc}^{abs}<8$\%, respectively.

We report Atacama Large Millimetre/submillimeter Array (ALMA) Cycle-5 follow-up observations of two candidate [CII] emitters at z ~ 6 in the ALMA deep field in SSA22 (ADF22). The candidates were detected blindly in a Cycle-2 ALMA survey covering ~ 5 square arcmins, with a single tuning, along with two CO lines associated with galaxies at lower redshifts. Various tests suggested at least one of the two > 6-sigma [CII] candidates should be robust (Hayatsu et al. 2017). Nevertheless, our new, deeper observations recover neither candidate, demonstrating a higher contamination rate than expected. The cause of the spurious detections is under investigation but at present it remains unclarified.

State-of-the-art rest-frame UV and FIR photometric and spectroscopic observations are now pushing the redshift frontiers of galaxy formation studies up to $z\sim9-11$ and beyond. Recent HST observations unveiled the presence of a star-forming galaxy exhibiting the Lyman break at $\lambda_{\rm obs}=1.47\pm0.01$ $\mu$m, i.e., a $z=11.09^{+0.08}_{-0.12}$ galaxy with a stellar mass of $\sim10^9 M_\odot$, demonstrating that galaxy build-up was well underway early in the epoch of reionization (EoR) at $z>10$. Targeted spectroscopy of a lensed Lyman break galaxy uncovers the earliest metals known to date up to $z=9.1096\pm0.0006$ by detecting the bright [OIII] 88~$\mu$m nebular line, indicating the onset of star formation 250 million years after the Big Bang, i.e., corresponding to a redshift of $z\sim15$. These latest findings lead us to a number of key questions: How and when metal enrichment happened in the EoR? What was the nature of the earliest-epoch star-forming galaxies at $z=10-15$? What was the spatial distribution of such galaxies, and what was the relation to the putative large-scale ionization bubbles during the EoR? What were the dark-halo masses of such earliest-epoch star-forming galaxies? To address all these questions, we need to uncover a statistically large number of $z=10-15$ galaxies in the pre-reionization era. Here we argue two possible pathways: (1) a wide-area, sensitive blind spectroscopic survey of [OIII] 88 $\mu$m line-emitting galaxies at submillimeter wavelengths, and (2) an ultra-wide-area, high-cadence photometric survey of transient sources at radio-to-(sub)millimeter wavelengths, together with the immediate follow-up spectroscopy with an ultra-wide-band spectrograph, to catch the pop-III $\gamma$-ray bursts.

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

The scatter in the relationship between the strength of [CII] 158$\mu$m emission and the star formation rate at high-redshift has been the source of much recent interest. Although the relationship is well-established locally, several intensely star-forming galaxies have been found whose [CII] 158$\mu$m emission is either weak, absent or spatially offset from the young stars. Here we present new ALMA data for the two most distant, gravitationally-lensed and spectroscopically-confirmed galaxies, A2744\_YD4 at $z=$8.38 and MACS1149\_JD1 at $z=$9.11, both of which reveal intense [OIII] 88$\mu$m emission. In both cases we provide stringent upper limits on the presence of [CII] 158$\mu$m with respect to [OIII] 88$\mu$m. We review possible explanations for this apparent redshift-dependent [CII] deficit in the context of our recent hydrodynamical simulations. Our results highlight the importance of using several emission line diagnostics with ALMA to investigate the nature of the interstellar medium in early galaxies.

We present results of the cross-correlation Lya intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Lya emitters (LAEs) at $z=5.7$ and $6.6$ in a total area of $4$ deg$^2$. Although overwhelming amount of data quality controls have been performed for the narrowband images and the LAE samples, we further conduct extensive analysis evaluating systematics of large-scale point spread function wings, sky subtractions, and unknown errors on the basis of physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Lya intensity of the narrowband images and the LAEs. We identify very diffuse Lya emission with the $3\sigma$ ($2\sigma$) significance at $>150$ comoving kpc (ckpc) far from the LAEs at $z = 5.7$ ($6.6$), beyond a virial radius of star-forming galaxies with $M_\mathrm{h} \sim 10^{11} M_\odot$. The diffuse Lya emission possibly extends up to $1,000$ ckpc with the surface brightness of $10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We confirm that the small-scale ($<150$ ckpc) Lya radial profiles of LAEs in our Lya intensity maps are consistent with those obtained by recent MUSE observations (Leclercq et al. 2017). Comparisons with numerical simulations suggest that the large-scale ($\sim 150-1,000$ ckpc) Lya emission are not explained by unresolved faint sources of neighboring galaxies including satellites, but by a combination of Lya photons emitted from the central LAE and other unknown sources, such as a cold-gas stream and galactic outflow. We find no evolution in the Lya radial profiles of our LAEs from $z=5.7$ to $6.6$, where theoretical models predict a flattening of the profile slope made by cosmic reionization, albeit with our moderately large observational errors.

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

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 present results from the LymAn Continuum Escape Survey (LACES), a Hubble Space Telescope (HST) program designed to characterize the ionizing radiation emerging from a sample of Lyman alpha emitting galaxies at redshift $z\simeq 3.1$. As many show intense [O III] emission characteristic of $z>6.5$ star-forming galaxies, they may represent valuable low redshift analogs of galaxies in the reionization era. Using HST Wide Field Camera 3 / UVIS $F336W$ to image Lyman continuum emission, we investigate the escape fraction of ionizing photons in this sample. For 61 sources, of which 77% are spectroscopically confirmed and 53 have measures of [O III] emission, we detect Lyman continuum leakage in 20%, a rate significantly higher than is seen in individual continuum-selected Lyman break galaxies. We estimate there is a 98% probability that $\leq 2$ of our detections could be affected by foreground contamination. Fitting multi-band spectral energy distributions (SEDs) to take account of the varying stellar populations, dust extinctions and metallicities, we derive individual Lyman continuum escape fractions corrected for foreground intergalactic absorption. We find escape fractions of 15 to 60% for individual objects, and infer an average 20% escape fraction by fitting composite SEDs for our detected samples. Surprisingly however, even a deep stack of those sources with no individual $F336W$ detections provides a stringent upper limit on the average escape fraction of less than 0.5%. We examine various correlations with source properties and discuss the implications in the context of the popular picture that cosmic reionization is driven by such compact, low metallicity star-forming galaxies.

The origin of dust in galaxy halos or in the circum-galactic medium (CGM) is still a mystery. We investigate if the radiation pressure in high-redshift ($z\sim 10$) galaxies can efficiently transport dust to halos. To clarify the first dust enrichment of galaxy halos in the early Universe, we solve the motion of a dust grain considering radiation pressure, gas drag, and gravity in the vertical direction of the galactic disc. Radiation pressure is estimated in a consistent manner with the stellar spectra and dust extinction. As a consequence, we find that dust grains with radii $a\sim 0.1~\mu$m successfully escape from the galactic disc if the ongoing star formation episode converts more than 15 per cent of the baryon content into stars and lasts $\gtrsim 30$ Myr, while larger and smaller grains are trapped in the disc because of gravity and gas drag, respectively. We also show that grain charge significantly enhances gas drag at a few--10 scale heights of the galactic disc, where the grain velocities are suppressed to $\sim 1$ km s$^{-1}$. There is an optimum dust-to-gas ratio ($\sim 10^{-3}$) in the galactic disc and an optimum virial mass $\sim 10^{10}$--$10^{11}$ M$_{\odot}$ for the transport of $a\sim 0.1~\mu$m grains to the halo. We conclude that early dust enrichment of galaxy halos at $z\gtrsim 10$ is important for the origin of dust in the CGM.

We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $\mu$m line and rest-frame 90 $\mu$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $\mu$m flux density of $137 \pm 26$ $\mu$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $\beta = 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($\tau_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $\tau_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $\tau_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass.

In the SSA22 field which exhibits a large-scale proto-cluster at $z=3.1$, we carried out a spectroscopic survey for Lyman Break Galaxies (LBGs) with the VLT/VIMOS and identified 78 confident LBGs at $z=2.5$--4. We stacked their spectra in the observer's frame by using a sophisticated method. Analyzing the composite spectrum, we have revealed that the large-scale proto-cluster at $z=3.1$ has a strong HI absorption dip of rest-frame equivalent width of $-1.7$ A. Another strong absorption dip found at $z=3.28$ is associated with a modestly high-density LBG peak, similar to that at $z=3.1$. We have also detected an HI transparency peak at $z=2.98$ in the composite spectrum, coincident with a void in the LBG distribution. In this paper, we also investigated the relation between LBGs, HI gas and AGNs at $z=3$--4 in the SSA22 field. Two AGNs at $z=3.353$ and 3.801 are, respectively, associated with the LBG concentration of an overdensity factor $\delta_{\rm LBG}\simeq2$ in the present statistics. Another structure at $z=3.453$ is remarkable: 20 comoving Mpc scale dense HI gas which is not associated with any apparent LBG overdensity but involving a pair of AGNs. Such structure may be a new type of the AGN-matter correlation. If the inhomogeneous structures over a comoving Gpc scale found in this paper are confirmed with sufficient statistics in the future, the SSA22 field will become a key region to test the standard cold dark matter structure formation scenario.

We present an analysis to disentangle the connection between physical quantities that characterize the conditions of ionized HII regions -- metallicity ($Z$), ionization parameter ($U$), and electron density ($n_\mathrm{e}$) -- and the global stellar mass ($M_\ast$) and specific star formation rate ($\mathrm{sSFR}=\mathrm{SFR}/M_\ast$) of the host galaxies. We construct composite spectra of galaxies at $0.027 \le z \le 0.25$ from Sloan Digital Sky Survey, separating the sample into bins of $M_\ast$ and sSFR, and estimate the nebular conditions from the emission line flux ratios. Specially, metallicity is estimated from the direct method based on the faint auroral lines [OIII]$\lambda$4363 and [OII]$\lambda\lambda$7320,7330. The metallicity estimates cover a wide range from $12+\log\mathrm{O/H}\sim7.6\textrm{--}8.9$. It is found that these three nebular parameters all are tightly correlated with the location in the $M_\ast$--sSFR plane. With simple physically-motivated ans\"atze, we derive scaling relations between these physical quantities by performing multi regression analysis. In particular, we find that $U$ is primarily controlled by sSFR, as $U \propto \mathrm{sSFR}^{0.43}$, but also depends significantly on both $Z$ and $n_\mathrm{e}$. The derived partial dependence of $U \propto Z^{-0.36}$ is weaker than the apparent correlation ($U\propto Z^{-1.52}$). The remaining negative dependence of $U$ on $n_\mathrm{e}$ is found to be $U \propto n_\mathrm{e}^{-0.29}$. The scaling relations we derived are in agreement with predictions from theoretical models and observations of each aspect of the link between these quantities. Our results provide a useful set of equations to predict the nebular conditions and emission-line fluxes of galaxies in semi-analytic models.

Context. Measuring how the physical properties of galaxies change across cosmic times is essential to understand galaxy formation and evolution. With the advent of numerous ground-based and space-borne instruments launched over the past few decades we now have exquisite multi-wavelength observations of galaxies from the FUV to the radio domain. To tap into this mine of data and obtain new insight into the formation and evolution of galaxies, it is essential that we are able to extract information from their SED. Aims. We present a completely new implementation of CIGALE. Written in python, its main aims are to easily and efficiently model the FUV to radio spectrum of galaxies and estimate their physical properties such as star formation rate, attenuation, dust luminosity, stellar mass, and many other physical quantities. Methods. To compute the spectral models, CIGALE builds composite stellar populations from simple stellar populations combined with highly flexible star formation histories, calculates the emission from gas ionised by massive stars, and attenuates both the stars and the ionised gas with a highly flexible attenuation curve. Based on an energy balance principle, the absorbed energy is then re-emitted by the dust in the mid- and far-infrared domains while thermal and non-thermal components are also included, extending the spectrum far into the radio range. A large grid of models is then fitted to the data and the physical properties are estimated through the analysis of the likelihood distribution. Results. CIGALE is a versatile and easy-to-use tool that makes full use of the architecture of multi-core computers, building grids of millions of models and analysing samples of thousands of galaxies, both at high speed. Beyond fitting the SEDs of galaxies and parameter estimations, it can also be used as a model-generation tool or serve as a library to build new applications.

We present the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch.

Recent observations with the Atacama Large Millimeter/submillimeter Array (ALMA) detected far-infrared emission lines such as the [OIII] 88 \mu m line from galaxies at $z \sim 7 - 9$. Far-infrared lines can be used to probe the structure and kinematics of such high-redshift galaxies as well as to accurately determine their spectroscopic redshifts. We use a cosmological simulation of galaxy formation to study the physical properties of [OIII] 88 \mu m emitters. In a comoving volume of 50 $h^{-1}$ Mpc on a side, we locate 34 galaxies with stellar masses greater than $10^8\ {\rm M_\odot}$ at $z = 9$, and more than 270 such galaxies at $z = 7$. We calculate the [OIII] 88 \mu m luminosities ($L_{\rm OIII}$) by combining a physical model of HII regions with emission line calculations using the photoionization code CLOUDY. We show that the resulting $L_{\rm OIII}$, for a given star formation rate, is slightly higher than predicted from the empirical relation for local galaxies, and is consistent with recent observations of galaxies at redshifts 7 - 9. Bright [OIII] emitters with $L_{\rm OIII} > 10^8 {\rm L_\odot}$ have stellar masses greater than $10^9\ {\rm M_\odot}$, star formation rates higher than $3\ {\rm M_\odot\ yr}^{-1}$, and the typical metallicity is $\sim 0.1\ {\rm Z_\odot}$. The galaxies are hosted by dark matter halos with masses greater than $10^{10.5}\ {\rm M_\odot}$. Massive galaxies show characteristic structure where the [OIII] emitting gas largely overlaps with young stars, but the emission peak is separated from the main stellar population, suggesting the stochastic and localized nature of star formation in the first galaxies. We propose to use the [OIII] 5007 \AA\ line, to be detected by James Webb Space Telescope (JWST), to study the properties of galaxies whose [OIII] 88 \mu m line emission has been already detected with ALMA.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

We investigate Lya, [OIII]5007, Ha, and [CII]158um emission from 1124 galaxies at z=4.9-7.0. Our sample is composed of 1092 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 34 galaxies at z=5.148-7.508 with deep ALMA [CII]158um data in the literature. Fluxes of strong rest-frame optical lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, we find that the rest-frame Ha equivalent width and the Lya escape fraction f_Lya positively correlate with the rest-frame Lya equivalent width EW^0_Lya. The f_Lya-EW^0_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5 significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW^0_Lya~0-30 A, and then decreases out to EW^0_Lya~130 A. We also identify an anti-correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW^0_Lya at the >99% confidence level. We carefully investigate physical origins of the correlations with stellar-synthesis and photoionization models, and find that a simple anti-correlation between EW_Lya^0 and metallicity explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study, indicating detections of metal-poor (~0.03 Zo) galaxies with EW^0_Lya~200 A.

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

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We report fourteen and twenty-eight protocluster candidates at z=5.7 and 6.6 over 14 and 19 deg^2 areas, respectively, selected from 2,230 (259) Lya emitters (LAEs) photometrically (spectroscopically) identified with Subaru/Hyper Suprime-Cam (HSC) deep images (Keck, Subaru, and Magellan spectra and the literature data). Six out of the 42 protocluster candidates include 1-12 spectroscopically confirmed LAEs at redshifts up to z=6.574. By the comparisons with the cosmological Lya radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates are progenitors of the present-day clusters with a mass of > 10^14 M_sun. We then investigate the correlation between LAE overdensity delta and Lya rest-frame equivalent width EW_Lya^rest, because the cosmological Lya RT model suggests that a slope of EW_Lya^rest-delta relation is steepened towards the epoch of cosmic reionization (EoR), due to the existence of the ionized bubbles around galaxy overdensities easing the escape of Lya emission from the partly neutral intergalactic medium (IGM). The available HSC data suggest that the slope of the EW_Lya^rest-delta correlation does not evolve from the post-reionization epoch z=5.7 to the EoR z=6.6 beyond the moderately large statistical errors. There is a possibility that we would detect the evolution of the EW_Lya^rest - delta relation from z=5.7 to 7.3 by the upcoming HSC observations providing large samples of LAEs at z=6.6-7.3.

We present observations of escaping Lyman Continuum (LyC) radiation from 34 massive star-forming galaxies and 12 weak AGN with reliably measured spectroscopic redshifts at $z$$\simeq$2.3-4.1. We analyzed Hubble Space Telescope Wide Field Camera 3 (WFC3) mosaics of the Early Release Science field in three UVIS filters to sample the rest-frame LyC over this redshift range. With our best current assessment of the WFC3 systematics, we provide 1$\sigma$ upper limits for the average LyC emission of galaxies at $z$=2.35, 2.75, and 3.60 to $\sim$28.5, 28.1, and 30.7 mag in image stacks of 11-15 galaxies in the WFC3/UVIS F225W, F275W, and F336W, respectively. The LyC flux of weak AGN at $z$=2.62 and 3.32 are detected at 28.3 and 27.4 mag with SNRs of $\sim$2.7 and 2.5 in F275W and F336W for stacks of 7 and 3 AGN, respectively, while AGN at $z$=2.37 are constrained to $\gtrsim$27.9 mag at 1$\sigma$ in a stack of 2 AGN. The stacked AGN LyC light profiles are flatter than their corresponding non-ionizing UV continuum profiles out to radii of r$\lesssim$0."9, which may indicate a radial dependence of porosity in the ISM. With synthetic stellar SEDs fit to UV continuum measurements longwards of Ly$\alpha$ and IGM transmission models, we constrain the absolute LyC escape fractions to $f_{\rm esc}^{\rm abs}$$\simeq$$22^{+44}_{-22}$% at $z$=2.35 and $\lesssim$55% at $z$=2.75 and 3.60, respectively. All available data for galaxies, including published work, suggests a more sudden increase of $f_{\rm esc}$ with redshift at $z$$\simeq$2. Dust accumulating in (massive) galaxies over cosmic time correlates with increased HI column density, which may lead to reducing $f_{\rm esc}$ more suddenly at $z$$\lesssim$2. This may suggest that star-forming galaxies collectively contributed to maintaining cosmic reionization at redshifts $z$$\gtrsim$2-4, while AGN likely dominated reionization at $z$$\lesssim$2.

We investigate the X-ray active galactic nucleus (AGN) properties of millimeter galaxies in the Great Observatories Origins Deep Survey South (GOODS-S) field detected with the Atacama Large Millimeter/submillimeter Array (ALMA), by utilizing the Chandra 7-Ms data, the deepest X-ray survey to date. Our millimeter galaxy sample comes from the ASAGAO survey covering 26 arcmin$^2$ (12 sources at a 1.2-mm flux-density limit of $\approx$0.6 mJy), supplemented by the deeper but narrower 1.3-mm survey of a part of the ASAGAO field by Dunlop et al.\ (2017). Fourteen out of the total 25 millimeter galaxies have Chandra counterparts. The observed AGN fractions at $z=1.5-3$ is found to be 90$^{+8}_{-19}$\% and $57^{+23}_{-25}$\% for the ultra/luminous infrared galaxies with $\log L_{\rm IR}/L_{\odot} = 12-12.8$ and $\log L_{\rm IR}/L_{\odot} = 11.5-12$, respectively. The majority ($\sim$2/3) of the ALMA and/or Herschel detected X-ray AGNs at $z=1.5-3$ appear to be star-formation dominant populations, having $L_{\rm X}$/ $L_{\rm IR}$ ratios smaller than the "simultaneous evolution" value expected from the local black-hole mass to stellar mass ($M_{\rm BH}$-$M_*$) relation. On the basis of the $L_{\rm X}$ and stellar mass relation, we infer that a large fraction of star-forming galaxies at $z=1.5-3$ have black hole masses smaller than those expected from the local $M_{\rm BH}$-$M_*$ relation. This is opposite to previous reports on luminous AGN at same redshifts detected in wider and shallower surveys, which are subject to selection biases against lower luminosity AGN. Our results are consistent with an evolutionary scenario that star formation occurs first, and an AGN-dominant phase follows later, in objects finally evolving into galaxies with classical bulges.

We investigate the impact of radiation pressure on spatial dust distribution inside H$_\mathrm{II}$ regions using one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. In order to investigate grain size effects as well, we introduce two additional fluid components describing large and small dust grains in the simulations. Relative velocity between dust and gas strongly depends on the drag force. We include collisional drag force and coulomb drag force. We find that, in a compact H$_\mathrm{II}$ region, a dust cavity region is formed by radiation pressure. Resulting dust cavity sizes (~0.2 pc) agree with observational estimates reasonably well. Since dust inside an H$_\mathrm{II}$ region is strongly charged, relative velocity between dust and gas is mainly determined by the coulomb drag force. Strength of the coulomb drag force is about 2-order of magnitude larger than that of the collisional drag force. In addition, in a cloud of mass $10^5$ $M_{\odot}$, we find that the radiation pressure changes the grain size distribution inside H$_\mathrm{II}$ regions. Since large (0.1 $\mathrm{\mu m}$) dust grains are accelerated more efficiently than small (0.01 $\mathrm{\mu m}$) grains, the large to small grain mass ratio becomes smaller by an order of magnitude compared with the initial one. Resulting dust size distributions depend on the luminosity of the radiation source. The large and small grain segregation becomes weaker when we assume stronger radiation source, since dust grain charges become larger under stronger radiation and hence coulomb drag force becomes stronger.

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.

We report the identification of four millimeter line emitting galaxies with the Atacama Large Milli/submillimeter Array (ALMA) in SSA22 Field (ADF22). We analyze the ALMA 1.1 mm survey data, with an effective survey area of 5 arcmin$^2$, a frequency range of 253.1--256.8 and 269.1--272.8 GHz, angular resolution of 0".7 and RMS noise of 0.8 mJy beam$^{-1}$ at 36 km s$^{-1}$ velocity resolution. We detect four line emitter candidates with significance levels above $6 \sigma$. We identify one of the four sources as a CO(9-8) emitter at $z = 3.1$ in a member of the proto-cluster known in this field. Another line emitter with an optical counterpart is likely a CO(4-3) emitter at $z = 0.7$. The other two sources without any millimeter continuum or optical/near-infrared counterpart are likely to be [CII] emitter candidates at $z = 6.0$ and $6.5$. The equivalent widths of the [CII] candidates are consistent with those of confirmed high-redshift [CII] emitters and candidates, and are a factor of 10 times larger than that of the CO(9-8) emitter detected in this search. The [CII] luminosity of the candidates are $4-7 \times 10^8~\rm L_\odot$. The star formation rates (SFRs) of these sources are estimated to be $10-20~\rm M_{\odot}~yr^{-1}$ if we adopt an empirical [CII] luminosity - SFR relation. One of them has a relatively low-S/N ratio, but shows features characteristic of emission lines. Assuming that at least one of the two candidates is a [CII] emitter, we derive a lower limit of [CII]-based star formation rate density (SFRD) at $z~\sim~6$. The resulting value of $> 10^{-2}$ $\rm M_\odot yr^{-1} Mpc^{-3}$ is consistent with the dust-uncorrected UV-based SFRD. Future millimeter/submillimeter surveys can be used to detect a number of high redshift line emitters, with which to study the star formation history in the early Universe.

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {\AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

Using galaxies as background light sources to map intervening Lya absorption is a novel approach to study the interplay among galaxies, the circum-galactic medium (CGM), and the intergalactic medium (IGM). Introducing a new measure of z = 3.1$ HI Lya absorption relative to the cosmic mean, Delta_NB497, estimated from photometric data of star-forming galaxies at 3.3 < z < 3.5, we have made two-dimensional Delta_NB497 maps in the z = 3.1 SSA22 proto-cluster region and two control fields (SXDS and GOODS-N fields) with a spatial resolution of ~ 5 comoving Mpc. The Delta_NB497 measurements in the SSA22 field are systematically larger than those in the control fields, and this HI absorption enhancement extends more than 50 comoving Mpc. The field-averaged (i.e., ~50 comoving Mpc scale) Delta_NB497 and the overdensity of Lya emitters (LAEs) seem to be correlated, while there is no clear dependency of the Delta_NB497 on the local LAE overdensity in a few comoving Mpc scale. These results suggest that diffuse HI gas spreads out in/around the SSA22 proto-cluster. We have also found an enhancement of Delta_NB497 at a projected distance < 100 physical kpc from the nearest z = 3.1 galaxies at least in the SSA22 field, which is probably due to HI gas associated with the CGM of individual galaxies. The HI absorption enhancement in the CGM-scale tends to be weaker around galaxies with stronger Lya emission, which suggests that the Lya escape fraction from galaxies depends on hydrogen neutrality in the CGM.

We present spatial correlations of galaxies and IGM HI in the COSMOS/UltraVISTA 1.62 deg$^2$ field. Our data consist of 13,415 photo-$z$ galaxies at $z\sim2-3$ with $K_s<23.4$ and the Ly$\alpha$ forest absorptions in the background quasar spectra selected from SDSS data with no signature of damped Ly$\alpha$ system contamination. We estimate a galaxy overdensity $\delta_{gal}$ in an impact parameter of 2.5 pMpc, and calculate the Ly$\alpha$ forest fluctuations $\delta_{\langle F\rangle}$ whose negative values correspond to the strong Ly$\alpha$ forest absorptions. We identify weak evidence of an anti-correlation between $\delta_{gal}$ and $\delta_{\langle F\rangle}$ with a Spearman's rank correlation coefficient of $-0.39$ suggesting that the galaxy overdensities and the Ly$\alpha$ forest absorptions positively correlate in space at the $\sim90\%$ confidence level. This positive correlation indicates that high-$z$ galaxies exist around an excess of HI gas in the Ly$\alpha$ forest. We find four cosmic volumes, dubbed $A_{obs}$-$D_{obs}$, that have extremely large (small) values of $\delta_{gal} \simeq0.8$ ($-1$) and $\delta_{\langle F\rangle}$ $\simeq0.1$ ($-0.4$), three out of which, $B_{obs}$-$D_{obs}$, significantly depart from the correlation, and weaken the correlation signal. We perform cosmological hydrodynamical simulations, and compare with our observational results. Our simulations reproduce the correlation, agreeing with the observational results. Moreover, our simulations have model counterparts of $A_{obs}$-$D_{obs}$, and suggest that the observations pinpoint, by chance, a galaxy overdensity like a proto-cluster, gas filaments lying on the sightline, a large void, and orthogonal low-density filaments. Our simulations indicate that the significant departures of $B_{obs}$-$D_{obs}$ are produced by the filamentary large-scale structures and the observation sightline effects.

We present the largest to date sample of hydrogen Lyman continuum (LyC) emitting galaxy candidates at any redshift, with 18 Ly$\alpha$ Emitters (LAEs) and 7 Lyman Break Galaxies (LBGs), obtained from the SSA22 field with Subaru/Suprime-Cam. The sample is based on the 159 LAEs and 136 LBGs observed in the field, all with spectroscopically confirmed redshifts, and these LyC candidates are selected as galaxies with counterpart in a narrow-band filter image which traces LyC at $z\geq 3.06$. Many LyC candidates show a spatial offset between the rest-frame non-ionizing ultraviolet (UV) detection and the LyC-emitting substructure or between the Ly$\alpha$ emission and LyC. The possibility of foreground contamination complicates the analysis of the nature of LyC emitters, although statistically it is highly unlikely that all candidates in our sample are contaminated by foreground sources. Many viable LyC LAE candidates have flux density ratios inconsistent with standard models, while also having too blue UV slopes to be foreground contaminants. Stacking reveals no significant LyC detection, suggesting that there is a dearth of objects with marginal LyC signal strength, perhaps due to a bimodality in the LyC emission. The foreground contamination-corrected $3\sigma$ upper limits of the observed average flux density ratios are $f_{LyC}/f_{UV}<0.08$ from stacking LAEs and $f_{LyC}/f_{UV}<0.02$ from stacking LBGs. There is a sign of a positive correlation between LyC and Ly$\alpha$, suggesting that both types of photons escape via a similar mechanism. The LyC detection rate among proto-cluster LBGs is seemingly lower compared to the field.

We present Keck/MOSFIRE spectra of the diagnostic nebular emission lines [OIII]5007,4959, [OII]3727, and H-beta for a sample of 15 redshift z=3.1-3.7 Ly-alpha emitters (LAEs) and Lyman break galaxies (LBGs). In conjunction with spectra from other surveys, we confirm earlier indications that LAEs have a much higher [OIII]/[OII] line ratio than is seen in similar redshift LBGs. By comparing their distributions on a [OIII]/[OII] versus R23 diagram, we demonstrate that this difference cannot arise solely because of their lower metallicities but most likely is due to a harder ionizing spectrum. Using measures of H-beta and recombination theory, we demonstrate, for a subset of our LAEs, that xi_ion - the number of Lyman continuum photons per UV luminosity - is indeed 0.2-0.5 dex larger than for typical LBGs at similar redshifts. Using photoionization models we estimate the effect this would have on both [OIII]/[OII] and R23 and conclude such a hard spectrum can only partially explain such intense line emission. The additional possibility is that such a large [OIII]/[OII] ratio is in part due to density rather than ionization bound nebular regions, which would imply a high escape fraction of ionizing photons. We discuss how further observations could confirm this possibility. Clearly LAEs with intense [OIII] emission represent a promising analog of those z>7 sources with similarly strong lines which are thought to be an important contributor to cosmic reionization.

We present new analysis of multi-epoch, H-band, scattered light images of the AB Aur system. We used a Monte Carlo, radiative transfer code to simultaneously model the system's SED and H-band polarized intensity imagery. We find that a disk-dominated model, as opposed to one that is envelope dominated, can plausibly reproduce AB Aur's SED and near-IR imagery. This is consistent with previous modeling attempts presented in the literature and supports the idea that at least a subset of AB Aur's spirals originate within the disk. In light of this, we also analyzed the movement of spiral structures in multi-epoch H-band total light and polarized intensity imagery of the disk. We detect no significant rotation or change in spatial location of the spiral structures in these data, which span a 5.8 year baseline. If such structures are caused by disk-planet interactions, the lack of observed rotation constrains the location of the orbit of planetary perturbers to be >47 AU.

Subaru/SuprimeCam narrowband photometry of the SSA22 field reveals the presence of four Lyman continuum (LyC) candidates among a sample of 14 AGN. Two show offsets and likely have stellar LyC in nature or are foreground contaminants. The remaining two LyC candidates are Type I AGN. We argue that the average LyC escape fraction of high redshift low luminosity AGN is not likely to be unity, as often assumed in the literature. From direct measurement we obtain the average LyC-to-UV flux density ratio and ionizing emissivity for a number of AGN classes and find it at least a factor of two lower than values obtained assuming f_esc = 1. Comparing to recent Ly{\alpha} forest measurements, AGNs at redshift z~3 make up at most ~12% and as little as ~5% of the total ionizing budget. Our results suggest that AGNs are unlikely to dominate the ionization budget of the Universe at high redshifts.

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

We have performed very large and high resolution cosmological hydrodynamic simulations in order to investigate detectability of nebular lines in the rest-frame ultraviolet (UV) to optical wavelength range from galaxies at $z > 7$. We find that the expected line fluxes are very well correlated with apparent UV magnitudes. The C IV $1549 \rm \AA$ and the C III] $1909 \rm \AA$ lines of galaxies brighter than $26~ \rm AB$ are detectable with current facilities such as the VLT/XShooter and the Keck/MOSFIRE. Metal lines such as C IV $1549 \rm \AA$, C III] $1909 \rm \AA$, [O II] $3727 \rm \AA$ and [O III] $4959/5007 \rm \AA$ can be good targets for the spectroscopic observation with Thirty Meter Telescope (TMT), European Extremely Large Telescope (E-ELT), Giant Magellan Telescope (GMT) and James Webb Space Telescope (JWST). H$\alpha$ and H$\beta$ lines are also expected to be detectable with these telescopes. We also predict detectability of the nebular lines for $z > 10$ galaxies that will be found with JWST, Wide-Field Infrared Survey Telescope (WFIRST) and First Light And Reionization Explorer (FLARE) $(11 \leq z \leq 15$). We conclude that the C IV $1549 \rm \AA$, C III] $1909 \rm \AA$, [O III] $4959/5007 \rm \AA$ and H$\beta$ lines even from $z \sim 15$ galaxies can be strong targets for TMT, ELT and JWST. We also find that the magnification by gravitational lensing is of great help to detect such high-$z$ galaxies. According to our model, the C III] $1909 \rm \AA$ line in $z > 9$ galaxy candidates is detectable even using the current facilities.

We made near infrared multicolor imaging observations of a disk around Herbig Be star HD100546 using Gemini/NICI. K (2.2\,$\mu$m), H$_2$O ice (3.06\,$\mu$m), and L'(3.8\,$\mu$m) disk images were obtained and we found the 3.1\,$\mu$m absorption feature in the scattered light spectrum, likely due to water ice grains at the disk surface. We compared the observed depth of the ice absorption feature with the disk model based on \cite{Oka2012} including water ice photodesorption effect by stellar UV photons. The observed absorption depth can be explained by the both disk models with/without photodesorption effect within the measurement accuracy, but slightly favors the model with photodesorption effects, implying that the UV photons play an important role on the survival/destruction of ice grains at the Herbig Ae/Be disk surface. Further improvement on the accuracy of the observations of the water ice absorption depth is needed to constrain the disk models.

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.

Using galaxies as background light sources to map the Lya absorption lines is a novel approach to study Damped Lya Absorbers (DLAs). We report the discovery of an intervening z = 3.335 +- 0.007 DLA along a galaxy sight-line identified among 80 Lyman Break Galaxy (LBG) spectra obtained with our Very Large Telescope/Visible Multi-Object Spectrograph survey in the SSA22 field. The measured DLA neutral hydrogen (HI) column density is log (NHI/cm^{-2}) = 21.68 +- 0.17. The DLA covering fraction over the extended background LBG is > 70% (2 sigma), yielding a conservative constraint on the DLA area of > 1 kpc^2. Our search for a counterpart galaxy hosting this DLA concludes that there is no counterpart galaxy with star formation rate (SFR) larger than a few Msun yr^{-1}, ruling out an unobscured violent star formation in the DLA gas cloud. We also rule out the possibility that the host galaxy of the DLA is a passive galaxy with Mstar > 5 x 10^{10} Msun or a heavily dust-obscured galaxy with E(B-V) > 2. The DLA may coincide with a large-scale overdensity of the spectroscopic LBGs. The occurrence rate of the DLA is compatible with that of DLAs found in QSO sight-lines.

Gravitationally lensed galaxies with magnification ~10-100 are routinely detected at high redshifts, but magnifications significantly higher than this are hampered by a combination of low probability and large source sizes. Magnifications of ~1000 may nonetheless be relevant in the case of intrinsically small, high-redshift objects with very high number densities. Here, we explore the prospects of detecting compact (< 10 pc), high-redshift (z > 7) Population III star clusters at such extreme magnifications in large-area surveys with planned telescopes like Euclid, WFIRST and WISH. We find that the planned WISH 100 sq. deg ultradeep survey may be able to detect a small number of such objects, provided that the total stellar mass of these star clusters is > 10000 solar masses. If candidates for such lensed Population III star clusters are found, follow-up spectroscopy of the surrounding nebula with the James Webb Space Telescope or groundbased Extremely Large Telescopes should be able to confirm the Population III nature of these objects. Multiband photometry of these objects with the James Webb Space Telescope also has the potential to confirm that the stellar initial mass function in these Population III star clusters is top-heavy, as supported by current simulations.

We aim to investigate the effect of the escaping ionizing radiation on the color selection of high redshift galaxies and identify candidate Lyman continuum (LyC) emitters. The intergalactic medium prescription of Inoue et al.(2014) and galaxy synthesis models of Bruzual&Charlot (2003) have been used to properly treat the ultraviolet stellar emission, the stochasticity of the intergalactic transmission and mean free path in the ionizing regime. Color tracks are computed by turning on/off the escape fraction of ionizing radiation. At variance with recent studies, a careful treatment of IGM transmission leads to no significant effects on the high-redshift broad-band color selection. The decreasing mean free path of ionizing photons with increasing redshift further diminishes the contribution of the LyC to broad-band colors. We also demonstrate that prominent LyC sources can be selected under suitable conditions by calculating the probability of a null escaping ionizing radiation. The method is applied to a sample of galaxies extracted from the GOODS-S field. A known LyC source at z=3.795 is successfully recovered as a LyC emitter candidate and another convincing candidate at z=3.212 is reported. A detailed analysis of the two sources (including their variability and morphology) suggests a possible mixture of stellar and non-stellar (AGN) contribution in the ultraviolet. Conclusions: Classical broad-band color selection of 2.5<z<4.5 galaxies does not prevent the inclusion of LyC emitters in the selected samples. Large fesc in relatively bright galaxies (L>0.1L*) could be favored by the presence of a faint AGN not easily detected at any wavelength. A hybrid stellar and non-stellar (AGN) ionizing emission could coexist in these systems and explain the tensions found among the UV excess and the stellar population synthesis models reported in literature.

We exploit wide-field Ly$alpha$ imaging with Subaru to probe the environment around TN J1338-1942, a powerful radio galaxy with a >100 kpc Ly$\alpha$ halo at z=4.11. We used a sample of Ly$\alpha$ emitters (LAEs) down to $\log(L_{\rm Ly\alpha} [erg\, s^{-1}])\sim 42.8$ to measure the galaxy density around TNJ1338, compared to a control sample from a blank field taken with the same instrument. We found that TNJ1338 resides in a region with a peak overdensity of $\delta_{\rm LAE}=2.8\pm 0.5$ on scales of $8\, h^{-1}\rm Mpc$ (on the sky) and $112\, h^{-1}\rm Mpc$ (line of sight) in comoving coordinates. Adjacent to this overdensity, we found a strong underdensity where virtually no LAEs are detected. We used a semi-analytical model of LAEs derived from the Millennium Simulation to compare our results with theoretical predictions. While the theoretical density distribution is consistent with the blank field, overdense regions such as that around TNJ1338 are very rare, with a number density of $6.4\times 10^{-8}\rm Mpc^{-3}$ (comoving), corresponding to the densest < 0.4 percentile at $z\simeq 4.1$. We also found that the Ly$\alpha$ luminosity function in the TNJ1338 field differs from that in the blank field: the number of bright LAEs ($\log(L_{\rm Ly\alpha}[erg\,s^{-1}]) \gtrsim 43.3$) is enhanced, while the number of fainter LAEs is relatively suppressed. These results suggest that some powerful radio galaxies associated with Ly$\alpha$ nebulae reside in extreme overdensities on $\sim 3$--$6\, \rm Mpc$ scales, where star-formation and AGN activity may be enhanced via frequent galaxy mergers or high rates of gas accretion from the surroundings.

We present near-IR and far-UV observations of the pre-transitional (gapped) disk in HD 169142 using NASA's Infrared Telescope Facility and Hubble Space Telescope. The combination of our data along with existing data sets into the broadband spectral energy distribution reveals variability of up to 45% between ~1.5-10 {\mu}m over a maximum timescale of 10 years. All observations known to us separate into two distinct states corresponding to a high near-IR state in the pre-2000 epoch and a low state in the post-2000 epoch, indicating activity within the <1 AU region of the disk. Through analysis of the Pa {\beta} and Br {\gamma} lines in our data we derive a mass accretion rate in May 2013 of (1.5 - 2.7) x 10^-9 Msun/yr. We present a theoretical modeling analysis of the disk in HD 169142 using Monte-Carlo radiative transfer simulation software to explore the conditions and perhaps signs of planetary formation in our collection of 24 years of observations. We find that shifting the outer edge (r = 0.3 AU) of the inner disk by 0.05 AU toward the star (in simulation of accretion and/or sculpting by forming planets) successfully reproduces the shift in NIR flux. We establish that the ~40-70 AU dark ring imaged in the NIR by Quanz et al. (2013) and Momose et al. (2013) and at 7 mm by Osorio et al. (2014) may be reproduced with a 30% scaled density profile throughout the region, strengthening the link to this structure being dynamically cleared by one or more planetary mass bodies.

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

We present an updated version of the so-called Madau model for the attenuation by the intergalactic neutral hydrogen against the radiation from distant objects. First, we derive a distribution function of the intergalactic absorbers from the latest observational statistics of the Ly$\alpha$ forest, Lyman limit systems, and damped Ly$\alpha$ systems. The distribution function excellently reproduces the observed redshift evolutions of the Ly$\alpha$ depression and of the mean-free-path of the Lyman continuum simultaneously. Then, we derive a set of the analytic functions which describe the mean intergalactic attenuation curve for objects at $z>0.5$. The new model predicts less (or more) Ly$\alpha$ attenuation for $z\simeq3$--5 ($z>6$) sources through usual broad-band filters relative to the original Madau model. This may cause a systematic difference in the photometric redshift estimates, which is, however, still as small as about 0.05. Finally, we find a more than 0.5 mag overestimation of the Lyman continuum attenuation in the original Madau model at $z>3$, which causes a significant overcorrection against direct observations of the Lyman continuum of galaxies.

We have performed a large cosmological hydrodynamics simulation tailored to the deep survey with the Hubble Space Telescope made in 2012, the so-called UDF12 campaign. After making a light-cone output, we have applied the same color selection criteria as the UDF12 campaign to select galaxies from our simulation, and then, have examined the physical properties of them as a proxy of the real observed UDF12 galaxies at $z > 7$. As a result, we find that the halo mass is almost linearly proportional to the observed ultraviolet (UV) luminosity ($4 \times 10^{11}~{\rm M_{\odot}}$ at $M_{\rm UV} = -21$). The dust attenuation and UV slope $\beta$ well correlates with the observed UV luminosity, which is consistent with observations quantitatively. The star formation rate (SFR) is also linearly proportional to the stellar mass and the specific SFR shows only a weak dependency on the mass. We also find an increasing star formation history with a time-scale of $\sim100$ Myr in the high-$z$ galaxies. An average metallicity weighted by the Lyman continuum luminosity reaches up to $>0.1$ Solar even at $z \sim 10$, suggesting a rapid metal enrichment. We also expect $\geq 0.1$ mJy at 350 GHz of the dust thermal emission from the galaxies with $H_{160} \leq 27$, which can be detectable with the Atacama Large Milimetre-submilimetre Array. The galaxies selected by the UDF12 survey contribute to only $52--12\%$ of the cosmic SFR density from $z \sim 7$ to $z \sim 10$, respectively. The James Webb Space Telescope will push the detection fraction up to $77--72\%$.

We investigate the potential use of nebular emission lines in the rest-frame far-infrared (FIR) for determining spectroscopic redshift of z>8 galaxies with the Atacama Large Millimeter/sub-millimeter Array (ALMA). After making a line emissivity model as a function of metallicity, especially for the [O III] 88 micron line which is likely to be the strongest FIR line from H II regions, we predict the line fluxes from high-z galaxies based on a cosmological hydrodynamics simulation of galaxy formation. Since the metallicity of galaxies reaches at ~0.2 Zsun even at z>8 in our simulation, we expect the [O III] 88 micron line as strong as 1.3 mJy for 27 AB objects, which is detectable at a high significance by <1 hour integration with ALMA. Therefore, the [O III] 88 micron line would be the best tool to confirm the spectroscopic redshifts beyond z=8.

Current data indicate that the reionization of the Universe was complete by redshift z~6-7, and while the sources responsible for this process have yet to be identified, star-forming galaxies are often considered the most likely candidates. However, the contribution from z>6 galaxies to cosmic reionization critically depends on the fraction of ionizing (Lyman continuum, LyC) photons escaping from these objects and into the intergalactic medium. At z<4, the escaping LyC flux can be measured directly, but the opacity of the neutral intergalactic medium precludes such measurements at higher redshifts. In a recent paper, we argue that since the LyC escape fraction regulates the contribution of nebular emission to the rest-frame optical/UV spectra of galaxies, the James Webb Space Telescope should be able to indirectly assess the LyC escape fraction for galaxies at z~6-9. JWST can, on the other hand, not constrain the fraction of LyC photons directly absorbed by dust, and this is where SPICA comes in. The dust continuum emission from gravitationally lensed LyC-leakers at z=6 may in principle be detectable with SPICA, thereby constraining the level of LyC extinction in these objects.

The definition of the galactic stellar mass estimated from the spectral energy distribution is ambiguous in the literature; whether the stellar mass includes the mass of the stellar remnants, i.e. white dwarfs, neutron stars, and black holes, is not well described. The remnant mass fraction in the total (living+remnant) stellar mass of a simple stellar population monotonically increases with the age of the population, and the initial mass function and metallicity affect the increasing rate. Since galaxies are composed of a number of stellar populations, the remnant mass fraction may depend on the total stellar mass of galaxies in a complex way. As a result, the shape of the stellar mass function of galaxies may change, depending on the definition of the stellar mass. In order to explore this issue, we ran a cosmological hydrodynamical simulation, and then, we have found that the remnant mass fraction indeed correlates with the total stellar mass of galaxies. However, this correlation is weak and the remnant fraction can be regarded as a constant which depends only on the redshift. Therefore, the shape of the stellar mass function is almost unchanged, but it simply shifts horizontally if the remnant mass includes or not. The shift is larger at lower redshift and it reaches 0.2-dex at z=0 for a Chabrier IMF. Since this causes a systematic difference, we should take care of the definition of the `stellar' mass, when comparing one's result with others.

The fraction of ionizing photons (fesc) that escape from z>6 galaxies is an important parameter when assessing the role of these objects in the reionization of the Universe, but the opacity of the intergalactic medium precludes a direct measurement of fesc for individual galaxies at these epochs. We argue, that since fesc regulates the impact of nebular emission on the spectra of galaxies, it should nonetheless be possible to indirectly probe fesc well into the reionization epoch. As a first step, we demonstrate that by combining measurements of the rest-frame UV slope beta with the equivalent width of the Hb emission line, galaxies with very high Lyman continuum escape fractions (fesc>0.5) should be identifiable up to z~9 through spectroscopy with the upcoming James Webb Space Telescope (JWST). By targeting strongly lensed galaxies behind low-redshift galaxy clusters, JWST spectra of sufficiently good quality can be obtained for M(1500)<-16.0 galaxies at z~7 and for M(1500)<-17.5 galaxies at z~9. Dust-obscured star formation may complicate the analysis, but supporting observations with ALMA or the planned SPICA mission may provide useful constraints on the dust properties of these galaxies.

During the course of our Lyman continuum imaging survey, we found that the spectroscopically confirmed Lya emitter LAE 221724+001716 at z = 3.10 in the SSA 22 field shows strong Lyman continuum emission (lambda_rest ~ 900 A) that escapes from this galaxy. However, another recent spectroscopic survey revealed that the supposed Lyman continuum emission could arise from a foreground galaxy at z = 1.76 if the emission line newly detected from the galaxy at lambda_obs ~ 3360 A is Lya. If this is the case, as the angular separation between these two galaxies is very small (~ 0.6"), LAE 221724+001716 at z = 3.10 could be amplified by the gravitational lensing caused by this intervening galaxy. Here we present a possible gravitational lens model for the system of LAE 221724+001716. First, we estimate the stellar mass of the intervening galaxy as Mstar ~ 3.5x10^9 Msun from its UV luminosity and ~ 3.0x10^7--2.4x10^9 Msun through the SED fitting. Then, we find that the gravitational magnification factor ranges from 1.01 to 1.16 using the so-called singular isothermal sphere model for strong lensing. While LAE 221724+001716 is the first system of an LAE-LAE lensing reported so far, the estimated magnification factor is not so significant because the stellar mass of the intervening galaxy is small.

The cosmic star formation rate density (CSFRD) has been observationally investigated out to redshift z~10. However, most of theoretical models for galaxy formation underpredict the CSFRD at z>1. Since the theoretical models reproduce the observed luminosity functions (LFs), luminosity densities (LDs), and stellar mass density at each redshift, this inconsistency does not simply imply that theoretical models should incorporate some missing unknown physical processes in galaxy formation. Here, we examine the cause of this inconsistency in UV wavelengths by using a mock catalog of galaxies generated by a semi-analytic model of galaxy formation. We find that this inconsistency is due to two observational uncertainties: dust obscuration correction and conversion from UV luminosity to star formation rate (SFR). The methods for correction of obscuration and SFR conversion used in observational studies result in the overestimation of CSFRD by ~ 0.1-0.3 dex and ~ 0.1-0.2 dex, respectively, compared to the results obtained directly from our mock catalog. We present new empirical calibrations for dust attenuation and conversion from observed UV LFs and LDs into CSFRD.

The timing and duration of the reionization epoch is crucial to the emergence and evolution of structure in the universe. The relative roles that star-forming galaxies, active galactic nuclei and quasars play in contributing to the metagalactic ionizing background across cosmic time remains uncertain. Deep quasar counts provide insights into their role, but the potentially crucial contribution from star-formation is highly uncertain due to our poor understanding of the processes that allow ionizing radiation to escape into the intergalactic medium (IGM). The fraction of ionizing photons that escape from star-forming galaxies is a fundamental free parameter used in models to "fine-tune" the timing and duration of the reionization epoch that occurred somewhere between 13.4 and 12.7 Gyrs ago (redshifts between 12 > z > 6). However, direct observation of Lyman continuum (LyC) photons emitted below the rest frame \ion{H}{1} ionization edge at 912 \AA\ is increasingly improbable at redshifts z > 3, due to the steady increase of intervening Lyman limit systems towards high z. Thus UV and U-band optical bandpasses provide the only hope for direct, up close and in depth, observations of the types of environment that favor LyC escape. By quantifying the evolution over the past 11 billion years (z < 3) of the relationships between LyC escape and local and global parameters ..., we can provide definitive information on the LyC escape fraction that is so crucial to answering the question of, how did the universe come to be ionized? Here we provide estimates of the ionizing continuum flux emitted by "characteristic" (L_{uv}^*) star-forming galaxies as a function of look back time and escape fraction, finding that at z = 1 (7.6 Gyrs ago) L_{uv}^* galaxies with an escape fraction of 1% have a flux of 10^{-19} ergs cm^{-2} s^{-1} \AA^{-1}.

Observations with Subaru telescope have detected surprisingly strong Lyman continuum (LyC; $\sim900$ \AA\ in the rest-frame) from some Lyman $\alpha$ emitters (LAEs) at $z=3.1$. We have examined the stellar population which simultaneously accounts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet of the LAEs. As a result, we have found that stellar populations with metallicity $Z\geq1/50 Z_\odot$ can explain the observed LyC strength only with a very top-heavy initial mass function (IMF; $< m > \sim 50 M_\odot$). However, the critical metallicity for such an IMF is expected to be much lower. A very young ($\sim1$ Myr) and massive ($\sim100$ $M_\odot$) extremely metal-poor ($Z\leq5\times10^{-4}Z_\odot$) or metal-free (so-called Population III) stellar population can also reproduce the observed LyC strength if the mass fraction of such 'primordial' stellar population is $\sim1$% in total stellar mass of the LAEs.

The relation between the ratio of infrared (IR) and ultraviolet (UV) flux densities (the infrared excess: IRX) and the slope of the UV spectrum (\beta) of galaxies plays a fundamental role in the evaluation of the dust attenuation of star forming galaxies especially at high redshifts. Many authors, however, pointed out that there is a significant dispersion and/or deviation from the originally proposed IRX-\beta relation depending on sample selection. We reexamined the IRX-\beta relation by measuring the far- and near-UV flux densities of the original sample galaxies with GALEX and AKARI imaging data, and constructed a revised formula. We found that the newly obtained IRX values were lower than the original relation because of the significant underestimation of the UV flux densities of the galaxies, caused by the small aperture of IUE, Further, since the original relation was based on IRAS data which covered a wavelength range of \lambda = 42--122\mum, using the data from AKARI which has wider wavelength coverage toward longer wavelengths, we obtained an appropriate IRX-\beta relation with total dust emission (TIR): \log(L_{\rm TIR}/L_{\rm FUV}) = \log [10^{0.4(3.06+1.58\beta)}-1] +0.22. This new relation is consistent with most of the preceding results for samples selected at optical and UV, though there is a significant scatter around it. We also found that even the quiescent class of IR galaxies follows this new relation, though luminous and ultraluminous IR galaxies distribute completely differently as well known before.

This paper investigate what is the main driver of the dust mass growth in the interstellar medium (ISM) by using a chemical evolution model of galaxy with metals (elements heavier than helium) in dust phase in addition to the total amount of metals. We consider asymptotic giant branch (AGB) stars, type II supernovae (SNe II) and the dust mass growth in the ISM as the sources of dust, and SN shocks as the destruction mechanism of dust. Further, to describe the dust evolution precisely, our model takes into account the age and metallicity (the ratio of metal mass to ISM mass) dependence of the sources of dust. We particularly focused on the dust mass growth, and found that the dust mass growth in the ISM is regulated by the metallicity. To quantify this aspect, we introduce a "critical metallicity", which is a metallicity at which the contribution of stars (AGB stars and SNe II) equals that of the dust mass growth in the ISM. If the star formation timescale is shorter, the value of the critical metallicity is higher, but the galactic age at which the metallicity reaches the critical metallicity is shorter. From observations, it was expected that the dust mass growth was the dominant source of dust in the Milky Way and dusty QSOs at high redshifts. By introducing the critical metallicity, it is clearly shown that the dust mass growth is the main source of dust in such galaxies with various star formation timescales and ages. The dust mass growth in the ISM is regulated by metallicity, and we stress that the critical metallicity works as an indicator to judge whether the grain growth in the ISM is dominant source of dust in a galaxy, especially because of the strong and nonlinear dependence on the metallicity.

The origin of cosmic dust is a fundamental issue in planetary science. This paper revisits the origin of dust in galaxies, in particular, in the Milky Way, by using a chemical evolution model of a galaxy composed of stars, interstellar medium, metals (elements heavier than helium), and dust. We start from a review of time-evolutionary equations of the four components, and then, we present simple recipes for the stellar remnant mass and yields of metal and dust based on models of stellar nucleosynthesis and dust formation. After calibrating some model parameters with the data from the solar neighborhood, we have confirmed a shortage of the stellar dust production rate relative to the dust destruction rate by supernovae if the destruction efficiency suggested by theoretical works is correct. If the dust mass growth by material accretion in molecular clouds is active, the observed dust amount in the solar neighborhood is reproduced. We present a clear analytic explanation of the mechanism for determining dust content in galaxies after the activation of accretion growth: a balance between accretion growth and supernova destruction. Thus, the dust content is independent of the uncertainty of the stellar dust yield after the growth activation. The timing of the activation is determined by a critical metal mass fraction which depends on the growth and destruction efficiencies. The solar system formation seems to have occurred well after the activation and plenty of dust would have existed in the proto-solar nebula.

We investigate the effect of photodesorption on the snow line position at the surface of a protoplanetary disk around a Herbig Ae/Be star, motivated by the detection of water ice particles at the surface of the disk around HD142527 by Honda et al. For this aim, we obtain the density and temperature structure in the disk with a 1+1D radiative transfer and determine the distribution of water ice particles in the disk by the balance between condensation, sublimation, and photodesorption. We find that photodesorption induced by the far-ultraviolet radiation from the central star depresses the ice-condensation front toward the mid-plane and pushes the surface snow line outward significantly when the stellar effective temperature exceeds a certain critical value. This critical effective temperature depends on the stellar luminosity and mass, the water abundance in the disk, and the yield of photodesorption. We present an approximate analytic formula for the critical temperature. We separate Herbig Ae/Be stars into two groups on the HR diagram according to the critical temperature; one is the disks where photodesorption is effective and from which we may not find ice particles at the surface, and the other is the disks where photodesorption is not effective. We estimate the snow line position at the surface of the disk around HD142527 to be 100--300 AU which is consistent with the water ice detection at >140 AU in the disk. All results depend on the dust grain size by a complex way and this point requires more work in future.

Population III galaxies, made partly or exclusively of metal-free stars, are predicted to exist at high redshifts and may produce very strong Lya emission. A substantial fraction of these Lya photons are likely absorbed in the intergalactic medium at z>6, but recent simulations suggest that significant Lya emission may be detectable up to z~8.5, i.e. well into the reionization epoch. Here, we argue that high-redshift population III galaxies with strong Lya emission can be identified in Hubble Space Telescope imaging data because of their unusual colours. We quantify this effect in some of the filters used in Y-band dropout searches for galaxies at z~8 and find that population III galaxies with high Lya fluxes may exhibit much bluer J-H colours at z=8-10 than any normal type of galaxy at these redshifts. This colour signature can arise even if pop III stars account for as little as ~1e-3 to ~1e-2 of the stellar mass in these galaxies. Some of the anomalously blue objects reported in current Y-band dropout samples do in fact meet the colour criteria for Lya-emitting population III galaxies.

Finding the first generation of galaxies in the early Universe is the greatest step forward for understanding galaxy formation and evolution. For strategic survey of such galaxies and interpretation of the obtained data, this paper presents an ultraviolet-to-optical spectral model of galaxies with a great care of the nebular emission. In particular, we present a machine-readable table of intensities of 119 nebular emission lines from Ly$\alpha$ to the rest-frame 1 $\mu$m as a function of metallicity from zero to the Solar one. Based on the spectral model, we present criteria of equivalent widths of Ly$\alpha$, He {\sc ii} $\lambda1640$, H$\alpha$, H$\beta$, [O {\sc iii}] $\lambda5007$ to select extremely metal-poor and metal-free galaxies although these criteria have uncertainty caused by the Lyman continuum escape fraction and the star formation duration. We also present criteria of broad-band colours which will be useful to select candidates for spectroscopic follow-up from drop-out galaxies. We propose the line intensity ratio of [O {\sc iii}] $\lambda5007$ to H$\beta$ $<0.1$ as the most robust criterion for $<1/1000$ of the Solar metallicity. This ratio of a galaxy with a few $M_\odot$ yr$^{-1}$ at $z\sim8$ is detectable by spectroscopy with the James Webb Space Telescope within a reasonable exposure time.

We present results of a survey of Lyman $\alpha$ emitters (LAEs) at $z=6.5$ which is thought to be the final epoch of the cosmic reionization. In a $\approx530$ arcmin$^2$ deep image of the SSA22 field taken through a narrowband filter NB912 installed in the Subaru/Suprime-Cam, we have found only 14 LAE candidates with $L_{\rm Ly\alpha}\ga3\times10^{42}$ erg s$^{-1}$. Even applying the same colour selection criteria, the number density of the LAE candidates is a factor of 3 smaller than that found at the same redshift in the Subaru Deep field (SDF). Assuming the number density in the SDF is a cosmic average, the probability to have a number density equal to or smaller than that found in the SSA22 field is only 7% if we consider fluctuation by the large-scale structure (i.e. cosmic variance) and Poisson error. Therefore, the SSA22 field may be a rare void at $z=6.5$. On the other hand, we have found that the number density of $i'$-drop galaxies with $25.5<z'<26.0$ in the SSA22 field agrees well with that in the SDF. If we consider a scenario that a larger neutral fraction of intergalactic hydrogen, $x_{\rm HI}$, in the SSA22 field obscures a part of Ly$\alpha$ emission, $x_{\rm HI}$ in the SSA22 field should be about 2 times larger than that in the SDF. This can be translated into $x_{\rm HI}<0.9$ at $z=6.5$ in the SSA22 field. A much larger survey area than previous ones is required to overcome a large fluctuation reported here and to obtain a robust constraint on $x_{\rm HI}$ at the end of the reionization from LAEs.

(Abridged) Deep NB359 imaging with Subaru by Iwata et al. have detected surprisingly strong Lyman continuum (LyC; ~900A in the rest-frame) from some LAEs at z=3.1. However, the redshifts might be misidentified due to a narrow wavelength coverage in previous spectroscopy. We here present new deep spectroscopy covering the observed 4,000-7,000A with VLT/VIMOS and Subaru/FOCAS of 8 LAEs detected in NB359. All the 8 objects have only one detectable emission line around 4,970A which is most likely to be Ly-A at z=3.1, and thus, the objects are certainly LAEs at the redshift. However, 5 of them show a ~0.''8 spatial offset between the Ly-A emission and the source detected in NB359. No indications of the redshifts of the NB359 sources are found although it is statistically difficult that all the 5 LAEs have a foreground object accounting for the NB359 flux. The rest 3 LAEs show no significant offset from the NB359 position. Therefore, they are truly LyC emitting LAEs at z=3.1. We also examine the stellar population which simultaneously accounts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet of the LAEs. We consider the latest statistics of Lyman limit systems to estimate the LyC optical depth in the IGM and an additional contribution of the bound-free LyC from photo-ionized nebulae to the LyC emissivity. As a result, we find that stellar populations with metallicity Z>=1/50Z_sun can explain the observed LyC strength only with a very top-heavy initial mass function (IMF; <m>~50 M_sun). However, the critical metallicity for such an IMF is expected to be much lower. A very young (~1 Myr) and massive (~100 M_sun) extremely metal-poor (Z<=5e-4Z_sun) or metal-free (so-called Population III) stellar population can reproduce the observed LyC strength. The required mass fraction of such `primordial' stellar population is ~1--10% in total stellar mass of the LAEs.

We use ultra-deep ultraviolet VLT/VIMOS intermediate-band and VLT/FORS1 narrow-band imaging in the GOODS Southern field to derive limits on the distribution of the escape fraction (f_esc) of ionizing radiation for L >~ L*(z=3) Lyman Break Galaxies (LBGs) at redshift 3.4--4.5. Only one LBG, at redshift z=3.795, is detected in its Lyman continuum (LyC; S/N~5.5), the highest redshift galaxy currently known with a direct detection. Its ultraviolet morphology is quite compact (R_eff=0.8, kpc physical). Three out of seven AGN are also detected in their LyC, including one at redshift z=3.951 and z850 = 26.1. From stacked data (LBGs) we set an upper limit to the average f_esc in the range 5%--20%, depending on the how the data are selected (e.g., by magnitude and/or redshift). We undertake extensive Monte Carlo simulations that take into account intergalactic attenuation, stellar population synthesis models, dust extinction and photometric noise in order to explore the moments of the distribution of the escaping radiation. Various distributions (exponential, log-normal and Gaussian) are explored. We find that the median f_esc is lower than ~6% with an 84% percentile limit not larger than 20%. If this result remains valid for fainter LBGs down to current observational limits, then the LBG population might be not sufficient to account for the entire photoionization budget at the redshifts considered here, with the exact details dependent upon the assumed ionizing background and QSO contribution thereto. It is possible that f_esc depends on the UV luminosity of the galaxies, with fainter galaxies having higher f_esc, and estimates of f_esc from a sample of faint LBG from the HUDF (i775<28.5) are in broad quantitative agreement with such a scenario.

The disk around AB Aur was imaged and resolved at 24.6\,$\mu$m using the Cooled Mid-Infrared Camera and Spectrometer on the 8.2m Subaru Telescope. The gaussian full-width at half-maximum of the source size is estimated to be 90 $\pm$ 6 AU, indicating that the disk extends further out at 24.6\,$\mu$m than at shorter wavelengths. In order to interpret the extended 24.6\,$\mu$m image, we consider a disk with a reduced surface density within a boundary radius $R_c$, which is motivated by radio observations that suggest a reduced inner region within about 100 AU from the star. Introducing the surface density reduction factor $f_c$ for the inner disk, we determine that the best match with the observed radial intensity profile at 24.6\,$\mu$m is achieved with $R_c$=88 AU and $f_c$=0.01. We suggest that the extended emission at 24.6\,$\mu$m is due to the enhanced emission from a wall-like structure at the boundary radius (the inner edge of the outer disk), which is caused by a jump in the surface density at $R_c$. Such reduced inner disk and geometrically thick outer disk structure can also explain the more point-like nature at shorter wavelengths. We also note that this disk geometry is qualitatively similar to a pre-transitional disk, suggesting that the AB Aur disk is in a pre-transitional disk phase.

We present the discovery of a candidate of giant radio-quiet Lyman-alpha (Lya) blob (RQLAB) in a large-scale structure around a high-redshift radio galaxy (HzRG) lying in a giant Lya halo, B3 J2330+3927 at redshift z=3.087. We obtained Lya imaging around B3 J2330+3927 with Subaru/Suprime-Cam to search for Lya emitters (LAEs) and absorbers (LAAs) at redshift z=3.09+-0.03. We detected candidate 127 LAEs and 26 LAAs in the field of view of 31' x 24'. We found that B3 J2330+3927 is surrounded by a 130 kpc Lya halo and a large-scale (60 x 20 comoving Mpc) filamentary structure. The large-scale structure contains one prominent local density peak with an overdensity of greater than 5, which is 8' (15 comoving Mpc) away from B3 J2330+3927. In this peak, we discovered a candidate 100 kpc RQLAB. The existence of both types of Lya nebulae in the same large-scale structure suggests that giant Lya nebulae need special large-scale environments to form. On smaller scales, however, the location of B3 J2330+3927 is not a significant local density peak in this structure, in contrast to the RQLAB. There are two possible interpretations of the difference of the local environments of these two Lya nebulae. Firstly, RQLAB may need a prominent (delta ~ 5) density peak of galaxies to form through intense star-bursts due to frequent galaxy interactions/mergers and/or continuous gas accretion in an overdense environment. On the other hand, Lya halo around HzRG may not always need a prominent density peak to form if the surrounding Lya halo is mainly powered by its radio and AGN activities. Alternatively, both RQLAB and Lya halo around HzRG may need prominent density peaks to form but we could not completely trace the density of galaxies because we missed evolved and dusty galaxies in this survey.

It is essential to know galactic emissivity and spectrum of Lyman continuum (LyC) in order to understand the cosmic reionization. Here we consider an escape of nebular LyC from galaxies and examine the consequent spectral energy distribution. It is usually assumed that hydrogen nebular LyC mostly produced by bound-free transitions is consumed within photo-ionized nebulae (so-called "on-the-spot" approximation). However, an escape of the continuum should be taken into account if stellar LyC escapes from galaxies through "matter-bounded" nebulae. We show that the escaping hydrogen bound-free LyC makes a strong bump just below the Lyman limit. Such a galaxy would be observed as a Lyman "bump" galaxy. This bump results from the radiation energy re-distribution of stellar LyC by nebulae. The strength of the bump depends on electron temperature in nebulae, escape fraction of stellar and nebular LyC, hardness of stellar LyC (i.e. metallicity, initial mass function, age and star formation history), and IGM attenuation. We can use the bump to find very young ($\sim1$ Myr), massive ($\sim100$ $M_\odot$), and extremely metal-poor (or metal-free) stellar populations at $z<4$. Because of the bump, 900 \AA-to-1500 \AA luminosity density ratio (per Hz) becomes maximum (2--3 times larger than the stellar intrinsic ratio) when about 40% of the stellar LyC is absorbed by nebulae. The total number of escaping LyC photons increases due to the escape of nebular LyC but does not exceed the stellar intrinsic one. The radiation energy re-distribution by nebulae reduces the mean energy of escaping LyC only by $\approx10$% relative to that of stellar LyC. Therefore, the effect of the escape of nebular LyC on the reionization process may be small.

Knowing the amount of ionizing photons from young star-forming galaxies is of particular importance to understanding the reionization process. Here we report initial results of Subaru/Suprime-Cam deep imaging observation of the SSA22 proto-cluster region at z=3.09, using a special narrow-band filter to optimally trace Lyman continuum (LyC) from galaxies at z~3. The unique wide field-of-view of Suprime-Cam enabled us to search for ionizing photons from 198 galaxies (73 Lyman break galaxies (LBGs) and 125 Ly-alpha emitters (LAEs)) with spectroscopically measured redshifts z~3.1. We detected LyC from 7 LBGs, as well as from 10 LAE candidates. Some of the detected galaxies show significant spatial offsets of LyC from non-ionizing UV emission. For some LBGs the observed non-ionizing UV to LyC flux density ratios are smaller than values expected from population synthesis models with a standard Salpeter initial mass function (IMF) with moderate dust attenuation (which is suggested from the observed UV slopes), even if we assume very transparent IGM along the sightlines of these objects. This implies an intrinsically bluer spectral energy distribution, e.g, that produced by a top-heavy IMF, for these LBGs. The observed flux desity ratios of non-ionizing UV to LyC of 7 detected LBGs range from 2.4 to 23.8 and the median is 6.6. The observed flux density ratios of the detected LAEs are even smaller than LBGs, if they are truly at z~3.1. We find that the median value of the flux density ratio for the deteced LBGs suggest that their escape fractions is likely to be higher than 4%, if the LyC escape is isotropic. The results imply that some of the LBGs in the proto-cluster at z~3 have the escape fraction significantly higher than that of galaxies (in a general field) at z~1 studied previously.

The temperature in the optically thick interior of protoplanetary discs is essential for the interpretation of millimeter observations of the discs, for the vertical structure of the discs, for models of the disc evolution and the planet formation, and for the chemistry in the discs. Since large icy grains have a large albedo even in the infrared, the effect of scattering of the diffuse radiation in the discs on the interior temperature should be examined. We have performed a series of numerical radiation transfer simulations including isotropic scattering by grains with various typical sizes for the diffuse radiation as well as for the incident stellar radiation. We also have developed an analytic model including isotropic scattering to understand the physics concealed in the numerical results. With the analytic model, we have shown that the standard two-layer approach is valid only for grey opacity (i.e. grain size $\ga10$ \micron) even without scattering. A three-layer interpretation is required for grain size $\la10$ \micron. When the grain size is 0.1--10 \micron, the numerical simulations show that isotropic scattering reduces the temperature of the disc interior. This reduction is nicely explained by the analytic three-layer model as a result of the energy loss by scatterings of the incident stellar radiation and of the warm diffuse radiation in the disc atmosphere. For grain size $\ga10$ \micron (i.e. grey scattering), the numerical simulations show that isotropic scattering does not affect the interior temperature. This is nicely explained by the analytic two-layer model; the energy loss by scattering in the disc atmosphere is exactly offset by the "green-house effect" due to scattering of the cold diffuse radiation in the interior.

We have examined the dust photoelectric heating in the intergalactic medium (IGM). The heating rate in a typical radiation field of the IGM is represented by $\Gamma_{\rm pe} = 1.2\times10^{-34}$ erg s$^{-1}$ cm$^{-3}$ $({\cal D}/10^{-4})(n_{\rm H}/10^{-5} {\rm cm^{-3}})^{4/3} (J_{\rm L}/10^{-21} {\rm erg s^{-1} cm^{-2} Hz^{-1} sr^{-1}})^{2/3} (T/10^4 {\rm K})^{-1/6}$, where ${\cal D}$ is the dust-to-gas mass ratio, $n_{\rm H}$ is the hydrogen number density, $J_{\rm L}$ is the mean intensity at the hydrogen Lyman limit of the background radiation, and $T$ is the gas temperature, if we assume the new X-ray photoelectric yield model by Weingartner et al. (2006) and the dust size distribution in the Milky Way by Mathis, Rumpl, & Nordsieck (1977). This heating rate dominates the HI and HeII photoionization heating rates when the hydrogen number density is less than $\sim10^{-6}$ cm$^{-3}$ if ${\cal D}=10^{-4}$ which is 1% of that in the Milky Way, although the heating rate is a factor of 2--4 smaller than that with the old yield model by Weingartner & Draine (2001). The grain size distribution is very important. If only large ($\ge0.1$ $\mu$m) grains exist in the IGM, the heating rate is reduced by a factor of $\simeq5$. Since the dust heating is more efficient in a lower density medium relative to the photoionization heating, it may cause an inverted temperature--density relation in the low density IGM suggested by Bolton et al. (2008). Finally, we have found that the dust heating is not very important in the mean IGM before the cosmic reionization.

We have made a Monte Carlo simulation of the intergalactic absorption in order to model the Lyman continuum absorption, which is required to estimate the escape fraction of the Lyman continuum from distant galaxies. To input into the simulation, we derive an empirical distribution function of the intergalactic absorbers which reproduces recent observational statistics of the Lyman $\alpha$ forest, Lyman limit systems (LLSs), and damped Lyman $\alpha$ systems (DLAs) simultaneously. In particular, we assume a common functional form of the number evolution along the redshift for all types of absorbers. The Lyman series transmissions in our simulation reproduce the observed redshift evolution of the transmissions excellently, and the Lyman continuum transmission also agrees with an observed estimation which is still quite rare in the literature. The probability distribution of the Lyman $\alpha$ opacity in our simulation is log-normal with a tail towards a large opacity. This tail is produced by DLAs. The probability distribution of the Lyman continuum opacity in our simulation also show a broad tail towards a large opacity. This tail is produced by LLSs. Because of the rarity of LLSs, we have a chance to have a clean line of sight in the Lyman continuum even for $z\sim4$ with a probability of about 20%. Our simulation expects a good correlation between the Lyman continuum opacity and the Lyman $\alpha$ opacity, which may be useful to estimate the former from the latter for an individual line of sight.

We discuss how we obtain the spatial distribution of ice on the surface of the circumstellar disk around young stars. Ice in the disks plays a very important role in various issues, for instance, on the disk structure, on the planet formation, on the isotopic anomaly in meteorites, and on the origin of the sea on the Earth. Therefore, the spatially resolved observation of the condensation/sublimation front of ice, so-called ``snow line'' is strongly required. Here, we propose a new method for obtaining the spatially resolved ``snow line'' on the circumstellar disks by observing 3 \micron H$_2$O ice feature in the scattered light. Based on radiative transfer considerations, we show that the feature is clearly imprinted in the spectrum of the scattered light from both optically thick and thin circumstellar disks. We also show that the scattered light and the H$_2$O ice feature from protoplanetary disks are detectable and spatially resolvable with the current instruments through a $H_2O$ narrowband filter around 3 \micron. Finally, we present a diagnostics of disk dust properties on the $K-H_2O$ and $K-L'$ two color diagram.

We analyse the attenuation properties of a sample of UV selected galaxies, with the use of the spectrophotometric model Grasil. In particular, we focus on the relation between dust attenuation and the reddening in the UV spectral region. We show that a realistic modelling of geometrical distribution of dust and of the different population of stars can explain the UV reddening of normal spiral galaxies also with a standard Milky Way dust. Our results clearly underline that it is fundamental to take into account that younger stars suffer a higher attenuation than older stars (the age-dependent extinction) because stars are born in more-than-average dusty environments. In this work we also find that the concentration of young stars on the galactic plane of spirals has a relevant impact on the expected UV colours, impact that has not been explored before this paper. Finally, we discuss the role of IMF in shaping the relation between UV reddening and dust attenuation, and we show that a Kroupa IMF is more consistent with observed data than the classical Salpeter IMF.

The escape fraction of ionizing photons from galaxies is a crucial quantity controlling the cosmic ionizing background radiation and the reionization. Various estimates of this parameter can be obtained in the redshift range, z=0--6, either from direct observations or from the observed ionizing background intensities. We compare them homogeneously in terms of the observed flux density ratio of ionizing ($\sim900$ \AA rest-frame) to non-ionizing ultraviolet ($\sim1500$ \AA rest-frame) corrected for the intergalactic absorption. The escape fraction is found to increase by an order of magnitude, from a value less than 0.01 at $z\la1$ to about 0.1 at $z\ga4$.

We discuss dust properties in the interstellar medium (ISM) of nearby normal galaxies, by comparing observations in the ultraviolet (UV) with simulations by a radiative transfer model. The observed UV colours of nearby galaxies show a reddening relative to their expected intrinsic colours. Some authors argued that the Milky Way dust cannot reproduce the reddening because of the prominent 2175 \AA absorption bump. Other authors proposed a reduction mechanism of the bump strength in an {\it attenuation law} derived from the ratio of the observed intensity to the intrinsic one through an age-selective attenuation (i.e., young stars are more attenuated selectively). We newly find that the wavelength dependence of the scattering albedo also has a strong effect on the UV colour; an albedo decreasing toward shorter wavelengths (except for the absorption bump range) produces a significant UV reddening. After comparing the observed UV colours of nearby normal galaxies with those expected from radiative transfer simulations assumed several dust models, we find two sorts of dust suitable for these galaxies: (1) dust with a bump and a smaller albedo for a shorter wavelength (except for the bump range), and (2) dust without any bump but with an almost constant albedo. If very small carbonaceous grains responsible for the common unidentified infrared emission band are also the bump carrier, the former dust is favorable. Finally, we derive mean attenuation laws of various dust models as a function of the UV attenuation, and derive some relations between the UV attenuation and observable/theoretical quantities.

Observed ultraviolet (UV) colours of nearby disc galaxies show a reddening relative to their expected intrinsic colours. Since the 2175 \AA bump found in the Milky Way's dust extinction law blues the UV colours, it might suggest that dust in nearby disc galaxies does not have the bump and that the Milky Way is exceptional. However, this conclusion can be modified by the effect of scatterings. If the scattering albedo decreases towards shorter wavelengths, observable UV colours redden. An extensive comparison between observed UV colours and those expected from radiative transfer simulations shows two types of dust suitable for nearby disc galaxies: (1) dust with a bump and a smaller albedo for a shorter wavelength (except for the bump range), and (2) dust without any bump but with an almost constant albedo. If very small carbonaceous grains responsible for the common unidentified infrared emission band are also the bump carrier, the former dust is favorable.

The thermal and fragmentation properties of star-forming clouds have important consequences on the corresponding characteristic stellar mass. The initial composition of the gas within these clouds is a record of the nucleosynthetic products of previous stellar generations. In this paper we present a model for the evolution of star-forming clouds enriched by metals and dust from the first supernovae, resulting from the explosions of metal-free progenitors with masses in the range 12 - 30 Msun and 140 - 260 Msun. Using a self-consistent approach, we show that: (i) metals depleted onto dust grains play a fundamental role, enabling fragmentation to solar or sub-solar mass scales already at metallicities Zcr = 10^{-6} Zsun; (ii) even at metallicities as high as 10^{-2} Zsun, metals diffused in the gas-phase lead to fragment mass scales which are ~ 100 Msun; (iii) C atoms are strongly depleted onto amorphous carbon grains and CO molecules so that CII plays a minor role in gas cooling, leaving OI as the main gas-phase cooling agent in low-metallicity clouds. These conclusions hold independently of the assumed supernova progenitors and suggest that the onset of low-mass star formation is conditioned to the presence of dust in the parent clouds.

We present results of VLT/FORS2 spectroscopy of galaxies at z~3 in the Hubble Deep Field-South (HDF-S). A sample of galaxies was drawn from the photo-z catalogue based on the HST/WFPC2 optical images and the deep near-infrared images obtained with VLT/ISAAC as a part of the FIRES project. We selected galaxies with photometric redshift between 2.5 and 4. Most of the selected galaxies are bright in rest-frame UV wavelengths and satisfy color selection criteria of Lyman break galaxies (LBGs) at z~3. The number of target galaxies with I(AB)<25.0 was 15. We identified new 5 firm and 2 probable redshifts in addition to confirmations of previously known 6 galaxies at z~3. We found 6 among these 13 galaxies lie at a quite narrow redshift range at z = 2.80+-0.01. We examined stellar populations of the galaxies with spectroscopic redshifts through comparisons of their optical and near-IR photometry data with template spectra generated by a population synthesis code. The ages from the onset of star formation for these star-forming galaxies with I<=25.0 are typically 50-200 Myr, and their stellar masses are between (0.5-5) x 10^10 M_sun, consistent with previous studies. We also compared these SED fitting results with those for distant red galaxies (DRGs) at z>2 discovered by FIRES. DRGs have larger stellar masses, larger dust attenuation than our UV-luminous LBG sample, and their star formation rates are often comparable to LBGs. These trends suggest that majority of DRGs are indeed the most massive systems at the redshift and are still in the active star-forming phase. Unless the number density of DRGs is much smaller than LBGs, estimates based on UV selected sample could miss substantial part of stellar mass density at z~3.

We have performed narrow-band imaging observations with the Very Large Telescope, aimed at detecting the Lyman continuum (LC) flux escaping from galaxies at z~3. We do not find any significant LC flux from our sample of two galaxies in the Hubble Deep Field South, at z=3.170 and 3.275. The corresponding lower limits on the F_{1400}/F_{900} flux density (per Hz) ratio are 15.6 and 10.2 (3-sigma confidence level). After correction for the intergalactic hydrogen absorption, the resulting limits on the relative escape fraction of the LC are compared with those obtained by different approaches, at similar or lower redshifts. One of our two objects has a relative escape fraction lower than the detection reported by Steidel et al. in a composite spectrum of z~3 galaxies. A larger number of objects is required to reach a significant conclusion. Our comparison shows the potential of narrow-band imaging for obtaining the best limit on the relative escape fraction at z~3. Stacking a significant number of galaxies observed through a narrow-band filter would provide constraint on the galactic contribution to the cosmic reionization.

We investigate the attenuation law seen through an interstellar medium (ISM) with clumpy spatial distributions of stars and dust. The clumpiness of the dust distribution is introduced by a multi-phase ISM model. We solve a set of radiative transfer equations with multiple anisotropic scatterings through the clumpy ISM in a 1-D plane-parallel geometry by using the mega-grain approximation, in which dusty clumps are regarded as very large particles (i.e. mega-grains). The clumpiness of the stellar distribution is introduced by the youngest stars embedded in the clumps. We assume a smooth spatial distribution for older stars. The youngest stars are surrounded by denser dusty gas and suffer stronger attenuation than diffuse older stars (i.e. age-selective attenuation). The apparent attenuation law is a composite of the attenuation laws for the clumpy younger stars and for the diffuse older stars with a luminosity weight. In general, the stellar population dominating the luminosity changes from older stars to younger stars as the wavelength decreases. This makes the attenuation law steep; the composite attenuation rapidly increases from small attenuation for older stars at a long wavelength to large attenuation for younger stars at a short wavelength. The resultant attenuation law of normal disc galaxies is expected to be much steeper than that of starburst galaxies observed by Calzetti et al. Finally, the Calzetti's attenuation law is regarded as a special case with a large density in our framework.

This paper examines the allowed amount of IG (intergalactic) dust, which is constrained by extinction and reddening of distant SNe and thermal history of IGM (intergalactic medium) affected by dust photoelectric heating. Based on the observational cosmic star formation history, we find an upper bound of $\chi$, the mass ratio of the IG dust to the total metal in the Universe, as $\chi\la 0.1$ for $10 {\rm \AA} \la a \la 0.1 \micron$ and $\chi\la 0.1(a/0.1 \micron)^{-1}$ for $0.1 \micron\la a\la1\micron$, where $a$ is a characteristic grain size of the IG dust. This upper bound of $\chi\sim0.1$ suggests that the dust-to-metal ratio in the IGM is smaller than the current Galactic value. The corresponding allowed density of the IG dust increases from $\sim10^{-34}$ g cm$^{-3}$ at $z=0$ to $\sim10^{-33}$ g cm$^{-3}$ at $z\sim1$, and keeps almost the value toward higher redshift. This causes IG extinction of $\la 0.2$ mag at the observer's $B$-band for $z\sim 1$ sources and that of $\la 1$ mag for higher redshift sources. Furthermore, if $E(B-V)\sim 0.1$ mag at the observer's frame against $z\ga1$ sources is detected, we can conclude that a typical size of the IG dust is $\la 100$ \AA. The 2175 \AA absorption feature of small graphite may be found as a local minimum at $z\sim2.5$ in a plot of the observed $E(B-V)$ as a function of the source redshift. Finally, the IGM mean temperature at $z\la1$ can be still higher than $10^{4}$ K, provided the size of the IG dust is $\la100$ \AA.

Using the near-infrared (NIR) observations of the afterglows of the high redshift ($5 \la z \la 25$) gamma-ray bursts (GRBs) which will be detected by the {\it Swift} satellite, we discuss a way to study the cosmic reionization history. In principle the details of the cosmic reionization history are well imprinted in the NIR spectra of the GRB afterglows. However the spectroscopy with a space telescope is required to obtain such an information for a very high redshift ($z\ga 15$) unless the neutral fraction of the high-$z$ universe is less than $10^{-6}$. The broad-band photometry has the higher sensitivity than that of the spectroscopy, so that the NIR photometric follow-up of the GRB afterglows is very promising to examine the cosmic reionization history. A few minutes exposure with a 8-m class ground-based telescope of the afterglow of the high-$z$ GRBs will reveal how many times the reionization occurred in the universe.

We investigate the GRB formation rate based on the relation between the spectral peak energy ($E_{p}$) and the isotropic luminosity. The $E_{p}$--luminosity relation covers the energy range of 50 -- 2000 keV and the luminosity range of $10^{50}$--$10^{54} erg/s, respectively. We find that the relation is considerably tighter than similar relations suggested previously. Using $E_{p}$--luminosity relation, we estimate the luminosity and the redshift of 684 GRBs with the unknown distances and derive the GRB formation rate as a function of the redshift. For $0 \le z \le 2$, the GRB formation rate is well correlated with the star formation rate while it increases monotonously from $z\sim 2$ out to $z \sim 12$. This behavior is consistent with the results of previous works using the lag--luminosity relation or the variability--luminosity relation.

Infrared (IR) luminosity of galaxies originating from dust emission can be used as an indicator of the star formation rate (SFR). Inoue et al. (2000, IHK) have derived a formula for the conversion from IR luminosity to SFR by using the following three quantities: the fraction of Lyman continuum luminosity absorbed by gas (f), the fraction of UV luminosity absorbed by dust (epsilon), and the fraction of dust heating from old (>10^8 yr) stellar populations (eta). We develop a method to estimate those three quantities based on the idea that the various way of SFR estimates should return the same SFR. After applying our method to samples of galaxies, the following results are obtained. First, our method is applied to star-forming galaxies, finding that f~0.6, epsilon~0.5, and eta~0.4 as representative values. Next, we apply the method to a starburst sample, which shows larger extinction than the star-forming galaxy sample. With the aid of f, epsilon, and eta, we estimate reliable SFRs. Moreover, the H-alpha luminosity, if the H-alpha extinction is corrected by using the Balmer decrement, is suitable for a statistical analysis of SFR, because the same correction factor for the Lyman continuum extinction is applicable to both normal and starburst galaxies over all the range of SFR. The metallicity dependence of f and epsilon is also tested: Only the latter proves to have a correlation with metallicity. As an extension of our result, we show that all UV, H-alpha, and IR comoving luminosity densities at z=0 give a consistent SFR (~ 3x10^{-2}h M_sun/Mpc^3). Useful formulae for SFR estimate are listed.

This paper investigates the evolution of the dust-to-metal ratio in galaxies based on a simple evolution model for the amount of metal and dust with infall. We take into account grain formation in stellar mass-loss gas, grain growth by the accretion of metallic atoms in a cold dense cloud, and grain destruction by SNe shocks. Especially, we propose that the accretion efficiency is independent of the star-formation history. This predicts various evolutionary tracks in the metallicity ($Z$)--dust-to-gas ratio ($\cal D$) plane depending on the star-formation history. In this framework, the observed linear $Z$--$\cal D$ relation of nearby spiral galaxies can be interpreted as a sequence of a constant galactic age. We emphasize that an observational study of the $Z$--$\cal D$ relation of galaxies at $z\sim 1$ is very useful to constrain the efficiencies of dust growth and destruction. We also suggest that the Lyman break galaxies at $z\sim 3$ have a very low dust-to-metal ratio, typically $\ltsim 0.1$. Although the effect of infall on the evolutionary tracks in the $Z$--$\cal D$ plane is quite small, the dispersion of the infall rate can disturb the $Z$--$\cal D$ relation with a constant galactic age.

This Letter investigates the amount of dust in the intergalactic medium (IGM). The dust photoelectric heating can be the most efficient heating mechanism in the IGM where the density is very small and there are a lot of hard ultraviolet photons. Comparing the observational thermal history of IGM with a theoretical one taking into account the dust photoelectric heating, we can put an upper limit on the dust-to-gas ratio, ${\cal D}$, in the IGM. Since the rate of the dust photoelectric heating depends on the size of dust, we find the following results: If the grain size is $\ga 100$ \AA, ${\cal D}$ at $z \sim 3$ is $\la 1/100$ Galactic value corresponding to $\Omega_{\rm dust}^{\rm IGM}\la 10^{-5}$. On the other hand, if the grain size is as small as $\sim 10$ \AA, ${\cal D}$ is $\la 1/1000$ Galactic value corresponding to $\Omega_{\rm dust}^{\rm IGM}\la 10^{-6}$.

Before the end of 2002 will be launched the GALEX satellite (a NASA/SMEX project) which will observe all the sky in Ultraviolet (UV) through filters at 1500 and 2300 A down to m(AB) 21. In 2004 will be launched the ASTRO-F satellite which will perform an all sky survey at Far-Infrared (FIR) wavelengths. The cross-correlation of both surveys will lead to very large samples of galaxies for which FIR and UV fluxes will be available. Using the FIR to UV flux ratio as a quantitative tracer of the dust extinction we will be able to measure the extinction in the nearby universe (z<0.2) and to perform a statistically significant analysis of the extinction as a function of galactic properties. Of particular interest is the construction of pure FIR and UV selected samples for which the extinction will be measured as templates for the observation of high redshift galaxies.

We examine what types of galaxies the conversion formula from dust infrared (IR) luminosity into the star formation rate (SFR) derived by Kennicutt (1998) is applicable to. The ratio of the observed IR luminosity, $L_{\rm IR}$, to the intrinsic bolometric luminosity of the newly ($\la$ 10 Myr) formed stars, $L_{\rm SF}$, of a galaxy can be determined by a mean dust opacity in the interstellar medium and the activity of the current star formation. We find that these parameters area being $0.5 \le L_{\rm IR}/L_{\rm SF} \le 2.0$ is very large, and many nearby normal and active star-forming galaxies really fall in this area. It results from offsetting two effects of a small dust opacity and a large cirrus contribution of normal galaxies relative to starburst galaxies on the conversion of the stellar emission into the dust IR emission. In conclusion, the SFR determined from the IR luminosity under the assumption of $L_{\rm IR}=L_{\rm SF}$ like Kennicutt (1998) is reliable within a factor of 2 for all galaxies except for dust rich but quiescent galaxies and extremely dust poor galaxies.

We discuss the dust distribution within photoionized regions. Assuming a geometry with a central dust cavity, which is strongly suggested by the literature, we can estimate the cavity's radius from the ratio of the infrared and radio fluxes by using a simple transfer model of Lyman continuum photons. We apply the method to a sample of the Galactic H {\sc ii} regions. The estimated typical radius of the dust cavity of the Galactic compact H {\sc ii} regions is about 30% of the Str{\"o}mgren radius. Taking account of uncertainties both of the observational data and the model, we can reject a dust distribution lacking a central cavity. Therefore, the dust cavity model is supported independently of the previous works. We discuss the formation mechanism of such a dust cavity and its detectability by present and future infrared facilities.

We examine Lyman continuum extinction (LCE) in H II regions by comparing infrared fluxes of 49 H II regions in the Galaxy, M31, M33, and the LMC with estimated production rates of Lyman continuum photons. A typical fraction of Lyman continuum photons that contribute to hydrogen ionization in the H II regions of three spiral galaxies is $\la$ 50 %. The fraction may become smaller as the metallicity (or dust-to-gas ratio) increases. We examine the LCE effect on estimated star formation rates (SFR) of galaxies. The correction factor for the Galactic dust-to-gas ratio is 2-5.

We re-examine the effect of Lyman continuum ($\lambda \leq 912$ \AA) extinction (LCE) by dust in H {\sc ii} regions in detail and discuss how it affects the estimation of the global star formation rate (SFR) of galaxies. To clarify the first issue, we establish two independent methods for estimating a parameter of LCE ($f$), which is defined as the fraction of Lyman continuum photons contributing to hydrogen ionization in an H {\sc ii} region. One of those methods determines $f$ from the set of Lyman continuum flux, electron density and metallicity. In the framework of this method, as the metallicity and/or the Lyman photon flux increase, $f$ is found to decrease. The other method determines $f$ from the ratio of infrared flux to Lyman continuum flux. Importantly, we show that $f \la 0.5$ via both methods in many H {\sc ii} regions of the Galaxy. Thus, it establishes that dust in such H {\sc ii} regions absorbs significant amount of Lyman continuum photons directly. To examine the second issue, we approximate $f$ to a function of only the dust-to-gas mass ratio (i.e., metallicity), assuming a parameter fit for the Galactic H {\sc ii} regions. We find that a characteristic $\hat{f}$, which is defined as $f$ averaged over a galaxy-wide scale, is 0.3 for the nearby spiral galaxies. This relatively small $\hat{f}$ indicates that a typical increment factor due to LCE for estimating the global SFR ($1/\hat{f}$) is large ($\sim 3$) for the nearby spiral galaxies. Therefore, we conclude that the effect of LCE is not negligible relative to other uncertainties of estimating the SFR of galaxies.

Infrared (IR) dust emission from galaxies is frequently used as an indicator of star formation rate (SFR). However, the effect of the dust-to-gas ratio (i.e., amount of the dust) on the conversion law from IR luminosity to SFR has not so far been considered. Then, in this paper, we present a convenient analytical formula including this effect. In order to obtain the dependence on the dust-to-gas ratio, we extend the formula derived in our previous paper, in which a theoretical formula converting IR luminosity to SFR was derived. That formula was expressed as ${\rm SFR}/(M_\odot~{\rm yr}^{-1})=\{3.3\times 10^{-10}(1- \eta)/(0.4-0.2f+0.6\epsilon)\} (L_{\rm IR}/L_\odot)$, where f is the fraction of ionizing photons absorbed by hydrogen, $\epsilon$ is the efficiency of dust absorption for nonionizing photons, $\eta$ is the cirrus fraction of observed dust luminosity, and $L_{\rm IR}$ is the observed luminosity of dust emission in the 8-1000-$\mu$m range. Our formula explains the IR excess of the Galaxy and the Large Magellanic Cloud. In the current paper, especially, we present the metallicity dependence of our conversion law between SFR and $L_{\rm IR}$. This is possible since both f and $\epsilon$ can be estimated via the dust-to-gas ratio, which is related to metallicity. We have confirmed that the relation between the metallicity and the dust-to-gas ratio is applied to both giant and dwarf galaxies. Finally, we apply the result to the cosmic star formation history. We find that the comoving SFR at z=3 calculated from previous empirical formulae is underestimated by a factor of 4-5.

We investigate the structure of the molecular cloud complexes (MCCs) as a group of several giant molecular clouds (GMCs) in the Galaxy. Then, we find that the mass-size relation which has been reported for the GMCs establishes well even for the very large MCCs whose size is about 1 kpc. Since the horizontal size of the MCCs is more than the thickness of the Galactic disk, we can no longer consider the MCCs to be a sphere. Thus, we construct a structural model of the MCCs, adopting a rectangular-solid geometry. As a result, our model explains the observed mass-size relation of the MCCs very well. {}From the estimated external pressure around the MCCs, we find that they are in a rough pressure balance with the interstellar medium. Moreover, we find there is observational deficiency of the MCCs with a large size and surface density. Then, we suggest that the external pressure has a significant effect on the structure and evolution of the MCCs. We also discuss the effect of H {\sc ii} regions in the MCCs.

It has been reported recently that there are some early-type spiral (Sa--Sab) galaxies having evident star-forming regions which concentrate in their own central 1-kpc. In such central region, is the mechanism of the star formation distinct from that in disks of spiral galaxies? To reveal this, we estimate the star formation efficiency (SFE) in this central 1-kpc star-forming region of some early-type spiral galaxies, taking account of the condition for this 1-kpc region to be self-gravitating. Using two indicators of present star formation rate (H$\alpha$ and infrared luminosity), we estimate the SFE to be a few percents. This is equivalent to the observational SFE in the disks of late-type spiral (Sb--) galaxies. This coincidence may support the universality of the mean SFE of spiral galaxies reported in the recent studies. That is, we find no evidence of distinct mechanism of the star formation in the central 1-kpc region of early-type galaxies. Also, we examine the structure of the central star-forming region, and discuss a method for estimating the mass of star-forming regions.

We construct a new algorithm for estimating the star formation rate (SFR) of galaxies from their infrared (IR) luminosity by developing the theory of the IR emission from a dusty H II region. The derived formula is ${\rm SFR}/(\MO {\rm yr^{-1}})=\{3.3\times 10^{-10}(1-\eta)/ (0.4-0.2f+0.6 \epsilon)\}(L{\rm ^{obs}_{IR}}/\LO)$, where $f$ is the fraction of ionizing photons absorbed by hydrogen, $\epsilon$ is the efficiency of dust absorption for nonionizing photons from OB stars, and $\eta$ is the cirrus fraction of observed IR luminosity. The previous conversion formulae of SFR from the IR luminosity is applicable to only the case where the observed IR luminosity is nearly equal to the bolometric luminosity (starburst galaxies etc.), except for some empirical formulae. On the other hand, our theoretical SFR is applicable to galaxies even with a moderate star formation activity. That is, our simple and convenient formula is significantly useful for estimating the SFR of various morphologies and types of galaxies --- from early elliptical to late spiral and irregular galaxies, or from active starburst to quiescent galaxies --- as far as they have neither an abnormal dust-to-gas ratio nor an evident active galactic nucleus.