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

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.

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

Observations of high-redshift quasars imply the presence of supermassive black holes already at z~ 7.5. An appealing and promising pathway to their formation is the direct collapse scenario of a primordial gas in atomic-cooling haloes at z ~ 10 - 20, when the $\rm H_2$ formation is inhibited by a strong background radiation field, whose intensity exceeds a critical value, $J_{\rm crit}$. To estimate $J_{\rm crit}$, typically, studies have assumed idealized spectra, with a fixed ratio of $\rm H_{2}$ photo-dissociation rate $k_{\rm H_2}$ to the $\rm H^-$ photo-detachment rate $k_{\rm H^-}$. This assumption, however, could be too narrow in scope as the nature of the background radiation field is not known precisely. In this work we argue that the critical condition for suppressing the $\rm H_2$ cooling in the collapsing gas could be described in a more general way by a combination of $k_{\rm H_2}$ and $k_{\rm H^-}$ parameters. By performing a series of cosmological zoom-in simulations with an encompassing set of $k_{\rm H_2}$ and $k_{\rm H^-}$, we examine the gas flow by following evolution of basic parameters of the accretion flow. We test under what conditions the gas evolution is dominated by $\rm H_{2}$ and/or atomic cooling. We confirm the existence of a critical curve in the $k_{\rm H_2}-k_{\rm H^-}$ plane, and provide an analytical fit to it. This curve depends on the conditions in the direct collapse, and reveals domains where the atomic cooling dominates over the molecular cooling. Furthermore, we have considered the effect of $\rm H_{2}$ self-shielding on the critical curve, by adopting three methods for the effective column density approximation in $\rm H_{2}$. We find that the estimate of the characteristic length-scale for shielding can be improved by using $\lambda_{\rm Jeans25}$, which is 0.25 times that of the local Jeans length.

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

We compute the evolution of interstellar dust in a hydrodynamic simulation of an isolated disc galaxy. We newly implement the evolution of full grain size distribution by sampling 32 grid points on the axis of the grain radius. We solve it consistently with the chemical enrichment and hydrodynamic evolution of the galaxy. This enables us to theoretically investigate spatially resolved evolution of grain size distribution in a galaxy. The grain size distribution evolves from a large-grain-dominated ($\gtrsim 0.1~\mu$m) phase to a small-grain production phase, eventually converging to a power-law-like grain size distribution similar to the so-called MRN distribution. We find that the small-grain abundance is higher in the dense ISM in the early epoch ($t\lesssim 1$ Gyr) because of efficient dust growth by accretion, while coagulation makes the small-grain abundance less enhanced in the dense ISM later. This leads to steeper extinction curves in the dense ISM than in the diffuse ISM in the early phase, while they show the opposite trend later. The radial trend is also described by faster evolution in the inner part. We also confirm that the simulation reproduces the observed trend in the relation between dust-to-gas ratio and metallicity, and in the radial gradients of dust-to-gas ratio and dust-to-metal ratio. Since the above change in the grain size distribution occurs at $t\sim 1$ Gyr, the age and density dependence of grain size distribution has a significant impact on the extinction curves even at high redshift.

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

We search for host galaxy candidates of nearby fast radio bursts (FRBs), FRB 180729.J1316+55, FRB 171020, FRB 171213, FRB 180810.J1159+83, and FRB 180814.J0422+73 (the second repeating FRB). We compare the absolute magnitudes and the expected host dispersion measure $\rm DM_{host}$ of these candidates with that of the first repeating FRB, FRB 121102, as well as those of long gamma ray bursts (LGRBs) and superluminous supernovae (SLSNe), the proposed progenitor systems of FRB 121102. We find that while the FRB 121102 host is consistent with those of LGRBs and SLSNe, the nearby FRB host candidates, at least for FRB 180729.J1316+55, FRB 171020, and FRB180814.J0422+73, either have a smaller $\rm DM_{host}$ or are fainter than FRB121102 host, as well as the hosts of LGRBs and SLSNe. In order to avoid the uncertainty in estimating $\rm DM_{host}$ due to the line-of-sight effect, we propose a galaxy-group-based method to estimate the electron density in the inter-galactic regions, and hence, $\rm DM_{IGM}$. The result strengthens our conclusion. We conclude that the host galaxy of FRB 121102 is atypical, and LGRBs and SLSNe are likely not the progenitor systems of at least most nearby FRB sources. {The recently reported two FRB hosts differ from the host of FRB 121102 and also the host candidates suggested in this paper. This is consistent with the conclusion of our paper and suggests that the FRB hosts are very diverse. }

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

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

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

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.

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.

Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires us to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ~10^{-6} pc and rapidly expands outward. (2) A central core forms, with a mass of ~1 Mo, supported by thermal gas pressure gradients and rotation. (3) Growing thermal gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ~5 x 10^{37}-5 x 10^{38} erg/s, of order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ~10^{-4}-10^{-3} pc, and ~0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, mainly because the radiation in the FLD is capable of escaping due to anisotropy in the optical depth and associated gradients.

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

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

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

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

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.

We perform three dimensional smoothed particle hydrodynamics (SPH) simulations of gas accretion onto the seeds of binary stars to investigate their short-term evolution. Our simulation setup is more realistic compared to the previous works by taking into account of dynamically evolving envelope with non-uniform distribution of gas density and angular momentum of accreting flow. Our initial condition includes a seed binary and a surrounding gas envelope, modelling the phase of core collapse of gas cloud when the fragmentation has already occurred. We assume that the seed binary has no eccentricity and no growth by gas accretion. The envelope is assumed to be an isothermal gas with no self-gravity. We run multiple simulations with different values of initial mass ratio $q_0$ (the ratio of secondary over primary mass) and gas temperature, and find a critical value of $q_{\rm c} = 0.25$ which distinguishes the later evolution of mass ratio $q$ as a function of time. If $q_0$ >~ $q_{\rm c}$, the secondary seed grows faster and $q$ increases monotonically towards unity. If $q_0$ <~ $q_{\rm c}$, on the other hand, the primary seed grows faster and $q$ is lower than $q_0$ at the end of the simulation. Based on our numerical results, we analytically calculate the long-term evolution of the seed binary including the growth of binary by gas accretion. We find that the seed binary with $q_0$ >~ $q_{\rm c}$ evolves towards an equal-mass binary star, and that with $q_0$ <~ $q_{\rm c}$ evolves to a binary with an extreme value of $q$. Binary separation is a monotonically increasing function of time for any $q_0$, suggesting that the binary growth by accretion does not lead to the formation of close binaries.

We provide detailed comparison between the AMR code Enzo-2.4 and the SPH/N- body code GADGET-3 in the context of isolated or cosmological direct baryonic collapse within dark matter (DM) halos to form supermassive black holes. Gas flow is examined by following evolution of basic parameters of accretion flows. Both codes show an overall agreement in the general features of the collapse, however, many subtle differences exist. For isolated models, the codes increase their spatial and mass resolutions at different pace, which leads to substantially earlier collapse in SPH than in AMR cases due to higher gravitational resolution in GADGET-3. In cosmological runs, the AMR develops a slightly higher baryonic resolution than SPH during halo growth via cold accretion permeated by mergers. Still, both codes agree in the buildup of DM and baryonic structures. However, with the onset of collapse, this difference in mass and spatial resolution is amplified, so evolution of SPH models begins to lag behind. Such a delay can have effect on formation/destruction rate of H2 due to UV background, and on basic properties of host halos. Finally, isolated non-cosmological models in spinning halos, with spin parameter {\lambda} ~ 0.01 - 0.07, show delayed collapse for greater {\lambda}, but pace of this increase is faster for AMR. Within our simulation setup, GADGET-3 requires significantly larger computational resources than Enzo- 2.4 during collapse, and needs similar resources, during the pre-collapse, cosmological structure formation phase. Yet it benefits from substantially higher gravitational force and hydrodynamic resolutions, except at the end of collapse.

Using the cosmological smoothed particle hydrodynamical code GADGET-3 we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of $3\leq z \leq7$ one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within $\sim 0.20$ dex. Over the stellar mass range we probe ($10^{10.39}\leq M_s/M_\odot \leq 10^{10.75}$ at $z =$ 3 and $10^{8.48}\leq M_s/M_\odot \leq 10^{9.55}$ at $z =$ 7) one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find the cumulative number density evolution exhibits a trend towards higher values which can be quantified by simple linear formulations going as $-0.10\Delta z$ for descendants and $0.12\Delta z$ for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of $\sim2$. This result is in excellent agreement, within $0.10$ dex, with abundance matching results over the same redshift range. However, we find that our more realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion, mergers, and star formation of individual galaxies will lead to a larger scatter in other physical properties such as metallicity and star-formation rate.

We use cosmological adaptive mesh refinement (AMR) code Enzo zoom-in simulations to study the long term evolution of the collapsing gas within dark matter (DM) halos at high redshifts. This direct collapse process is a leading candidate for rapid formation of supermassive black hole (SMBH) seeds at high z. To circumvent the Courant condition at small radii, we have used the sink particle method, and focus on the evolution on scales ~0.01-10 pc. The collapse proceeds in two stages, with the secondary runaway happening within the central 10 pc, and with no detected fragmentation. The sink particles form when the collapsing gas requires additional refinement of the grid size at the highest refinement level. Their mass never exceeds ~10^3 Mo, with the sole exception of the central seed which grows dramatically to ~ 2 x 10^6 Mo in ~2 Myr, confirming the feasibility of this path to the SMBH. The time variability of angular momentum axis in the accreted gas results in the formation of two misaligned disks --- a small inner disk, and a more massive, outer disk which is inclined by ~45^o to the inner disk. The self-gravity of these disks is heavily diluted --- both disks lie within the Roche limit of the central seed. While the inner disk is geometrically thin and weakly asymmetric, the outer disk flares due to turbulent motions as a result of the massive inflow along a pair of penetrating filaments. The geometry of inflow via filaments determines the dominant and secondary Fourier modes in this disk --- these modes have a non-self-gravitational origin. We do not confirm that m=1 is a principal mode that drives the inflow in the presence of a central massive object. While the positions of the disks depend on the scale chosen to break the self-similar collapse, the overall configuration appears to be generic, and is expected to form when the central seed becomes sufficiently massive.

The Bullet Cluster has provided some of the best evidence for the $\Lambda$ cold dark matter ($\Lambda\rm{CDM}$) model via direct empirical proof of the existence of collisionless dark matter, while posing a serious challenge owing to the unusually high inferred pairwise velocities of its progenitor clusters. Here we investigate the probability of finding such a high-velocity pair in large-volume N-body simulations, particularly focusing on differences between halo finding algorithms. We find that algorithms that do not account for the kinematics of infalling groups yield vastly different statistics and probabilities. When employing the ROCKSTAR halo finder that considers particle velocities, we find numerous Bullet-like pair candidates that closely match not only the high pairwise velocity, but also the mass, mass ratio, separation distance, and collision angle of the initial conditions that have been shown to produce the Bullet Cluster in non-cosmological hydrodynamic simulations. The probability of finding a high pairwise velocity pair among haloes with $M_{\rm halo}\geq10^{14} M_{\odot}$ is $4.6\times 10^{-4}$ using ROCKSTAR, while it is $\approx 34\times$ lower using a friends-of-friends (FOF) based approach as in previous studies. This is because the typical spatial extent of Bullet progenitors is such that FOF tends to group them into a single halo despite clearly distinct kinematics. Further requiring an appropriately high average mass among the two progenitors, we find the comoving number density of potential Bullet-like candidates to be on the order of $\approx10^{-10} \,{\rm Mpc}^{-3}$. Our findings suggest that $\Lambda\rm{CDM}$ straightforwardly produces massive, high relative velocity halo pairs analogous to Bullet Cluster progenitors, and hence the Bullet Cluster does not present a challenge to the $\Lambda\rm{CDM}$ model.

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

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.

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

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

It has been shown observationally that star formation (SF) correlates tightly with the presence of molecular hydrogen (H$_2$). Therefore it would be important to investigate its implication on galaxy formation in a cosmological context. In the present work, we track the H$_2$ mass fraction within our cosmological smoothed particle hydrodynamics (SPH) code GADGET-3 using an equilibrium analytic model by Krumholz et al. This model allows us to regulate the star formation in our simulation by the local abundance of H$_2$ rather than the total cold gas density, and naturally introduce the dependence of star formation on metallicity. We investigate implications of the equilibrium H$_2$-based SF model on galaxy population properties, such as the stellar-to-halo mass ratio (SHMR), baryon fraction, cosmic star formation rate density (SFRD), galaxy specific SFR, galaxy stellar mass functions (GSMF), and Kennicutt-Schmidt (KS) relationship. The advantage of our work over the previous ones is having a large sample of simulated galaxies in a cosmological volume from high-redshift to $z=0$. We find that low-mass halos with M$_{DM}<10^{10.5}$ M$_\odot$ are less efficient in producing stars in the H$_2$-based SF model at $z\geq6$, which brings the simulations to a better agreement with observational estimates of SHMR and GSMF at the low-mass end. This is particularly evident by a reduction in the number of low-mass galaxies at M$_\star\leq10^{8}$ M$_\odot$ in the GSMF. The overall SFRD is also reduced at high-$z$ in the H$_2$ run, which results in slightly higher SFRD at low-redshift due to more abundant gas available for star formation at later times. This new H$_2$ model is able to reproduce the empirical KS relationship at $z=0$ naturally without the need for setting its normalization by hand, and overall it seems to have more advantages than the previous pressure-based SF model.

Estimating the mass of a supermassive black hole (SMBH) in an active galactic nucleus (AGN) usually relies on the assumption that the broad line region (BLR) is virialized. However, this assumption seems invalid in BLR models that consists of an accretion disk and its wind. The disk is likely Keplerian and therefore virialized. However, the wind material must, beyond a certain point, be dominated by an outward force that is stronger than gravity. Here, we analyze hydrodynamic simulations of four different disk winds: an isothermal wind, a thermal wind from an X-ray heated disk, and two line-driven winds, one with and the other without X-ray heating and cooling. For each model, we check whether gravity governs the flow properties, by computing and analyzing the volume-integrated quantities that appear in the virial theorem: internal, kinetic, and gravitational energies, We find that in the first two models, the winds are non-virialized whereas the two line-driven disk winds are virialized up to a relatively large distance. The line-driven winds are virialized because they accelerate slowly so that the rotational velocity is dominant and the wind base is very dense. For the two virialized winds, the so-called projected virial factor scales with inclination angle as $1/ \sin^2{i}$. Finally, we demonstrate that an outflow from a Keplerian disk becomes unvirialized more slowly when it conserves the gas specific angular momentum -- as in the models considered here, than when it conserves the angular velocity -- as in the so-called magneto-centrifugal winds.

We present the results of a numerical study examining the effect of H_2-based star formation (SF) model on the rest-frame UV luminosity function and star formation rate function (SFRF) of z>=6 galaxies, and the implications for reionization. Using cosmological hydrodynamical simulations outfitted with an H_2-SF model, we find good agreement with our previous results (non-H_2 SF model) and observations at Muv<=-18. However at Muv>-18, we find that the LF deviates from both our previous work and current observational extrapolations, producing significantly fewer low-luminosity galaxies and exhibiting additional turnover at the faint end. We constrain the redshift evolution of this turnover point using a modified Schechter function that includes additional terms to quantify the turnover magnitude (Muv^t) and subsequent slope (Beta). We find that Muv^t evolves from Muv^t=-17.33 (at z=8) to -15.38 (z=6), while Beta becomes shallower by \Delta Beta=0.22 during the same epoch. This occurs in an Muv range which will be observable by JWST. By integrating the SFRF, we determine that even though H_2-SF model significantly reduces the number density of low luminosity galaxies at Muv>-18, it does not suppress the total SFR density enough to affect the capability of SF to maintain reionization.

The passage of muons through matter is mostly affected by their Coulomb interactions with electrons and nuclei. The muon interactions with electrons lead to continuous energy loss and stopping of muons, while their scattering off nuclei lead to angular 'diffusion'. By measuring both the number of stopped muons and angular changes in muon trajectories we can estimate density and identify materials. Here we demonstrate the material identification using data taken at Los Alamos with the Mini Muon Tracker.

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

We examine the duty cycle and the history of star formation (SFH) for high-redshift galaxies at z>=6 using cosmological hydrodynamic simulations. We find that, even though individual galaxies have bursty SFH, the averaged SFH between z~15 to z=6 can be characterized well by either an exponentially increasing functional form with characteristic time-scales of 70 Myr to 200 Myr for galaxies with stellar masses Ms~10^6 Msun to >10^10 Msun respectively, or by a simple power-law form which exhibits a similar mass dependent time-scales. Using the SFH of individual galaxies, we measure the duty cycle of star formation (DC_SFH); i.e., the fraction of time a galaxy of a particular mass spends above a star formation rate (SFR) threshold which would make it observable to the Hubble Space Telescope (HST) during a given epoch. We also examine the fraction of galaxies at a given redshift that are brighter than a rest-frame UV magnitude (Muv ~ -18), which is sufficient enough to make them observable (DC_Muv). We find that both DC_SFH and DC_Muv make a sharp transition from zero (for galaxies with Ms <= 10^7 Msun) to unity (for Ms > 10^9 Msun). The measured duty cycle is also manifested in the intrinsic scatter in the Ms-SFR relationship (~ 1 dex) and Ms-Muv relationship (\Delta Muv ~ +-1 mag). We provide analytic fits to the DC as a function of Ms using a sigmoid function, which can be used to correct for catalogue incompleteness. We consider the effects of duty cycle to the observational estimate of galaxy stellar mass functions (GSMF) and the star formation rate density (SFRD), and find that it results in a much shallower low-mass end slopes of the GSMF and a reduction of >~ 70% of our intrinsic SFRD, making our simulation results more compatible with observational estimates.

(Abridged) We investigate non-spherical behavior of gas accreting onto a central supermassive black hole performing simulations using the SPH code GADGET-3 including radiative cooling and heating by the central X-ray source. As found in earlier 1D studies, our 3D simulations show that the accretion mode depends on the X-ray luminosity (L_X) for a fixed density at infinity and accretion efficiency. In the low L_X limit, gas accretes in a stable, spherically symmetric fashion. In the high L_X limit, the inner gas is significantly heated up and expands, reducing the central mass inflow rate. The expanding gas can turn into a strong enough outflow capable of expelling most of the gas at larger radii. For some intermediate L_X, the accretion flow becomes unstable developing prominent non-spherical features, the key reason for which is thermal instability (TI) as shown by our analyses. Small perturbations of the initially spherically symmetric accretion flow that is heated by the intermediate L_X quickly grow to form cold and dense clumps surrounded by overheated low density regions. The cold clumps continue their inward motion forming filamentary structures; while the hot infalling gas slows down because of buoyancy and can even start outflowing through the channels in between the filaments. We found that the ratio between the mass inflow rates of the cold and hot gas is a dynamical quantity depending on several factors: time, spatial location, and L_X; and ranges between 0 and 4. We briefly discuss astrophysical implications of such TI-driven fragmentation of accreting gas on the formation of clouds in narrow and broad line regions of AGN, the formation of stars, and the observed variability of the AGN luminiosity.

We study the gas accretion onto a supermassive black hole (SMBH) using the 3D SPH code GADGET-3 on scales of 0.1-200 pc. First we test our code with spherically symmetric, adiabatic Bondi accretion problem. We find that our simulation can reproduce the expected Bondi accretion flow very well for a limited amount of time until the effect of outer boundary starts to be visible. We also find artificial heating of gas near the inner accretion boundary due to the artificial viscosity of SPH. Second, we implement radiative cooling and heating due to X-rays, and examine the impact of thermal feedback by the central X-ray source. The accretion flow roughly follows the Bondi solution for low central X-ray luminosities, however, the flow starts to exhibit non-spherical fragmentation due to thermal instability for a certain range of central L_X, and a strong overall outflow develops for greater L_X. The cold gas develops filamentary structures that fall into the central SMBH, whereas the hot gas tries to escape through the channels in-between the cold filaments. Such fragmentation of accreting gas can assist in the formation of clouds around AGN, induce star-formation, and contribute to the observed variability of narrow-line regions.

We present the results of a numerical study comparing photometric and physical properties of simulated z=6-9 galaxies to the observations taken by the WFC3 instrument aboard the Hubble Space Telescope. Using cosmological hydrodynamical simulations we find good agreement with observations in color-color space at all studied redshifts. We also find good agreement between observations and our Schechter luminosity function fit in the observable range, Muv<= -18, provided that a moderate dust extinction effect exists for massive galaxies. However beyond what currently can be observed, simulations predict a very large number of low-mass galaxies and evolving steep faint-end slopes from alpha_L = -2.15 at z=6 to alpha_L = -2.64 at z=9, with a dependence of |alpha_L| \propto (1+z)^0.59. During the same epoch, the normalization phi* increases and the characteristic magnitude Muv* becomes moderately brighter with decreasing redshift. We find similar trends for galaxy stellar mass function with evolving low-mass end slope from alpha_M = - 2.26 at z=6 to alpha_M = -2.87 at z=9, with a dependence of |alpha_M| \propto (1+z)^0.65. Together with our recent result on the high escape fraction of ionizing photons for low-mass galaxies, our results suggest that the low-mass galaxies are important contributor of ionizing photons for the reionisation of the Universe at z>=6.

The existence of 1E0657-56 poses a challenge to the concordance Lambda cold dark matter model. Here we investigate the velocity distribution of dark matter halo pairs in large N-body simulations with differing box sizes (250Mpc/h-2 Gpc/h) and resolutions. We examine statistics such as the halo masses, pairwise halo velocities (v_{12}), and pair separation distances. We then compare our results to the initial conditions (ICs) required to reproduce the observational properties of 1E0657-56 in non-cosmological hydrodynamical simulations. We find that the high velocity tail of the v_{12} distribution extends to greater velocities as we increase the box size. We also find that the number of high-v_{12} pairs increases as we increase the particle count and resolution with a fixed box size, however, this increase is mostly due to lower mass halos which do not match the observed masses of 1E0657-56. We find that the redshift evolution is not very strong for the v_{12} distribution function between z=0.0 and z~0.5. We identify some pairs whose v_{12} resemble the required ICs, however, even the best candidates have either wrong halo mass ratios, or too large separations. Our simulations suggest that it is very difficult to produce such ICs at z=0.0, 0.296,& 0.489 in comoving volumes as large as (2Gpc/h)^3. Based on the extrapolation of our cumulative v_{12} function, we find that one needs a simulation with a comoving box size of (4.48Gpc/h)^3 and 2240^3 DM particles in order to produce at least one pair of halos that resembles the required v_{12} and observed masses of 1E0657-56. We find that the probability of finding a halo pair with v_{12}>=3000km/s and masses >=10^{14}Msun to be 2.76x10^{-8} at z=0.489. We conclude that either 1E0657-56 is incompatible with the concordance LCDM universe, or the ICs suggested by the non-cosmological simulations must be revised to give a lower value of v_{12}.

There are mainly two different approaches to measure the cosmic star formation history: direct star formation rate density (SFRD) and stellar mass density rhostar as functions of redshift. Compilations of current observations seem to show a disparity in the two quantities, in the sense that the integral of SFRD is higher than the observed rhostar (after considering gas recycling). Using cosmological smoothed particle hydrodynamics simulations based on the concordance Lambda cold dark matter model, we show that the two quantities become more consistent with each other when we consider the observed galaxy mass limit. The comparison between simulations and (dust corrected) observed cosmic SFRD shows a good agreement, while the observed rhostar is significantly lower than the simulation results. This can be reconciled if the current high-$z$ galaxy surveys are missing faint low-mass galaxies due to their flux limit. Our simulated GSMFs have steep low-mass end slopes of alpha < -2 at z>3, and when these numerous low-mass galaxies are included, the total rhostar matches with the integral of SFRD.

(Abridged) We test how accurately the smoothed particle hydrodynamics (SPH) numerical technique can follow spherically-symmetric Bondi accretion. Using the 3D SPH code GADGET-3, we perform simulations of gas accretion onto a central supermassive black hole (SMBH) of mass 10^8 M_sun within the radial range of 0.1 - 200 pc. We carry out simulations without and with radiative heating by a central X-ray corona and radiative cooling. For an adiabatic case, the radial profiles of hydrodynamical properties match the Bondi solution, except near the inner and outer radius of the computational domain. We find that adiabatic Bondi accretion can be reproduced for durations of a few dynamical times at the Bondi radius, and for longer times if the outer radius is increased. With radiative heating and cooling included, the spherically accreting gas takes a longer time to reach a steady-state than the adiabatic Bondi accretion runs, and in some cases does not reach a steady-state even within several hundred dynamical times. We find that artificial viscosity in the GADGET code causes excessive heating near the inner radius, making the thermal properties of the gas inconsistent with a physical solution. This overheating occurs typically only in the supersonic part of the flow, so that it does not affect the mass accretion rate. We see that increasing the X-ray luminosity produces a lower central mass inflow rate, implying that feedback due to radiative heating is operational in our simulations. With a sufficiently high X-ray luminosity, the inflowing gas is radiatively heated up, and an outflow develops. We conclude that the SPH simulations can capture the gas dynamics needed to study radiative feedback provided artificial viscosity alters only highly supersonic part of the inflow.

Prompted by various analyses of long (Type II) GRB rates and their relationship to the cosmic star-formation history, metallicity and luminosity function evolution, we systematically analyze these effects with a Monte Carlo code. We test various cosmic star-formation history models including analytical and empirical models as well as those derived from cosmological simulations. We also explore expressions for metallicity enhancement of the GRB rate with redshift, as presented in the literature, and discuss improvements to these analytic expressions from the point of view of galactic evolution. These are also compared to cosmological simulations on metal enrichment. Additionally we explore possible evolutionary effects of the GRB rate and luminosity function with redshift. The simulated results are tested with the observed Swift sample including the L, z, and peak flux (log N-log P) distributions. The observational data imply that an increase in the GRB rate is necessary to account for the observations at high redshift, although the form of this enhancement is unclear. A rate increase due to lower metallicity at higher redshift may not be the singular cause and is subject to a variety of uncertainties. Alternatively, evolution of the GRB luminosity function break with redshift shows promise as a possible alternative.

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

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.

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

We study the properties of Lyman-alpha emitters (LAEs) and Lyman-break galaxies (LBGs) at z=3-6 using cosmological SPH simulations. We investigate two simple scenarios for explaining the observed Ly-a and rest-frame UV luminosity functions (LFs) of LAEs: (i) the "escape fraction" scenario, in which the "effective" escape fraction (including the IGM attenuation) of Ly-a photons is f_Lya ~0.1 (0.15) at z=3 (6), and (ii) the "stochastic" scenario, in which the fraction of LAEs that are turned on at z=3 (6) is \Cstoc ~0.07 (0.2) after correcting for the IGM attenuation. Our comparisons with a number of different observations suggest that the stochastic scenario is preferred over the escape fraction scenario. We find that the mean values of stellar mass, metallicity and black hole mass hosted by LAEs are all smaller in the stochastic scenario than in the escape fraction scenario. In our simulations, the galaxy stellar mass function evolves rapidly, as expected in hierarchical structure formation. However, its evolution is largely compensated by a beginning decline in the specific star formation rate, resulting in little evolution of the rest-frame UV LF from z=6 to 3. The rest-frame UV LF of both LAEs and LBGs at z=3 & 6 can be described well by the stochastic scenario provided the extinction is moderate, E(B-V) ~0.15, for both populations, although our simulation might be overpredicting the number of bright LBGs at z=6. We also discuss the correlation function and bias of LAEs. The Ly-a LFs at z=6 in a field-of-view of 0.2 deg^2 show a significantly larger scatter owing to cosmic variance relative to that in a 1 deg^2 field, and the scatter seen in the current observational estimates of the Ly-a LF can be accounted for by cosmic variance.

We calculate the cross-correlation function (CCF) between damped Ly-a systems (DLAs) and Lyman break galaxies (LBGs) using cosmological hydrodynamic simulations at z=3. We compute the CCF with two different methods. First, we assume that there is one DLA in each dark matter halo if its DLA cross section is non-zero. In our second approach we weight the pair-count by the DLA cross section of each halo, yielding a cross-section-weighted CCF. We also compute the angular CCF for direct comparison with observations. Finally, we calculate the auto-correlation functions of LBGs and DLAs, and their bias against the dark matter distribution. For these different approaches, we consistently find that there is good agreement between our simulations and observational measurements by Cooke et al. and Adelberger et al. Our results thus confirm that the spatial distribution of LBGs and DLAs can be well described within the framework of the concordance Lambda CDM model. We find that the correlation strengths of LBGs and DLAs are consistent with the actual observations, and in the case of LBGs it is higher than would be predicted by low-mass galaxy merger models.

We develop a new ``Multicomponent and Variable Velocity'' (MVV) galactic outflow model for cosmological smoothed particle hydrodynamic (SPH) simulations. The MVV wind model reflects the fact that the wind material can arise from different phases in the interstellar medium (ISM), and the mass-loading factor in the MVV model is a function of galaxy stellar mass. We find that the simulation with the MVV outflow has the following characteristics: (i) the intergalactic medium (IGM) is hardly heated up, and the mean IGM temperature is almost the same as in the no-wind run; (ii) it has lower cosmic star formation rates (SFRs) compared to the no-wind run, but higher SFRs than the constant velocity wind run; (iii) it roughly agrees with the observed IGM metallicity, and roughly follows the observed evolution of Omega(Civ); (iv) the lower mass galaxies have larger mass-loading factors, and the low-mass end of galaxy stellar mass function is flatter than in the previous simulations. Therefore, the MVV outflow model mildly alleviates the problem of too steep galaxy stellar mass function seen in the previous SPH simulations. In summary, the new MVV outflow model shows reasonable agreement with observations, and gives better results than the constant velocity wind model.

We compute the infrared (IR) emission from high-redshift galaxies in cosmological smoothed particle hydrodynamics simulations by coupling the output of the simulation with the population synthesis code `GRASIL' by Silva et al. Based on the stellar mass, metallicity and formation time of each star particle, we estimate the full spectral energy distribution of each star particle from ultraviolet to IR, and compute the luminosity function of simulated galaxies in the Spitzer broadband filters for direct comparison with the available Spitzer observations.

We investigate the effects of the change of cosmological parameters and star formation (SF) models on the cosmic SF history using cosmological smoothed particle hydrodynamics (SPH) simulations based on the cold dark matter (CDM) model. We vary the cosmological parameters within 1-sigma error from the WMAP best-fit parameters, and find that such changes in cosmological parameters mostly affect the amplitude of the cosmic SF history. At high redshift (hereafter high-z), the star formation rate (SFR) is sensitive to the cosmological parameters that control the small-scale power of the primordial power spectrum, while the cosmic matter content becomes important at lower redshifts. We also test two new SF models: 1) the `Pressure' model based on the work by Schaye & Dalla Vecchia (2008), and 2) the `Blitz' model that takes the effect of molecular hydrogen formation into account, based on the work by Blitz & Rosolowsky (2006). Compared to the previous conventional SF model, the Pressure model reduces the SFR in low-density regions and shows better agreement with the observations of the Kennicutt-Schmidt law. This model also suppresses the early star formation and shifts the peak of the cosmic SF history toward lower redshift, more consistently with the recent observational estimates of cosmic SFR density. The simulations with the new SF model also predict lower global stellar mass densities at high-z, larger populations of low-mass galaxies and a higher gas fraction in high-z galaxies. Our results suggest that there is room left in the model uncertainties to reconcile the discrepancy that was found between the theory and observations of cosmic SF history and stellar mass density. Nevertheless, our simulations still predict higher stellar mass densities than most of the observational estimates.

We summarize the results from numerical simulations of mass outflows from AGN. We focus on simulations of outflows driven by radiation from large-scale inflows. We discuss the properties of these outflows in the context of the so-called AGN feedback problem. Our main conclusion is that this type of outflows are efficient in removing matter but inefficient in removing energy.

(Abridged) We measure and analyze the energy, momentum, and mass feedback efficiencies due to radiation from AGN in relatively large scale outflows. Our measurements are based on the axisymmetric and time-dependent radiation-hydrodynamical simulations recently presented in Kurosawa & Proga. These simulations follow dynamics of gas under the influences of the AGN radiation and the gravity of the central 10^8 Msun black hole on scales from ~0.01 to ~10pc. We compare our Mdot_a-rho_o relation with that predicted by the Bondi accretion model. For high luminosities comparable to the Eddington limit, the power-law fit (Mdot_a ~ rho_o^q) to our models yields q=~0.5 instead of q=1.0 which is predicted by the Bondi model. This difference is caused by the outflows which are important for the overall mass budget at high luminosities. The maximum momentum and mass feedback efficiencies found in our models are ~10^-2 and ~10^-1, respectively. However, the outflows are much less important energetically: their thermal and kinetic powers in units of the radiative luminosity are ~10^-5 and ~10^-4, respectively. The efficiencies do not increase monotonically with the accretion luminosity but rather peak around the Eddington limit beyond which a steady state disk-wind-like solution exists. Our energy feedback efficiencies are significantly lower than 0.05, which is required in some cosmological and galaxy merger simulations. The low feedback efficiencies found here could have significant implications on the mass growth of super massive black holes in the early universe. We stress however that we have not considered the innermost parts of the accretion and outflow where radiation and matter interact most strongly. The feedback from this region could have efficiencies significantly above the low values found here.

We examine the past and current work on the star formation (SF) histories of dwarf galaxies in cosmological hydrodynamic simulations. The results obtained from different numerical methods are still somewhat mixed, but the differences are understandable if we consider the numerical and resolution effects. It remains a challenge to simulate the episodic nature of SF history in dwarf galaxies at late times within the cosmological context of a cold dark matter model. More work is needed to solve the mysteries of SF history of dwarf galaxies employing large-scale hydrodynamic simulations on the next generation of supercomputers.

When, and how, did the first galaxies and supermassive black holes (SMBH) form, and how did they reionization the Universe? First galaxy formation and cosmic reionization are among the last frontiers in studies of cosmic structure formation. We delineate the detailed astrophysical probes of early galaxy and SMBH formation afforded by observations at centimeter through submillimeter wavelengths. These observations include studies of the molecular gas (= the fuel for star formation in galaxies), atomic fine structure lines (= the dominant ISM gas coolant), thermal dust continuum emission (= an ideal star formation rate estimator), and radio continuum emission from star formation and relativistic jets. High resolution spectroscopic imaging can be used to study galaxy dynamics and star formation on sub-kpc scales. These cm and mm observations are the necessary compliment to near-IR observations, which probe the stars and ionized gas, and X-ray observations, which reveal the AGN. Together, a suite of revolutionary observatories planned for the next decade from centimeter to X-ray wavelengths will provide the requisite panchromatic view of the complex processes involved in the formation of the first generation of galaxies and SMBHs, and cosmic reionization.

We present the results of a numerical study on the effects of metal enrichment and metal cooling on galaxy formation and cosmic star formation (SF) history using cosmological hydrodynamic simulations. We find following differences in the simulation with metal cooling when compared to the run without it: (1) the cosmic star formation rate (SFR) is enhanced by about 50 & 20% at z=1 & 3, respectively; (2) the gas mass fraction in galaxies is lower; (3) the total baryonic mass function (gas + star) at z=3 does not differ significantly, but shows an increase in the number of relatively massive galaxies at z=1; (4) the baryonic mass fraction of intergalactic medium (IGM) is reduced at z<3 due to more efficient cooling and gas accretion onto galaxies. Our results suggest that the metal cooling enhances the galaxy growth by two different mechanisms: (1) increase of SF efficiency in the local interstellar medium (ISM), and (2) increase of IGM accretion onto galaxies. The former process is effective throughout most of the cosmic history, while the latter is effective only at z<3 when the IGM is sufficiently enriched by metals owing to feedback.

We study the incidence rate of damped Ly-a systems associated with the host galaxies of gamma-ray bursts (GRB-host-DLAs) as functions of neutral hydrogen column density (N_HI) and projected star formation rate (SFR) using cosmological SPH simulations. Assuming that the occurrence of GRBs is correlated with the local SFR, we find that the median N_HI of GRB-host-DLAs progressively shifts to lower N_HI values with increasing redshift, and the incidence rate of GRB-host-DLAs with log N_HI > 21.0 decreases rapidly at z>=6. Our results suggest that the likelihood of observing the signature of IGM attenuation in GRB afterglows increases towards higher redshift, because it will not be blocked by the red damping wing of DLAs in the GRB host galaxies. This enhances the prospects of using high-redshift GRBs to probe the reionization history of the Universe. The overall incidence rate of GRB-host-DLAs decreases monotonically with increasing redshift, whereas that of QSO-DLAs increases up to z=6. A measurement of the difference between the two incidence rates would enable an estimation of the value of \eta_grb, which is the mass fraction of stars that become GRBs for a given amount of star formation. Our predictions can be tested by upcoming high-z GRB missions, including JANUS (Joint Astrophysics Nascent Universe Scout) and SVOM (Space multi-band Variable Object Monitor).

Damped Lyman-alpha systems (DLAs) are useful probes of star formation and galaxy formation at high redshift. We study the physical properties of DLAs and their relationship to Lyman-break galaxies using cosmological hydrodynamic simulations based on the concordance Lambda cold dark matter model. Fundamental statistics such as global neutral hydrogen (HI) mass density, HI column density distribution function, DLA rate-of-incidence and mean halo mass of DLAs are reproduced reasonably well by the simulations, but with some deviations that need to be understood better in the future. We discuss the feedback effects by supernovae and galactic winds on the DLA distribution. We also compute the [C_II] emission from neutral gas in high-z galaxies, and make predictions for the future observations by ALMA and SPICA. Agreement and disagreement between simulations and observations are discussed, as well as the future directions of our DLA research.

We present the results of a numerical study of a galactic wind model and its implications on the properties of damped Lyman-alpha absorbers (DLAs) using cosmological hydrodynamic simulations. We vary both the wind strength and the internal parameters of the the wind model in a series of cosmological SPH simulations that include radiative cooling and heating by a UV background, star formation, and feedback from supernovae and galactic winds. To test our simulations, we examine the DLA `rate-of-incidence' as a function of halo mass, galaxy apparent magnitude, and impact parameter. We find that the statistical distribution of DLAs does not depend on the exact values of internal numerical parameters that control the decoupling of hydrodynamic forces when the gas is ejected from starforming regions. The DLA rate-of-incidence in our simulations at z=3 is dominated by the faint galaxies with apparent magnitude R_AB < 25.5. However, interestingly in a `strong wind' run, the differential distribution of DLA sight-lines is peaked at Mhalo = 10^{12} Msun/h (R_AB~27), and the mean DLA halo mass is Mmean=10^{12.4} Msun/h (R_AB ~ 26). These mass-scales are much larger than those if we ignore winds, because galactic wind feedback suppresses the DLA cross section in low-mass halos and increases the relative contribution to the DLA incidence from more massive halos. The DLAs in our simulations are more compact than the present-day disk galaxies, and the impact parameter distribution is very narrow unless we limit the search for the host galaxy to only bright LBGs. The comoving number density of DLAs is higher than that of LBGs down to R_AB=30 mag if the physical radius of each DLA is smaller than 5 kpc/h_70. We discuss conflicts between current simulations and observations, and potential problems with simulations based on the CDM model.

We measure the three-dimensional topology of large-scale structure in the Sloan Digital Sky Survey (SDSS). This allows the genus statistic to be measured with unprecedented statistical accuracy. The sample size is now sufficiently large to allow the topology to be an important tool for testing galaxy formation models. For comparison, we make mock SDSS samples using several state-of-the-art N-body simulations: the Millennium run of Springel et al. (2005)(10 billion particles), Kim & Park (2006) CDM models (1.1 billion particles), and Cen & Ostriker (2006) hydrodynamic code models (8.6 billion cell hydro mesh). Each of these simulations uses a different method for modeling galaxy formation. The SDSS data show a genus curve that is broadly characteristic of that produced by Gaussian random phase initial conditions. Thus the data strongly support the standard model of inflation where Gaussian random phase initial conditions are produced by random quantum fluctuations in the early universe. But on top of this general shape there are measurable differences produced by non-linear gravitational effects (cf. Matsubara 1994), and biasing connected with galaxy formation. The N-body simulations have been tuned to reproduce the power spectrum and multiplicity function but not topology, so topology is an acid test for these models. The data show a ``meatball'' shift (only partly due to the Sloan Great Wall of Galaxies; this shift also appears in a sub-sample not containing the Wall) which differs at the 2.5\sigma level from the results of the Millennium run and the Kim & Park dark halo models, even including the effects of cosmic variance.

Taking three independent approaches, we investigate the simultaneous constraints set on the cosmic star formation history from various observations, including stellar mass density and extragalactic background light (EBL). We compare results based on: 1) direct observations of past light-cone, 2) a model using local fossil evidence constrained by SDSS observations at z~0 (the `Fossil' model), and 3) theoretical ab initio models from three calculations of cosmic star formation history: (a) new (1024)^3 Total Variation Diminishing (TVD) cosmological hydrodynamic simulation, (b) analytic expression of Hernquist & Springel based on cosmological Smoothed Particle Hydrodynamics (SPH) simulations, and (c) semi-analytic model of Cole et al. We find good agreement among the three independent approaches up to the order of observational errors, except that all the models predict bolometric EBL of I_tot ~= 37-52 nW m^-2 sr^-1, which is at the lower edge of the the observational estimate by Hauser & Dwek. We emphasize that the Fossil model that consists of two components -- spheroids and disks --, when normalized to the local observations, provides a surprisingly simple but accurate description of the cosmic star formation history and other observable quantities. Our analysis suggests that the consensus global parameters at z=0 are: Omega_* = 0.0023+-0.0004, I_EBL = 43+-7 nW m^-2 sr^-1 rho_SFR=(1.06+-0.22)e-2 Msun yr^-1 Mpc^-3, j_bol = (3.1+-0.2)e8 Lsun Mpc^-3.

We discuss the detectability of high-redshift galaxies via [CII] 158 micron line emission by coupling an analytic model with cosmological Smoothed Particle Hydrodynamics (SPH) simulations that are based on the concordance Lambda cold dark matter (CDM) model. Our analytic model describes a multiphase interstellar medium irradiated by the far ultra-violet radiation from local star-forming regions, and it calculates thermal and ionization equilibrium between cooling and heating. The model allows us to predict the mass fraction of a cold neutral medium (CNM) embedded in a warm neutral medium (WNM). Our cosmological SPH simulations include a treatment of radiative cooling/heating, star formation, and feedback effects from supernovae and galactic winds. Using our method, we make predictions for the [CII] luminosity from high-redshift galaxies which can be directly compared with upcoming observations by the Atacama Large Millimeter Array (ALMA) and the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We find that the number density of high-redshift galaxies detectable by ALMA and SPICA via [CII] emission depends significantly on the amount of neutral gas which is highly uncertain. Our calculations suggest that, in a CDM universe, most [CII] sources at z=3 are faint objects with \Snu < 0.01 mJy. Lyman-break galaxies (LBGs) brighter than R_AB=23.5 mag are expected to have flux densities \Snu = 1-3 mJy depending on the strength of galactic wind feedback. The recommended observing strategy for ALMA and SPICA is to aim at very bright LBGs or star-forming DRG/BzK galaxies.

Cosmological hydrodynamic simulations have significantly improved over the past several years, and we have already shown that the observed properties of Lyman-break galaxies (LBGs) at z=3 can be explained well by the massive galaxies in the simulations. Here we extend our study to z=6 and show that we obtain good agreement for the LBGs at the bright-end of the luminosity function (LF). Our simulations also suggest that the cosmic star formation rate density has a peak at z= 5-6, and that the current LBG surveys at z=6 are missing a significant number of faint galaxies that are dimmer than the current magnitude limit. Together, our results suggest that the universe could be reionized at z=6 by the Pop II stars in ordinary galaxies. We also estimate the LF of Lyman-alpha emitters (LAEs) at z=6 by relating the star formation rate in the simulation to the Ly-alpha luminosity. We find that the simulated LAE LFs agree with the observed data provided that the net escape fraction of Ly-alpha photon is f_{Ly-alpha} <= 0.1. We investigate two possible scenarios for this effect: (1) all sources in the simulation are uniformly dimmer by a factor of 10 through attenuation, and (2) one out of ten LAEs randomly lights up at a given moment. We show that the correlation strength of the LAE spatial distribution can possibly distinguish the two scenarios.

We perform a spectrophotometric analysis of galaxies at redshifts z = 4 - 6 in cosmological SPH simulations of a Lambda CDM universe. Our models include radiative cooling and heating by a uniform UV background, star formation, supernova feedback, and a phenomenological model for galactic winds. Analysing a series of simulations of varying boxsize and particle number allows us to isolate the impact of numerical resolution on our results. Specifically, we determine the luminosity functions in B, V, R, i', and z' filters, and compare the results with observed galaxy surveys done with the Subaru telescope and the Hubble Space Telescope. We find that the simulated galaxies have UV colours consistent with observations and fall in the expected region of the colour-colour diagrams used by the Subaru group. Assuming a uniform extinction of E(B-V) = 0.15, we also find reasonable agreement between simulations and observations in the space density of UV bright galaxies at z = 3 - 6, down to the magnitude limit of each survey. For the same moderate extinction level of E(B-V) ~ 0.15, the simulated luminosity functions match observational data, but have a steep faint-end slope with alpha ~ -2.0. We discuss the implications of the steep faint-end slope found in the simulations.

Recent observations have revealed a population of red massive galaxies at high redshift which are challenging to explain in hierarchical galaxy formation models. We analyze this ``massive galaxy problem'' with two different types of hydrodynamic simulations -- Eulerian TVD and SPH -- of a concordance Lambda CDM universe. We consider two separate but connected aspects of the problem posed by these extremely red objects (EROs): (1) the mass-scale, and (2) red colors. We perform spectrophotometric analyses of simulated galaxies in B, z, R, I, Js, Ks, K filters, and compare their near-IR properties with observations at redshift z=1-3. We find that the simulated galaxies brighter than the magnitude limit of K_vega=20 mag have stellar masses Mstar >= 10^{11} Msun/h and a number density of a few \times 10^{-4} h^3 Mpc^{-3} at z~2, in good agreement with the observed number density in the K20 survey. Therefore, our hydrodynamic simulations do not exhibit the ``mass-scale problem''. The answer to the ``redness problem'' is less clear because of our poor knowledge of the amount of dust extinction in EROs and the uncertain fraction of star-forming EROs. However, our simulations can account for the observed comoving number density of ~ 1e-4 Mpc^{-3} at z=1-2 if we assume a uniform extinction of E(B-V)=0.4 for the entire population of simulated galaxies. The critical question of the existence of a significant population of old, passively evolving massive EROs should continue to be tested with future hydrodynamic simulations. Upcoming observations of the thermal emission of dust in 24 micron by the Spitzer Space Telescope will help to better estimate the dust content of EROs at z=1-3, and thus to further constrain the star formation history of the Universe, and theoretical models of galaxy formation. (abridged)

We compare two cosmological hydrodynamic simulation codes in the context of hierarchical galaxy formation: The SPH code GADGET, and the Eulerian AMR code ENZO. Both codes represent dark matter with the N-body method, but use different gravity solvers and fundamentally different approaches to hydrodynamics. We compare the GADGET `entropy conserving' SPH formulation with two ENZO methods: The piecewise parabolic method (PPM), and the artificial viscosity-based scheme used in the ZEUS code. In this paper we focus on a comparison of cosmological simulations that follow either only dark matter, or also adiabatic baryonic gas. The dark matter-only runs agree generally quite well, provided ENZO is run with a comparatively fine root grid and a low overdensity threshold for mesh refinement, otherwise the abundance of low-mass halos is suppressed. This is due to the hierarchical particle-mesh method used to compute gravitational forces in ENZO, which tends to deliver lower force resolution than the tree algorithm of GADGET. At comparable force resolution, we find that the latter offers substantially better performance and lower memory consumption than the present gravity solver in ENZO. In simulations that include adiabatic gas dynamics, we find general agreement in the distribution functions of temperature, entropy, and density for gas of moderate to high overdensity, as found inside dark matter halos. However, there are some significant differences at lower overdensities. We argue that these discrepancies are presumably owing to differences in the shock-capturing abilities of the different methods. In particular, ZEUS hydro leads to some unphysical heating at early times in preshock regions. Overall, the GADGET hydro results are bracketed by those for ENZO/ZEUS and ENZO/PPM. (abridged)

We study the properties of galaxies at z=2 in a Lambda cold dark matter universe, using two different types of hydrodynamic simulation methods -- Eulerian TVD and smoothed particle hydrodynamics (SPH) -- and a spectrophotometric analysis in the U_n, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) and Steidel et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 1e-2 h^3 Mpc^-3 for E(B-V)=0.15, which is roughly twice that of the number density found by Erb et al. (2004) for the UV bright sample. This suggests that roughly half of the massive galaxies with M*>10^{10} Msun/h at z=2 are UV bright population, and the other half is bright in the infra-red wavelengths. The most massive galaxies at z=2 have stellar masses >= 10^{11-12} Msun. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, together with significant contribution by starbursts reaching up to 1000 Msun/yr which lie on top of the continuous component. TVD simulations indicate a more sporadic star formation history than the SPH simulations. Our results do not imply that hierarchical galaxy formation fails to account for the observed massive galaxies at z>=1. The global star formation rate density in our simulations peaks at z>=5, a much higher redshift than predicted by the semianalytic models. This star formation history suggests early build-up of the stellar mass density, and predicts that 70 (50, 30)% of the total stellar mass at z=0 had already been formed by z=1 (2, 3). Upcoming observations by Spitzer and Swift might help to better constrain the star formation history at high redshift.

We study the photometric properties of Lyman-break galaxies (LBGs) formed by redshift z=3 in a set of large cosmological smoothed-particle hydrodynamics simulations of the Lambda cold dark matter (CDM) model. Our numerical simulations include radiative cooling and heating with a uniform UV background, star formation, supernova feedback, and a phenomenological model for galactic winds. Analysing a series of simulations of varying boxsize and particle number allows us to isolate the impact of numerical resolution on our results. We compute spectra of simulated galaxies using a population synthesis model, and derive colours and luminosity functions of galaxies at z=3 after applying local dust extinction and absorption by the intergalactic medium (IGM). We find that the simulated galaxies have U-G and G-R colours consistent with observations, provided that intervening absorption by the IGM is applied. The observed properties of LBGs, including their number density, colours, and luminosity functions, can be explained if LBGs are identified with the most massive galaxies at z=3, having typical stellar mass of M_{star} ~ 1e10 Msun/h, a conclusion broadly consistent with earlier studies based on hydrodynamic simulations of the Lamda CDM model. We also find that most simulated LBGs were continuously forming stars at a high rate for more than one Gyr up until z=3, but with numerous starbursts lying on top of the continuous component. Interestingly, our simulations suggest that more than 50% of the total stellar mass and star formation rate in the Universe are accounted for by galaxies that are not detected in the current generation of LBG surveys.

We study the properties of galaxies at z=2 in a Lambda CDM universe, using two different types of hydrodynamic simulation methods (Eulerian TVD and SPH) and a spectrophotometric analysis in the Un, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 2e-2 h^3/Mpc^3, roughly one order of magnitude larger than that of Lyman break galaxies at z=3. The most massive galaxies at z=2 have stellar masses >~1e11 Msun, and their observed-frame G-R colors lie in the range 0.0<G-R<1.0. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, although the TVD simulation indicates a more sporadic star formation history than the SPH simulations. Of order half of their stellar mass was already assembled by z~4. The reddest massive galaxies at z=2 with G-R >= 1.0 and Mstar>1e10 Msun/h finished the build-up of their stellar mass by z~3. Interestingly, our study suggests that the majority of the most massive galaxies at z=2 should be detectable at rest-frame UV wavelengths, contrary to some recent claims made on the basis of near-IR studies of galaxies at the same epoch, provided the median extinction is less than E(B-V)<0.3. However, our results also suggest that the fraction of stellar mass contained in galaxies that pass the color-selection criteria could be as low as 50% of the total stellar mass in the Universe at z=2. Our simulations suggest that the missing stellar mass is contained in fainter (R>25.5) and intrinsically redder galaxies. Our results do not suggest that hierarchical galaxy formation fails to account for the massive galaxies at z>=1. (abridged)

We investigate possible signatures of feedback from galactic superwinds on the metallicity of the Lyman alpha forest, using a set of high resolution hydrodynamic simulations of a Lambda cold dark matter model. Simulations produce metals self-consistently, based on one single parameter, the metal yield, which in turn is constrained by metallicity in the intra-cluster gas. We follow metals as a separate density species. For the metallicity of Lyman alpha clouds with column density of N_{HI}~10^{14.5}-10^{15.5}/cm^2 at z=2-4 we find reasonable agreement between simulations, both with and without GSW, and observations (Schaye et al). A unique signature and sensitive test of GSW is, however, provided by lower density regions with gas density of \rho/<\rho> = 0.01-1.0 and a corresponding column density of 10^{12}-10^{14}/cm^2. Without GSW we predict that both the mean and median metallicity of \lya clouds in this column density range at z=2-4 should have Z<10^{3}Zsun, since these small systems support little star formation. GSW contaminate these regions, however, and also there is a significant fraction (~25%) of Lyman alpha clouds in this column density range which have a high metallicity excess of 10^{-2}Zsun, resulting in a mean metallicity of ~10^{-2}Zsun. In addition, we find that there is a minimum in the median metallicity for clouds of N_{HI}~10^{13}-10^{14}/cm^2 in the case with GSW, whereas without GSW the metallicity decreases monotonically and rapidly with decreasing column density. Finally, we predict that the ratio of secondary (e.g., N) to primary metals (e.g., O,C) is expected to be smaller by a factor of 10 in clouds of N_{HI}~10^{14.5}/cm^2 compared to that in large galaxies; this factor increases to >50 for N_{HI}< 10^{13.5}/cm^2.

Zero-field muon-spin-relaxation (muSR) and magnetic-susceptibility measurements have been carried out in La_2-x_Sr_x_Cu_1-y_Zn_y_O_4_ with x = 0.10, 0.115 and 0.13 changing the Zn-concentration y finely up to 0.10, with the aim to clarify effects of the nonmagnetic impurity Zn on the Cu-spin dynamics and superconductivity. The muSR measurements have revealed that, in each x, a magnetic order of Cu spins is observed at y ~ 0.0075, while it disappears and Cu spins turn into a fast fluctuating state for y > 0.03. From the magnetic-susceptibility measurements, on the other hand, it has been found that the volume fraction of the superconducting state rapidly decreases through the slight doping of Zn and that its y dependence corresponds to the y dependence of the volume fraction of the fast fluctuating region of Cu spins estimated from the muSR results. Both the rapid decrease of the volume fraction of the superconducting state and the formation of the magnetic order through the slight doping of Zn can be interpreted in terms of the development of the so-called "Swiss cheese" model. That is, it is concluded that Cu spins in a non-superconducting region around Zn exhibit slowing-down of the fluctuations or form an incoherent or coherent static magnetic order. Moreover, the formation of a non-superconducting region around Zn is considered to be due to the pinning of the dynamical spin correlation or dynamical stripe correlations.

We study the evolution of the global stellar mass density in a Lambda cold dark matter universe using two different types of hydrodynamical simulations (Eulerian TVD and SPH) and the analytical model of Hernquist & Springel (2003). We find that the theoretical calculations all predict both a higher stellar mass density at z~3 than indicated by current observations, and that the peak of the cosmic star formation rate history should lie at z~5. Such a star formation history implies that as much as (70%, 30%) of the total stellar mass density today must already have formed by z=(1, 3). Our results suggest that current observations at z~3 are missing as much as 50% of the total stellar mass density in the Universe, perhaps owing to an inadequate allowance for dust obscuration in star-forming galaxies, limited sample sizes, or cosmic variance. We also compare our results with some of the updated semi-analytic models of galaxy formation.

We simulate the evolution of the intergalactic medium (IGM) in a universe dominated by a cosmological constant. We find that within a few Hubble times from the present epoch, the baryons will have two primary phases: one phase composed of low-density, low-temperature, diffuse, ionized gas which cools exponentially with cosmic time due to adiabatic expansion, and a second phase of high-density, high-temperature gas in virialized dark matter halos which cools much more slowly by atomic processes. The mass fraction of gas in halos converges to ~40% at late times, about twice its calculated value at the present epoch. We find that in a few Hubble times, the large scale filaments in the present-day IGM will rarefy and fade away into the low-temperature IGM, and only islands of virialized gas will maintain their physical structure. We do not find evidence for fragmentation of the diffuse IGM at later times. More than 99% of the gas mass will maintain a steady ionization fraction above 80% within a few Hubble times. The diffuse IGM will get extremely cold and dilute but remain highly ionized, as its recombination time will dramatically exceed the age of the universe.

We use a large-scale Lambda-CDM hydrodynamical simulation to assess the dependence of the cosmic Star Formation Rate (SFR) on overdensity of luminosity. The `morphology-density relation' is qualitatively captured when the SFR is plotted as a function of redshift, divided into the quartiles of the overdensity distribution.

We use cosmological smoothed-particle hydrodynamics (SPH) simulations of the Lambda cold dark matter (CDM) model to study the abundance of damped Lyman-alpha absorbers (DLAs) in the redshift range z=0-4.5. We compute the cumulative DLA abundance by using the relation between DLA cross-section and the total halo mass inferred from the simulations. Our approach includes a phenomenological model of galactic wind. We employ the ``conservative entropy'' formulation of SPH, and utilise a series of simulations of varying box-size and particle number to isolate the impact of numerical resolution on our results. We show that the DLA abundance was overestimated in previous studies for three reasons: (1) the overcooling of gas occurring with non-conservative formulations of SPH, (2) a lack of numerical resolution, and (3) an inadequate treatment of feedback. Our new results for the total neutral hydrogen mass density, DLA abundance, and column density distribution function all agree reasonably well with observational estimates at redshift z=3, indicating that DLAs arise naturally from radiatively cooled gas in dark matter haloes that form in a Lambda CDM universe. Our simulations suggest a moderate decrease in DLA abundance by roughly a factor of two from z=4.5 to 3, consistent with observations. A significant decline in abundance from z=3 to z=1, followed by weak evolution from z=1 to z=0, is also indicated, but our low-redshift results need to be interpreted with caution because they are based on coarser simulations than the ones employed at high redshift. Our highest resolution simulation also suggests that the halo mass-scale below which DLAs do not exist is slightly above 10^8 Msun/h at z=3-4, somewhat lower than previously estimated. (abridged)

We study the distribution of the star formation rate and metallicity of damped Lyman-alpha absorbers using cosmological SPH simulations of the Lambda cold dark matter model in the redshift range z=0-4.5. Our approach includes a phenomenological model of galactic wind. We find that there is a positive correlation between the projected stellar mass density and the neutral hydrogen column density (NHI) of DLAs for high NHI systems, and that there is a good correspondence in the spatial distribution of stars and DLAs in the simulations. The evolution of typical star-to-gas mass ratios in DLAs can be characterised by an increase from about 2 at z=4.5 to 3 at z=3, to 10 at z=1, and finally to 20 at z=0. We also find that the projected SFR density in DLAs follows the Kennicutt law well at all redshifts, and the simulated values are consistent with the recent observational estimates of this quantity by Wolfe et al. (2003a,b). The rate of evolution in the mean metallicity of simulated DLAs as a function of redshift is mild, and is consistent with the rate estimated from observations. The predicted metallicity of DLAs is generally sub-solar in our simulations, and there is a significant scatter in the distribution of DLA metallicity for a given NHI. However, we find that the median metallicity of simulated DLAs is close to that of Lyman-break galaxies, which is higher than the values typically observed for DLAs by nearly an order of magnitude. This discrepancy with observations could be due to an inadequate treatment of SN feedback in our current simulations, perhaps indicating that metals are not expelled efficiently enough from DLAs by outflows. Alternatively, the current observations might be missing the majority of the high metallicity DLAs due to selection effects. (abridged)

We study the distribution of star formation rate (SFR) and metallicity of damped Lyman-alpha absorbers (DLAs) using cosmological smoothed particle hydrodynamics (SPH) simulations of the Lambda cold dark matter model. Our simulations include a phenomenological model for feedback by galactic winds which allows us to examine the effect of galactic outflows on the distribution of SFR and metallicity of DLAs. For models with strong galactic winds, we obtain good agreement with recent observations with respect to total neutral hydrogen mass density, N_HI column-density distribution, abundance of DLAs, and for the distribution of SFR in DLAs. However, we also find that the median metallicity of simulated DLAs is higher than the values typically observed by nearly an order of magnitude. This discrepancy with observations could be due to shortcomings in the treatment of the supernova feedback or the multiphase structure of the gas in our current simulations. Recent observations by Wolfe et al. (2003a,b) seem to point to the same problem; i.e. the observed DLA metallicities are much lower than those expected from the (either observed or simulated) DLA star formation rates, a puzzle that has been known as the ``missing metals'' problem for the globally averaged quantities.

While second and higher order correlations of the light distribution have received extensive study, the lowest order probability distribution function (PDF) -- the probability that a unit volume of space will emit a given amount of light -- has received very little attention. We estimate this function with the aid of hydrodynamic simulations of the Lambda CDM model, finding it significantly different from the mass density PDF, and not simply related to it by linear bias or any of the other prescriptions commonly adopted. If the optical light PDF is, in reality, similar to what we find in the simulations, then some measures of Omega_M based on mass-to-light ratio and the cosmic virial theorem will have significantly underestimated Omega_M. Basically, the problem is one of selection bias, with galaxy forming regions being unrepresentative of the dark matter distribution in a way not described by linear bias. Knowledge of the optical PDF and the plausible assumption of a log-normal distribution for the matter PDF will allow one to correct for these selection biases. We find that this correction (which amounts to 20-30%) brings the values of Omega_M estimated by using the mass-to-light ratio and the cosmic virial theorem to the range 0.2-0.3, in better agreement with the WMAP result than the uncorrected estimates. In addition, the relation between mass and light PDFs gives us insight concerning the nature of the void phenomenon. In particular our simulation indicates that 20% of mass is distributed in voids, which occupy 85% of volume in the universe.

We simulate the future evolution of the observed inhomogeneities in the local universe assuming that the global expansion rate is dominated by a cosmological constant. We find that within two Hubble times (~ 30 billion years) from the present epoch, large-scale structures will freeze in comoving coordinates and the mass distribution of bound objects will stop evolving. The Local Group will get somewhat closer to the Virgo cluster in comoving coordinates, but will be pulled away from the Virgo in physical coordinates due to the accelerated expansion of the Universe. In the distant future there will only be one massive galaxy within our event horizon, namely the merger product of the Andromeda and the Milky Way galaxies. All galaxies that are not gravitationally bound to the Local Group will recede away from us and eventually exit from our event horizon. More generally, we identify the critical interior overdensity above which a shell of matter around an object will remain bound to it at late times.

We study the evolution of Lyman Break Galaxies (LBGs) from z=5 to z=0 by tracing the merger trees of galaxies in a large-scale hydrodynamic simulation based on a Lambda cold dark matter model. In particular, we emphasize on the range of properties of the sample selected by the rest-frame V band luminosity, in accordance with recent near-IR observations. The predicted rest-frame V band luminosity function agrees well with the observed one when dust extinction is taken into account. The stellar content and the star formation histories of LBGs are also studied. We find that the LBGs intrinsically brighter than Mv=-21.0 at z=3 have stellar masses of at least 10^9\Msun, with a median of 10^{10}h^{-1}\Msun. The brightest LBGs (Mv<-23) at z=3 merge into clusters/groups of galaxies at z=0, as suggested from clustering studies of LBGs. Roughly one half of the galaxies with -23<Mv<-22 at z=3 fall into groups/clusters, and the other half become typical L* galaxies at z=0 with stellar mass of ~10^{11}\Msun. Descendants of LBGs at the present epoch have formed roughly 30% of their stellar mass by z=3, and the half of their current stellar population is 10 Gyr old, favoring the scenario that LBGs are the precursors of the present day spheroids. We find that the most luminous LBGs have experienced a starburst within 500 Myr prior to z=3, but also have formed stars continuously over a period of 1 Gyr prior to z=3 when all the star formation in progenitors is coadded. We also study the evolution of the mean stellar metallicity distribution of galaxies, and find that the entire distribution shifts to lower metallicity at higher redshift. The observed sub-solar metallicity of LBGs at z=3 is naturally predicted in our simulation.

The luminosity function of galaxies is derived from a cosmological hydrodynamic simulation of a Lambda cold dark matter (CDM) universe with the aid of a stellar population synthesis model. At z=0, the resulting B band luminosity function has a flat faint end slope of \alpha \approx -1.15 with the characteristic luminosity and the normalization in a fair agreement with observations, while the dark matter halo mass function is steep with a slope of \alpha \approx -2. The colour distribution of galaxies also agrees well with local observations. We also discuss the evolution of the luminosity function, and the colour distribution of galaxies from z=0 to 5. A large evolution of the characteristic mass in the stellar mass function due to number evolution is compensated by luminosity evolution; the characteristic luminosity increases only by 0.8 mag from z=0 to 2, and then declines towards higher redshift, while the B band luminosity density continues to increase from z=0 to 5 (but only slowly at z>3).

We study star formation history and stellar metallicity distribution in galaxies in a Lambda cold dark matter universe using a hydrodynamic cosmological simulation. Our model predicts star formation rate declining in time exponentially from an early epoch to the present with the time-scale of 6 Gyr, which is consistent with the empirical Madau plot with modest dust obscuration. Star formation in L^* galaxies continues intermittently to the present also with an exponentially declining rate of a similar time-scale, whereas in small galaxies star formation ceases at an early epoch. The mean age of the extant stars decreases only slowly with increasing redshift, and exceeds 1 Gyr at z=3. Normal galaxies contain stars with a wide range of metallicity and age: stars formed at z<1 have metallicity of 0.1-1.0 Zsun, while old stars take a wide range of values from 10^{-6} Zsun to 3.0 Zsun. The mean metallicity of normal galaxies is in the range 0.1-1.0 Zsun. Dwarf galaxies that contain only old stars have a wide range of mean metallicity (10^{-4}-1.0 Zsun), but on average they are metal deficient compared with normal galaxies.

We carry out an extensive study of the cosmic Mach number (\mach) on scales of R=5, 10 and 20h^-1Mpc using an LCDM hydrodynamical simulation. We particularly put emphasis on the environmental dependence of \mach on overdensity, galaxy mass, and galaxy age. We start by discussing the difference in the resulting \mach according to different definitions of \mach and different methods of calculation. The simulated Mach numbers are slightly lower than the linear theory predictions even when a non-linear power spectrum was used in the calculation, reflecting the non-linear evolution in the simulation. We find that the observed \mach is higher than the simulated mean by more than 2-standard deviations, which suggests either that the Local Group is in a relatively low-density region or that the true value of \Omega_m is ~ 0.2, significantly lower than the simulated value of 0.37. We show from our simulation that the Mach number is a weakly decreasing function of overdensity. We also investigate the correlations between galaxy age, overdensity and \mach for two different samples of galaxies --- DWARFs and GIANTs. Older systems cluster in higher density regions with lower \mach, while younger ones tend to reside in lower density regions with larger \mach, as expected from the hierarchical structure formation scenario. However, for DWARFs, the correlation is weakened by the fact that some of the oldest DWARFs are left over in low-density regions during the structure formation history. For giant systems, one expects blue-selected samples to have higher \mach than red-selected ones. We briefly comment on the effect of the warm dark matter on the expected Mach number.

We compute the cosmic star formation rate (SFR) and the rest-frame comoving luminosity density in various pass-bands as a function of redshift using large-scale \Lambda-CDM hydrodynamical simulations with the aim of understanding their behavior as a function of redshift. To calculate the luminosity density of galaxies, we use an updated isochrone synthesis model which takes metallicity variations into account. The computed SFR and the UV-luminosity density have a steep rise from z=0 to 1, a moderate plateau between z=1 - 3, and a gradual decrease beyond z=3. The raw calculated results are significantly above the observed luminosity density, which can be explained either by dust extinction or the possibly inappropriate input parameters of the simulation. We model the dust extinction by introducing a parameter f; the fraction of the total stellar luminosity (not galaxy population) that is heavily obscured and thus only appears in the far-infrared to sub-millimeter wavelength range. When we correct our input parameters, and apply dust extinction with f=0.65, the resulting luminosity density fits various observations reasonably well, including the present stellar mass density, the local B-band galaxy luminosity density, and the FIR-to-submm extragalactic background. Our result is consistent with the picture that \sim 2/3 of the total stellar emission is heavily obscured by dust and observed only in the FIR. The rest of the emission is only moderately obscured which can be observed in the optical to near-IR wavelength range. We also argue that the steep falloff of the SFR from z=1 to 0 is partly due to the shock-heating of the universe at late times, which produces gas which is too hot to easily condense into star-forming regions.