Dispersion measure (DM) of Fast Radio Bursts (FRBs) are commonly used as a indicator of distance assuming that DM in excess of the expected amount within the Milky Way in the direction of each FRB arise mostly from the inter-galactic medium. However, the assumption might not be true if, for example, most FRB progenitors are embedded in ionized circumstellar material (CSM, e.g. supernova remnant). In this study, we jointly analyze distributions of DM, flux density, and fluence of the FRB samples observed by the Parkes telescope and the Australian Square Kilometre Array Pathfinder (ASKAP) using analytical models of FRBs, to constrain fractions of various DM components that shape the overall DM distribution and emission properties of FRBs. Comparing the model predictions with the observations we find that the typical amount of DM in each FRB host galaxy is $\sim 120$ cm$^{-3}$pc which is naturally explained as a combination of interstellar medium (ISM) and halo of an ordinary galaxy, without additional contribution from ionized CSM that is directly associated with an FRB progenitor. Furthermore, we also find that observed flux densities of FRBs do not statistically suffer strong $K$-correction, i.e. the typical luminosity density of FRBs does not significantly change within the range of emitting frequency $\nu_{\rm rest} \sim$ 1--4 GHz.

We present the results of CO observations toward 14 host galaxies of long-duration gamma-ray bursts (GRBs) at z = 0.1-2.5 by using the Atacama Large Millimeter/submillimeter Array. We successfully detected CO(3-2) or CO(4-3) emission in eight hosts (z = 0.3-2), which more than doubles the sample size of GRB hosts with CO detection. The derived molecular gas mass is $M_{\rm gas} = (0.2-6) \times 10^{10}$ $M_{\odot}$ assuming metallicity-dependent CO-to-H$_2$ conversion factors. By using the largest sample of GRB hosts with molecular gas estimates (25 in total, of which 14 are CO-detected) including results from the literature, we compared molecular gas properties with those of other star-forming galaxies (SFGs). The GRB hosts tend to have a higher molecular gas mass fraction ($\mu_{\rm gas}$) and a shorter gas depletion timescale ($t_{\rm depl}$) as compared with other SFGs at similar redshifts especially at $z \lesssim 1$. This could be a common property of GRB hosts or an effect introduced by the selection of targets which are typically above the main-sequence line. To eliminate the effect of selection bias, we analyzed $\mu_{\rm gas}$ and $t_{\rm depl}$ as a function of the distance from the main-sequence line ($\delta$MS). We find that the GRB hosts follow the same scaling relations as other SFGs, where $\mu_{\rm gas}$ increases and $t_{\rm depl}$ decreases with increasing $\delta {\rm MS}$. No molecular gas deficit is observed when compared to other SFGs of similar SFR and stellar mass. These findings suggest that the same star-formation mechanism is expected to be happening in GRB hosts as in other SFGs.

We present the clustering properties of low-$z$ $(z\leq1.4)$ galaxies selected by the Hyper Suprime-Cam Subaru Strategic Program Wide layer over $145$ deg$^{2}$. The wide-field and multi-wavelength observation yields $5,064,770$ galaxies at $0.3\leq z\leq1.4$ with photometric redshifts and physical properties. This enables the accurate measurement of angular correlation functions and subsequent halo occupation distribution (HOD) analysis allows the connection between baryonic properties and dark halo properties. The fraction of less-massive satellite galaxies at $z\lesssim1$ is found to be almost constant at $\sim20\%$, but it gradually decreases beyond $M_{\star} \sim 10^{10.4}h^{-2}M_{\odot}$. However, the abundance of satellite galaxies at $z>1$ is quite small even for less-massive galaxies due to the rarity of massive centrals at high-$z$. This decreasing trend is connected to the small satellite fraction of Lyman break galaxies at $z>3$. The stellar-to-halo mass ratios at $0.3\leq z\leq1.4$ are almost consistent with the predictions obtained using the latest empirical model; however, we identify small excesses from the theoretical model at the massive end. The pivot halo mass is found to be unchanged at $10^{11.9-12.1}h^{-1}M_{\odot}$ at $0.3\leq z\leq1.4$, and we systematically show that $10^{12}h^{-1}M_{\odot}$ is a universal pivot halo mass up to $z\sim5$ that is derived using only the clustering/HOD analyses. Nevertheless, halo masses with peaked instantaneous baryon conversion efficiencies are much smaller than the pivot halo mass regardless of a redshift, and the most efficient stellar-mass assembly is thought to be in progress in $10^{11.0-11.5}h^{-1}M_{\odot}$ dark haloes.

We present the results of CO(1-0) and CO(4-3) observations of the host galaxy of a long-duration gamma-ray burst GRB080207 at z = 2.0858 by using the Karl G. Jansky Very Large Array and the Atacama Large Millimeter/submillimeter Array. The host is detected in CO(1-0) and CO(4-3), becoming the first case for a GRB host with more than two CO transitions detected combined with CO(2-1) and CO(3-2) in the literature. Adopting a metallicity-dependent CO-to-H2 conversion factor, we derive a molecular gas mass of Mgas = 8.7 x 10^10 Modot, which places the host in a sequence of normal star-forming galaxies in a Mgas-star-formation rate (SFR) plane. A modified blackbody fit to the far-infrared--millimeter photometry results in a dust temperature of 37 K and a dust mass of Mdust = 1.5 x 10^8 Modot. The spatially-resolving CO(4-3) observations allow us to examine the kinematics of the host. The CO velocity field shows a clear rotation and is reproduced by a rotation-dominated disk model with a rotation velocity of 350 km/s and a half-light radius of 2.4 kpc. The CO spectral line energy distribution derived from the four CO transitions is similar to that of starburst galaxies, suggesting a high excitation condition. Comparison of molecular gas properties between the host and normal (main-sequence) galaxies at similar redshifts shows that they share common properties such as gas mass fraction, gas depletion timescale, gas-to-dust ratio, location in the Mgas-SFR (or surface density) relation, and kinematics, suggesting that long-duration GRBs can occur in normal star-forming environments at z ~ 2.

We have searched for optical variability in the host galaxy of the radio variable source which is possibly associated with fast radio burst (FRB) 150418. We compare images of the galaxy taken 1 day after the burst using Subaru/Suprime-Cam with images taken $\sim$ 1 year after the burst using Gemini-South/GMOS. No optical variability is found between the two epochs with a limiting absolute magnitude $\gtrsim -19$ (AB). This limit applies to optical variability of the putative active galactic nucleus in the galaxy on a timescale of $\sim$ 1 year, and also to the luminosity of an optical counterpart of FRB~150418 one day after the burst should it have occurred in this galaxy.

The origin of fast radio bursts (FRBs), bright millisecond radio transients, is still somewhat of a mystery. Several theoretical models expect that the FRB accompanies an optical afterglow (e.g., Totani 2013; Kashiyama et al. 2013). In order to investigate the origin of FRBs, we perform $gri$-band follow-up observations of FRB~151230 (estimated $z \lesssim 0.8$) with Subaru/Hyper Suprime-Cam at $8$, $11$, and $14$~days after discovery. The follow-up observation reaches a $50\%$ completeness magnitude of $26.5$~mag for point sources, which is the deepest optical follow-up of FRBs to date. We find $13$ counterpart candidates with variabilities during the observation. We investigate their properties with multicolor and multi-wavelength observations and archival catalogs. Two candidates are excluded by the non-detection of FRB~151230 in the other radio feed horns that operated simultaneously to the detection, as well as the inconsistency between the photometric redshift and that derived from the dispersion measure of FRB~151230. Eight further candidates are consistent with optical variability seen in AGNs. Two more candidates are well fitted with transient templates (Type IIn supernovae), and the final candidate is poorly fitted with all of our transient templates and is located off-center of an extended source. It can only be reproduced with rapid transients with a faint peak and rapid decline and the probability of chance coincidence is $\sim3.6\%$. We also find that none of our candidates are consistent with Type Ia supernovae, which rules out the association of Type Ia supernovae to FRB~151230 at $z\leq0.6$ and limits the dispersion measure of the host galaxy to $\lesssim300$~pc~cm$^{-3}$ in a Type Ia supernova scenario.

We investigate how the statistical properties of dispersion measure (DM) and apparent flux density/fluence of (non-repeating) fast radio bursts (FRBs) are determined by unknown cosmic rate density history [$\rhoFRB (z)$] and luminosity function (LF) of the transient events. We predict the distributions of DMs, flux densities, and fluences of FRBs taking account of the variation of the receiver efficiency within its beam, using analytical models of $\rhoFRB (z)$ and LF. Comparing the predictions with the observations, we show that the cumulative distribution of apparent fluences suggests that FRBs originate at cosmological distances and $\rhoFRB$ increases with redshift resembling cosmic star formation history (CSFH). We also show that a LF model with a bright-end cutoff at log$_{10}L_\nu$ [erg s$^{-1}$Hz$^{-1}$] $\sim$ 34 are favored to reproduce the observed DM distribution if $\rhoFRB (z)\propto$ CSFH, although the statistical significance of the constraints obtained with the current size of the observed sample is not high. Finally, we find that the correlation between DM and flux density of FRBs is potentially a powerful tool to distinguish whether FRBs are at cosmological distances or in the local universe more robustly with future observations.

We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $\alpha=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($\alpha=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4\pi$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate.

We present the cross-correlation between 151 luminous quasars ($M_{ \mathrm{UV}} < -26$) and 179 protocluster candidates at $z \sim 3.8$, extracted from the Wide imaging survey ($ \sim 121~ $deg$^2$) performed with a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We find that only two out of 151 quasars reside in regions that are more overdense compared to the average field at $ > 4 \sigma $. The distributions of the distance between quasars and the nearest protoclusters and the significance of the overdensity at the position of quasars are statistically identical to those found for $g$-dropout galaxies, suggesting that quasars tend to reside in almost the same environment as star-forming galaxies at this redshift. Using stacking analysis, we find that the average density of $g$-dropout galaxies around quasars is slightly higher than that around $g$-dropout galaxies on $1.0 - 2.5$ pMpc scales, while at $ < 0.5$ pMpc that around quasars tends to be lower. We also find that quasars with higher UV-luminosity or with more massive black holes tend to avoid the most overdense regions, and that the quasar near zone sizes are anti-correlated with overdensity. These findings are consistent with a scenario in which the luminous quasar at $z \sim4 $ resides in structures that are less massive than those expected for the progenitors of today's rich clusters of galaxies, and possibly that luminous quasars may be suppressing star formation in their close vicinity.

We constrain the quasar contribution to cosmic reionization based on our deep optical survey of z~6 quasars down to z_R=24.15 using Subaru/Suprime-Cam in three UKIDSS-DXS fields covering 6.5 deg^2. In Kashikawa et al. (2015), we select 17 quasar candidates and report our initial discovery of two low-luminosity quasars (M_1450~ -23) from seven targets, one of which might be a Lyman alpha emitting galaxy. From an additional optical spectroscopy, none of the four candidates out of the remaining ten turn out to be genuine quasars. Moreover, the deeper optical photometry provided by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) shows that, unlike the two already-known quasars, the i-z and z-y colors of the last six candidates are consistent with M- or L-type brown dwarfs. Therefore, the quasar luminosity function (QLF) in the previous paper is confirmed. Compiling QLF measurements from the literature over a wide magnitude range, including an extremely faint AGN candidate from Parsa et al. (2017}, to fit them with a double power-law, we find that the best-fit faint-end slope is alpha=-2.04^+0.33_-0.18 (-1.98^+0.48_-0.21) and characteristic magnitude is M_1450^*=-25.8^+1.1_-1.9 (-25.7^+1.0_-1.8) in the case of two (one) quasar detection. Our result suggests that, if the QLF is integrated down to M_1450=-18, quasars produce ~1-12% of the ionizing photons required to ionize the whole universe at z~6 with 2sigma confidence level, assuming that the escape fraction is f_esc=1 and the IGM clumpy factor is C=3. Even when the systematic uncertainties are taken into account, our result supports the scenario that quasars are the minor contributors of reionization.

We present the results of clustering analyses of Lyman break galaxies (LBGs) at $z\sim3$, $4$, and $5$ using the final data release of the Canada--France--Hawaii Telescope Legacy Survey (CFHTLS). Deep- and wide-field images of the CFHTLS Deep Survey enable us to obtain sufficiently accurate two-point angular correlation functions to apply a halo occupation distribution analysis. Mean halo masses, calculated as $\langle M_{h} \rangle=10^{11.7}-10^{12.8}h^{-1}M_{\odot}$, increase with stellar-mass limit of LBGs. The threshold halo mass to have a central galaxy follows the same increasing trend with the low-$z$ results, whereas the threshold halo mass to have a satellite galaxy shows higher values at $z=3-5$ than $z=0.5-1.5$ over the entire stellar mass range. Satellite fractions of dropout galaxies, even at less massive haloes, are found to drop sharply from $z=2$ down to less than $0.04$ at $z=3-5$. These results suggest that satellite galaxies form inefficiently within dark haloes at $z=3-5$ even for less massive satellites with $M_{\star}<10^{10}M_{\odot}$. We compute stellar-to-halo mass ratios (SHMRs) assuming a main sequence of galaxies, which is found to provide consistent SHMRs with those derived from a spectral energy distribution fitting method. The observed SHMRs are in good agreement with the model predictions based on the abundance-matching method within $1\sigma$ confidence intervals. We derive observationally, for the first time, $M_{{\rm h}}^{{\rm pivot}}$, which is the halo mass at a peak in the star-formation efficiency, at $3<z<5$, and it shows a little increasing trend with cosmic time at $z>3$. In addition, $M_{{\rm h}}^{{\rm pivot}}$ and its normalization are found to be almost unchanged during $0<z<5$. Our study shows an observational evidence that galaxy formation is ubiquitously most efficient near a halo mass of $M_{{\rm h}}\sim10^{12}M_{\odot}$ over cosmic time.

We discuss the host galaxy metallicity distribution of all long gamma-ray bursts (GRBs) whose redshifts are known to be $< 0.4$, including newly obtained spectroscopic datasets of the host galaxies of GRB 060614, 090417B, and 130427A. We compare the metallicity distribution of the low-redshift sample to the model predictions, and constrain the relation between metallicity and GRB occurrence. We take account of spatial variation of metallicities among star forming regions within a galaxy. We found that the models, in which only low-metallicity stars produce GRBs with a sharp cutoff of GRB production efficiency around 12+log(O/H) $\sim$ 8.3, can well reproduce the observed distribution, while the models with no metallicity dependence are not consistent with the observations. We also discuss possible sampling biases we may suffer by collecting long GRBs whose redshifts are known, presenting the photometric observations of the host galaxy of GRB 111225A at $z = 0.297$ whose redshift has been undetermined until $\sim$ 2.3 years after the burst.

The cosmological evolution of the binary black hole (BH) merger rate and the energy density of the gravitational-wave (GW) background are investigated. To evaluate the redshift dependence of the BH formation rate, BHs are assumed to originate from low-metallicity stars, and the relations between the star formation rate, metallicity and stellar mass of galaxies are combined with the stellar mass function at each redshift. As a result, it is found that when the energy density of the GW background is scaled with the merger rate at the local Universe, the scaling factor does not depend on the critical metallicity for the formation of BHs. Also taking into account the merger of binary neutron stars, a simple formula to express the energy spectrum of the GW background is constructed for the inspiral phase. The relation between the local merger rate and the energy density of the GW background will be examined by future GW observations.

We present the discovery of three protoclusters at $z\sim3\mathrm{-}4$ with spectroscopic confirmation in the Canada-France-Hawaii Telescope (CFHT) Legacy Survey Deep Fields. In these fields, we investigate the large-scale projected sky distribution of $z\sim3\mathrm{-}6$ Lyman break galaxies and identify 21 protocluster candidates from regions that are overdense at more than $4\sigma$ overdensity significance. Based on cosmological simulations, it is expected that more than $76\%$ of these candidates will evolve into a galaxy cluster of at least a halo mass of $10^{14}\,\mathrm{M_\odot}$ at $z=0$. We perform follow-up spectroscopy for eight of the candidates using Subaru/FOCAS, KeckII/DEIMOS, and Gemini-N/GMOS. In total we target 462 dropout candidates and obtain 138 spectroscopic redshifts. We confirm three real protoclusters at $z=3\mathrm{-}4$ with more than five members spectroscopically identified, and find one to be an incidental overdense region by mere chance alignment. The other four candidate regions at $z\sim5\mathrm{-}6$ require more spectroscopic follow-up in order to be conclusive. A $z=3.67$ protocluster, which has eleven spectroscopically confirmed members, shows a remarkable core-like structure composed of a central small region ($<0.5\,\mathrm{physical\>Mpc}$) and an outskirts region ($\sim1.0\,\mathrm{physical\>Mpc}$). The Ly$\alpha$ equivalent widths of members of the protocluster are significantly smaller than those of field galaxies at the same redshift while there is no difference in the UV luminosity distributions. These results imply that some environmental effects start operating as early as at $z\sim4$ along with the growth of the protocluster structure.

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $\Omega_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

The galaxy population hosting long-duration GRBs provides a means to constrain the progenitor and an opportunity to use these violent explosions to characterize the nature of the high-redshift universe. Studies of GRB host galaxies in emission reveal a population of star-forming galaxies with great diversity, spanning a wide range of masses, metallicities, and redshifts. However, as a population GRB hosts are significantly less massive and poorer in metals than the hosts of other core-collapse transients, suggesting that GRB production is only efficient at metallicities significantly below Solar. GRBs may also prefer compact galaxies, and dense and/or central regions of galaxies, more than other types of core-collapse explosion. Meanwhile, studies of hosts in absorption against the luminous GRB optical afterglow provide a unique means of unveiling properties of the ISM in even the faintest and most distant galaxies; these observations are helping to constrain the chemical evolution of galaxies and the properties of interstellar dust out to very high redshifts. New ground- and space-based instrumentation, and the accumulation of larger and more carefully-selected samples, are continually enhancing our view of the GRB host population.

The spectrum of the supernova relic neutrino (SRN) background from past stellar collapses including black hole formation (failed supernovae) is calculated. The redshift dependence of the black hole formation rate is considered on the basis of the metallicity evolution of galaxies. Assuming the mass and metallicity ranges of failed supernova progenitors, their contribution to SRNs is quantitatively estimated for the first time. Using this model, the dependences of SRNs on the cosmic star formation rate density, shock revival time and equation of state are investigated. The shock revival time is introduced as a parameter that should depend on the still unknown explosion mechanism of core collapse supernovae. The dependence on equation of state is considered for failed supernovae, whose collapse dynamics and neutrino emission are certainly affected. It is found that the low-energy spectrum of SRNs is mainly determined by the cosmic star formation rate density. These low-energy events will be observed in the Super-Kamiokande experiment with gadolinium-loaded water.

We examine how the small-scale ($<$\,kpc) variation of metallicity within a galaxy, which is found in nearby galaxies, affect the observational estimates of metallicity in the explosion sites of transient events such as core-collapse supernovae (CC SNe) and gamma-ray bursts (GRBs). Assuming the same luminosity, metallicity, and spatial distributions of \ion{H}{ii}\ regions (hereafter HIIRs) as observed in M31, we compute the apparent metallicities that we would obtain when the spectrum of a target region is blended with those of surrounding HIIRs within the length scale of typical spatial resolution. When the spatial resolution of spectroscopy is $\lesssim$ 0.5\,kpc, which is typical for the existing studies of CC SN sites, we find that the apparent metallicities reflect the metallicities of target regions, but with significant systematic uncertainties in some cases. When the spatial resolution is $\gtrsim$ 1.0 kpc, regardless of the target regions (which has a wide range of metallicity that spans $\sim$ 0.6 dex for the M31 HIIRs), we always obtain the apparent metallicities similar to the average metallicity of the M31 HIIRs. Given that the apparent metallicities measured with $\gtrsim$\,kpc scale resolution do not necessarily reflect the immediate environment of the stellar explosions, the current observational estimates of high metallicities for some of the long GRB host galaxies do {\em not} rule out the hypothesis that the long GRBs are exclusively born in a low-metallicity environment.

We present the discovery of one or two extremely faint z~6 quasars in 6.5 deg^2 utilizing a unique capability of the wide-field imaging of the Subaru/Suprime-Cam. The quasar selection was made in (i'-z_B) and (z_B-z_R) colors, where z_B and z_R are bandpasses with central wavelengths of 8842A and 9841A, respectively. The color selection can effectively isolate quasars at z~6 from M/L/T dwarfs without the J-band photometry down to z_R<24.0, which is 3.5 mag. deeper than SDSS. We have selected 17 promising quasar candidates. The follow-up spectroscopy for seven targets identified one apparent quasar at z=6.156 with M_1450=-23.10. We also identified one possible quasar at z=6.041 with a faint continuum of M_1450=-22.58 and a narrow Lyman-alpha emission with HWHM=427 km/s, which cannot be distinguished from Lyman-alpha emitters. We derive the quasar luminosity function at z~6 by combining our faint quasar sample with the bright quasar samples by SDSS and CFHQS. Including our data points invokes a higher number density in the faintest bin of the quasar luminosity function than the previous estimate employed. This suggests a steeper faint-end slope than lower-z, though it is yet uncertain based on a small number of spectroscopically identified faint quasars and several quasar candidates are still remain to be diagnosed. The steepening of the quasar luminosity function at the faint-end does increase the expected emission rate of the ionizing photon, however, it only changes by a factor of ~2-6. This was found to be still insufficient for the required photon budget of reionization at z~6.

We present near-infrared spectroscopy of the host galaxy of dark GRB 080325 using Subaru/MOIRCS. The obtained spectrum provides a clear detection of H$\alpha$ emission and marginal [NII]$\lambda$6584. The host is a massive (M$_{*}\sim10^{11}$M$_{\odot}$), dusty ($A_{V}\sim 1.2$) star-forming galaxy at z=1.78. The star formation rate calculated from the H$\alpha$ luminosity (35.6-47.0 M$_{\odot}$ yr$^{-1}$) is typical among GRB host galaxies (and star-forming galaxies generally) at z $>$1; however, the specific star formation rate is lower than normal star-forming galaxies at redshift $\sim$ 1.6, in contrast to the high specific star formation rates measured for many of other GRB hosts. The metallicity of the host is estimated to be 12+log(O/H)$_{\rm KK04}$$=$8.88. We emphasize that this is one of the most massive distant host galaxies for which metallcity is measured with emission-line diagnostics. The metallicity is fairly high among GRB hosts. However, this is still lower than the metallicity of normal star-forming galaxies of the same mass at z$\sim$1.6. The metallicity offset from normal star-forming galaxies is close to a typical value of other GRB hosts and indicates that GRB host galaxies are uniformly biased toward low metalicity over a wide range of redshift and stellar mass. The low-metallicity nature of the GRB 080325 host is likely not attributable to the fundamental metallicity relation of star-forming galaxies beacuse it is a metal-poor outlier from the relation and has a low sSFR. Thus we conclude that metallicity is important to the mechanism that produced this GRB.

We discuss how we can detect and identify counterparts of fast radio bursts (FRBs) in future optical follow up observations of FRBs if real-time alert of FRBs becomes available. We consider kilonovae as candidates of FRB optical counterparts, as expected in the case that FRBs originate from mergers of double neutron star binaries. Although theoretical predictions on luminosities of kilonovae are still highly uncertain, recent models suggest that kilonovae can be detected at redshifts up to z $\sim$ 0.3 within the range of the uncertainties. We expect $\sim$ 1--5 unrelated supernovae (SNe) down to a similar variability magnitude in 5 days interval within the typical error radius of a FRB. We show that, however, a kilonova can be distinguished from these SNe by its rapid decay and/or color evolution, making it possible to verify the existence of a kilonova associated with a FRB. We also discuss the case that SNe Ia are FRB optical counterparts, as it might be if FRBs originate from double white dwarf binaries. Verification of this scenario is also possible, since the chance probability of finding a SNe Ia having consistent explosion time with that of a FRB within the FRB error region is small (typically $\lesssim$ 0.01).

We performed a systematic spectroscopic observation of a protocluster at $z=6.01$ in the Subaru Deep Field. We took spectroscopy for all 53 $i'$-dropout galaxies down to $z'=27.09\,\mathrm{mag}$ in/around the protocluster region. From these observations, we confirmed that 28 galaxies are at $z\sim6$, of which ten are clustered in a narrow redshift range of $\Delta z<0.06$. To trace the evolution of this primordial structure, we applied the same $i'$-dropout selection and the same overdensity measurements used in the observations to a semi-analytic model built upon the Millennium Simulation. We obtain a relation between the significance of overdensities observed at $z\sim6$ and the predicted dark matter halo mass at $z=0$. This protocluster with $6\sigma$ overdensity is expected to grow into a galaxy cluster with a mass of $\sim5\times10^{14}\,\mathrm{M_\odot}$ at $z=0$. Ten galaxies within $10\,\mathrm{comoving\>Mpc}$ of the overdense region can, with more than an 80% probability, merge into a single dark matter halo by $z=0$. No significant differences appeared in UV and Ly$\alpha$ luminosities between the protocluster and field galaxies, suggesting that this protocluster is still in the early phase of cluster formation before the onset of any obvious environmental effects. However, further observations are required to study other properties, such as stellar mass, dust, and age. We do find that galaxies tend to be in close pairs in this protocluster. These pair-like subgroups will coalesce into a single halo and grow into a more massive structure. We may witness an onset of cluster formation at $z\sim6$ toward a cluster as seen in local universe.

The unprecedentedly bright optical afterglow of GRB 130606A located by Swift at a redshift close to the reionization era (z = 5.913) provides a new opportunity to probe the ionization status of intergalactic medium (IGM). Here we present an analysis of the red Ly alpha damping wing of the afterglow spectrum taken by Subaru/FOCAS during 10.4-13.2 hr after the burst. We find that the minimal model including only the baseline power-law and HI absorption in the host galaxy does not give a good fit, leaving residuals showing concave curvature in 8400-8900 A with an amplitude of about 0.6% of the flux. Such a curvature in the short wavelength range cannot be explained either by extinction at the host with standard extinction curves, intrinsic curvature of afterglow spectra, or by the known systematic uncertainties in the observed spectrum. The red damping wing by intervening HI gas outside the host can reduce the residual by about 3 sigma statistical significance. We find that a damped Ly alpha system is not favored as the origin of this intervening HI absorption, from the observed Ly beta and metal absorption features. Therefore absorption by diffuse IGM remains as a plausible explanation. A fit by a simple uniform IGM model requires HI neutral fraction of f_HI ~ 0.1-0.5 depending on the distance to the GRB host, implying high f_HI IGM associated with the observed dark Gunn-Peterson (GP) troughs. This gives a new evidence that the reionization is not yet complete at z ~ 6.

We present a new model of the extragalactic background light (EBL) and corresponding gamma-gamma opacity for intergalactic gamma-ray absorption from z = 0 up to z = 10, based on a semi-analytical model of hierarchical galaxy formation that reproduces key observed properties of galaxies at various redshifts. Including the potential contribution from Population III stars and following the cosmic reionization history in a simplified way, the model is also broadly consistent with available data concerning reionization, particularly the Thomson scattering optical depth constraints from WMAP. In comparison with previous EBL studies up to z ~ 3-5, our predicted gamma-gamma opacity is in general agreement for observed gamma-ray energy below 400/(1 + z) GeV, whereas it is a factor of ~ 2 lower above this energy because of a correspondingly lower cosmic star formation rate, even though the observed UV luminosity is well reproduced by virtue of our improved treatment of dust obscuration and direct estimation of star formation rate. The horizon energy at which the gamma-ray opacity is unity does not evolve strongly beyond z ~ 4 and approaches ~ 20 GeV. The contribution of Population III stars is a minor fraction of the EBL at z = 0, and is also difficult to distinguish through gamma-ray absorption in high-z objects, even at the highest levels allowed by the WMAP constraints. Nevertheless, the attenuation due to Population II stars should be observable in high-z gamma-ray sources by telescopes such as Fermi or CTA and provide a valuable probe of the evolving EBL in the rest-frame UV. The detailed results of our model are publicly available in numerical form at the URL http://www.slac.stanford.edu/~yinoue/Download.html.

We investigate the relation between stellar mass ($M_\star$), star formation rate (SFR), and metallicity ($Z$) of galaxies, so called the fundamental metallicity relation, in the galaxy sample of the Sloan Digital Sky Survey Data Release 7. We separate the galaxies into narrow redshift bins and compare the relation at different redshifts, and find statistically significant ($> 99$%) evolution. We test various observational effects that might cause seeming $Z$ evolution, and find it difficult to explain the evolution of the relation only by the observational effects. In the current sample of low redshift galaxies, galaxies with different $M_\star$ and SFR are sampled from different redshifts, and there is degeneracy between $M_\star$/SFR and redshift. Hence it is not straightforward to distinguish a relation between $Z$ and SFR from a relation between $Z$ and redshift. The separation of the intrinsic relation from the redshift evolution effect is a crucial issue to understand evolution of galaxies.

We present results of a search for a supernova (SN) component associated with GRB 100418A at the redshift of 0.624. The field of GRB 100418A was observed with FOCAS on Subaru 8.2m telescope under a photometric condition (seeing 0.3"-0.4") on 2010 May 14 (UT). The date corresponds to 25.6 days after the burst trigger (15.8 days in the restframe). We did imaging observations in V, Rc, and Ic bands, and two hours of spectrophotometric observations. We got the resolved host galaxy image which elongated 1.6" (= 11 kpc) from north to south. No point source was detected on the host galaxy. The time variation of Rc-band magnitude shows that the afterglow of GRB 100418A has faded to Rc \sim > 24 without SN like rebrightening, when we compare our measurement to the reports in GCN circulars. We could not identify any SN feature such as broad emission-lines or bumps in our spectrum. Assuming the SN is fainter than the 3{\sigma} noise spectrum of our observation, we estimate the upper limit on the SN absolute magnitude MIc,obs > -17.2 in observer frame Ic-band. This magnitude is comparable to the faintest type Ic SNe. We also estimate host galaxy properties from the spectrum. The host galaxy of GRB 100418A is relatively massive (log M_{star}/M_{sun} = 9.54) compared to typical long GRB host galaxies, and has 12+log(O/H) = 8.75.

We predict the metallicity probability distribution function (PDF) of long gamma-ray burst (GRB) host galaxies at low-redshifts ($z \la 0.3$) when GRBs occur only in low-metallicity environment, assuming empirical formulations of galaxy properties. We discuss contribution of high-metallicity galaxies to the cosmic rate of low-metallicity GRBs, taking internal dispersion of metallicity within each galaxy into account. Assuming GRBs trace low-metallicity star formation $< \Zcrit$: 12+log$_{10}$(O/H) = 8.2, we find that GRB host galaxies may have systematically higher-metallicity than that of GRB progenitors. Furthermore, we expect $\ga$ 10% of the host galaxies to have 12+log$_{10}$(O/H) $> 8.8$, if galaxies have internal dispersion of metallicity comparable to that observed in the Milky Way. Our results show that the low-metallicity scenario of GRB progenitors can be reconciled with the recent discoveries of the high-metallicity host galaxies of GRBs. We also show possible bimodality in the host metallicity PDF that results from the single progenitor model of GRBs. If found in future observation, the bimodality can be a clue to constrain the nature of GRB progenitors.

We study the relationship between the metallicity of gamma-ray burst (GRB) progenitors and the probability distribution function (PDF) of GRB host galaxies as a function of luminosity using cosmological hydrodynamic simulations of galaxy formation. We impose a maximum limit to the gas metallicity in which GRBs can occur, and examine how the predicted luminosity PDF of GRB host galaxies changes in the simulation. We perform the Kolmogorov-Smirnov test, and show that the result from our simulation agrees with the observed luminosity PDF of core-collapse supernovae (SNe) host galaxies when we assume that the core-collapse SNe trace star formation. When we assume that GRBs occur only in a low-metallicity environment with $Z\lesssim 0.1 \Zsun$, GRBs occur in lower luminosity galaxies, and the simulated luminosity PDF becomes quantitatively consistent with the observed luminosity PDF. The observational bias against the host galaxies of optically dark GRBs owing to dust extinction may be another reason for the lower luminosities of GRB host galaxies, but the observed luminosity PDF of GRB host galaxies cannot be reproduced solely by the dust bias in our simulation.

We present optical and near infrared observations of GRB 080325 classified as a "Dark GRB". Near-infrared observations with Subaru/MOIRCS provided a clear detection of afterglow in Ks band, although no optical counterpart was reported. The flux ratio of rest-wavelength optical to X-ray bands of the afterglow indicates that the dust extinction along the line of sight to the afterglow is Av = 2.7 - 10 mag. This large extinction is probably the major reason for optical faintness of GRB 080325. The J - Ks color of the host galaxy, (J - Ks = 1.3 in AB magnitude), is significantly redder than those for typical GRB hosts previously identified. In addition to J and Ks bands, optical images in B, Rc, i', and z' bands with Subaru/Suprime-Cam were obtained at about one year after the burst, and a photometric redshift of the host is estimated to be z_{photo} = 1.9. The host luminosity is comparable to L^{*} at z \sim 2 in contrast to the sub-L^{*} property of typical GRB hosts at lower redshifts. The best-fit stellar population synthesis model for the host shows that a large dust extinction (Av = 0.8 mag) attributes to the red nature of the host and that the host galaxy is massive (M_{*} = 7.0 \times 10^{10} Msun) which is one of the most massive GRB hosts previously identified. By assuming that the mass-metallicity relation for star-forming galaxies at z \sim 2 is applicable for the GRB host, this large stellar mass suggests the high metallicity environment around GRB 080325, consistent with inferred large extinction.

Long gamma-ray bursts (GRBs) have been suggested to occur preferentially in low metallicity environment. We discuss the possibility and theoretical aspects of using Lyman alpha emission properties of long GRB host galaxies as a metallicity indicator of high redshift GRB environments, where direct metallicity measurements are not easy. We propose to use the fraction of Ly-alpha emitters (LAEs) in long GRB host galaxies as a function of UV luminosity, which can be compared with star-formation-rate weighted LAE fraction of Lyman-break galaxies as the standard in the case of no metallicity dependence. There are two important effects of metallicity dependence of long GRB rate to change the LAE fraction of host galaxies. One is the enhancement of intrinsic Ly-alpha equivalent width (EW) by stronger ionizing UV luminosity of low metallicity stellar population, and the other is extinction by interstellar dust to change the observable EW. Based on a latest theoretical model of LAEs that reproduce observations, we argue that the latter is likely to work in the opposite direction to the former, i.e., to decrease LAE fraction if GRBs preferentially occur in low-metallicity environments, because of the clumpy inter-stellar medium effect. The high LAE fraction of GRB host galaxies indicated by observations is quantitatively be explained by the LAE model if GRBs occur when Z < 0.1Z_sun, although this result is still indicative because of the limited statistics and theoretical uncertainties. This result demonstrates that the LAE statistics of GRB hosts may give us useful information in the future.

One of the possible origins of short gamma-ray bursts (SGRBs) is merging of compact binaries, and the effect of large kick velocity is a signature that can be used as an observational test for this hypothesis. Intracluster SGRBs that escaped from a host galaxy in a galaxy cluster are interesting in this context, since they would escape more easily by cluster tidal force, and would have brighter afterglow luminosity by dense intracluster gas, than those in general field galaxies. Here we calculate the escape fraction of compact binaries from their host galaxies in a galaxy cluster, and discuss some observational implications. We found that the escape fraction strongly depends on the nature of dark matter subhalos associated with member galaxies. If the amount of dark matter around member galaxies is not large and the gravitational potential for an escaping binary is determined mostly by stellar mass, most of SGRBs should escape and be observed as hostless, which is a much higher fraction than those in the field. Hence, statistics of intracluster SGRBs could give important information about the dark matter distribution in galaxy clusters, as well as hints for the origin of SGRBs.