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.

This is the second paper of a series to shed light on how galactic properties depend on the intergalactic medium (IGM) environment traced by the Ly$\alpha$ forest. In this paper, we observationally investigate the IGM--galaxy connection using the publicly available 3D IGM tomography data (CLAMATO: Lee et al. 2016, 2018) and several galaxy catalogs in the COSMOS field. We measure the cross-correlation function (CCF) for $570$ galaxies with spec-$z$ measurements, and detect a correlation with the IGM up to $50$ $h^{-1}$ comoving Mpc. We show that galaxies with stellar masses of $10^9-10^{10}$ M$_\odot$ are the dominant contributor to the total CCF signal. We also investigate CCFs for several galaxy populations: Ly$\alpha$ emitters (LAEs), H$\alpha$ emitters (HAEs), [OIII] emitters (O3Es), active galactic nuclei (AGNs), and submillimeter galaxies (SMGs), and detect the highest signal in AGNs and SMGs at large scales ($r\geq5$ $h^{-1}$ Mpc), but in LAEs at small scales ($r<5$ $h^{-1}$ Mpc). We find that they live in various IGM environments --HAEs trace the IGM in a similar manner to the continuum-selected galaxies, but LAEs and O3Es tend to reside in higher density regions. Additionally, LAEs' CCF is flat up to $r\sim3$ $h^{-1}$ Mpc, indicating that they tend to avoid the highest density regions. For AGNs and SMGs, the CCF peak at $r=5-6$ $h^{-1}$ Mpc implies that they tend to be in locally lower density regions. We suspect that it is due to the IGM HI photoionization by AGNs, i.e., the proximity effect.

This is the first paper in a series to systematically investigate the environmental dependence of galactic properties traced by Ly$\alpha$ forest absorption in the interagalactic medium (IGM). Using our cosmological hydrodynamic simulations, we investigate the IGM--galaxy connection at $z=2$ by two methods: (I) cross-correlation analysis between galaxies and the fluctuation of Ly$\alpha$ forest transmission ($\delta_\text{F}$); and (II) comparing the overdensity of neutral hydrogen (HI) and galaxies. Our simulations reproduce observed cross-correlation functions (CCF) between the Ly$\alpha$ forest and Lyman-break galaxies. We further investigate the variation of the CCF using subsamples divided by dark matter halo mass ($M_\text{DH}$), galaxy stellar mass ($M_*$), and star-formation rate (SFR), and find that the CCF signal becomes stronger with increasing $M_\text{DH}$, $M_*$, and SFR. The CCFs between galaxies and gas-density fluctuation are also found to have similar trends. Therefore, the variation of the $\delta_\text{F}$--CCF depending on $M_\text{DH}$, $M_*$, and SFR is due to varying gas density around galaxies. We conclude that $M_\text{DH}$ is the most sensitive parameter to characterize the correlation between galaxies and the IGM distribution. Method (II) also finds correlations between galaxies and HI overdensities. Overall, our analyses confirm the spatial correlation between galaxies and IGM HI, with more massive galaxies being clustered in higher-density regions as expected from the $\Lambda$CDM paradigm.

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

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.

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.

We measure the scale dependence and redshift dependence of 21 cm line emitted from the neutral hydrogen gas at redshift 1<z<5 using full cosmological hydrodynamic simulations by taking the ratios between the power spectra of HI-dark matter cross correlation and dark matter auto-correlation. The neutral hydrogen distribution is computed in full cosmological hydrodynamic simulations including star formation and supernova feedback under a uniform ultra-violet background radiation. We find a significant scale dependence of HI bias at z>3 on scales of k>1h/Mpc, but it is roughly constant at lower redshift z<3. The redshift evolution of HI bias is relatively slow compared to that of QSOs at similar redshift range. We also measure a redshift space distortion (RSD) of HI gas to explore the properties of HI clustering. Fitting to a widely applied theoretical prediction, we find that the constant bias is consistent with that measured directly from the real-space power spectra, and the velocity dispersion is marginally consistent with the linear perturbation prediction. Finally we compare the results obtained from our simulation and the Illustris simulation, and conclude that the detailed astrophysical effects do not affect the scale dependence of HI bias very much, which implies that the cosmological analysis using 21 cm line of HI will be robust against the uncertainties arising from small-scale astrophysical processes such as star formation and supernova feedback.

Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of infrared (IR)-emitting region, IR luminosity function and IRX--$\beta$ relation. After confirming that our models reproduce the local galaxy properties, we examine the STUDIES sample at $z\approx 1-4$, finding that the simulation reproduces the aforementioned quantities except for the $z\gtrsim 2$ IR luminosity function at the massive end ($\sim 10^{13}$ L$_{\odot}$). This means that the current simulation correctly reproduces the overall scaling between the size and luminosity (or star formation rate) of dusty region, but lacks extreme starburst phenomena at $z\gtrsim 2$. We also discuss extinction curves and possible AGN contribution.

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

To study the dust evolution in the cosmological structure formation history, we perform a smoothed particle hydrodynamic simulation with a dust enrichment model in a cosmological volume. We adopt the dust evolution model that represents the grain size distribution by two sizes and takes into account stellar dust production and interstellar dust processing. We examine the dust mass function and the scaling properties of dust in terms of the characteristics of galaxies. The simulation broadly reproduces the observed dust mass functions at redshift $z = 0$, except that it overproduces the massive end at dust mass $M_\mathrm{d} \gtrsim 10^{8}$ ${\rm M}_\odot$. This overabundance is due to overproducing massive gas/metal-rich systems, but we also note that the relation between stellar mass and gas-phase metallicity is reproduced fairly well by our recipe. The relation between dust-to-gas ratio and metallicity shows a good agreement with the observed one at $z=0$, which indicates successful implementation of dust evolution in our cosmological simulation. Star formation consumes not only gas but also dust, causing a decreasing trend of the dust-to-stellar mass ratio at the high-mass end of galaxies. We also examine the redshift evolution up to $z \sim~ 5$, and find that the galaxies have on average the highest dust mass at $z = 1-2$. For the grain size distribution, we find that galaxies with metallicity $\sim 0.3~ Z_\odot$ tend to have the highest small-to-large grain abundance ratio; consequently, the extinction curves in those galaxies have the steepest ultraviolet slopes.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

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.

In this paper, we study the effects of a recent drop in star formation rate (SFR) on the spectra of epoch of reionization (EoR) galaxies, and the resulting degeneracy with the spectral features produced by extreme Lyman continuum leakage. In order to study these effects in the wavelength range relevant for the upcoming James Webb Space Telescope (JWST), we utilize synthetic spectra of simulated EoR galaxies from cosmological simulations together with synthetic spectra of partially quenched mock galaxies. We find that rapid declines in the SFR of EoR galaxies could seriously affect the applicability of methods that utilize the equivalent width of Balmer lines and the ultraviolet spectral slope to assess the escape fraction of EoR galaxies. In order to determine if the aforementioned degeneracy can be avoided by using the overall shape of the spectrum, we generate mock NIRCam observations and utilize a classification algorithm to identify galaxies that have undergone quenching. We find that while there are problematic cases, JWST/NIRCam or NIRSpec should be able to reliably identify galaxies with redshifts $z\sim 7$ that have experienced a significant decrease in the SFR (by a factor 10-100) in the past 50-100 Myr with a success rate $\gtrsim 85\%$. We also find that uncertainties in the dust-reddening effects on EoR galaxies significantly affect the performance of the results of the classification algorithm. We argue that studies that aim to characterize the dust extinction law most representative in the EoR would be extremely useful.

To investigate the evolution of dust in a cosmological volume, we perform hydrodynamic simulations, in which the enrichment of metals and dust is treated self-consistently with star formation and stellar feedback. We consider dust evolution driven by dust production in stellar ejecta, dust destruction by sputtering, grain growth by accretion and coagulation, and grain disruption by shattering, and treat small and large grains separately to trace the grain size distribution. After confirming that our model nicely reproduces the observed relation between dust-to-gas ratio and metallicity for nearby galaxies, we concentrate on the dust abundance over the cosmological volume in this paper. The comoving dust mass density has a peak at redshift $z\sim 1$--2, coincident with the observationally suggested dustiest epoch in the Universe. {In the local Universe}, roughly 10 per cent of the dust is contained in the intergalactic medium (IGM), where only 1/3--1/4 of the dust survives against dust destruction by sputtering. We also show that the dust mass function is roughly reproduced at $\lesssim 10^8$ M$_\odot$, while the massive end still has a discrepancy, which indicates {the necessity of stronger feedback in massive galaxies}. %%The relation showed that accretion is essential for dusty galaxies. In addition, our model broadly reproduces the observed radial profile of dust surface density in the circum-galactic medium (CGM). While our model satisfies the observational constraints for the dust extinction {on cosmological scales}, it predicts that the dust in the CGM and IGM is dominated by large ($> 0.03~\mu$m) grains, which is in tension with the steep reddening curves {observed} in the CGM.

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.

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

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

There are two major theoretical issues for the star formation law (the relation between the surface densities of molecular gas and star formation rate on a galaxy scale): (i) At low metallicity, it is not obvious that star-forming regions are rich in H$_2$ because the H$_2$ formation rate depends on the dust abundance; and (ii) whether or not CO really traces H$_2$ is uncertain, especially at low metallicity. To clarify these issues, we use a hydrodynamic simulation of an isolated disc galaxy with a spatial resolution of a few tens parsecs. The evolution of dust abundance and grain size distribution is treated consistently with the metal enrichment and the physical state of the interstellar medium. We compute the H$_2$ and CO abundances using a subgrid post-processing model based on the dust abundance and the dissociating radiation field calculated in the simulation. We find that when the metallicity is $\lesssim 0.4$ Z$_\odot$ ($t<1$ Gyr), H$_2$ is not a good tracer of star formation rate because H$_2$-rich regions are limited to dense compact regions. At $Z\gtrsim 0.8$ Z$_\odot$, a tight star formation law is established for both H$_2$ and CO. At old ($t \sim 10$ Gyr) ages, we also find that adopting the so-called MRN grain size distribution with an appropriate dust-to-metal ratio over the entire disc gives reasonable estimates for the H$_2$ and CO abundances. For CO, improving the spatial resolution of the simulation is important while the H$_2$ abundance is not sensitive to sub-resolution structures at $Z\gtrsim 0.4$ Z$_\odot$.

To understand the evolution of extinction curve, we calculate the dust evolution in a galaxy using smoothed particle hydrodynamics simulations incorporating stellar dust production, dust destruction in supernova shocks, grain growth by accretion and coagulation, and grain disruption by shattering. The dust species are separated into carbonaceous dust and silicate. The evolution of grain size distribution is considered by dividing grain population into large and small gains, which allows us to estimate extinction curves. We examine the dependence of extinction curves on the position, gas density, and metallicity in the galaxy, and find that extinction curves are flat at $t \lesssim 0.3$ Gyr because stellar dust production dominates the total dust abundance. The 2175 \AA\ bump and far-ultraviolet (FUV) rise become prominent after dust growth by accretion. At $t \gtrsim 3$ Gyr, shattering works efficiently in the outer disc and low density regions, so extinction curves show a very strong 2175 \AA\ bump and steep FUV rise. The extinction curves at $t\gtrsim 3$ Gyr are consistent with the Milky Way extinction curve, which implies that we successfully included the necessary dust processes in the model. The outer disc component caused by stellar feedback has an extinction curves with a weaker 2175 \AA\ bump and flatter FUV slope. The strong contribution of carbonaceous dust tends to underproduce the FUV rise in the Small Magellanic Cloud extinction curve, which supports selective loss of small carbonaceous dust in the galaxy. The snapshot at young ages also explain the extinction curves in high-redshift quasars.

We perform smoothed particle hydrodynamics (SPH) simulations of an isolated galaxy with a new treatment for dust formation and destruction. To this aim, we treat dust and metal production self-consistently with star formation and supernova feedback. For dust, we consider a simplified model of grain size distribution by representing the entire range of grain sizes with large and small grains. We include dust production in stellar ejecta, dust destruction by supernova (SN) shocks, grain growth by accretion and coagulation, and grain disruption by shattering. We find that the assumption of fixed dust-to-metal mass ratio becomes no longer valid when the galaxy is older than 0.2 Gyr, at which point the grain growth by accretion starts to contribute to the nonlinear rise of dust-to-gas ratio. As expected in our previous one-zone model, shattering triggers grain growth by accretion since it increases the total surface area of grains. Coagulation becomes significant when the galaxy age is greater than $\sim$ 1 Gyr: at this epoch the abundance of small grains becomes high enough to raise the coagulation rate of small grains. We further compare the radial profiles of dust-to-gas ratio $(\mathcal{D})$ and dust-to-metal ratio $(\mathcal{D}/Z)$ (i.e., depletion) at various ages with observational data. We find that our simulations broadly reproduce the radial gradients of dust-to-gas ratio and depletion. In the early epoch ($\lesssim 0.3$ Gyr), the radial gradient of $\mathcal{D}$ follows the metallicity gradient with $\mathcal{D}/Z$ determined by the dust condensation efficiency in stellar ejecta, while the $\mathcal{D}$ gradient is steeper than the $Z$ gradient at the later epochs because of grain growth by accretion. The framework developed in this paper is applicable to any SPH-based galaxy evolution simulations including cosmological ones.

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.

We present a search for low energy antineutrino events coincident with the gravitational wave events GW150914 and GW151226, and the candidate event LVT151012 using KamLAND, a kiloton-scale antineutrino detector. We find no inverse beta-decay neutrino events within $\pm 500$ seconds of either gravitational wave signal. This non-detection is used to constrain the electron antineutrino fluence and the luminosity of the astrophysical sources.

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 present basic properties of $\sim$3,300 emission line galaxies detected by the FastSound survey, which are mostly H$\alpha$ emitters at $z \sim$ 1.2-1.5 in the total area of about 20 deg$^2$, with the H$\alpha$ flux sensitivity limit of $\sim 1.6 \times 10^{-16} \rm erg \ cm^{-2} s^{-1}$ at 4.5 sigma. This paper presents the catalogs of the FastSound emission lines and galaxies, which will be open to the public in the near future. We also present basic properties of typical FastSound H$\alpha$ emitters, which have H$\alpha$ luminosities of $10^{41.8}$-$10^{43.3}$ erg/s, SFRs of 20--500 $M_\odot$/yr, and stellar masses of $10^{10.0}$--$10^{11.3}$ $M_\odot$. The 3D distribution maps for the four fields of CFHTLS W1--4 are presented, clearly showing large scale clustering of galaxies at the scale of $\sim$ 100--600 comoving Mpc. Based on 1,105 galaxies with detections of multiple emission lines, we estimate that contamination of non-H$\alpha$ lines is about 4% in the single-line emission galaxies, which are mostly [OIII]$\lambda$5007. This contamination fraction is also confirmed by the stacked spectrum of all the FastSound spectra, in which H$\alpha$, [NII]$\lambda \lambda$6548,6583, [SII]$\lambda \lambda$6717, 6731, and [OI]$\lambda \lambda$6300,6364 are seen.

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$ detection of RSD. Our constraint is consistent with the prediction of general relativity $f\sigma_8\sim 0.392$ within the $1-\sigma$ confidence level. When we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $f\sigma_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $f\sigma_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

In the late stages of nuclear burning for massive stars ($M>8~M_{\sun}$), the production of neutrino-antineutrino pairs through various processes becomes the dominant stellar cooling mechanism. As the star evolves, the energy of these neutrinos increases and in the days preceding the supernova a significant fraction of emitted electron anti-neutrinos exceeds the energy threshold for inverse beta decay on free hydrogen. This is the golden channel for liquid scintillator detectors because the coincidence signature allows for significant reductions in background signals. We find that the kiloton-scale liquid scintillator detector KamLAND can detect these pre-supernova neutrinos from a star with a mass of $25~M_{\sun}$ at a distance less than 690~pc with 3$\sigma$ significance before the supernova. This limit is dependent on the neutrino mass ordering and background levels. KamLAND takes data continuously and can provide a supernova alert to the community.

We report a measurement of the neutrino-electron elastic scattering rate of 862 keV 7Be solar neutrinos based on a 165.4 kton-day exposure of KamLAND. The observed rate is 582 +/- 90 (kton-day)^-1, which corresponds to a 862 keV 7Be solar neutrino flux of (3.26 +/- 0.50) x 10^9 cm^-2s^-1, assuming a pure electron flavor flux. Comparing this flux with the standard solar model prediction and further assuming three flavor mixing, a nu_e survival probability of 0.66 +/- 0.14 is determined from the KamLAND data. Utilizing a global three flavor oscillation analysis, we obtain a total 7Be solar neutrino flux of (5.82 +/- 0.98) x 10^9 cm^-2s^-1, which is consistent with the standard solar model predictions.

We present the results from a large near-infrared spectroscopic survey with Subaru/FMOS (\textit{FastSound}) consisting of $\sim$ 4,000 galaxies at $z\sim1.4$ with significant H$\alpha$ detection. We measure the gas-phase metallicity from the [N~{\sc ii}]$\lambda$6583/H$\alpha$ emission line ratio of the composite spectra in various stellar mass and star-formation rate bins. The resulting mass-metallicity relation generally agrees with previous studies obtained in a similar redshift range to that of our sample. No clear dependence of the mass-metallicity relation with star-formation rate is found. Our result at $z\sim1.4$ is roughly in agreement with the fundamental metallicity relation at $z\sim0.1$ with fiber aperture corrected star-formation rate. We detect significant [S~{\sc ii}]$\lambda\lambda$6716,6731 emission lines from the composite spectra. The electron density estimated from the [S~{\sc ii}]$\lambda\lambda$6716,6731 line ratio ranges from 10 -- 500 cm$^{-3}$, which generally agrees with that of local galaxies. On the other hand, the distribution of our sample on [N~{\sc ii}]$\lambda$6583/H$\alpha$ vs. [S~{\sc ii}]$\lambda\lambda$6716,6731/H$\alpha$ is different from that found locally. We estimate the nitrogen-to-oxygen abundance ratio (N/O) from the N2S2 index, and find that the N/O in galaxies at $z\sim1.4$ is significantly higher than the local values at a fixed metallicity and stellar mass. The metallicity at $z\sim1.4$ recalculated with this N/O enhancement taken into account decreases by 0.1 -- 0.2 dex. The resulting metallicity is lower than the local fundamental metallicity relation.

We search for electron anti-neutrinos ($\overline{\nu}_e$) from long and short-duration gamma-ray bursts~(GRBs) using data taken by the KamLAND detector from August 2002 to June 2013. No statistically significant excess over the background level is found. We place the tightest upper limits on $\overline{\nu}_e$ fluence from GRBs below 7 MeV and place first constraints on the relation between $\overline{\nu}_e$ luminosity and effective temperature.

We present a result of a blind search for [CII] 158 $\mu$m emitters at $z\sim 4.5$ using ALMA Cycle~0 archival data. We collected extra-galactic data covering at 330-360 GHz (band~7) from 8 Cycle~0 projects from which initial results have been already published. The total number of fields is 243 and the total on-source exposure time is 19.2 hours. We searched for line emitters in continuum-subtracted data cubes with spectral resolutions of 50, 100, 300 and 500 km/s. We could not detect any new line emitters above a 6-$\sigma$ significance level. This result provides upper limits to the [CII] luminosity function at $z\sim 4.5$ over $L_{\rm [CII]} \sim 10^8 - 10^{10} L_{\odot}$ or star formation rate, SFR $\sim$ 10-1000 M$_{^\odot}$/yr. These limits are at least 2 orders of magnitude larger than the [CII] luminosity functions expected from the $z \sim 4$ UV luminosity function or from numerical simulation. However, this study demonstrates that we would be able to better constrain the [CII] luminosity function and to investigate possible contributions from dusty galaxies to the cosmic star-formation rate density by collecting Cycle~1+2 archival data as the ALMA Patchy Deep Survey.

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

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 present cosmological hydrodynamic simulations performed to study evolution of galaxy population. The simulations follow timed release of mass, energy, and metals by stellar evolution and employ phenomenological treatments of supernova feedback, pre-supernova feedback modeled as feedback by radiation pressure from massive stars, and quenching of gas cooling in large halos. We construct the fiducial model so that it reproduces the observationally estimated galaxy stellar mass functions and the relationships between the galaxy stellar mass and the host halo mass from $z = 4$ to 0. We find that the fiducial model constructed this way naturally explains the cosmic star formation history, the galaxy downsizing, and the star formation rate and metallicity of the star-forming galaxies. The simulations without the quenching of the gas cooling in large halos overproduce massive galaxies at $z < 2$ and fail to reproduce galaxy downsizing. The simulations that do not employ the radiation pressure feedback from young stars predict too strong redshift evolution of the mass-metallicity relation. Furthermore, the slope of the relation becomes too steep at low redshift without the radiation pressure feedback. The metallicity dependence in the radiation pressure feedback is a key to explain the observed mass-metallicity relation. These facts indicate that these two processes in addition to supernova feedback are essential for galaxy evolution. Our simple phenomenological model is suitable to construct a mock galaxy sample to study physical properties of observed galaxy populations.

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.

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.

We use large cosmological Smoothed-Particle-Hydrodynamics simulations to study the formation and evolution of sub-millimetre galaxies (SMGs). In our previous work, we studied the statistical properties of ultra-violet selected star-forming galaxies at high redshifts. We populate the same cosmological simulations with SMGs by calculating the reprocess of stellar light by dust grains into far-infrared to millimetre wavebands in a self-consistent manner. We generate light-cone outputs to compare directly the statistical properties of the simulated SMGs with available observations. Our model reproduces the submm source number counts and the clustering amplitude. We show that bright SMGs with flux $S > 1$ mJy reside in halos with mass of $\sim 10^{13} M_{\odot}$ and have stellar masses greater than $10^{11}\sim \rm M_{\odot}$. The angular cross-correlation between the SMGs and Lyman-$\alpha$ emitters is significantly weaker than that between the SMGs and Lyman-break galaxies. The cross-correlation is also weaker than the auto-correlation of the SMGs. The redshift distribution of the SMGs shows a broad peak at $z \sim 2$, where Bright SMGs contribute significantly to the global cosmic star formation rate density. Our model predicts that there are hundreds of SMGs with $S > 0.1$ mJy at $z > 5$ per 1 square degree field. Such SMGs can be detected by ALMA.

We report a measurement of the neutrino-electron elastic scattering rate from 8B solar neutrinos based on a 123 kton-day exposure of KamLAND. The background-subtracted electron recoil rate, above a 5.5 MeV analysis threshold is 1.49+/-0.14(stat)+/-0.17(syst) events per kton-day. Interpreted as due to a pure electron flavor flux with a 8B neutrino spectrum, this corresponds to a spectrum integrated flux of 2.77+/-0.26(stat)+/-0.32(syst) x 10^6 cm^-2s^-1. The analysis threshold is driven by 208Tl present in the liquid scintillator, and the main source of systematic uncertainty is due to background from cosmogenic 11Be. The measured rate is consistent with existing measurements and with Standard Solar Model predictions which include matter enhanced neutrino oscillation.

We use very large cosmological Smoothed-Particle-Hydrodynamics simulations to study the properties of high redshift Lyman-alpha emitters (LAEs). We identify star-forming galaxies at z=3.1 in a cosmological volume of 100 h^-1Mpc on a side. We develop a phenomenological model of absorption, scattering and escape of Lyman-alpha photons on the assumption that the clumpiness of the inter-stellar medium in a galaxy is correlated with the larger scale substructure richness. The radiative transfer effect proposed by Neufeld (1991) allows a large fraction of Lyman-alpha photons to escape from a clumpy galaxy even if it contains a substantial amount of dust. Our model reproduces, for the first time, all of the following observed properties of LAEs at z=3.1: the angular correlation function, ultra-violet and Lyman-alpha luminosity functions, and the equivalent width distribution. A simple model that takes only dust absorption into account fails in matching the observational data, suggesting that the kind of effect we consider is needed. Our model also predicts a bimodal age distribution for LAEs. There are old, massive and dusty LAEs, similar to recently found high redshift LAEs. The large LAEs have escape fractions of Lyman-alpha photons of f_esc ~ 0.05-0.1.

From multi-wavelength observations of LAEs,we know that while many LAEs appear to be young and less massive,a noticeable fraction of LAEs possess much older populations of stars and larger stellar mass.How these two classes of LAEs are concordant with the hierarchical galaxy formation scenario has not been understood clearly so far.In this paper,we model LAEs by three-dimensional cosmological simulations of dark halo merger in a CDM universe.As a result,it is shown that the age of simulated LAEs can spread over a wide range from 2*10^6yr to 9*10^8yr.Also,we find that there are two types of LAEs, in one of which the young half-mass age is comparable to the mean age of stellar component,and in the other of which the young half-mass age is appreciably shorter than the mean age.We define the former as Type 1 LAEs and the latter as Type 2 LAEs.A Type 1 corresponds to early starburst in a young galaxy,whereas a Type 2 does to delayed starburst in an evolved galaxy,as a consequence of delayed accretion of a subhalo onto a larger parent halo.Thus,the same halo can experience a Type 2 LAE-phase as well as a Type 1 LAE-phase in the merger history.Type 1s are expected to be younger than 1.5*10^8yr,less dusty,and less massive with stellar mass M*<10^8 Msun,while Type 2s are older than 1.5*10^8yr,even dustier,and as massive as M*~10^8-10^10Msun.The fraction of Type 2s in all LAEs is a function of redshift.Type 2s discriminated clearly from Type 1s in two color diagram of z'-H vs J-K.We find that the brightness distribution of LyA in Type 2s is more extended than the main stellar component,in contrast to Type 1s.This is not only because delayed starbursts tend to occur in the outskirts of a parent galaxy,but also because LyA photons are effectively absorbed by dust in an evolved galaxy.Hence,the extent of LyA emission may be an additional measure to distinguish Type 2s from Type 1s

Recently discovered large-scale structure of Ly\alpha Emitters (LAEs) raises a novel challenge to the cold dark matter (CDM) cosmology. The structure is extended over more than 50 Mpc at redshift z=3.1, and exhibits a considerably weak angular correlation. Such properties of LAE distributions appear to be incompatible with the standard biased galaxy formation scenario in the CDM cosmology. In this paper, by considering the possibility that LAEs are short-lived events, we attempt to build up the picture of LAEs concordant with the CDM cosmology. We find that if the lifetime of LAEs is as short as (6.7 \pm 0.6) \times 10^7 yr, the distributions of simulated galaxies successfully match the extension and morphology of large-scale structure of LAEs at z=3.1, and also the weak angular correlation function. This result implies that LAEs at z=3.1 do not necessarily reside in high density peaks, but tends to be located in less dense regions, in a different way from the expectation by the standard biased galaxy formation scenario. In addition, we make a prediction for the angular correlation function of LAEs at redshifts higher than 3. It is found that the prediction deviates from that by the standard biased galaxy formation scenario even at redshifts 4 < z < 6.