The impact of the streaming between baryons and dark matter on the first structures has been actively explored by recent studies. We investigate how the key results are affected by two popular approximations. One is to implement the streaming by accounting for only the relative motion while assuming "baryons trace dark matter" spatially at the initialization of simulation. This neglects the smoothing on the gas density taking place before the initialization. In our simulation initialized at $z_i=200$, it overestimates the gas density power spectrum up to 40\% at $k\approx10^2~h~\mbox{Mpc}^{-1}$ at $z=20$. Halo mass ($M_h$) and baryonic fraction in halos ($f_{b,h}$) are also overestimated, but the relation between the two remains unchanged. The other approximation tested is to artificially amplify the density/velocity fluctuations in the cosmic mean density to simulate the first minihalos that form in overdense regions. This gives a head start to the halo growth while the subsequent growth is similar to that in the mean density. The growth in a true overdense region, on the other hand, is accelerated gradually in time. For example, raising $\sigma_8$ by 50\% effectively transforms $z\rightarrow\sqrt{1.5}z$ in the halo mass growth history while in 2-$\sigma$ overdensity, the growth is accelerated by a constant in redshift: $z\rightarrow{z+4.8}$. As a result, halos are more grown in the former than in the latter before $z\approx27$ and vice versa after. The $f_{b,h}$-$M_h$ relation is unchanged in those cases as well, suggesting that the Pop III formation rate for a given $M_h$ is insensitive to the tested approximations.

Line intensity mapping (LIM) is an emerging observational method to study the large-scale structure of the Universe and its evolution. LIM does not resolve individual sources but probes the fluctuations of integrated line emissions. A serious limitation with LIM is that contributions of different emission lines from sources at different redshifts are all confused at an observed wavelength. We propose a deep learning application to solve this problem. We use conditional generative adversarial networks to extract designated information from LIM. We consider a simple case with two populations of emission line galaxies; H$\rm\alpha$ emitting galaxies at $z = 1.3$ are confused with [OIII] emitters at $z = 2.0$ in a single observed waveband at 1.5 $\rm\mu$m. Our networks trained with 30,000 mock observation maps are able to extract the total intensity and the spatial distribution of H$\rm\alpha$ emitting galaxies at $z = 1.3$. The intensity peaks are successfully located with 74% precision. The precision increases to 91% when we combine the results of 5 networks. The mean intensity and the power spectrum are reconstructed with an accuracy of $\sim$10%. The extracted galaxy distributions at a wider range of redshift can be used for studies on cosmology and on galaxy formation and evolution.

Atmospheric escape is an important process that controls the long-term evolution of close-in planets. We perform radiation hydrodynamics simulations of photo-evaporation of exoplanets' atmospheres to study the effect of photoelectric heating by far-ultraviolet (FUV) radiation. Specifically, we consider a close-in hot Jupiter around a hot A-star. Hot main-sequence stars emit not only extreme ultraviolet radiation but also FUV radiation, and thus can drive strong atmospheric escape by photoelectric heating. We show that the planetary atmosphere escapes at a rate of $\dot{M}\sim10^{14}\,\mathrm{g/s}$ if the atmosphere contains heavy elements and dust grains with the level of ten percent of the solar metallicity. Close-in planets around hot stars can lose a significant fraction of the atmosphere during the long-term evolution. We also explore the metallicity dependence of the FUV driven escape. The mass-loss rate increases with increasing the atmosphere's metallicity because of the enhanced photoelectric heating, but the stellar FUV flux decreases with increasing stellar metallicity. The net dependence nearly cancels if we assume the same or similar metallicities for the planet and the host star. We derive an accurate estimate for the mass-loss rate as a function of FUV flux and metallicity, and of the planet's characteristics. The FUV driven atmospheric escape may be a key process to understand and explain the so-called sub-Jovian desert.

Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zel'dovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $z\sim0.1-0.2$. We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic ray protons at shocks in and around galaxy clusters. We find the acceleration efficiency must be below 5\% with a $2\sigma$ confidence level when the hydrostatic mass bias of clusters is assumed to be 30\%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the non-thermal cosmic-ray pressure in the ICM can introduce only a $\le 3\%$ level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array.

We study the formation of very metal-poor stars under protostellar radiative feedback effect. We use cosmological simulations to identify low-mass dark matter halos and star-forming gas clouds within them. We then follow protostar formation and the subsequent long-term mass accretion phase of over one million years using two-dimensional radiation-hydrodynamics simulations. We show that the critical physical process that sets the final mass is formation and expansion of a bipolar HII region. The process is similar to the formation of massive primordial stars, but radiation pressure exerted on dust grains also contributes to halting the accretion flow in the low-metallicity case. We find that the net feedback effect in the case with metallicity $Z = 10^{-2}~Z_{\odot}$ is stronger than in the case with $Z \sim 1~Z_{\odot}$. With decreasing metallicity, the radiation pressure effect becomes weaker, but photoionization heating of the circumstellar gas is more efficient owing to the reduced dust attenuation. In the case with $Z = 10^{-2}~Z_{\odot}$, the central star grows as massive as 200 solar-masses, similarly to the case of primordial star formation. We conclude that metal-poor stars with a few hundred solar masses can be formed by gas accretion despite the strong radiative feedback.

We study the emission properties of thermonuclear explosions in a helium white dwarf (WD) tidal disruption event (TDE). We consider a TDE where a 0.2 $M_{\odot}$ helium WD is disrupted by a $10^{2.5}\,M_{\odot}$ intermediate-mass black hole (IMBH). The helium WD is not only tidally disrupted but is also detonated by the tidal compression and by succeeding shocks. We focus on the emission powered by radioactive nuclei in the unbound TDE ejecta. We perform hydrodynamic simulations coupled with nuclear reactions, post-process detailed nucleosynthesis calculations, and then radiative transfer simulations. We thus derive multi-band light curves and spectra. The helium WD TDE shows rapid ($\Delta t_{1\mathrm{mag}}\simeq5\text{--}10$ days) and relatively faint ($L_{\mathrm{peak}}\simeq10^{42}\,\mathrm{erg}\,\mathrm{s}^{-1}$) light curves, because the ejecta mass and $^{56}$Ni mass are low ($0.12\,M_{\odot}$ and $0.03\,M_{\odot}$, respectively). The spectra show strong calcium and Fe-peak features and very weak silicon features, reflecting the peculiar elemental abundance. The key feature is the Doppler shift of the spectral lines up to $\simeq\pm12,000\,\mathrm{km}\,\mathrm{s}^{-1}$, depending on the viewing angle, due to the bulk motion of the ejecta. Our model matches well with two rapid and faint transients reported in Pursiainen et al. (2018). The particular model presented here does not match with observed SNe Iax, calcium-rich transients, or .Ia explosion candidates, either in the spectra or light curves. However, we expect a large variety of the observational signatures once a wide range of the WD/black hole masses and orbital parameters are considered. This study contributes to the search for WD TDEs with current and upcoming surveys, and to the identification of IMBHs as disrupters in the TDEs.

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

We study the enrichment and mixing of r-process elements in ultrafaint dwarf galaxies (UFDs). We assume that r-process elements are produced by neutron-star mergers (NSMs), and examine multiple models with different natal kick velocities and explosion energies. To this end, we perform cosmological simulations of galaxy formation to follow mixing of the dispersed r-process elements driven by star formation and the associated stellar feedback in progenitors of UFDs. We show that the observed europium abundance in Reticulum II is reproduced by our inner explosion model where a NSM is triggered at the centre of the galaxy, whereas the relatively low abundance in Tucana III is reproduced if a NSM occurs near the virial radius of the progenitor galaxy. The latter case is realised only if the neutron-star binary has a large natal kick velocity and travels over a long distance of a kilo-parsec before merger. In both the inner and outer explosion cases, it is necessary for the progenitor galaxy to sustain prolonged star formation over a few hundred million years after the NSM, so that the dispersed r-process elements are well mixed within the inter-stellar medium. Short-duration star formation results in inefficient mixing, and then a large variation is imprinted in the stellar europium abundances, which is inconsistent with the observations of Reticulum II and Tucana III.

Cross-correlation analysis of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intra-cluster medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ auto-power spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ auto-power spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} \%$ ($68\%$ C.L.), which is consistent with recent measurements by mass calibration techniques.

Galaxy imaging surveys enable us to map the cosmic matter density field through weak gravitational lensing analysis. The density reconstruction is compromised by a variety of noise originating from observational conditions, galaxy number density fluctuations, and intrinsic galaxy properties. We propose a deep-learning approach based on generative adversarial networks (GANs) to reduce the noise in the weak lensing map under realistic conditions. We perform image-to-image translation using conditional GANs in order to produce noiseless lensing maps using the first-year data of the Subaru Hyper Suprime-Cam (HSC) survey. We train the conditional GANs by using 30000 sets of mock HSC catalogs that directly incorporate observational effects. We show that an ensemble learning method with GANs can reproduce the one-point probability distribution function (PDF) of the lensing convergence map within a $0.5-1\sigma$ level. We use the reconstructed PDFs to estimate a cosmological parameter $S_{8} = \sigma_{8}\sqrt{\Omega_{\rm m0}/0.3}$, where $\Omega_{\rm m0}$ and $\sigma_{8}$ represent the mean and the scatter in the cosmic matter density. The reconstructed PDFs place tighter constraint, with the statistical uncertainty in $S_8$ reduced by a factor of $2$ compared to the noisy PDF. This is equivalent to increasing the survey area by $4$ without denoising by GANs. Finally, we apply our denoising method to the first-year HSC data, to place $2\sigma$-level cosmological constraints of $S_{8} < 0.777 \, ({\rm stat}) + 0.105 \, ({\rm sys})$ and $S_{8} < 0.633 \, ({\rm stat}) + 0.114 \, ({\rm sys})$ for the noisy and denoised data, respectively.

We study the gravitational instability and fragmentation of primordial protostellar discs by using high-resolution cosmological hydrodynamics simulations. We follow the formation and evolution of spiral arms in protostellar discs, examine the dynamical stability, and identify a physical mechanism of secondary protostar formation. We use linear perturbation theory based on the spiral-arm instability (SAI) analysis in our previous studies. We improve the analysis by incorporating the effects of finite thickness and shearing motion of arms, and derive the physical conditions for SAI in protostellar discs. Our analysis predicts accurately the stability and the onset of arm fragmentation that is determined by the balance between self-gravity and gas pressure plus the Coriolis force. Formation of secondary and multiple protostars in the discs is explained by the SAI, which is driven by self-gravity and thus can operate without rapid gas cooling. We can also predict the typical mass of the fragments, which is found to be in good agreement with the actual masses of secondary protostars formed in the simulation.

We study statistics of the 21-cm signal from the epoch of reionization. We propose to use [OIII] line emitting galaxies to cross-correlate with the 21-cm signal from $z = 7 - 10$. To this aim, we employ simulations of reionization obtained post-processing the high-resolution cosmological hydrodynamic simulation Massive Black-II with the three-dimensional (3D) radiative transfer code CRASH to follow the propagation of ionizing photons from a variety of sources. We show that, during the early phases of reionization, the 21-cm signal is positively correlated with the spatial distribution of the [OIII] emitters on large scales $(k < 1h~\rm cMpc^{-1})$. This positive correlation is generated by the temperature-galaxy correlation and it is a few times larger than when we assume that the heating is saturated. As the reionized regions expand, the correlation changes its sign to negative from $z = 10$ to 8. The signals at this epoch can be detected by combining the Square Kilometre Array (SKA) and a wide-field [OIII] emitter survey. We also calculate the cross-power spectrum with a 3D [OIII] intensity field, aiming at exploiting future intensity mapping observations. We conclude that high-redshift [OIII] line emitters can be used to probe the reionization process when the intergalactic medium is largely neutral.

Weak gravitational lensing is a powerful probe of the large-scale cosmic matter distribution. Wide-field galaxy surveys allow us to generate the so-called weak lensing maps, but actual observations suffer from noise due to imperfect measurement of galaxy shape distortions and to the limited number density of the source galaxies. In this paper, we explore a deep-learning approach to reduce the noise. We develop an image-to-image translation method with conditional adversarial networks (CANs), which learn efficient mapping from an input noisy weak lensing map to the underlying noise field. We train the CANs using $30000$ image pairs obtained from $1000$ ray-tracing simulations of weak gravitational lensing. We show that the trained CANs reproduce the true one-point probability distribution function (PDF) of the noiseless lensing map with a bias less than $1\sigma$ on average, where $\sigma$ is the statistical error. We perform a Fisher analysis to make forecast for cosmological parameter inference with the one-point lensing PDF. By our denoising method using CANs, the first derivative of the PDF with respect to the cosmic mean matter density and the amplitude of the primordial curvature perturbations becomes larger by $\sim50\%$. This allows us to improve the cosmological constraints by $\sim30-40\%$ with using observational data from ongoing and upcoming galaxy imaging surveys.

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.

We study influence by models of inter-stellar medium (ISM) on properties of galaxies in cosmological simulations. We examine three models widely used in previous studies. The ISM models impose different equations of state on dense gas. Using zoom-in simulations, we demonstrate that switching the ISM models can control formation of giant clumps in massive discs at redshifts $z\sim1$--$2$ while their initial conditions and the other settings such as stellar feedback are unchanged. Thus, not only feedback but ISM models can also be responsible for clumpy morphologies of simulated galaxies. We find, however, that changing the ISM models hardly affects global properties of galaxies, such as the total stellar and gas masses, star formation rate, metallicity and stellar angular momentum, irrespective of the significant difference of clumpiness; namely the ISM models only change clumpiness of discs. In addition, our approach provides a test to investigate impact by clump formation on the evolution of disc galaxies using the same initial conditions and feedback. We find that clump formation does not significantly alter the properties of galaxies and therefore could not be the causes of starburst or quenching.

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

One of the major goals of next-generation space-borne and ground-based telescopes is to detect and characterize the first galaxies that were in place in the first few hundred million years after the big bang. We study the spectral energy distribution (SED) of the first galaxies and discuss the prospects for detection and identification. We consider very young star-forming galaxies at $z=15$ and incorporate the contribution from pre-main-sequence (PMS) stars. Unlike in the present-day galaxies, primordial protostars are not embedded in dusty gas clouds, and hence the light from them can be visible at a wide range of wavelengths. We use MESA code to follow the PMS evolution and use the BT-Settl model to calculate the SED of individual PMS stars. We show that PMS stars contribute to boost the flux in the mid-infrared, and that the galaxy SED at very early evolutionary phases is overall redder than at later phases. The infrared flux contribution is comparable to that caused by emission lines powered by massive stars. We argue that the contribution from PMS stars is important for characterizing young galaxies in the early Universe and also for target selection with future deep galaxy surveys.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

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

We study the stellar dynamics of the first star clusters after intermediate-mass black holes (IMBHs) are formed via runaway stellar collisions. We use the outputs of cosmological simulations of Sakurai et al. (2017) to follow the star cluster evolution in a live dark matter (DM) halo. Mass segregation within a cluster promotes massive stars to be captured by the central IMBH occasionally, causing tidal disruption events (TDEs). We find that the TDE rate scales with the IMBH mass as $\dot{N}_{\rm TDE}\sim0.3\,{\rm Myr}^{-1}(M_{\rm IMBH}/1000\,{\rm M}_{\odot})^2$. The DM component affects the star cluster evolution by stripping stars from the outer part. When the DM density within the cluster increases, the velocity dispersion of the stars increases, and then the TDE rate decreases. By the TDEs, the central IMBHs grow to as massive as $700-2500\,{\rm M}_{\odot}$ in 15 million years. The IMBHs are possible seeds for the formation of supermassive BHs observed at $z\gtrsim 6-7$, if a large amount of gas is supplied through galaxy mergers and/or large-scale gas accretion, or they might remain as IMBHs from the early epochs to the present-day Universe.

We present a novel scenario for the formation of carbon-enhanced metal-poor (CEMP) stars. Carbon enhancement at low stellar metallicities is usually considered a consequence of faint or other exotic supernovae. An analytical estimate of cooling times in low-metallicity gas demonstrates a natural bias, which favours the formation of CEMP stars as a consequence of inhomogeneous metal mixing: carbon-rich gas has a shorter cooling time and can form stars prior to a potential nearby pocket of carbon-normal gas, in which star formation is then suppressed due to energetic photons from the carbon-enhanced protostars. We demonstrate that this scenario provides a natural formation mechanism for CEMP stars from carbon-normal supernovae, if inhomogeneous metal mixing provides carbonicity differences of at least one order of magnitude separated by >10pc. In our fiducial (optimistic) model, 8% (83%) of observed CEMP-no stars ([Ba/Fe]<0) can be explained by this formation channel. This new scenario may change our understanding of the first supernovae and thereby our concept of the first stars. Future 3D simulations are required to assess the likelihood of this mechanism to occur in typical high-redshift galaxies.

Fragmentation of spiral arms can drive the formation of giant clumps and induce intense star formation in disc galaxies. Based on the spiral-arm instability analysis of our Paper I, we present linear perturbation theory of dynamical instability of self-gravitating spiral arms of magnetised gas, focusing on the effect of toroidal magnetic fields. Spiral arms can be destabilised by the toroidal fields which cancel Coriolis force, i.e. magneto-Jeans instability. Our analysis can be applied to multi-component systems that consist of gas and stars. To test our analysis, we perform ideal magneto-hydrodynamics simulations of isolated disc galaxies and examine the simulation results. We find that our analysis can characterise dynamical instability leading arms to fragment and form clumps if magnetic fields are nearly toroidal. We propose that dimensionless growth rate of the most unstable perturbation, which is computed from our analysis, can be used to predict fragmentation of spiral arms within an orbital time-scale. Our analysis is applicable as long as magnetic fields are nearly toroidal. Using our analytic model, we estimate a typical mass of clumps forming from spiral-arm fragmentation to be consistent with observed giant clumps $\sim10^{7-8}~{\rm M_\odot}$. Furthermore, we find that, although the magnetic destabilisation can cause low-density spiral arms to fragment, the estimated mass of resultant clumps is almost independent from strength of magnetic fields since marginal instability occurs at long wavelengths which compensate the low densities of magnetically destabilised arms.

We perform a suite of 3D radiation hydrodynamics simulations of photoevaporation of molecular gas clumps illuminated by external massive stars. We study the fate of solar-mass clumps and derive their lifetimes with varying the gas metallicity over a range of $10^{-3} \,Z_\odot \leq Z \leq Z_\odot $. Our simulations incorporate radiation transfer of far ultraviolet (FUV) and extreme ultraviolet (EUV) photons, and follow atomic/molecular line cooling and dust-gas collisional cooling. Nonequilibrium chemistry is coupled with the radiative transfer and hydrodynamics in a self-consistent manner. We show that radiation-driven shocks compress gas clumps to have a volume that is set by the pressure-equilibrium with the hot ambient gas. Radiative cooling enables metal-rich clumps to condense and to have small surface areas, where photoevaporative flows are launched. For our fiducial set-up with an O-type star at a distance of 0.1 parsec, the resulting photoevaporation rate is as small as $\sim 10^{-5} M_{\odot}/{\rm yr}$ for metal-rich clumps, but is larger for metal-poor clumps that have larger surface areas. The clumps are continuously accelerated away from the radiation source by the so-called rocket effect, and can travel over $\sim$1 parsec within the lifetime. We also study photoevaporation of clumps in a photo-dissociation region. Photoelectric heating is inefficient for metal-poor clumps that contain a smaller amount of grains, and thus they survive for over $10^5$ years. We conclude that the gas metallicity strongly affects the clump lifetime and thus determines the strength of feedback from massive stars in star-forming regions.

Extremely metal poor (EMP) stars in the Milky Way inherited the chemical composition of the gas out of which they formed. They therefore carry the chemical fingerprint of the first stars in their spectral lines. It is commonly assumed that EMP stars form from gas that was enriched by only one progenitor supernova ('mono-enriched'). However, recent numerical simulations show that the first stars form in small clusters. Consequently, we expect several supernovae to contribute to the abundances of an EMP star ('multi-enriched'). We analyse seven recently observed EMP stars from the TOPoS survey by applying the divergence of the chemical displacement and find that J1035+0641 is mono-enriched ($p_{mono}=53\%$) and J1507+0051 is multi-enriched ($p_{mono}=4\%$). For the remaining five stars we can not make a distinct prediction ($p_{mono} \lesssim 50\%$) due to theoretical and observational uncertainties. Further observations in the near-UV will help to improve our diagnostic and therefore contribute to constrain the nature of the first stars.

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

We study gravitational collapse of low-metallicity gas clouds and the formation of protostars by three-dimensional hydrodynamic simulations. Grain growth, non-equilibrium chemistry, molecular cooling, and chemical heating are solved in a self-consistent manner for the first time. We employ the realistic initial conditions for the abundances of metal and dust, and the dust size distribution obtained from recent Population III supernova calculations. We also introduce the state-of-the-art particle splitting method based on the Voronoi tessellation and achieve an extremely high mass resolution of ~10^{-5} Msun (10 earth masses) in the central region. We follow the thermal evolution of several early clouds with various metallicities. We show that the condition for cloud fragmentation depends not only on the gas metallicity but also on the collapse timescale. In many cases, the cloud fragmentation is prevented by the chemical heating owing to molecular hydrogen formation even though dust cooling becomes effective. Meanwhile, in several cases, efficient OH and H2O cooling promotes the cloud elongation, and then cloud "filamentation" is driven by dust thermal emission as a precursor of eventual fragmentation. While the filament fragmentation is driven by rapid gas cooling with >10^{-5} Zsun, fragmentation occurs in a different manner by the self-gravity of a circumstellar disk with <10^{-5} Zsun. We use a semi-analytic model to estimate the number fraction of the clouds which undergo the filament fragmentation to be a several percents with 10^{-5}--10^{-4} Zsun. Overall, our simulations show a viable formation path of the recently discovered Galactic low-mass stars with extremely small metallicities.

We perform a suite of radiation hydrodynamics simulations of photoevaporating disks with varying the metallicity in a wide range of $10^{-3} \, Z_\odot \leq Z \leq 10^{0.5} \, Z_\odot $. We follow the disk evolution for over $\sim 5000$ years by solving hydrodynamics, radiative transfer, and non-equilibrium chemistry. Our chemistry model is updated from the first paper of this series by adding X-ray ionization and heating. We study the metallicity dependence of the disk photoevaporation rate and examine the importance of X-ray radiation. In the fiducial case with solar metallicity, including the X-ray effects does not significantly increase the photoevaporation rate when compared to the case with ultra-violet (UV) radiation only. At sub-solar metallicities in the range of $Z \gtrsim 10^{-1.5} \, Z_\odot $, the photoevaporation rate increases as metallicity decreases owing to the reduced opacity of the disk medium. The result is consistent with the observational trend that disk lifetimes are shorter in low metallicity environments. Contrastingly, the photoevaporation rate decreases at even lower metallicities of $Z \lesssim 10^{-1.5} \, Z_\odot $, because dust-gas collisional cooling remains efficient compared to far UV photoelectric heating whose efficiency depends on metallicity. The net cooling in the interior of the disk suppresses the photoevaporation. However, adding X-ray radiation significantly increases the photoevaporation rate, especially at $Z \sim 10^{-2}\, Z_\odot$. Although the X-ray radiation itself does not drive strong photoevaporative flows, X-rays penetrate deep into the neutral region in the disk, increase the ionization degree there, and reduce positive charges of grains. Consequently, the effect of photoelectric heating by far UV radiation is strengthened by the X-rays and enhances the disk photoevaporation.

Extremely metal-poor (EMP) stars in the Milky Way (MW) allow us to infer the properties of their progenitors by comparing their chemical composition to the metal yields of the first supernovae. This method is most powerful when applied to mono-enriched stars, i.e. stars that formed from gas that was enriched by only one previous supernova. We present a novel diagnostic to identify this subclass of EMP stars. We model the first generations of star formation semi-analytically, based on dark matter halo merger trees that yield MW-like halos at the present day. Radiative and chemical feedback are included self-consistently and we trace all elements up to zinc. Mono-enriched stars account for only $\sim 1\%$ of second generation stars in our fiducial model and we provide an analytical formula for this probability. We also present a novel analytical diagnostic to identify mono-enriched stars, based on the metal yields of the first supernovae. This new diagnostic allows us to derive our main results independently from the specific assumptions made regarding Pop III star formation, and we apply it to a set of observed EMP stars to demonstrate its strengths and limitations. Our results may provide selection criteria for current and future surveys and therefore contribute to a deeper understanding of EMP stars and their progenitors.

The extragalactic $\gamma$-ray background (EGB) arises from the accumulation of $\gamma$-ray emissions from resolved and unresolved extragalactic sources as well as diffuse processes. It is important to study the statistical properties of the EGB in the context of cosmological structure formation. Known astrophysical $\gamma$-ray sources such as blazars, star-forming galaxies, and radio galaxies are expected to trace the underlying cosmic matter density distribution. We explore the correlation of the EGB from Fermi-LAT data with the large-scale matter density distribution from the Subaru Hyper Suprime-Cam (HSC) SSP survey. We reconstruct an unbiased surface matter density distribution $\kappa$ at $z<1$ by applying weak-gravitational lensing analysis to the first-year HSC data. We then calculate the $\gamma - \kappa$ cross-correlation. Our measurements are consistent with a null detection, but a weak correlation is found at angular scales of 30-60 arcmin, especially when distant source galaxies at $z > 1$ are used for the lensing $\kappa$ reconstruction. The large-scale correlation suggests strong clustering of high-redshift $\gamma$-ray sources such as blazars. However, the inferred bias factor of $4-5$ is larger by about a factor of two than results from other clustering analyses. The final HSC data covering 1,400 squared degrees will play an essential role to determine accurately the blazar bias at $z > 0.5$.

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

We run a suite of hydrodynamics simulations of tidal disruption events (TDEs) of a white dwarf (WD) by a black hole (BH) with a wide range of WD/BH masses and orbital parameters. We implement nuclear reactions to study nucleosynthesis and its dynamical effect through release of nuclear energy. The released nuclear energy effectively increases the fraction of unbound ejecta. This effect is weaker for a heavy WD with 1.2 $\mathrm{M}_{\odot}$, because the specific orbital energy distribution of the debris is predominantly determined by the tidal force, rather than by the explosive reactions. The elemental yield of a TDE depends critically on the initial composition of a WD, while the BH mass and the orbital parameters also affect the total amount of synthesized elements. Tanikawa et al. (2017) find that simulations of WD-BH TDEs with low resolution suffer from spurious heating and inaccurate nuclear reaction results. In order to examine the validity of our calculations, we compare the amounts of the synthesized elements with the upper limits of them derived in a way where we can avoid uncertainties due to low resolution. The results are largely consistent, and thus support our findings. We find particular TDEs where early self-intersection of a WD occurs during the first pericentre passage, promoting formation of an accretion disc. We expect that relativistic jets and/or winds would form in these cases because accretion rates would be super-Eddington. The WD-BH TDEs result in a variety of events depending on the WD/BH mass and pericentre radius of the orbit.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

Protoplanetary disks are thought to have lifetimes of several million years in the solar neighborhood, but recent observations suggest that the disk lifetimes are shorter in a low metallicity environment. We perform a suite of radiation hydrodynamics simulations of photoevaporation of protoplanetary disks to study the disk structure and its long-term evolution of $\sim 10000$ years, and the metallicity dependence of mass-loss rate. Our simulations follow hydrodynamics, extreme and far ultra-violet radiative transfer, and non-equilibrium chemistry in a self-consistent manner. Dust grain temperatures are also calculated consistently by solving the radiative transfer of the stellar irradiation and grain (re-)emission. We vary the disk gas metallicity over a wide range of $10^{-4}~ Z_\odot \leq Z \leq 10 ~Z_\odot$. The photoevaporation rate is lower with higher metallicity in the range of $10^{-1} \,Z_\odot \lesssim Z \lesssim 10 \,Z_\odot$, because dust shielding effectively prevents far-ultra violet (FUV) photons from penetrating into and heating the dense regions of the disk. The photoevaporation rate sharply declines at even lower metallicities in $10^{-2} \,Z_\odot \lesssim Z \lesssim 10^{-1}\,Z_\odot$, because FUV photoelectric heating becomes less effective than dust-gas collisional cooling. The temperature in the neutral region decreases, and photoevaporative flows are excited only in an outer region of the disk. At $10^{-4}\,Z_\odot \leq Z \lesssim 10^{-2}\,Z_\odot$, HI photoionization heating acts as a dominant gas heating process and drives photoevaporative flows with roughly a constant rate. The typical disk lifetime is shorter at $Z=0.3~Z_\odot$ than at $Z = Z_\odot$, being consistent with recent observations of the extreme outer galaxy.

We perform a set of cosmological simulations of early structure formation with incorporating baryonic streaming motions. We present a case where a significantly elongated gas cloud with $\sim\!10^4\,$solar masses (${\rm M_\odot}$) is formed in a pre-galactic ($\sim\!10^7\,{\rm M_\odot}$) dark halo. The gas streaming into the halo compresses and heats the massive filamentary cloud to a temperature of $\sim\!10,000\,$Kelvin. The gas cloud cools rapidly by atomic hydrogen cooling, and then by molecular hydrogen cooling down to $\sim\!400\,$Kelvin. The rapid decrease of the temperature and hence of the Jeans mass triggers fragmentation of the filament to yield multiple gas clumps with a few hundred solar masses. We estimate the mass of the primordial star formed in each fragment by adopting an analytic model based on a large set of radiation hydrodynamics simulations of protostellar evolution. The resulting stellar masses are in the range of $\sim\!50$-$120\,{\rm M_\odot}$. The massive stars gravitationally attract each other and form a compact star cluster. We follow the dynamics of the star cluster using a hybrid $N$-body simulation. We show that massive star binaries are formed in a few million years through multi-body interactions at the cluster center. The eventual formation of the remnant black holes will leave a massive black hole binary, which can be a progenitor of strong gravitational wave sources similar to those recently detected by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO).

Recent detections of the cross-correlation of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) enable unique studies of cluster astrophysics and cosmology. In this work, we present constraints on the amplitude of the non-thermal pressure fraction in galaxy clusters, $\alpha_0$, and the amplitude of the matter power spectrum, $\sigma_8$, using measurements of the tSZ power spectrum from Planck, and the tSZ-WL cross-correlation from Planck and the Red Cluster Sequence Lensing Survey. We fit the data to a semi-analytic model with the covariance matrix using $N$-body simulations. We find that the tSZ power spectrum alone prefers $\sigma_8 \sim 0.85$ and a large fraction of non-thermal pressure ($\alpha_0 \sim 0.2$-$0.3$). The tSZ-WL cross-correlation on the other hand prefers a significantly lower $\sigma_8 \sim 0.6$, and low $\alpha_0 \sim 0.05$. We show that this tension can be mitigated by allowing for a steep slope in the stellar-mass-halo-mass relation, which would cause a reduction of the gas in low-mass halos. In such a model, the combined data prefer $\sigma_8 \sim 0.7$ and $\alpha_0 \sim 0.2$, consistent with predictions from hydrodynamical simulations.

Weak lensing three-point statistics are powerful probes of the structure of dark matter halos. We propose to use the correlation of the positions of galaxies with the shapes of background galaxy pairs, known as the halo-shear-shear correlation (HSSC), to measure the mean halo ellipticity and the abundance of subhalos in a statistical manner. We run high-resolution cosmological $N$-body simulations and use the outputs to measure the HSSC for galaxy halos and cluster halos. Non-spherical halos cause a characteristic azimuthal variation of the HSSC, and massive subhalos in the outer region near the virial radius contribute to $\sim10\%$ of the HSSC amplitude. Using the HSSC and its covariance estimated from our $N$-body simulations, we make forecast for constraining the internal structure of dark matter halos with future galaxy surveys. With 1000 galaxy groups with mass greater than $10^{13.5}\, h^{-1}M_{\odot}$, the average halo ellipticity can be measured with an accuracy of ten percent. A spherical, smooth mass distribution can be ruled out at a $\sim5\sigma$ significance level. The existence of subhalos whose masses are in 1-10 percent of the main halo mass can be detected with $\sim10^4$ galaxies/clusters. We conclude that the HSSC provides valuable information on the structure of dark halos and hence on the nature of dark matter.

Supernovae Type-Ia (SNeIa) play a significant role in exploring the history of the expansion of the Universe, since they are the best-known standard candles with which we can accurately measure the distance to the objects. Finding large samples of SNeIa and investigating their detailed characteristics have become an important issue in cosmology and astronomy. Existing methods relied on a photometric approach that first measures the luminance of supernova candidates precisely and then fits the results to a parametric function of temporal changes in luminance. However, it inevitably requires multi-epoch observations and complex luminance measurements. In this work, we present a novel method for classifying SNeIa simply from single-epoch observation images without any complex measurements, by effectively integrating the state-of-the-art computer vision methodology into the standard photometric approach. Our method first builds a convolutional neural network for estimating the luminance of supernovae from telescope images, and then constructs another neural network for the classification, where the estimated luminance and observation dates are used as features for classification. Both of the neural networks are integrated into a single deep neural network to classify SNeIa directly from observation images. Experimental results show the effectiveness of the proposed method and reveal classification performance comparable to existing photometric methods with multi-epoch observations.

Formation of supermassive stars (SMSs) with mass ~10^4 Msun is a promising pathway to seed the formation of supermassive black holes in the early universe. The so-called direct-collapse (DC) model postulates that such an SMS forms in a hot gas cloud irradiated by a nearby star-forming galaxy. We study the DC SMS formation in a fully cosmological context using three-dimensional radiation hydrodynamics simulations. We initialize our simulations using the outputs of the cosmological simulation of Chon et al. (2016), where two DC gas clouds are identified. The long-term evolution over a hundred thousand years is followed from the formation of embryo protostars through their growth to SMSs. We show that the strength of the tidal force by a nearby galaxy determines the multiplicity of the formed stars and affects the protostellar growth. In one case, where a collapsing cloud is significantly stretched by strong tidal force, multiple star-disk systems are formed via filament fragmentation. Small-scale fragmentation occurs in each circumstellar disk, and more than 10 stars with masses of a few times 10^3 Msun are finally formed. Interestingly, about a half of them are found as massive binary stars. In the other case, the gas cloud collapses nearly spherically under a relatively weak tidal field, and a single star-disk system is formed. Only a few SMSs with masses ~ 10^4 Msun are found already after evolution of a hundred thousand years, and the SMSs are expected to grow further by gas accretion and to leave massive blackholes at the end of their lives.

We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD Runge-Kutta integration. We apply our spatially fifth- and seventh-order schemes to a suite of simulations of collisionless self-gravitating systems and electrostatic plasma simulations, including linear and nonlinear Landau damping in one dimension and Vlasov--Poisson simulations in a six-dimensional phase space. The high-order schemes achieve a significantly improved accuracy in comparison with the third-order positive-flux-conserved scheme adopted in our previous study. With the semi-Lagrangian time integration, the computational cost of our high-order schemes does not significantly increase, but remains roughly the same as that of the third-order scheme. Vlasov--Poisson simulations on $128^3 \times 128^3$ mesh grids have been successfully performed on a massively parallel computer.

The origin of super-massive black holes in the early universe remains poorly understood.Gravitational collapse of a massive primordial gas cloud is a promising initial process,but theoretical studies have difficulty growing the black hole fast enough.We report numerical simulations of early black hole formation starting from realistic cosmological conditions.Supersonic gas motions left over from the Big Bang prevent early gas cloud formation until rapid gas condensation is triggered in a proto-galactic halo. A protostar is formed in the dense, turbulent gas cloud, and it grows by sporadic mass accretion until it acquires 34,000 solar masses.The massive star ends its life with a catastrophic collapse to leave a black hole -- a promising seed for the formation of a monstrous black hole.

Fragmentation of a spiral arm is thought to drive the formation of giant clumps in galaxies. Using linear perturbation analysis for self-gravitating spiral arms, we derive an instability parameter and define the conditions for clump formation. We extend our analysis to multi-component systems that consist of gas and stars in an external potential. We then perform numerical simulations of isolated disc galaxies with isothermal gas, and compare the results with the prediction of our analytic model. Our model describes accurately the evolution of the spiral arms in our simulations, even when spiral arms dynamically interact with one another. We show that most of the giant clumps formed in the simulated disc galaxies satisfy the instability condition. The clump masses predicted by our model are in agreement with the simulation results, but the growth time-scale of unstable perturbations is overestimated by a factor of a few. We also apply our instability analysis to derive scaling relations of clump properties. The expected scaling relation between the clump size, velocity dispersion, and circular velocity is slightly different from that given by the Toomre instability analysis, but neither is inconsistent with currently available observations. We argue that the spiral-arm instability is a viable formation mechanism of giant clumps in gas-rich disc galaxies.

We study the formation of massive black holes in the first star clusters. We first locate star-forming gas clouds in proto-galactic haloes of $\gtrsim \!10^7\,{\rm M}_{\odot}$ in cosmological hydrodynamics simulations and use them to generate the initial conditions for star clusters with masses of $\sim \!10^5\,{\rm M}_{\odot}$. We then perform a series of direct-tree hybrid $N$-body simulations to follow runaway stellar collisions in the dense star clusters. In all the cluster models except one, runaway collisions occur within a few million years, and the mass of the central, most massive star reaches $\sim \!400-1900\,{\rm M}_{\odot}$. Such very massive stars collapse to leave intermediate-mass black holes (IMBHs). The diversity of the final masses may be attributed to the differences in a few basic properties of the host haloes such as mass, central gas velocity dispersion, and mean gas density of the central core. Finally, we derive the IMBH mass to cluster mass ratios, and compare them with the observed black hole to bulge mass ratios in the present-day Universe.

We propose to use degree-scale angular clustering of fast radio bursts (FRBs) to identify their origin and the host galaxy population. We study the information content in autocorrelation of the angular positions and dispersion measures (DM) and in cross-correlation with galaxies. We show that the cross-correlation with Sloan Digital Sky Survey (SDSS) galaxies will place stringent constraints on the mean physical quantities associated with FRBs. If $\sim$10,000 FRBs are detected with $\lesssim \rm deg$ resolution in the SDSS field, the clustering analysis with the intrinsic DM scatter of $100\, {\rm pc}/{\rm cm}^3$ can constrain the global abundance of free electrons at $z\lt1$ and the large-scale bias of FRB host galaxies (the statistical relation between the distribution of host galaxies and cosmic matter density field) with fractional errors (with a $68\%$ confidence level) of $\sim10\%$ and $\sim20\%$, respectively. The mean near-source dispersion measure and the delay time distribution of FRB rates relative to the global star forming rate can be also determined by combining the clustering and the probability distribution function of DM. Our approach will be complementary to high-resolution ($\ll {\rm deg}$) event localization using e.g., VLA and VLBI for identifying the origin of FRBs and the source environment. We strongly encourage future observational programs such as CHIME, UTMOST, and HIRAX to survey FRBs in the SDSS field.

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

We examine the possibility that dark matter consists of charged massive particles (CHAMPs) in view of the cosmic microwave background (CMB) anisotropies. The evolution of cosmological perturbations of CHAMP with other components is followed in a self-consistent manner, without assuming that CHAMP and baryons are tightly coupled. We incorporate for the first time the "kinetic re-coupling" of the Coulomb scattering, which is characteristic of heavy CHAMPs. By a direct comparison of the predicted CMB temperature/polarization auto-correlations in CHAMP models and the observed spectra in the Planck mission, we show that CHAMPs leave sizable effects on CMB spectra if they are lighter than $10^{11}\,{\rm GeV}$. Our result can be applicable to any CHAMP as long as its lifetime is much longer than the cosmic time at the recombination ($\sim 4 \times 10^{5}\, {\rm yr}$). An application to millicharged particles is also discussed.

Dark Matter (DM) models providing possible alternative solutions to the small- scale crisis of standard cosmology are nowadays of growing interest. We consider DM interacting with light hidden fermions via well motivated fundamental operators showing the resultant matter power spectrum is suppressed on subgalactic scales within a plausible parameter region. Our basic description of the evolution of cosmological perturbations relies on a fully consistent first principles derivation of a perturbed Fokker-Planck type equation, generalizing existing literature. The cosmological perturbation of the Fokker-Planck equation is presented for the first time in two different gauges, where the results transform into each other according to the rules of gauge transformation. Furthermore, our focus lies on a derivation of a broadly applicable and easily computable collision term showing important phenomenological differences to other existing approximations. As one of the main results and concerning the small-scale crisis, we show the equal importance of vector and scalar boson mediated interactions between DM and light fermions.

Gravitational collapse of a massive primordial gas cloud is thought to be a promising path for the formation of supermassive blackholes in the early universe. We study conditions for the so-called direct collapse (DC) blackhole formation in a fully cosmological context. We combine a semianalytic model of early galaxy formation with halo merger trees constructed from dark matter $N$-body simulations. We locate a total of 68 possible DC sites in a volume of $20\;h^{-1}\;\mathrm{Mpc}$ on a side. We then perform hydrodynamics simulations for 42 selected halos to study in detail the evolution of the massive clouds within them. We find only two successful cases where the gas clouds rapidly collapse to form stars. In the other cases, gravitational collapse is prevented by the tidal force exerted by a nearby massive halo, which otherwise should serve as a radiation source necessary for DC. Ram pressure stripping disturbs the cloud approaching the source. In many cases, a DC halo and its nearby light source halo merge before the onset of cloud collapse. Only when the DC halo is assembled through major mergers, the gas density increases rapidly to trigger gravitational instability. Based on our cosmological simulations, we conclude that the event rate of DC is an order of magnitude smaller than reported in previous studies, although the absolute rate is still poorly constrained. It is necessary to follow the dynamical evolution of a DC cloud and its nearby halo(s) in order to determine the critical radiation flux for DC.

Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reionization on it. Towards that end, we perform a set of fully-coupled, radiation-hydrodynamics simulations from cosmological initial conditions, extending the mass resolution of previous work to the scale of minihalos. Pre-reionization structure is evolved until a redshift $z_i$ at which the ionizing radiation from external sources arrives to sweep an R-type ionization front supersonically across the volume in a few Myr, until it is trapped on the surfaces of minihalos and converted to D-type, after which the minihalo gas is removed by photoevaporative winds. Small-scale density structures during this time lead to a high ($>$10) clumping factor for ionized gas, which hugely boosts the recombination rate until the structures are disrupted by the hydrodynamic feedback after $\sim 10-100~\rm{Myr}$. For incoming stellar radiation with intensity $J_{21}$ in a $200~h^{-1}~\rm{kpc}$ box with the mean density contrast $\bar\delta$, the number of extra recombinations per H atom, on top of what is expected from homogeneously distributed gas, is given by $0.32[J_{21}]^{0.12}[(1+z_i)/11]^{-1.7}[1+\bar\delta]^{2.5}$. In models in which most of the volume is ionized toward the end of reionization, this can add more than one recombination per H atom to the ionizing photon budget to achieve reionization.

The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have less concentrated structure than fully-convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit $\gtrsim 10^5~M_\odot$ derived for the fully-convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with $\dot M \lesssim 0.1~M_\odot~{\rm yr}^{-1}$. With $\dot{M} \simeq 0.3 - 1~M_\odot~{\rm yr}^{-1}$, the star becomes GR-unstable during the helium-burning stage at $M \simeq 2 - 3.5~\times 10^5~M_\odot$. In an extreme case with $10~M_\odot~{\rm yr}^{-1}$, the star does not collapse until the mass reaches $\simeq 8.0\times 10^5~M_\odot$, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.

We derive constraints on dark matter (DM) annihilation cross section and decay lifetime from cross-correlation analyses of the data from Fermi-LAT and weak lensing surveys that cover a wide area of $\sim660$ squared degrees in total. We improve upon our previous analyses by using an updated extragalactic $\gamma$-ray background data reprocessed with the Fermi Pass 8 pipeline, and by using well-calibrated shape measurements of about twelve million galaxies in the Canada-France-Hawaii Lensing Survey (CFHTLenS) and Red-Cluster-Sequence Lensing Survey (RCSLenS). We generate a large set of full-sky mock catalogs from cosmological $N$-body simulations and use them to estimate statistical errors accurately. The measured cross correlation is consistent with null detection, which is then used to place strong cosmological constraints on annihilating and decaying DM. For leptophilic DM, the constraints are improved by a factor of $\sim100$ in the mass range of O(1) TeV when including contributions from secondary $\gamma$ rays due to the inverse-Compton upscattering of background photons. Annihilation cross-sections of $\langle \sigma v \rangle \sim 10^{-23}\, {\rm cm}^3/{\rm s}$ are excluded for TeV-scale DM depending on channel. Lifetimes of $\sim 10^{25}$ sec are also excluded for the decaying TeV-scale DM. Finally, we apply this analysis to wino DM and exclude the wino mass around 200 GeV. These constraints will be further tightened, and all the interesting wino DM parameter region can be tested, by using data from future wide-field cosmology surveys.

We present an application of machine-learning (ML) techniques to source selection in the optical transient survey data with Hyper Suprime-Cam (HSC) on the Subaru telescope. Our goal is to select real transient events accurately and in a timely manner out of a large number of false candidates, obtained with the standard difference-imaging method. We have developed the transient selector which is based on majority voting of three ML machines of AUC Boosting, Random Forest, and Deep Neural Network. We applied it to our observing runs of Subaru-HSC in 2015 May and August, and proved it to be efficient in selecting optical transients. The false positive rate was 1.0% at the true positive rate of 90% in the magnitude range of 22.0--25.0 mag for the former data. For the latter run, we successfully detected and reported ten candidates of supernovae within the same day as the observation. From these runs, we learned the following lessons: (1) the training using artificial objects is effective in filtering out false candidates, especially for faint objects, and (2) combination of ML by majority voting is advantageous.

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

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

Light gravitinos of mass $\lesssim \mathcal{O} (10)$ eV are of particular interest in cosmology, offering various baryogenesis scenarios without suffering from the cosmological gravitino problem. The gravitino may contribute considerably to the total matter content of the Universe and affect structure formation from early to present epochs. After the gravitinos decouple from other particles in the early Universe, they free-stream and consequently suppress density fluctuations of (sub-)galactic length scales. Observations of structure at the relevant length-scales can be used to infer or constrain the mass and the abundance of light gravitinos. We derive constraints on the light gravitino mass using the data of cosmic microwave background (CMB) lensing from Planck and of cosmic shear from the Canada France Hawaii Lensing Survey, combined with analyses of the primary CMB anisotropies and the signature of baryon acoustic oscillations in galaxy distributions. The obtained constraint on the gravitino mass is $m_{3/2} < 4.7$ eV (95% C.L.), which is substantially tighter than the previous constraint from clustering analysis of Ly-$\alpha$ forests.

We present coupled stellar evolution (SE) and 3D radiation-hydrodynamic (RHD) simulations of the evolution of primordial protostars, their immediate environment, and the dynamic accretion history under the influence of stellar ionizing and dissociating UV feedback. Our coupled SE-RHD calculations result in a wide diversity of final stellar masses covering 10 Msun $\lesssim M_* \lesssim$ 1000 Msun. The formation of very massive ($\gtrsim$ 250 Msun) stars is possible under weak UV feedback, whereas ordinary massive (a few x 10 Msun) stars form when UV feedback can efficiently halt the accretion. This may explain the peculiar abundance pattern of a Galactic metal-poor star recently reported by Aoki et al. (2014), possibly the observational signature of very massive precursor primordial stars. Weak UV feedback occurs in cases of variable accretion, in particular when repeated short accretion bursts temporarily exceed 0.01 Msun/yr, causing the protostar to inflate. In the bloated state, the protostar has low surface temperature and UV feedback is suppressed until the star eventually contracts, on a thermal adjustment timescale, to create an HII region. If the delay time between successive accretion bursts is sufficiently short, the protostar remains bloated for extended periods, initiating at most only short periods of UV feedback. Disk fragmentation does not necessarily reduce the final stellar mass. Quite the contrary, we find that disk fragmentation enhances episodic accretion as many fragments migrate inward and are accreted onto the star, thus allowing continued stellar mass growth under conditions of intermittent UV feedback. This trend becomes more prominent as we improve the resolution of our simulations. We argue that simulations with significantly higher resolution than reported previously are needed to derive accurate gas mass accretion rates onto primordial protostars.

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

Measurements of the cross-correlation between the extragalactic gamma-ray background (EGB) and large-scale structure provide a novel probe of dark matter on extragalactic scales. We focus on luminous red galaxies (LRGs) as optimal targets to search for the signal of dark matter annihilation. We measure the cross-correlation function of the EGB taken from the Fermi Large Area Telescope with the LRGs from the Sloan Digital Sky Survey. Statistical errors are calculated using a large set of realistic mock LRG catalogs. The amplitude of the measured cross-correlation is consistent with null detection. Based on an accurate theoretical model of the distribution of dark matter associated with LRGs, we exclude dark matter annihilation cross-sections over $\langle \sigma v\rangle =3\times10^{-25}-10^{-26}\, {\rm cm}^3 \,{\rm s}^{-1}$ for a 10 GeV dark matter. We further investigate systematic effects due to uncertainties in the Galactic gamma-ray foreground emission, which we find to be an order of magnitude smaller than the current statistical uncertainty. We also estimate the contamination from astrophysical sources in the LRGs by using known scaling relations between gamma-ray luminosity and star-formation rate, finding them to be negligibly small. Based on these results, we suggest that LRGs remain ideal targets for probing dark matter annihilation with future EGB measurement and galaxy surveys. Increasing the number of LRGs in upcoming galaxy surveys such as LSST would lead to big improvements of factors of several in sensitivity.

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

Ongoing and future wide-field galaxy surveys can be used to locate a number of clusters of galaxies with cosmic shear measurement alone. We study constraints on cosmological models using statistics of weak lensing selected galaxy clusters. We extend our previous theoretical framework to model the statistical properties of clusters in variants of cosmological models as well as in the standard $\Lambda$CDM model. Weak lensing selection of clusters does not rely on the conventional assumption such as the relation between luminosity and mass and/or hydrostatic equilibrium, but a number of observational effects compromise robust identification. We use a large set of realistic mock weak-lensing catalogs as well as analytic models to perform a Fisher analysis and make forecast for constraining two competing cosmological models, $w$CDM model and $f(R)$ model proposed by Hu & Sawicki, with our lensing statistics. We show that weak lensing selected clusters are excellent probe of cosmology when combined with cosmic shear power spectrum even in presence of galaxy shape noise and masked regions. With the information of weak lensing selected clusters, the precision of cosmological parameter estimate can be improved by a factor of $\sim1.6$ and $\sim8$ for $w$CDM model and $f(R)$ model, respectively. Hyper Suprime-Cam survey with sky coverage of $1250$ squared degrees can constrain the equation of state of dark energy $w_{0}$ with a level of $\Delta w_0 \sim0.1$. It can also constrain the additional scalar degree of freedom in $f(R)$ model with a level of $|f_{R0}| \sim5\times10^{-6}$, when constraints from cosmic microwave background measurements are incorporated. Future weak lensing surveys with sky coverage of $20,000$ squared degrees will place tighter constraints on $w_{0}$ and $|f_{R0}|$ even without cosmic microwave background measurements.

While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10~Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the standard scale-invariant PPS is modified at small length scales as $P(k) \sim k^{m_{\rm s}}$ with $m_{\rm s} > 1$. We run a series of cosmological hydrodynamic simulations to examine the dependence of the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models with $m_{\rm s} > 1$, star-forming gas clouds are formed at $z > 100$, when formation of hydrogen molecules is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars formed in such "hot" clouds grow very rapidly by accretion to become extremely massive stars that may leave massive black holes with a few hundred solar-masses at $z > 100$. The shape of the PPS critically affects the properties and the formation epoch of the first generation of stars. Future experiments of the CMB polarization and the spectrum distortion may provide important information on the nature of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum.

We explore a variety of statistics of clusters selected with cosmic shear measurement by utilizing both analytic models and large numerical simulations. We first develop a halo model to predict the abundance and the clustering of weak lensing selected clusters. Observational effects such as galaxy shape noise are included in our model. We then generate realistic mock weak lensing catalogs to test the accuracy of our analytic model. To this end, we perform full-sky ray-tracing simulations that allow us to have multiple realizations of a large continuous area. We model the masked regions on the sky using the actual positions of bright stars, and generate 200 mock weak lensing catalogs with sky coverage of ~1000 squared degrees. We show that our theoretical model agrees well with the ensemble average of statistics and their covariances calculated directly from the mock catalogues. With a typical selection threshold, ignoring shape noise correction causes overestimation of the clustering of weak lensing selected clusters with a level of about 10%, and shape noise correction boosts the cluster abundance by a factor of a few. We calculate the cross-covariances using the halo model with accounting for the effective reduction of the survey area due to masks. The covariance of the cosmic shear auto power spectrum is affected by the mode-coupling effect that originates from sky masking. Our model and the results can be readily used for cosmological analysis with ongoing and future weak lensing surveys.

We present a novel method for particle splitting in smoothed particle hydrodynamics simulations. Our method utilizes the Voronoi diagram for a given particle set to determine the position of fine daughter particles. We perform several test simulations to compare our method with a conventional splitting method in which the daughter particles are placed isotropically over the local smoothing length. We show that, with our method, the density deviation after splitting is reduced by a factor of about two compared with the conventional method. Splitting would smooth out the anisotropic density structure if the daughters are distributed isotropically, but our scheme allows the daughter particles to trace the original density distribution with length scales of the mean separation of their parent. We apply the particle splitting to simulations of the primordial gas cloud collapse. The thermal evolution is accurately followed to the hydrogen number density of 10^12 /cc. With the effective mass resolution of ~10^-4 Msun after the multi-step particle splitting, the protostellar disk structure is well resolved. We conclude that the method offers an efficient way to simulate the evolution of an interstellar gas and the formation of stars.

We study the impact of baryonic physics on cosmological parameter estimation with weak lensing surveys. We run a set of cosmological hydrodynamics simulations with different galaxy formation models. We then perform ray-tracing simulations through the total matter density field to generate 100 independent convergence maps of 25 ${\rm deg}^2$ field-of-view, and use them to examine the ability of the following three lensing statistics as cosmological probes; power spectrum, peak counts, and Minkowski functionals. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of 1400 ${\rm deg}^2$, these three statistics provide tight constraints on the matter density, density fluctuation amplitude, and dark energy equation of state, but parameter bias is induced by the baryonic processes such as gas cooling and stellar feedback. The bias can be significant when the statistical errors become small in future observations with a much larger survey area. We find the bias is induced in different directions in the parameter space depending on the statistics employed. While the two-point statistic, i.e. power spectrum, yields robust results against baryonic effects, the overall constraining power is weak compared with peak counts and Minkowski functionals. On the other hand, using one of peak counts or Minkowski functionals, or combined analysis with multiple statistics, results in biased parameter estimate. The bias can be as large as $1\sigma$ for HSC survey, and will be more significant for upcoming wider area surveys. We suggest to use an optimized combination so that the baryonic effects on parameter estimation are mitigated. Such `calibrated' combination can place stringent and robust constraints on cosmological parameters.

A promising formation channel of SMBHs at redshift 6 is the so-called DC model, which posits that a massive seed BH forms through gravitational collapse of a $\sim 10^5~M_\odot$ SMS. We study the evolution of such a SMS growing by rapid mass accretion. In particular, we examine the impact of time-dependent mass accretion of repeating burst and quiescent phases that are expected to occur with a self-gravitating circumstellar disk. We show that the stellar evolution with such episodic accretion differs qualitatively from that expected with a constant accretion rate, even if the mean accretion rate is the same. Unlike the case of constant mass accretion, whereby the star expands roughly following $R_* \simeq 2.6 \times 10^3 R_\odot (M_*/100~M_\odot)^{1/2}$, the protostar can substantially contract during the quiescent phases between accretion bursts. The stellar effective temperature and ionizing photon emissivity increase accordingly as the star contracts, which can cause strong ionizing feedback and halt the mass accretion onto the star. With a fixed duration of the quiescent phase $\Delta t_{\rm q}$, such contraction occurs in early evolutionary phases, i.e. for $M_* \lesssim 10^3~M_\odot$ with $\Delta t_{\rm q} \simeq 10^3$ yr. For later epochs and larger masses but the same $\Delta t_{\rm q}$, contraction is negligible even during quiescent phases. With larger quiescent times $\Delta t_{\rm q}$, however, the star continues to contract during quiescent phases even for the higher stellar masses. We show that such behavior is well understood by comparing the interval time and the thermal relaxation time for a bloated surface layer. We conclude that the UV radiative feedback becomes effective if the quiescent phase associated by the burst accretion is longer than $\sim 10^3$ yr, which is possible in an accretion disk forming in the direct collapse model.

Recent observations of the cosmic microwave background (CMB) anisotropies and the distribution of galaxies, galaxy clusters, and the Lyman Alpha forest have constrained the shape of the power spectrum of matter fluctuations on large scales k < few h/Mpc. We explore a new technique to constrain the matter power spectrum on smaller scales, assuming the dark matter is a Weakly Interacting Massive Particle (WIMP) that annihilates at early epochs. Energy released by dark matter annihilation can modify the spectrum of CMB temperature fluctuations and thus CMB experiments such as Planck have been able to constrain the quantity f <sigma v> /m < 1/88 picobarn c / GeV, where f is the fraction of energy absorbed by gas, <sigma v> is the annihilation rate assumed constant, and m is the particle mass. We assume the standard scale-invariant primordial matter power spectrum of P_prim(k) ~ k^{n_s} at large scales k < k_p, while we adopt the modified power law of P_prim(k) ~ k_p^{n_s} (k/k_p)^{m_s} at small scales. We then aim at deriving constraints on m_s. For m_s > n_s, the excess small-scale power results in a much larger number of nonlinear small mass halos, particularly at high redshifts. Dark matter annihilation in these halos releases sufficient energy to partially ionize the gas, and consequently modify the spectrum of CMB fluctuations. We show that the recent Planck data can already be used to constrain the power spectrum on small scales. For a simple model with an NFW profile with halo concentration parameter c_200 = 5 and f <sigma v> / m = 1/100 picobarn c / GeV, we can limit the mass variance sigma_{max} < 100 at the 95% confidence level, corresponding to a power law index m_s < 1.43 (1.63) for k_p = 100 (1000) h/Mpc. Our results are also relevant to theories that feature a running spectral index.

We perform a large set of cosmological simulations of early structure formation and follow the formation and evolution of 1540 star-forming gas clouds to derive the mass distribution of primordial stars. The star formation in our cosmological simulations is characterized by two distinct populations, the so-called Population III.1 stars and primordial stars formed under the influence of far ultraviolet (FUV) radiation (Population III.2D stars). In this work, we determine the stellar masses by using the dependences on the physical properties of star-forming cloud and/or the external photodissociating intensity from nearby primordial stars, which are derived from the results of two-dimensional radiation hydrodynamic simulations of protostellar feedback. The characteristic mass of the Pop III stars is found to be a few hundred solar masses at z ~ 25, and it gradually shifts to lower masses with decreasing redshift. At high redshifts z > 20, about half of the star-forming gas clouds are exposed to intense FUV radiation and thus give birth to massive Pop III.2D stars. However, the local FUV radiation by nearby Pop III stars becomes weaker at lower redshifts, when typical Pop III stars have smaller masses and the mean physical separation between the stars becomes large owing to cosmic expansion. Therefore, at z < 20, a large fraction of the primordial gas clouds host Pop III.1 stars. At z =< 15, the Pop III.1 stars are formed in relatively cool gas clouds due to efficient radiative cooling by H_2 and HD molecules; such stars have masses of a few x 10 Msun. Since the stellar evolution and the final fate are determined by the stellar mass, Pop III stars formed at different epochs play different roles in the early universe.

We investigate the condition for the formation of low-mass second-generation stars in the early universe. It has been proposed that gas cooling by dust thermal emission can trigger fragmentation of a low-metallicity star-forming gas cloud. In order to determine the critical condition in which dust cooling induces the formation of low-mass stars, we follow the thermal evolution of a collapsing cloud by a one-zone semi-analytic collapse model. Earlier studies assume the dust amount in the local universe, where all refractory elements are depleted onto grains, and/or assume the constant dust amount during gas collapse. In this paper, we employ the models of dust formation and destruction in early supernovae to derive the realistic dust compositions and size distributions for multiple species as the initial conditions of our collapse calculations. We also follow accretion of heavy elements in the gas phase onto dust grains, i.e., grain growth, during gas contraction. We find that grain growth well alters the fragmentation property of the clouds, and that this still does not approach to the value in the local universe. The critical conditions can be written by the gas metallicity Zcr and the initial depletion efficiency fdep,0 of gas-phase metal onto grains, or dust-to-metal mass ratio, as (Zcr/10^{-5.5} Zsun) = (fdep,0/0.18)^{-0.44} with small scatters in the range of Zcr = [0.06--3.2]x10^{-5} Zsun. We also show that the initial dust composition and size distribution are important to determine Zcr.

We investigate the origin of carbon-enhanced metal-poor (CEMP) stars starting from the recently discovered $\rm [Fe/H]<-7.1$ star SMSS J031300 (Keller et al. 2014). We show that the elemental abundances observed on the surface of SMSS J031300 can be well fit by the yields of faint, metal free, supernovae. Using properly calibrated faint supernova explosion models, we study, for the first time, the formation of dust grains in such carbon-rich, iron-poor supernova ejecta. Calculations are performed assuming both unmixed and uniformly mixed ejecta and taking into account the partial destruction by the supernova reverse shock. We find that, due to the paucity of refractory elements beside carbon, amorphous carbon is the only grain species to form, with carbon condensation efficiencies that range between (0.15-0.84), resulting in dust yields in the range (0.025-2.25)M$_{\odot}$. We follow the collapse and fragmentation of a star forming cloud enriched by the products of these faint supernova explosions and we explore the role played by fine structure line cooling and dust cooling. We show that even if grain growth during the collapse has a minor effect of the dust-to-gas ratio, due to C depletion into CO molecules at an early stage of the collapse, the formation of CEMP low-mass stars, such as SMSS J031300, could be triggered by dust cooling and fragmentation. A comparison between model predictions and observations of a sample of C-normal and C-rich metal-poor stars supports the idea that a single common pathway may be responsible for the formation of the first low-mass stars.

We present the first measurement of the cross-correlation of weak gravitational lensing and the extragalactic gamma-ray background emission using data from the Canada-France-Hawaii Lensing Survey and the Fermi Large Area Telescope. The cross-correlation is a powerful probe of signatures of dark matter annihilation, because both cosmic shear and gamma-ray emission originate directly from the same DM distribution in the universe, and it can be used to derive constraints on dark matter annihilation cross-section. We show that the measured lensing-gamma correlation is consistent with a null signal. Comparing the result to theoretical predictions, we exclude dark matter annihilation cross sections of <sigma v> =10^{-24}-10^{-25} cm^3 s^-1 for a 100 GeV dark matter. If dark matter halos exist down to the mass scale of 10^-6 M_sun, we are able to place constraints on the thermal cross sections <sigma v> ~ 3 x 10^{-26} cm^3 s^-1 for a 10 GeV dark matter annihilation into tau^{+} tau^{-}. Future gravitational lensing surveys will increase sensitivity to probe annihilation cross sections of <sigma v> ~ 5 x 10^{-26} cm^3 s^-1 even for a 100 GeV dark matter. Detailed modeling of the contributions from astrophysical sources to the cross correlation signal could further improve the constraints by ~ 40-70 %.

Redshift space distortion (RSD) observed in galaxy redshift surveys is a powerful tool to test gravity theories on cosmological scales, but the systematic uncertainties must carefully be examined for future surveys with large statistics. Here we employ various analytic models of RSD and estimate the systematic errors on measurements of the structure growth-rate parameter, $f\sigma_8$, induced by non-linear effects and the halo bias with respect to the dark matter distribution, by using halo catalogues from 40 realisations of $3.4 \times 10^8$ comoving $h^{-3}$Mpc$^3$ cosmological N-body simulations. We consider hypothetical redshift surveys at redshifts z=0.5, 1.35 and 2, and different minimum halo mass thresholds in the range of $5.0 \times 10^{11}$ -- $2.0 \times 10^{13} h^{-1} M_\odot$. We find that the systematic error of $f\sigma_8$ is greatly reduced to ~5 per cent level, when a recently proposed analytical formula of RSD that takes into account the higher-order coupling between the density and velocity fields is adopted, with a scale-dependent parametric bias model. Dependence of the systematic error on the halo mass, the redshift, and the maximum wavenumber used in the analysis is discussed. We also find that the Wilson-Hilferty transformation is useful to improve the accuracy of likelihood analysis when only a small number of modes are available in power spectrum measurements.

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

One of the milestones in the cosmic history is the formation of the first luminous objects and Hydrogen reionization. The standard theory of cosmic structure formation predicts that the first generation of stars were born about a few hundred million years after the Big Bang. The dark Universe was then lit up once again, and eventually filled with ultraviolet photons emitted from stars, galaxies, and quasars. The exact epoch of the cosmic reionization and the details of the process, even the dominant sources, are not known except the fact that the universe was reionized early on. Signatures of reionization are expected to be imprinted in the cosmic microwave background radiation, especially in its large scale polarization. Future CMB experiments, together with other probes such as 21 cm surveys, will provide rich information on the process of reionization. We review recent studies on reionization. The implications from available observations in a wide range of wavelengths are discussed. Results from state-of-the-art computer simulations are presented. Finally, we discuss prospects for exploring the first few hundred million years of the cosmic history.

The measurement of cosmic shear using weak gravitational lensing is a challenging task that involves a number of complicated procedures. We study in detail the systematic errors in the measurement of weak lensing Minkowski Functionals (MFs). Specifically, we focus on systematics associated with galaxy shape measurements, photometric redshift errors, and shear calibration correction. We first generate mock weak-lensing catalogs that directly incorporate the actual observational characteristics of the Canada-France-Hawaii Lensing Survey (CFHTLenS). We then perform a Fisher analysis using the large set of mock catalogs for various cosmological models. We find that the statistical error associated with the observational effects degrades the cosmological parameter constraints by a factor of a few. The Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of ~1400 deg2 will constrain the dark energy equation of the state parameter with an error of Delta w_0 ~ 0.25 by the lensing MFs alone, but biases induced by the systematics can be comparable to the 1-sigma error. We conclude that the lensing MFs are powerful statistics beyond the two-point statistics, only if well-calibrated measurement of both the redshifts and the shapes of source galaxies is performed. Finally, we analyze the CFHTLenS data to explore the ability of the MFs to break degeneracies between a few cosmological parameters. Using a combined analysis of the MFs and the shear correlation function, we derive the matter density Omega_m0= 0.256+0.054-0.046.

We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in each gas cloud and the subsequent protostellar evolution until the gas mass accretion onto the protostar is halted by stellar radiative feedback. The accretion rates differ significantly among the primordial gas clouds which largely determine the final stellar masses. For low accretion rates the growth of a protostar is self-regulated by radiative feedback effects and the final mass is limited to several tens of solar masses. At high accretion rates the protostar's outer envelope continues to expand and the effective surface temperature remains low; such protostars do not exert strong radiative feedback and can grow in excess to one hundred solar masses. The obtained wide mass range suggests that the first stars play a variety of roles in the early universe, by triggering both core-collapse supernovae and pair-instability supernovae as well as by leaving stellar mass black holes. We find certain correlations between the final stellar mass and the physical properties of the star-forming cloud. These correlations can be used to estimate the mass of the first star from the properties of the parent cloud or of the host halo, without following the detailed protostellar evolution.

Dust grains in low-metallicity star-forming regions may be responsible for the formation of the first low-mass stars. The minimal conditions to activate dust-induced fragmentation require the gas to be pre-enriched above a critical dust-to-gas mass ratio Dcr=[2.6--6.3]x10^-9 with the spread reflecting the dependence on the grain properties. The recently discovered Galactic halo star SDSS J102915+172927 has a stellar mass of 0.8 Msun and a metallicity of Z=4.5x10^-5 Zsun and represents an optimal candidate for the dust-induced low-mass star formation. Indeed, for the two most plausible Population III supernova progenitors, with 20 Msun and 35 Msun, the critical dust-to-gas mass ratio can be overcome provided that at least 0.4 Msun of dust condenses in the ejecta, allowing for moderate destruction by the reverse shock. Here we show that even if dust formation in the first supernovae is less efficient or strong dust destruction does occur, grain growth during the collapse of the parent gas cloud is sufficiently rapid to activate dust cooling and likely fragmentation into low-mass and long-lived stars. Silicates and magnetite grains can experience significant grain growth in the density range 10^9 /cc < nH<10^12 /cc by accreting gas-phase species (SiO, SiO2, and Fe) until their gas-phase abundance drops to zero, reaching condensation efficiencies =1. The corresponding increase in the dust-to-gas mass ratio allows dust-induced cooling and fragmentation to be activated at 10^12 /cc < nH < 10^14 /cc, before the collapsing cloud becomes optically thick to continuum radiation. We show that for all the initial conditions that apply to the parent cloud of SDSS J102915+172927, dust-driven fragmentation is able to account for the formation of the star.

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

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

Starting from the geometrical concept of a 4-dimensional de-Sitter configuration of spheres in Euclidean 3-space and modelling voids in the Universe as spheres, we show that a uniform distribution over this configuration space implies a power-law for the void number density which is consistent with results from the excursion set formalism and with data, for an intermediate range of void volumes. The scaling dimension of the large scale structure can be estimated as well. We also discuss the effect of restricting the survey geometry on the void statistics. This work is a new application of de-Sitter geometry to cosmology and also provides a new geometrical perspective on self-similarity in cosmology.

We explore the discovery potential of light gravitino mass m_{3/2} by combining future cosmology surveys and collider experiments. The former probe the imprint of light gravitinos in the cosmic matter density field, whereas the latter search signatures of a supersymmetry breaking mechanism. Free-streaming of light gravitinos suppresses the density fluctuations at galactic and sub-galactic length scales, where weak gravitational lensing can be used as a powerful probe. We perform numerical simulations of structure formation to quantify the effect. We then run realistic ray-tracing simulations of gravitational lensing to measure the cosmic shear in models with light gravitino. We forecast the possible reach of future wide-field surveys by Fisher analysis; the light gravitino mass can be determined with an accuracy of m_{3/2}=4\pm 1 eV by a combination of the Hyper Suprime Cam survey and cosmic microwave background anisotropy data obtained by Planck satellite. The corresponding accuracy to be obtained by the future Large Synoptic Survey Telescope is \delta m_{3/2}=0.6 eV. Data from experiments at Large Hadron Collider at 14 TeV will provide constraint at m_{3/2} \simeq 5 eV in the minimal framework of gauge-mediated supersymmetry breaking (GMSB) model. We conclude that a large class of the GMSB model can be tested by combining the cosmological observations and the collider experiments.

Supermassive stars (SMSs) forming via very rapid mass accretion (Mdot >~ 0.1 Msun/yr) could be precursors of supermassive black holes observed beyond redshift of about 6. Extending our previous work, we here study the evolution of primordial stars growing under such rapid mass accretion until the stellar mass reaches 10^{4 - 5} Msun. Our stellar evolution calculations show that a star becomes supermassive while passing through the "supergiant protostar" stage, whereby the star has a very bloated envelope and a contracting inner core. The stellar radius increases monotonically with the stellar mass, until =~ 100 AU for M_* >~ 10^4 Msun, after which the star begins to slowly contract. Because of the large radius the effective temperature is always less than 10^4 K during rapid accretion. The accreting material is thus almost completely transparent to the stellar radiation. Only for M_* >~ 10^5 Msun can stellar UV feedback operate and disturb the mass accretion flow. We also examine the pulsation stability of accreting SMSs, showing that the pulsation-driven mass loss does not prevent stellar mass growth. Observational signatures of bloated SMSs should be detectable with future observational facilities such as the James Webb Space Telescope. Our results predict that an inner core of the accreting SMS should suffer from the general relativistic instability soon after the stellar mass exceeds 10^5 Msun. An extremely massive black hole should form after the collapse of the inner core.

We extend the subhalo abundance matching method to assign galaxy color to subhalos. We separate a luminosity-binned subhalo sample into two groups by a secondary subhalo property which is presumed to be correlated with galaxy color. The two subsamples then represent red and blue galaxy populations. We explore two models for the secondary property, namely subhalo assembly time and local dark matter density around each subhalo. The model predictions for the galaxy two-point correlation functions are compared with the recent results from the Sloan Digital Sky Survey. We show that the observed color dependence of galaxy clustering can be reproduced well by our method applied to cosmological N-body simulations without baryonic processes. We then compare the model predictions for the color-dependent galaxy-mass cross correlation functions with the results from gravitational lensing observations. The comparison allows us to distinguish the models, and also to discuss what subhalo property should be used to assign color to subhalos accurately. We show that the extended abundance matching method using the local dark matter density as a color proxy provides an accurate description of the galaxy populations in the local universe. We also study impacts of scatter in the local dark matter density-color relations. Introducing scatter improves agreements of our model predictions with the observed red and blue galaxy clustering and is needed to explain observed correlation functions in finer color bins. Finally, we study red galaxy fraction profiles in galaxy group- and cluster-sized halos and find the red fraction profiles have a relatively strong dependence on our model parameters. We argue that the red fraction profiles can be an important observational clue, in addition to galaxy clustering and galaxy-galaxy lensing, to explore the galaxy-(sub)halo connections.

Sky masking is unavoidable in wide-field weak lensing observations. We study how masks affect the measurement of statistics of matter distribution probed by weak gravitational lensing. We first use 1000 cosmological ray-tracing simulations to examine in detail the impact of masked regions on the weak lensing Minkowski Functionals (MFs). We consider actual sky masks used for a Subaru Suprime-Cam imaging survey. The masks increase the variance of the convergence field and the expected values of the MFs are biased. The bias then affects the non-Gaussian signals induced by the gravitational growth of structure. We then explore how masks affect cosmological parameter estimation. We calculate the cumulative signal-to-noise ratio S/N for masked maps to study the information content of lensing MFs. We show that the degradation of S/N for masked maps is mainly determined by the effective survey area. We also perform simple chi^2 analysis to show the impact of lensing MF bias due to masked regions. Finally, we compare ray-tracing simulations with data from a Subaru 2 deg^2 survey in order to address if the observed lensing MFs are consistent with those of the standard cosmology. The resulting chi^2/n_dof = 29.6/30 for three combined MFs, obtained with the mask effects taken into account, suggests that the observational data are indeed consistent with the standard LambdaCDM model. We conclude that the lensing MFs are powerful probe of cosmology only if mask effects are correctly taken into account.

Observational identification of the first stars is one of the great challenges in the modern astronomy. Although a single first star is too faint to be detected, supernova explosions of the first stars can be bright enough. An important question is whether such supernovae can be detected in the limited observational area with realistic observational resources. We perform detailed simulations to study the detectability of superluminous supernovae (SLSNe) at high redshifts, using the observationally-calibrated star formation rate density and supernova occurrence rate. We show that a 100 deg^2 survey with the limiting magnitude of 26 mag in near-infrared wavelengths will be able to discover about 10 SLSNe at z > 10. If the survey is extended to 200 deg^2 with 27 mag depth, about 10 SLSNe can be discovered at z > 15. We emphasize that the observations at > 3 um are crucial to detect and select SLSNe at z > 10. Our simulations are also applied to the planned survey with Euclid, WFIRST, and WISH. These surveys will be able to detect about 1000, 400, and 3000 SLSNe up to z ~ 5, 7, and 12, respectively. We conclude that detection of SLSNe at z > 10 is in fact achievable in the near future.

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^\alpha$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $\alpha$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

In a low-metallicity gas, rapid cooling by dust thermal emission is considered to induce cloud fragmentation and play a vital role in the formation of low-mass stars (<~ 1 M_sun) in metal-poor environments. We investigate how the growth of dust grains through accretion of heavy elements in the gas phase onto grain surfaces alters the thermal evolution and fragmentation properties of a collapsing gas cloud. We calculate directly grain growth and dust emission cooling in a self-consistent manner. We show that MgSiO3 grains grow sufficiently at gas densities nH = 10^{10}, 10^{12}, and 10^{14} /cc for metallicities Z = 10^{-4}, 10^{-5}, and 10^{-6} Zsun, respectively, where the cooling of the collapsing gas cloud is enhanced. The condition for efficient dust cooling is insensitive to the initial condensation factor of pre-existing grains within the realistic range of 0.001--0.1, but sensitive to metallicity. The critical metallicity is Zcrit ~ 10^{-5.5} Zsun for the initial grain radius r_{MgSiO3,0} <~ 0.01 um and Zcrit ~ 10^{-4.5} Zsun for r_{MgSiO3,0} >~ 0.1 um. The formation of a recently discovered low-mass star with extremely low metallicity (<= 4.5x10^{-5} Zsun) could have been triggered by grain growth.

We study the formation of non-linear structures in Warm Dark Matter (WDM) models and in a Long-Lived Charged Massive Particle (CHAMP) model. CHAMPs with a decay lifetime of about 1 yr induce characteristic suppression in the matter power spectrum at subgalactic scales through acoustic oscillations in the thermal background. We explore structure formation in such a model. We also study three WDM models, where the dark matter particles are produced through the following mechanisms: i) WDM particles are produced in the thermal background and then kinematically decoupled; ii) WDM particles are fermions produced by the decay of thermal heavy bosons; and iii) WDM particles are produced by the decay of non-relativistic heavy particles. We show that the linear matter power spectra for the three models are all characterised by the comoving Jeans scale at the matter-radiation equality. Furthermore, we can also describe the linear matter power spectrum for the Long-Lived CHAMP model in terms of a suitably defined characteristic cut-off scale k_{Ch}, similarly to the WDM models. We perform large cosmological N-body simulations to study the non-linear growth of structures in these four models. We compare the halo mass functions, the subhalo mass functions, and the radial distributions of subhalos in simulated Milky Way-size halos. We study the models with k_{cut}=51, 410, 820 h/Mpc, and confirm that these statistics are indeed similar between the different WDM models and the Long-Lived CHAMP model. The result suggests that the cut-off scale k_{cut} not only characterises the linear power spectra but also can be used to predict the non-linear clustering properties. The radial distribution of subhalos in Milky Way-size halos is consistent with the observed distribution for k_{cut}~50-800 h/Mpc; such models resolve the so-called "missing satellite problem".

Understanding the gas content of high redshift halos is crucial for studying the formation of the first generation of galaxies and reionization. Recently, Tseliakhovich & Hirata showed that the relative "stream" velocity between the dark matter and baryons at the time of recombination - formally a second order effect, but an unusually large one - can influence the later structure formation history of the Universe. We quantify the effect of the stream velocity on the so-called "characteristic mass" - the minimum mass of a dark matter halo capable of retaining most of its baryons throughout its formation epoch - using three different high-resolution sets of cosmological simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between the dark matter and baryons. In order to understand this effect theoretically, we generalize the linear theory filtering mass to properly account for the difference between the dark matter and baryonic density fluctuation evolution induced by the stream velocity. We show that the new filtering mass provides an accurate estimate for the characteristic mass, while other theoretical ansatzes for the characteristic mass are substantially less precise.

We present the first 3D simulations to include the effects of dark matter annihilation feedback during the collapse of primordial mini-halos. We begin our simulations from cosmological initial conditions and account for dark matter annihilation in our treatment of the chemical and thermal evolution of the gas. The dark matter is modelled using an analytical density profile that responds to changes in the peak gas density. We find that the gas can collapse to high densities despite the additional energy input from the dark matter. No objects supported purely by dark matter annihilation heating are formed in our simulations. However, we find that the dark matter annihilation heating has a large effect on the evolution of the gas following the formation of the first protostar. Previous simulations without dark matter annihilation found that protostellar discs around Population III stars rapidly fragmented, forming multiple protostars that underwent mergers or ejections. When dark matter annihilation is included, however, these discs become stable to radii of 1000 AU or more. In the cases where fragmentation does occur, it is a wide binary that is formed.

We study the thermal evolution of primordial star-forming gas clouds using three-dimensional cosmological simulations. We critically examine how assumptions and approximations made in calculating radiative cooling rates affect the dynamics of the collapsing gas clouds. We consider two important molecular hydrogen cooling processes that operate in a dense primordial gas; H_2 line cooling and continuum cooling by H_2 collision-induced emission. To calculate the optically thick cooling rates, we follow the Sobolev method for the former, whereas we perform ray-tracing for the latter. We also run the same set of simulations using simplified fitting functions for the net cooling rates. We compare the simulation results in detail. We show that the time- and direction-dependence of hydrodynamic quantities such as gas temperature and local velocity gradients significantly affects the optically thick cooling rates. Gravitational collapse of the cloud core is accelerated when the cooling rates are calculated by using the fitting functions. The structure and evolution of the central pre-stellar disk are also affected. We conclude that physically motivated implementations of radiative transfer are necessary to follow accurately the thermal and chemical evolution of a primordial gas to high densities.

We study the imprint of dark matter (DM) annihilation on the global 21 cm signal from the Dark Ages to Cosmic Reionization. Motivated by recent observations, we focus on three DM candidates: (i) a 10 GeV Bino-like neutralino (ii) a 200 GeV Wino and (iii) a 1 TeV heavier particle annihilating into leptons. For each DM candidate we assume two values for the thermally averaged annihilation cross section \sigma v, the standard thermal value \sigma v_th = 3 x 10^-26 cm^3 s^-1 and the maximum value allowed by WMAP7 data, \sigma v_max. We include the enhancement of DM annihilations due to collapsed structures, detailed estimates of energy deposition into the intergalactic medium (IGM), as well realistic prescriptions for astrophysical sources of UV and X-ray radiation. In these models, the additional heat input from DM annihilation suppresses the mean 21cm brightness temperature offset by \delta Tb ~ few-100 mK. In particular, the very deep \delta Tb ~ -150 mK absorption feature at ~20 < z < 25 predicted by popular models of the first galaxies is considerably reduced or totally erased by some of the considered DM candidates. Such an enhancement in IGM heating could come from either DM annihilations or a stronger-than-expected astrophysical component (i.e. abundant early X-ray sources). However, we find that the two signatures are not degenerate, since the DM heating is dominated by halos several orders of magnitude smaller than those hosting galaxies, whose fractional abundance evolves more slowly resulting in a smaller gradient: d \delta Tb/d \nu < 4 mK/MHz in the range \nu ~60-80 MHz. The detection of such signals by future radio telescopes would be clear evidence of DM energy injection at high-redshifts.

We study the formation of low-mass and extremely metal-poor stars in the early universe. Our study is motivated by the recent discovery of a low-mass (M < 0.8 Msun) and extremely metal-poor (Z <= 4.5 x 10^{-5} Zsun) star in the Galactic halo by Caffau et al. We propose a model that early supernova (SN) explosions trigger the formation of low-mass stars via shell fragmentation. We first perform one-dimensional hydrodynamic simulations of the evolution of an early SN remnant. We show that the shocked shell undergoes efficient radiative cooling and then becomes gravitationally unstable to fragment and collapse in about ten million years. We then follow the thermal evolution of the collapsing fragments using a one-zone code. Our one-zone calculation treats chemistry and radiative cooling self-consistently in low-metallicity gas. The collapsing gas cloud evolves roughly isothermally, until it cools rapidly by dust continuum emission at the density 10^{13}-10^{14} /cc. The cloud core then becomes thermally and gravitationally unstable and fragments. We argue that early SNe can trigger the formation of low-mass stars in the extremely metal-poor environment as Caffau et al. discovered recently.

The first stars in the universe ionized the ambient primordial gas through various feedback processes. "Second-generation" primordial stars potentially form from this disturbed gas after its recombination. In this Letter, we study the late formation stage of such second-generation stars, where a large amount of gas accretes onto the protostar and the final stellar mass is determined when the accretion terminates. We directly compute the complex interplay between the accretion flow and stellar ultraviolet (UV) radiation, performing radiation-hydrodynamic simulations coupled with stellar evolution calculations. Because of more efficient H2 and HD cooling in the pre-stellar stage, the accretion rates onto the star are ten times lower than in the case of the formation of the first stars. The lower accretion rates and envelope density result in the occurrence of an expanding bipolar HII region at a lower protostellar mass M_* \simeq 10Msun, which blows out the circumstellar material, thereby quenching the mass supply from the envelope to the accretion disk. At the same time the disk loses mass due to photoevaporation by the growing star. In our fiducial case the stellar UV feedback terminates mass accretion onto the star at M_* \simeq 17Msun. Although the derived masses of the second-generation primordial stars are systematically lower than those of the first generation, the difference is within a factor of only a few. Our results suggest a new scenario, whereby the majority of the primordial stars are born as massive stars with tens of solar masses, regardless of their generations.

We study the cosmological information contained in the Minkowski Functionals (MFs) of weak gravitational lensing convergence maps. We show that the MFs provide strong constraints on the local type primordial non-Gaussianity parameter f_NL. We run a set of cosmological N-body simulations and perform ray-tracing simulations of weak lensing, to generate 100 independent convergence maps of 25 deg^2 field-of-view for f_NL = -100, 0 and 100. We perform a Fisher analysis to study the degeneracy among other cosmological parameters such as the dark energy equation of state parameter w and the fluctuation amplitude sigma_8. We use fully nonlinear covariance matrices evaluated from 1000 ray-tracing simulations. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam survey with the proposed survey area of 1500 deg^2, the primordial non-Gaussianity can be constrained with a level of f_NL ~ 80 and w ~ 0.036 by weak lensing MFs. If simply scaled by the effective survey area, a 20000 deg^2 lensing survey using Large Synoptic Survey Telescope will give constraints of f_NL ~ 25 and w ~ 0.013. We show that these constraints can be further improved by a tomographic method using source galaxies in multiple redshift bins.

We show model light curves of superluminous supernova 2006gy on the assumption that the supernova is powered by the collision of supernova ejecta and its dense circumstellar medium. The initial conditions are constructed based on the shock breakout condition, assuming that the circumstellar medium is dense enough to cause the shock breakout within it. We perform a set of numerical light curve calculations by using a one-dimensional multigroup radiation hydrodynamics code STELLA. We succeeded in reproducing the overall features of the early light curve of SN 2006gy with the circumstellar medium whose mass is about 15 Msun (the average mass-loss rate ~ 0.1 Msun/yr). Thus, the progenitor of SN 2006gy is likely a very massive star. The density profile of the circumstellar medium is not well constrained by the light curve modeling only, but our modeling disfavors the circumstellar medium formed by steady mass loss. The ejecta mass is estimated to be comparable to or less than 15 Msun and the explosion energy is expected to be more than 4e51 erg. No 56Ni is required to explain the early light curve. We find that the multidimensional effect, e.g., the Rayleigh-Taylor instability, which is expected to take place in the cool dense shell between the supernova ejecta and the dense circumstellar medium, is important in understanding supernovae powered by the shock interaction. We also show the evolution of the optical and near-infrared model light curves of high-redshift superluminous supernovae. They can be potentially used to identify SN 2006gy-like superluminous supernovae in the future optical and near-infrared transient surveys.

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

We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations including the stability of King spheres, the gravitational instability and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64^6 grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

We investigate the formation of the first stars at the end of the cosmic dark ages with a suite of three-dimensional, moving mesh simulations that directly resolve the collapse of the gas beyond the formation of the first protostar at the centre of a dark matter minihalo. The simulations cover more than 25 orders of magnitude in density and have a maximum spatial resolution of 0.05 R_sun, which extends well below the radius of individual protostars and captures their interaction with the surrounding gas. In analogy to previous studies that employed sink particles, we find that the Keplerian disc around the primary protostar fragments into a number of secondary protostars, which is facilitated by H2 collisional dissociation cooling and collision-induced emission. The further evolution of the protostellar system is characterized by strong gravitational torques that transfer angular momentum between the secondary protostars formed in the disc and the surrounding gas. This leads to the migration of about half of the secondary protostars to the centre of the cloud in a free-fall time, where they merge with the primary protostar and enhance its growth to about five times the mass of the second most massive protostar. By the same token, a fraction of the protostars obtain angular momentum from other protostars via N-body interactions and migrate to higher orbits. On average, only every third protostar survives until the end of the simulation. However, the number of protostars present at any given time increases monotonically, suggesting that the system will continue to grow beyond the limited period of time simulated here.

We develop an analytic halo model for the distribution of dust around galaxies. The model results are compared with the observed surface dust density profile measured through reddening of background quasars in the Sloan Digital Sky Survey (SDSS) reported by Menard et al.(2010). We assume that the dust distribution around a galaxy is described by a simple power law, similarly to the mass distribution, but with a sharp cut-off at $\alpha R_{\rm vir}$ where $R_{\rm vir}$ is the galaxy's virial radius and $\alpha$ is a model parameter. Our model reproduces the observed dust distribution profile very well over a wide range of radial distance of $10 - 10^{4} h^{-1}$kpc. For the characteristic galaxy halo mass of $2\times 10^{12} h^{-1}M_{\odot}$ estimated for the SDSS galaxies, the best fit model is obtained if $\alpha$ is greater than unity, which suggests that dust is distributed to over a few hundred kilo-parsecs from the galaxies. The observed large-scale dust distribution profile is reproduced if we assume the total amount of dust is equal to that estimated from the integrated stellar evolution over the cosmic time.

We study the lensing convergence power spectrum and its covariance for a standard LCDM cosmology. We run 400 cosmological N-body simulations and use the outputs to perform a total of 1000 independent ray-tracing simulations. We compare the simulation results with analytic model predictions. The semi-analytic model based on Smith et al.(2003) fitting formula underestimates the convergence power by ~30% at arc-minute angular scales. For the convergence power spectrum covariance, the halo model reproduces the simulation results remarkably well over a wide range of angular scales and source redshifts. The dominant contribution at small angular scales comes from the sample variance due to the number fluctuations of halos in a finite survey volume. The signal-to-noise ratio for the convergence power spectrum is degraded by the non-Gaussian covariances by up to a factor 5 for a weak lensing survey to z_s ~1. The probability distribution of the convergence power spectrum estimators, among the realizations, is well approximated by a chi-square distribution with broadened variance given by the non-Gaussian covariance, but has a larger positive tail. The skewness and kurtosis have non-negligible values especially for a shallow survey. We argue that a prior knowledge on the full distribution may be needed to obtain an unbiased estimate on the ensemble averaged band power at each angular scale from a finite volume survey.

Observations of high-redshift supernovae (SNe) open a novel opportunity to study the massive star population in the early Universe. We study the detectability of superluminous SNe with upcoming optical and near-infrared (NIR) surveys. Our calculations are based on the cosmic star formation history, the SN occurence rate, the characteristic colour and the light curve of the SNe that are all calibrated by available observations. We show that 15-150 SNe up to z ~ 4 will be discovered by the proposed Subaru/Hyper Suprime-Cam deep survey: 30 deg^2 survey with 24.5 AB mag depth in z-band for 3 months. With its ultra-deep layer (3.5 deg^2 with 25.6 AB mag depth in z-band for 4 months), the highest redshift can be extended to z ~ 5. We further explore the detectability by upcoming NIR survey utilizing future satellites such as Euclid, WFIRST, and WISH. The wide-field NIR surveys are very efficient to detect high-redshift SNe. With a hypothetical deep NIR survey for 100 deg^2 with 26 AB mag depth at 1-4 um, at least ~ 50 SNe will be discovered at z>3 in half a year. The number of the detected SNe can place a strong constraint on the stellar initial mass function or its slope especially at the high-mass end. Superluminous SNe at high redshifts can be distinguished from other types of SNe by the long time-scale of their light curves in the observer's frame, the optical colours redder than other core-collapse SNe and the NIR colours redder than any other types of SNe.

It has been recently shown that the relative velocity between the dark matter and the baryons (vbc) at the time of recombination can affect the structure formation in the early universe (Tseliakhovich & Hirata 2010). We statistically quantify this effect using large cosmological simulations. We use three different high resolution sets of simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between dark matter and baryons. We show that the total number density of halos is suppressed by ~ 20% at z = 25 for vbc = 1{\sigma}(vbc), where {\sigma}(vbc) is the variance of the relative velocity, while for vbc = 3.4{\sigma}(vbc) the relative suppression at the same redshift reaches 50%, remaining at or above the 30% level all the way to z = 11. We also find high abundance of "empty halos", i.e., halos that have gas fraction below half of the cosmic mean baryonic fraction fb. Specifically we find that for vbc = 1{\sigma}(vbc) all halos below 10^5M\odot are empty at z \geq 19. The high abundance of empty halos results in significant delay in the formation of gas rich mini-halos and the first galaxies.

We explore the mass distribution of material associated with galaxies from the observation of gravitational weak lensing for the galaxy mass correlation function with the aid of $N$-body simulations of dark matter. The latter is employed to unfold various contributions that contribute to the integrated line of sight mass density. We conclude that galaxies have no definite edges of the matter distribution, extending to the middle to neighbouring galaxies with the density profile roughly $r^{-2.4}$ beyond the virial radius. The mass distributed beyond the virial radius (gravitationally bound radius) explains the gap seen in the mass density estimates, the global value $\Omega_m\sim 0.27$ and typically $\Omega_{\rm gal} \sim 0.15$ from the luminosity density multiplied by the mass to light ratio. We suggest to use a physical method of gravitational lensing to characterise galaxy samples rather than characterise them with photometric means.

The first stars fundamentally transformed the early universe by emitting the first light and by producing the first heavy elements. These effects were predetermined by the mass distribution of the first stars, which is thought to have been fixed by a complex interplay of gas accretion and protostellar radiation. We performed radiation-hydrodynamics simulations that followed the growth of a primordial protostar through to the early stages as a star with thermo-nuclear burning. The circumstellar accretion disk was evaporated by ultraviolet radiation from the star when its mass was 43 times that of the Sun. Such massive primordial stars, in contrast to the often postulated extremely massive stars, may help explain the fact that there are no signatures of the pair-instability supernovae in abundance patterns of metal-poor stars in our galaxy.

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation (Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to $E_{\rm iso} \gtrsim 10^{55-57}$ ergs, providing a unique probe of the high-redshift Universe. We consider both the so-called Pop III.1 stars (primordial) and Pop III.2 stars (primordial but affected by radiation from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical evolution of the intergalactic medium. We show that Pop III.2 GRBs occur more than 100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained by the current radio transient searches. We expect $\sim10-10^4$ radio afterglows above $\sim$ 0.3 mJy on the sky with $\sim 1$ year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while we expect to observe maximum of $N < 20$ GRBs per year integrated over at $z>6$ for Pop III.2 and $N < 0.08$ per year integrated over at $z>10$ for Pop III.1 with EXIST, and $N < 0.2$ for Pop III.2 GRBs per year integrated over at $z > 6$ with \textit{Swift}.

We run very large cosmological N-body hydrodynamical simulations in order to study statistically the baryon fractions in early dark matter halos. We critically examine how differences in the initial conditions affect the gas fraction in the redshift range z = 11--21. We test three different linear power spectra for the initial conditions: (1) A complete heating model, which is our fiducial model; this model follows the evolution of overdensities correctly, according to Naoz & Barkana (2005), in particular including the spatial variation of the speed of sound of the gas due to Compton heating from the CMB. (2) An equal-{\delta} model, which assumes that the initial baryon fluctuations are equal to those of the dark matter, while conserving sigma8 of the total matter. (3) A mean cs model, which assumes a uniform speed of sound of the gas. The latter two models are often used in the literature. We calculate the baryon fractions for a large sample of halos in our simulations. Our fiducial model implies that before reionization and significant stellar heating took place, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was ~ 3 * 10^4 Msun. However, the alternative models yield a wrong (higher by about 50%) minimum mass, since the system retains a memory of the initial conditions. We also demonstrate this using the "filtering mass" from linear theory, which accurately describes the evolution of the baryon fraction throughout the simulated redshift range.

Primordial stars formed in the early universe are thought to be hosted by compact dark matter (DM) halos. If DM consists of Weakly Interacting Massive Particles (WIMPs), such stars may be powered by DM annihilation during the early phases of their evolutions. We study the pre-main sequence evolutions of the primordial star using a detailed stellar evolution code under the assumption that the annihilation of adiabatically contracted WIMPs DM within the star provides a sufficient energy to sustain the stellar equilibrium. We follow the evolution of accreting stars using several gas mass accretion rates derived from cosmological simulations. We show that the stellar mass becomes very large, up to 900 - 1000 M_sun when the star reaches the main-sequence phase for a reasonable set of model parameters such as DM particle mass and the annihilation cross section. During the dark star phase, the star expands over a thousand solar-radii, while the surface temperature remains below 10^4 K. The energy generated by nuclear reactions is not dominant during this phase. We also study models with different gas mass accretion rates and the DM particle masses. All our models for different DM particle masses pass the dark star phase. The final mass of the dark stars is essentially unchanged for DM mass of m_DM <= 10 GeV. Gravitational collapse of the massive dark stars will leave massive black holes with mass as large as 1000 M_sun in the early universe.

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

We review our current understanding of how the first galaxies formed at the end of the cosmic dark ages, a few 100 million years after the Big Bang. Modern large telescopes discovered galaxies at redshifts greater than seven, whereas theoretical studies have just reached the degree of sophistication necessary to make meaningful predictions. A crucial ingredient is the feedback exerted by the first generation of stars, through UV radiation, supernova blast waves, and chemical enrichment. The key goal is to derive the signature of the first galaxies to be observed with upcoming or planned next-generation facilities, such as the James Webb Space Telescope or Atacama Large Millimeter Array. From the observational side, ongoing deep-field searches for very high-redshift galaxies begin to provide us with empirical constraints on the nature of the first galaxies.

We study how primordial non-Gaussianity affects the pairwise velocity probability density function (PDF) using an analytical model and cosmological N-body simulations. We adopt the local type non-Gaussian models characterized by f_{nl}, and examine both the linear velocity difference PDF and the linear pairwise velocity PDF. We show explicitly how f_{nl} induces correlations between originally independent velocities along the parallel and the perpendicular to the line of separation directions. We compare the model results with measurements from N-body simulations of the non-Gaussian models. Linear theory fails to predict the PDF in the f_{nl} models. Therefore we develop an analytic model based on the Zeldovich approximation to describe the evolution of the velocity PDF. Our analytical model and simulation results show remarkably good agreement in both the parallel and the perpendicular directions for the PDF profiles, as well as the change in the PDF due to primordial non-Gaussianity. The agreement is particularly good for relatively small separations (< 10 h^{-1} Mpc). The inclusion of the evolution of the velocity PDF is important to obtain a good description on the signature of primordial non-Gaussianity in the PDF. Our model provides the foundation to constrain f_{nl} using the peculiar velocity in future surveys.

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

We study the sample variance of the matter power spectrum for the standard Lambda Cold Dark Matter universe. We use a total of 5000 cosmological N-body simulations to study in detail the distribution of best-fit cosmological parameters and the baryon acoustic peak positions. The obtained distribution is compared with the results from the Fisher matrix analysis with and without including non-Gaussian errors. For the Fisher matrix analysis, we compute the derivatives of the matter power spectrum with respect to cosmological parameters using directly full nonlinear simulations. We show that the non-Gaussian errors increase the unmarginalized errors by up to a factor 5 for k_{max}=0.4h/Mpc if there is only one free parameter provided other parameters are well determined by external information. On the other hand, for multi-parameter fitting, the impact of the non-Gaussian errors is significantly mitigated due to severe parameter degeneracies in the power spectrum. The distribution of the acoustic peak positions is well described by a Gaussian distribution, with its width being consistent with the statistical interval predicted from the Fisher matrix. We also examine systematic bias in the best-fit parameter due to the non-Gaussian errors. The bias is found to be smaller than the 1 sigma statistical error for both the cosmological parameters and the acoustic scale positions.

In this paper we extend our previous works on spherically symmetric accretion onto black holes and super-spinars to the case in which the fluid has a finite angular momentum initially. We run 2.5D and 3D general relativistic hydrodynamic simulations of the accretion of a fat disk. We study how the accretion process changes by changing the values of the parameters of our model. We show that the value of the fluid angular momentum critically determines turn-on and off the production of powerful equatorial outflows around super-spinars. For corotating disks, equatorial outflows are efficiently generated, even for relatively low spin parameters or relatively large super-spinar radii. For counterrotating disks, equatorial outflows are instead significantly suppressed, and they are possible only in limited cases. We also study accretion around a tilted disk.

The collapse of dense cores with different metallicities is studied by hydrodynamical calculations coupled with detailed chemical and radiative processes. For this purpose, we construct a simple chemical network with non-equilibrium reactions among 15 chemical species, which reproduces the abundance of important molecular coolants by more detailed network very well. The evolution is followed until the formation of a hydrostatic protostar at the center. In a lower-metallicity gas cloud, the temperature during the collapse remains high owing to less efficient cooling. Using the temperature evolution at the center as a function the density, we discuss the possibility of fragmentation during the dust-cooling phase. The critical metallicity for the fragmentation is 10^{-5}Z_sun assuming moderate elongation of the cloud cores at the onset of this phase. From the density and velocity distributions at the time of protostar formation, we evaluate the mass accretion rate in the subsequent accretion phase. Using these accretion rates, we also calculate the evolution of the protostars under the assumption of stationary accretion flow. Finally, we discuss possible suppression of fragmentation by heating of the ambient gas by protostellar radiation, which is considered important in the contemporary star formation. We argue that it is negligible for <10^{-2}Zsun.

We present the results of three-dimensional general relativistic hydrodynamic simulations of adiabatic and spherically symmetric accretion in Kerr space-time. We consider compact objects with spin parameter $|a_*| \le 1$ (black holes) and with $|a_*| > 1$ (super-spinars). Our full three-dimensional simulations confirm the formation of equatorial outflows for high values of $|a_*|$, as found in our previous work in 2.5 dimensions. We show that the critical value of $|a_*|$ determining the onset of powerful outflows depends mainly on the radius of the compact object. The phenomenon of equatorial outflows can hardly occur around a black hole and may thus be used to test the bound $|a_*| \le 1$ for astrophysical black hole candidates.

Within 5-10 years, very long baseline interferometry facilities will be able to observe the "shadow" of super-massive black hole candidates. This will allow, for the first time, to test gravity in the strong field regime. In this paper, we study numerically the photon orbits in the $\delta = 2$ Tomimatsu-Sato space-time. The $\delta = 2$ Tomimatsu-Sato space-time is a stationary, axisymmetric, and asymptotically flat exact solution of the vacuum Einstein equations. We compare the associated shadow with the one of Kerr black holes. The shape of the shadow in the $\delta = 2$ Tomimatsu-Sato space-time is oblate and the difference between the two axes can be as high as 6% when viewed on the equatorial plane. We argue that future space sub-mm interferometers (e.g. VSOP-3) may distinguish the two cases, and thus are able to test the Cosmic Censorship Conjecture.

In this paper we continue our study on the accretion process onto super-spinning Kerr objects with no event horizon (super-spinars). We discuss the counterpart of the Bondi accretion onto black holes. We first report the results of our numerical simulations. We found a quasi steady-state configuration for any choice of the parameters of our model. The most interesting feature is the presence of hot outflows. Unlike jets and outflows produced around black holes, which are thought to be powered by magnetic fields and emitted from the poles, here the outflows are produced by the repulsive gravitational force at a small distance from the super-spinar and are ejected around the equatorial plane. In some circumstances, the amount of matter in the outflow is considerable, which can indeed significantly reduce the gas mass accretion rate. Finally, we discuss a possible scenario of the accretion process in more realistic situations, which cannot be simulated by our code.

The accretion process onto spinning objects in Kerr spacetimes is studied with numerical simulations. Our results show that accretion onto compact objects with Kerr parameter (characterizing the spin) $|a| < M$ and $|a| > M$ is very different. In the super-spinning case, for $|a|$ moderately larger than $M$, the accretion onto the central object is extremely suppressed due to a repulsive force at short distance. The accreting matter cannot reach the central object, but instead is accumulated around it, forming a high density cloud that continues to grow. The radiation emitted in the accretion process will be harder and more intense than the one coming from standard black holes; e.g. $\gamma$-rays could be produced as seen in some observations. Gravitational collapse of this cloud might even give rise to violent bursts. As $|a|$ increases, a larger amount of accreting matter reaches the central object and the growth of the cloud becomes less efficient. Our simulations find that a quasi-steady state of the accretion process exists for $|a|/M \gtrsim 1.4$, independently of the mass accretion rate at large radii. For such high values of the Kerr parameter, the accreting matter forms a thin disk at very small radii. We provide some analytical arguments to strengthen the numerical results; in particular, we estimate the radius where the gravitational force changes from attractive to repulsive and the critical value $|a|/M \approx 1.4$ separating the two qualitatively different regimes of accretion. We briefly discuss the observational signatures which could be used to look for such exotic objects in the Galaxy and/or in the Universe.

We propose a model in which intermediate-mass black holes (IMBHs) with mass of ~10000 Msun are formed in early dark matter halos. We carry out detailed stellar evolution calculations for accreting primordial stars including annihilation energy of dark matter particles. We follow the stellar core evolution consistently up to gravitational collapse. We show that very massive stars, as massive as 10000 Msun, can be formed in an early dark matter halo. Such stars are extremely bright with Log L/Lsun > 8.2. They gravitationally collapse to form IMBHs. These black holes could have seeded the formation of early super-massive blackholes.

Star formation remains an unsolved problem in astrophysics. Numerical studies of large-scale structure simulations cannot resolve the whole process and their approach usually assumes that only gas denser than a typical threshold can host and form stars. We investigate the onset of cosmological star formation and compare several very-high-resolution, three-dimensional, N-body/SPH simulations that include non-equilibrium, atomic and molecular chemistry, star formation prescriptions, and feedback effects. We study how primordial star formation depends on gas density thresholds, cosmological parameters and initial set-ups. For mean-density initial conditions, we find that standard low-density star-formation threshold (0.2 h^2/cm3) models predict the onset of star formation at z~25-31, depending on the adopted cosmology. In these models stars are formed regardless of the time between the moment when the threshold is reached and the effective runaway collapse. At high redshift, this time interval represents a significant fraction of the Hubble time and thus this assumption can induce large artificial off-sets to the onset of star formation. Choosing higher density thresholds (135 h^2/cm3) allows the entire cooling process to be followed, and the onset of star formation is then estimated to be at redshift z~12-16. When isolated, rare, high-density peaks are considered, the chemical evolution is much faster and the first star formation episodes occur at z > 40, almost regardless of the choice for the density threshold. These results could have implications for the formation redshift of the first cosmological objects, as inferred from direct numerical simulations of mean-density environments, and on the studies of the reionization history of the universe.

The standard theory of cosmic structure formation posits that the present-day rich structure of the Universe developed through gravitational amplification of tiny matter density fluctuations generated in its very early history. Recent observations of the cosmic microwave background, large-scale structure, and distant supernovae determined the energy content of the Universe and the basic statistics of the initial density field with great accuracy. It has become possible to make accurate predictions for the formation and nonlinear growth of structure from early to the present epochs. We review recent progress in the theory of structure formation in the universe. We focus on the formation of the first cosmological objects. Results from state-of-the-art numerical simulations are presented. Finally, we discuss prospects for future observations of the first generation of stars and galaxies.

Observations made using large ground-based and space-borne telescopes have probed cosmic history all the way from the present-day to a time when the Universe was less than a tenth of its present age. Earlier on lies the remaining frontier, where the first stars, galaxies, and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

We use 5000 cosmological N-body simulations of 1(Gpc/h)^3 box for the concordance LCDM model in order to study the sampling variances of nonlinear matter power spectrum. We show that the non-Gaussian errors can be important even on large length scales relevant for baryon acoustic oscillations (BAO). Our findings are (1) the non-Gaussian errors degrade the cumulative signal-to-noise ratios (S/N) for the power spectrum amplitude by up to a factor of 2 and 4 for redshifts z=1 and 0, respectively. (2) There is little information on the power spectrum amplitudes in the quasi-nonlinear regime, confirming the previous results. (3) The distribution of power spectrum estimators at BAO scales, among the realizations, is well approximated by a Gaussian distribution with variance that is given by the diagonal covariance component. (4) For the redshift-space power spectrum, the degradation in S/N by non-Gaussian errors is mitigated due to nonlinear redshift distortions. (5) For an actual galaxy survey, the additional shot noise contamination compromises the cosmological information inherent in the galaxy power spectrum, but also mitigates the impact of non-Gaussian errors. The S/N is degraded by up to 30% for a WFMOS-type survey. (6) The finite survey volume causes additional non-Gaussian errors via the correlations of long-wavelength fluctuations with the fluctuations we want to measure, further degrading the S/N values by about 30% even at high redshift z=3.

The cosmic dark ages are the mysterious epoch during which the pristine gas began to condense and ultimately form the first stars. Although these beginnings have long been a topic of theoretical interest, technology has only recently allowed the beginnings of observational insight into this epoch. Many questions surround the formation of stars in metal-free gas and the history of the build-up of metals in the intergalactic medium: (1) What were the properties of the first stellar and galactic sources to form in pristine (metal-free) gas? (2) When did the epoch of Population III (metal-free) star formation take place and how long did it last? (3) Was the stellar initial mass function dramatically different for the first stars and galaxies? These questions are all active areas of theoretical research. However, new observational constraints via the direct detection of Population III star formation are vital to making progress in answering the broader questions surrounding how galaxies formed and how the cosmological properties of the universe have affected the objects it contains.

We calculate the evolution of zero-metallicity Population III (Pop III) stars whose mass grows from the initial mass of $\sim 1M_{\odot}$ by accreting the surrounding gases. Our calculations cover a whole evolutionary stages from the pre-main sequence, via various nuclear burning stages, through the final core collapse or pair-creation instability phases. We adopt the following stellar mass-dependent accretion rates which are derived from cosmological simulations of early structure formation based on the low mass dark matter halos at redshifts $z \sim 20$: (1) the accretion rates for the first generation (Pop III.1) stars and (2) the rates for zero-metallicity but the second generation (Pop III.2) stars which are affected by radiation from the Pop III.1 stars. For comparison, we also study the evolution with the mass-dependent accretion rates which are affected by radiatibe feedback. We show that the final mass of Pop III.1 stars can be as large as $\sim 1000M_{\odot}$, beyond the mass range ($140 - 300M_{\odot}$) for the pair-instability supernovae. Such massive stars undergo core-collapse to form intermediate-mass black holes, which may be the seeds for merger trees to supermassive black holes. On the other hand, Pop III.2 stars become less massive ($\lsim 40 - 60M_{\odot}$), being in the mass range of ordinary iron core-collapse stars. Such stars explode and eject heavy elements to contribute to chemical enrichment of the early universe as observed in the abundance patterns of extremely metal-poor stars in the Galactic halo.

We calculate the spin temperature and 21 cm brightness of early \HII regions around the first stars. We use outputs from cosmological radiation-hydrodynamics simulations of the formation and evolution of early \HII regions. In the pre-reionization era, \HII regions around massive primordial stars have diameters of a few kpc. The gas within the \HII regions is almost fully ionized, but begins recombining after the central stars die off. The relic \HII regions are then seen as bright {\it emission} sources in hydrogen 21 cm. We make brightness temperature maps of the \HII regions, accounting for radiative coupling with Lyman-$\alpha$ photons in a simplified manner. The spin temperature in the relic \HII region is close to the gas kinetic temperature, generally several hundred to several thousand degrees. We show that the relic \HII region can be as bright as $\delta T_{\rm b} \sim 100$ mK in differential temperature against the cosmic microwave background for an angular resolution of sub-arcseconds. While individual early \HII patches will not be identified by currently planned radio telescopes, the collective fluctuations from early \HII regions might imprint signatures in the 21 cm background.

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the following contribution present the story of Dark Stars. In this first part, we describe the conditions under which dark stars form in the early universe: 1) high dark matter densities, 2) the annihilation products get stuck inside the star, and 3) dark matter heating wins over all other cooling or heating mechanisms.

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the previous contribution present the story of Dark Stars. In this second part, we describe the structure of Dark Stars and predict that they are very massive ($\sim 800 M_\odot$), cool (6000 K), bright ($\sim 10^6 L_\odot$), long-lived ($\sim 10^6$ years), and probable precursors to (otherwise unexplained) supermassive black holes. Later, once the initial dark matter fuel runs out and fusion sets in, dark matter annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again.

We investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a \Lambda Cold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum.

The nature of the first generation of stars in the Universe remains largely unknown. Observations imply the existence of massive primordial stars early in the history of the universe, and the standard theory for the growth of cosmic structure predicts that structures grow hierarchically through gravitational instability. We have developed an ab initio computer simulation of the formation of primordial stars that follows the relevant atomic and molecular processes in a primordial gas in an expanding universe. The results show that primeval density fluctuations left over from the Big Bang can drive the formation of a tiny protostar with a mass of just one percent that of the Sun. The protostar is a seed for the subsequent formation of a massive primordial star.

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the the universe. This talk presents the story of these Dark Stars. We make predictions that the first stars are very massive ($\sim 800 M_\odot$), cool (6000 K), bright ($\sim 10^6 L_\odot$), long-lived ($\sim 10^6$ years), and probable precursors to (otherwise unexplained) supermassive black holes. Later, once the initial DM fuel runs out and fusion sets in, DM annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again.

We critically examine how well the evolution of large-scale density perturbations is followed in cosmological $N$-body simulations. We first run a large volume simulation and perform a mode-by-mode analysis in three-dimensional Fourier space. We show that the growth of large-scale fluctuations significantly deviates from linear theory predictions. The deviations are caused by {\it nonlinear} coupling with a small number of modes at largest scales owing to finiteness of the simulation volume. We then develop an analytic model based on second-order perturbation theory to quantify the effect. Our model accurately reproduces the simulation results. For a single realization, the second-order effect appears typically as ``zig-zag'' patterns around the linear-theory prediction, which imprints artificial ``oscillations'' that lie on the real baryon-acoustic oscillations. Although an ensemble average of a number of realizations approaches the linear theory prediction, the dispersions of the realizations remain large even for a large simulation volume of several hundred megaparsecs on a side. For the standard $\Lambda$CDM model, the deviations from linear growth rate are as large as 10 percent for a simulation volume with $L = 500h^{-1}$Mpc and for a bin width in wavenumber of $\Delta k = 0.005h$Mpc$^{-1}$, which are comparable to the intrinsic variance of Gaussian random realizations. We find that the dispersions scales as $\propto L^{-3/2} \Delta k^{-1/2}$ and that the mean dispersion amplitude can be made smaller than a percent only if we use a very large volume of $L > 2h^{-1}$Gpc. The finite box size effect needs to be appropriately taken into account when interpreting results from large-scale structure simulations for future dark energy surveys using baryon acoustic oscillations.

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

The detection of large quantities of dust in z ~ 6 quasars by infrared and radio surveys presents puzzles for the formation and evolution of dust in these early systems. Previously (Li et al. 2007), we showed that luminous quasars at z > 6 can form through hierarchical mergers of gas-rich galaxies. Here, we calculate the dust properties of simulated quasars and their progenitors using a three-dimensional Monte Carlo radiative transfer code, ART^2 -- All-wavelength Radiative Transfer with Adaptive Refinement Tree. ART^2 incorporates a radiative equilibrium algorithm for dust emission, an adaptive grid for inhomogeneous density, a multiphase model for the ISM, and a supernova-origin dust model. We reproduce the SED and dust properties of SDSS J1148+5251, and find that the infrared emission are closely associated with the formation and evolution of the quasar host. The system evolves from a cold to a warm ULIRG owing to heating and feedback from stars and AGN. Furthermore, the AGN has significant implications for the interpretation of observation of the hosts. Our results suggest that vigorous star formation in merging progenitors is necessary to reproduce the observed dust properties of z~6 quasars, supporting a merger-driven origin for luminous quasars at high redshifts and the starburst-to-quasar evolutionary hypothesis. (Abridged)

Non-linear evolution of density fluctuations yields secondary anisotropies in the cosmic microwave background (CMB), which are correlated with the same density fluctuations that can be measured by weak lensing (WL) surveys. We study the CMB-WL correlation using analytical models as well as $N$-body simulations. We show that an analytical model based upon the time derivative of non-linear matter power spectrum agrees with simulations. All-sky cosmic-variance limited CMB and WL surveys allow us to measure the correlation from non-linearity with high significance ($50\sigma$) for $l_{\rm max}=10^4$, whereas the forthcoming missions such as Planck and LSST are expected to yield only marginal detections. The CMB-WL correlation is sensitive to the time derivative of structure growth. We study how this property may be used to constrain the nature of dark energy. While the expected constraints are not very strong, they may provide a cross check of results from other observations.

We use cosmological hydrodynamic simulations with unprecedented resolution to study the formation of primordial stars in an ionized gas at high redshifts. Our approach includes all the relevant atomic and molecular physics to follow the thermal evolution of a prestellar gas cloud to very high densities of ~10^{18} cm^{-3}. We locate a star-forming gas cloud within a reionized region in our cosmological simulation. The first run-away collapse is triggered when the gas cloud's mass is ~40 Msun. We show that the cloud core remains stable against chemo-thermal instability and also against gravitational deformation throughout its evolution. Consequently, a single proto-stellar seed is formed, which accretes the surrounding hot gas at the rate ~10^{-3} Msun/year. We carry out proto-stellar evolution calculations using the inferred accretion rate. The resulting mass of the star when it reaches the zero-age main sequence is M_ZAMS ~40 Msun. We argue that, since the obtained M_ZAMS is as large as the mass of the collapsing parent cloud, the final stellar mass should be close to this value. Such massive, rather than exceptionally massive, primordial stars are expected to cause early chemical enrichment of the Universe by exploding as black hole-forming super/hypernovae, and may also be progenitors of high redshift gamma-ray bursts. The elemental abundance patterns of recently discovered hyper metal-poor stars suggest that they might have been born from the interstellar medium that was metal-enriched by supernovae of these massive primordial stars.

The discovery of luminous quasars at redshift z ~ 6 indicates the presence of supermassive black holes (SMBHs) of mass ~10^9 Msun when the Universe was less than one billion years old. This finding presents several challenges for theoretical models. Here, we present the first multi-scale simulations that, together with a self-regulated model for the SMBH growth, produce a luminous quasar at z ~ 6.5 in the LCDM paradigm. We follow the hierarchical assembly history of the most massive halo in a ~ 3 Gpc^3 volume, and find that this halo of ~ 8x 10^{12} Msun forming at z ~ 6.5 after several major mergers is able to reproduce a number of observed properties of SDSS J1148+5251, the most distant quasar detected at z =6.42 (Fan et al. 2003). Moreover, the SMBHs grow through gas accretion below the Eddington limit in a self-regulated manner owing to feedback. We find that the progenitors experience significant star formation (up to 10^4 Msun/yr) preceding the major quasar phase such that the stellar mass of the quasar host reaches 10^{12} Msun at z ~ 6.5, consistent with observations of significant metal enrichment in SDSS J1148+5251. Our results provide a viable formation mechanism for z ~ 6 quasars in the standard LCDM cosmology, and demonstrate a common, merger-driven origin for the rarest quasars and the fundamental SMBH-host correlation in a hierarchical Universe.(Abridged)

We explore observational constraints on possible deviations from Newtonian gravity by means of large-scale clustering of galaxies. We measure the power spectrum and the bispectrum of Sloan Digital Sky Survey galaxies and compare the result with predictions in an empirical model of modified gravity. Our model assumes an additional Yukawa-like term with two parameters that characterize the amplitude and the length scale of the modified gravity. The model predictions are calculated using two methods; the second-order perturbation theory and direct N-body simulations. These methods allow us to study non-linear evolution of large-scale structure. Using the simulation results, we find that perturbation theory provides reliable estimates for the power spectrum and the bispectrum in the modified Newtonian model. We also construct mock galaxy catalogues from the simulations, and derive constraints on the amplitude and the length scale of deviations from Newtonian gravity. The resulting constraints from power spectrum are consistent with those obtained in our earlier work, indicating the validity of the previous empirical modeling of gravitational nonlinearity in the modified Newtonian model. If linear biasing is adopted, the bispectrum of the SDSS galaxies yields constraints very similar to those from the power spectrum. If we allow for the nonlinear biasing instead, we find that the ratio of the quadratic to linear biasing coefficients, b_2/b_1, should satisfy -0.4 < b_2/b_1<0.3 in the modified Newtonian model.

We have performed a large set of high-resolution cosmological simulations using smoothed particle hydrodynamics to study the formation of the first luminous objects in the LCDM cosmology. We follow the collapse of primordial gas clouds in eight early structures and document the scatter in the properties of the first star-forming clouds. Our first objects span formation redshifts from z~10 to z~50 and cover an order of magnitude in halo mas. We find that the physical properties of the central star-forming clouds are very similar in all of the simulated objects despite significant differences in formation redshift and environment. The physical properties of the clouds have little correlation with spin, mass, or assembly history of the host halo. The collapse of protostellar objects at higher redshifts progresses much more rapidly due to the higher densities, which accelerates the formation of molecular hydrogen, enhances initial cooling and shortens the dynamical timescales. The mass of the star-forming clouds cover a broad range, from a few hundred to a few thousand solar masses, and exhibit various morphologies: some of have disk-like structures nearly rotational supported; others form flattened spheroids; still others form bars. All of them develop a single proto-stellar `seed' which does not fragment into multiple objects up to the moment that the central gas becomes optically thick to H2 cooling lines. At this time, the instantaneous mass accretion rate onto the centre varies significantly from object to object, with disk-like structuress have the smallest mass accretion rates. The formation epoch and properties of the star-forming clouds are sensitive to the values of cosmological parameters.

(Abridged) We present the results of three-dimensional radiation-hydrodynamics simulations of the formation and evolution of early HII/HeIII regions around the first stars. Cooling (by H2 and HD) and recollapse of the gas in the relic HII region is also followed in a full cosmological context, until second-generation stars are formed. A large HII region with a few kiloparsec diameter is formed, within which a smaller HeIII region is embedded. Radiative feedback effect quenches further star-formation within the halo for a hundred million years. Accretion onto remnant blackholes will be inefficient. Recombination radiation within the HII region is weak, but persists for 50 million years. We also follow the thermal and chemical evolution of the photo-ionized gas in the relic HII region. The gas cools by HD line cooling down to a few tens Kelvin. At high redshifts (z>10), the minimum gas temperature is limited by T_CMB. Because of its low temperature, the characteristic mass of a Jeans-unstable gas clump is ~ 40 Msun, and is significantly smaller than a typical clump mass for early primordial gas clouds. We find no evidence of fragmentation by this epoch. Together with the small cloud mass, this result indicates that massive, rather than very massive, primordial stars may form in the relic HII region. Such stars might be responsible for early metal-enrichment of the interstellar medium from which recently discovered hyper metal-poor stars were born.

We present the first hydrodynamic N-body simulations of primordial gas clouds responsible for the reionisation process in dark energy cosmologies. We compare the cosmological constant scenario with a SUGRA quintessence model with marked dynamics in order to highlight effects due to the different acceleration histories imposed by the dark energy. We show that both the number density of gas clouds and their clumpiness keep a record of the expansion rate during evolution, similar to the non-linear dark matter profile at virialisation, as was recently demonstrated by Dolag et al. (2004). Varying the shape of the primordial power spectrum, we show how this effect is mitigated by a running spectral index decreasing the power at small scales. Our results demonstrate that, in order to constrain the dark energy from large scale structures, one must track its effects down to the distribution of luminous matter.

We study the formation of the first generation of stars in the standard cold dark matter model, using a very high-resolution hydordynamic simulations. Our simulation achieves a dynamic range of 10^{10} in length scale. With accurate treatment of atomic and molecular physics, it allows us to study the chemo-thermal evolution of primordial gas clouds to densities up to n = 10^{16}/cc without assuming any a priori equation of state; a six orders of magnitudes improvement over previous three-dimensional calculations. All the relevant atomic and molecular cooling and heating processes, including cooling by collision-induced continuum emission, are implemented. For calculating optically thick H2 cooling at high densities, we use the Sobolev method. To examine possible gas fragmentation owing to thermal instability, we compute explicitly the growth rate of isobaric perturbations. We show that the cloud core does not fragment in either the low-density or high-density regimes. We also show that the core remains stable against gravitational deformation and fragmentation. We obtain an accurate gas mass accretion rate within a 10 Msun innermost region around the protostar. The protostar is accreting the surrounding hot gas at a rate of 0.001-0.01 Msun/yr. From these findings we conclude that primordial stars formed in early minihalos are massive. We carry out proto-stellar evolution calculations using the obtained accretion rate. The resulting mass of the first star is M_ZAMS = 60-100 Msun, with the exact mass dependent on the actual accretion rate.

We simulate the growth of galaxies and their central supermassive black holes by implementing a suite of semi-analytic models on the output of the Millennium Run, a very large simulation of the concordance LCDM cosmogony. Our procedures follow the detailed assembly history of each object and are able to track the evolution of all galaxies more massive than the Small Magellanic Cloud throughout a volume comparable to that of large modern redshift surveys. In this first paper we supplement previous treatments of the growth and activity of central black holes with a new model for `radio' feedback from those AGN that lie at the centre of a quasistatic X-ray emitting atmosphere in a galaxy group or cluster. We show that for energetically and observationally plausible parameters such a model can simultaneously explain: (i) the low observed mass drop-out rate in cooling flows; (ii) the exponential cut-off at the bright end of the galaxy luminosity function; and (iii) the fact that the most massive galaxies tend to be bulge-dominated systems in clusters and to contain systematically older stars than lower mass galaxies. This success occurs because static hot atmospheres form only in the most massive structures, and radio feedback (in contrast, for example, to supernova or starburst feedback) can suppress further cooling and star formation without itself requiring star formation. We discuss possible physical models which might explain the accretion rate scalings required for our phenomenological `radio mode' model to be successful.

(Abridged) The violent hierarchical nature of the LCDM cosmology poses serious difficulties for the formation of disk galaxies. To help resolve these issues, we describe a new, merger-driven scenario for the cosmological formation of disk galaxies at high redshifts that supplements the standard model based on dissipational collapse.In this picture, large gaseous disks may be produced from high-angular momentum mergers of systems that are gas-dominated, i.e. M_gas/(M_gas +M_star > 0.5 at the height of the merger. Pressurization from the multiphase structure of the interstellar medium prevents the complete conversion of gas into stars during the merger, and if enough gas remains to form a disk, the remnant eventually resembles a disk galaxy. We perform numerical simulations of galaxy mergers to study how supernovae feedback strength, supermassive black hole growth and feedback, progenitor gas fraction, merger mass-ratio, and orbital geometry impact the formation of remnant disks. We find that disks can build angular momentum through mergers and the degree of rotational support of the baryons in the merger remnant is primarily related to feedback processes associated with star formation. Disk-dominated remnants are restricted to form in mergers that are gas-dominated at the time of final coalescence and gas-dominated mergers typically require extreme progenitor gas fractions (>80%). We also show that the formation of rotationally-supported stellar systems in mergers is not restricted to idealized orbits, or major or minor mergers. We suggest that the hierarchical nature of the LCDM cosmology and the physics of the interstellar gas may act together to form spiral galaxies by building the angular momentum of disks through early, gas-dominated mergers.

In Figure 6 we inadvertently labeled the proxy circular velocity as the virial velocity of the dark matter halo instead of what is actually plotted, the maximum circular velocity of the dark matter halo. The maximum halo circular velocity is a much better estimate of the disk V_c than is V_vir. This confusion influenced the discussion of the Tully-Fisher relation in our paper. In fact, Figure 6 demonstrates that it is possible to simultaneously reproduce both the local Tully-Fisher relation and luminosity function using semi-analytic techniques applied to the standard LCDM cosmology, thus contradicting previous studies of this issue and our own discussion in Section 3.6.

The currently standard theory of cosmic structure formation posits that the present-day clumpy appearance of the universe developed through gravitational amplification of the matter density fluctuations that are generated in the very early universe. The energy content of the univese and the basic statistics of the initial density field have been determined with a reasonable accuracy from recent observations of the cosmic microwave background, large-scale structure, and distant supernovae. It has become possible to make accurate predictions from the standard model. We review the latest observations and the recent progress in the theory of structure formation at low and high redshifts. Two promising methods to probe large-scale matter distribution are introduced and the future prospects are discussed. Results from state-of-the-art cosmological simulations are also presented.

We study the evolution of supernova (SN) remnants of the first stars, taking proper account of the radiative feedback of the progenitor stars on the surroundings. We carry out a series of one-dimensional hydrodynamic simulations with radiative cooling, starting from initial configurations that are drawn from the results of our earlier radiation hydrodynamic simulations of the first HII regions. In low-mass (< 10^6 M_sun) halos, the stellar radiation significantly reduces the ambient gas density prior to the SN explosion. The blastwave quickly propagates over the halo's virial radius, leading to complete evacuation of the gas even with the input energy of 10^50 erg. We find that a large fraction of the remnant's thermal energy is lost in 0.1-10 Myr by line cooling, whereas, for larger explosion energies, the remnant expands even more rapidly with decreasing interior density, and cools predominantly via inverse Compton process. In higher mass halos, the gas density near the explosion site remains high and the SN shock is heavily confined; the thermal energy of the remnant is quickly radiated away by free-free emission, even if the total input energy exceeds the binding energy of halos by two orders of magnitude. We show that the efficiency of halo destruction is determined not only by the explosion energy but also by the gas density profile, and thus controlled by radiative feedback prior to the explosion. Several implications of our results for the formation of first quasars and second-generation stars in the universe are also discussed.

The cold dark matter model has become the leading theoretical paradigm for the formation of structure in the Universe. Together with the theory of cosmic inflation, this model makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a novel framework for the quantitative physical interpretation of such surveys. This combines the largest simulation of the growth of dark matter structure ever carried out with new techniques for following the formation and evolution of the visible components. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with next generation surveys.

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular $\Lambda$ Cold Dark Matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of -0.8 or less. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement. Voids very clearly correspond to minima in the smoothed initial density field. We find a universal void mass profile of the form $\rho(<r)/\rho(r_{\rm eff}) \propto \exp[(r/r_{\rm eff})^\alpha]$ where $r_{\rm eff}$ is the effective radius of a void and $\alpha\sim 2$. The mass function of haloes in voids is steeper than that of haloes that populate denser regions. The abundances of void haloes seem to evolve somewhat more strongly between redshifts $\sim 1$ and 0 than the global abundances of haloes. (abridged)

In spite of the growing observational evidence for dark matter and dark energy in the universe, their physical nature is largely unknown. In fact, several authors have proposed modifications of Newton's law of gravity at cosmological scales to account for the apparent acceleration of the cosmic expansion. Inspired by such suggestions, we attempt to constrain possible deviations from Newtonian gravity by means of the clustering of SDSS (Sloan Digital Sky Survey) galaxies. To be specific, we assume a simple modification of Newton's law with an additional Yukawa-type term characterized by the amplitude \alpha and the length scale \lambda. Adopting spatially-flat universes dominated by cold dark matter and/or dark energy, we solve a linear perturbation equation for the growth of density fluctuations. In particular, we find an exact analytic solution for the Einstein -- de Sitter case. Following the Peacock-Dodds prescription, we compute the nonlinear power spectra of mass fluctuations, perform a statistical comparison with the SDSS galaxy data, and derive constraints in the \alpha-\lambda plane; for instance, we obtain the constraints of -0.5<\alpha <0.6 and -0.8<\alpha<0.9 (99.7% confidence level) for \lambda=5h^{-1}Mpc and 10h^{-1}Mpc, respectively. We also discuss several future possibilities for improving our analysis.

We present a simple model that accurately describes various statistical properties of peculiar velocities of dark matter halos. We pay particular attention to the following two effects; first, the evolution of the halo peculiar velocity depends on the local matter density, instead of the global density. Second, dark matter halos are biased tracers of the underlying mass distribution, thus halos tend to be located preferentially at high density regions. For the former, we develop an empirical model calibrated with N-body simulations, while for the latter, we use a conventional halo bias models based on the extended Press-Schechter model combined with an empirical log-normal probability distribution function of the mass density distribution. We find that compared with linear theory, the present model significantly improves the accuracy of predictions of statistical properties of the halo peculiar velocity field including the velocity dispersion, the probability distribution function, and the pairwise velocity dispersion at large separations. Thus our model predictions may be useful in analyzing future observations of the peculiar velocities of galaxy clusters.

We study the temperature structure of the intergalactic medium (IGM) using a large cosmological N-body/SPH simulation. We employ a two-temperature model for the thermal evolution of the ionized gas, in which the relaxation process between electrons and ions is explicitly included. In the diffuse, hot IGM, the relaxation time is comparable to the age of the Universe and hence the electron temperature in post-shock regions remains significantly smaller than the ion temperature. We show that, at the present epoch, a large fraction of the warm/hot intergalactic medium (WHIM) has a well-developed two temperature structure, with typical temperature differences of order a factor of a few. Consequently, the fraction of metals in various ionization states such as OVI, OVII, and OVIII, as well as their line emissivities, can differ locally by more than an order of magnitude from those computed with a single temperature model. It is thus necessary to follow the evolution of the electron temperature explicitly to determine absorption and emission by the WHIM. Although equipartition is nearly achieved in the denser intracluster medium (ICM), we find an appreciable systematic deviation between the gas-mass weighted electron temperature and the mean temperature even at half the virial radii of clusters. There is thus a reservoir of warm (Te < 1keV) gas in and around massive clusters. Our results imply that relaxation processes need to be considered in describing and interpreting observational data from existing X-ray telescopes as well as from future missions designed to detect the WHIM, such as the Diffuse Intergalactic Oxygen Surveyor and the Missing Baryon Explorer.

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, QSOs and also the CMB temperature map), is known to consist of two components; a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a Lambda-CDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating those two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive halos with M>10^14h^-1Msun at z<1 and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arc minute, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at sub-arcminute scales.

We study the formation and evolution of HII regions around the first stars formed at z=10-30. We use a one-dimensional Lagrangian hydrodynamics code which self-consistently incorporates radiative transfer and non-equilibrium primordial gas chemistry. The star-forming region is defined as a spherical gas cloud with a Population III star embedded at the center. We explore a large parameter space of host halo mass, gas density profile, and stellar luminosity. The formation of the HII region is characterized by initial slow expansion of a weak D-type ionization front near the center, followed by rapid propagation of an R-type front throughout the outer gas envelope. We find that the transition between the two front types is indeed a critical condition for the complete ionization of halos of cosmological interest. In small mass (< 10^6 M_sun) halos, the transition takes place within a few 10^5 yr, yielding high escape fractions (>80%) of both ionizing and photodissociating photons. In larger mass (> 10^7 M_sun) halos, the ionization front remains to be of D-type over the lifetime of the massive star, the HII region is confined well inside the virial radius, and the escape fractions are essentially zero. We derive an analytic formula that reproduces well the results of our simulations. We discuss immediate implications of the present results for the star formation history and early reionization of the Universe.

The cosmic infrared background (IRB) at wavelengths between 1 micron and 3 microns provides a useful probe of early star-formation prior to and during reionization. To explain the high optical depth to electron scattering, as measured by the Wilkinson Microwave Anisotropy Probe (WMAP), one requires significant star-formation activity at redshifts 10 and higher. In addition to massive stars, the IRB flux may be contributed by a population of early miniquasars. We study the relative contributions from first stars, supernovae and quasars to the IRB for reasonable star formation rates at high redshift. If miniquasars radiate efficiently at the Eddington-limit, current background measurements limit the fraction of mass in first stars that is converted to seed black holes to be roughly less than 10%. In the case of supernovae, though an individual supernova is much brighter than the progenitor star, due to the shorter lifetime of order few months, the fractional contribution to the IRB remains at a level of 10% and below when compared to the same contribution from stars. The bright supernovae may, however, be directly detectable by future large ground-based and space telescopes.

We analyze the dynamical and thermal evolution of dark matter and ICM in hydrodynamical Tree-SPH simulations of galaxy clusters. Starting from a sample of 17 high-resolution objects, with virial mass ranging from 3 10^{14} to 1.7 10^{15} M_sun, we follow the build-up of the systems in dark matter and hot gas through the repeated merging of satellites along their merging history trees. We measure the self-bound mass fraction of subhaloes as a function of time after the merging, estimate the satellite mean orbital properties as a function of the mass ratio with the main cluster at merging time, and study the evolution of their internal velocity dispersion, gas temperature and entropy as the substructure is disrupted by various dynamical processes, reaching eventually thermo-dynamic equilibrium in the gravitational potential of the main cluster. We model some relevant properties of subhalo orbits, as the time of the first pericentric and apocentric passages, and the typical distances and velocities at the corresponding times. This survival study can be used to interpret the dynamics of observed merging clusters; as an example we apply our results to the system 1E0657-56. We show that, in the light of our results, the most likely interpretation of the data for this cluster points to the merger of a small group with mass M\approx 1 10^{13} M_sun with a massive cluster with M\approx 1.3 10^{15} M_sun.

We explore the ability of weak lensing surveys to locate massive clusters. We use both analytic models of dark matter halos and mock weak lensing surveys generated from a large cosmological N-body simulation. The analytic models describe average properties of weak lensing halos and predict the number counts, enabling us to compute an effective survey selection function. We test the model prediction for the peak number counts in weak lensing mass maps against the mock numerical data, and find that the noise due to intrinsic galaxy ellipticities causes a systematic effect which increases the peak counts. We develop a correction scheme for the systematic effect in an empirical manner, and show that, after the correction, the model prediction agrees well with the mock data. The mock data is also used to examine the completeness and efficiency of the weak lensing halo search with fully taking into account the noise and the projection effect by large-scale structures. We show that the detection threshold of S/N=4-5 gives an optimal balance between completeness and efficiency. Our results suggest that, for a weak lensing survey with a galaxy number density of ng=30/arcmin^2 with a mean redshift z=1, the mean number of peaks in the 10sq deg area is N_peak=62 for a detection threshold S/N=4. The contamination rate is 42%, and thus, on average, 36 out of 62 peaks (at least) are signals from real halos. Weak lensing surveys thus provide a reasonably efficient way to searching for massive clusters.

If the dark matter in the Universe is made of weakly self-interacting particles, they may self-annihilate and emit gamma-rays. We use high resolution numerical simulations to estimate directly the annihilation flux from the central regions of the Milky Way and from dark matter substructures in its halo. Although such estimates remain uncertain because of their strong dependence on the structure of the densest regions, our numerical experiments suggest that less direct calculations have overestimated the emission both from the centre and from halo substructure. We estimate a maximal enhancement of at most a factor of a few with respect to a smooth spherical halo of standard Navarro-Frenk-White (NFW) structure. We discuss detection strategies for the next generation of gamma-ray detectors and find that the annihilation flux may be detectable, regardless of uncertainties about the densest regions, for the annihilation cross-sections predicted by currently popular elementary particle models for the dark matter.

We combine fast radiative transfer calculations with high resolution hydrodynamical simulations to study an epoch of early hydrogen reionisation by primordial stellar sources at redshifts 15<z<30. With relatively conservative assumptions, population III star formation proceeds in a self-regulated manner both locally and globally and, for a conventional LCDM cosmology, can significantly reionise the intergalactic medium between 15<z<20 as long as a large fraction of ionising photons can escape from these earliest galaxies. We then combine these results with our earlier work focusing on the role of population II stars in galaxies with virial temperatures >10^4K at redshifts 5<z< 20. Hence, we construct a complete reionisation history of the Universe which matches the Thomson optical depths as measured by the WMAP satellite as well as the evolution of the Gunn Peterson optical depth as seen in the asborption spectra of the higest redshift quasars. We find that even with conservative estimates for the impact of negative feedback mechanisms, primordial stellar sources contribute significantly to early reionisation. Future observations of a Thomson optical depth of tau_e>~0.13 would bolster the claim for the existence of population III stars similar to the ones studied here.

We describe hydrodynamical simulations of galaxy formation in a Lambda cold dark matter (CDM) cosmology performed using a subresolution model for star formation and feedback in a multiphase interstellar medium (ISM). In particular, we demonstrate the formation of a well-resolved disk galaxy. The surface brightness profile of the galaxy is exponential, with a B-band central surface brightness of 21.0 mag arcsec^-2 and a scale-length of R_d = 2.0 h^-1 kpc. We find no evidence for a significant bulge component. The simulated galaxy falls within the I-band Tully-Fisher relation, with an absolute magnitude of I = -21.2 and a peak stellar rotation velocity of V_rot=121.3 km s^-1. While the total specific angular momentum of the stars in the galaxy agrees with observations, the angular momentum in the inner regions appears to be low by a factor of ~2. The star formation rate of the galaxy peaks at ~7 M_sun yr^-1 between redshifts z=2-4, with the mean stellar age decreasing from \~10 Gyrs in the outer regions of the disk to ~7.5 Gyrs in the center, indicating that the disk did not simply form inside-out. The stars exhibit a metallicity gradient from 0.7 Z_sun at the edge of the disk to 1.3 Z_sun in the center. Using a suite of idealized galaxy formation simulations with different models for the ISM, we show that the effective pressure support provided by star formation and feedback in our multiphase model is instrumental in allowing the formation of large, stable disk galaxies. If ISM gas is instead modeled with an isothermal equation of state, or if star formation is suppressed entirely, growing gaseous disks quickly violate the Toomre stability criterion and undergo catastrophic fragmentation.

The birth and death of the first generation of stars have important implications for the thermal state and chemical properties of the intergalactic medium (IGM) in the early universe. Sometime after recombination, the neutral, chemically pristine gas was reionized by ultraviolet photons emitted from the first stars, but also enriched with heavy elements when these stars ended their lives as energetic supernovae. Using the results from previous high-resolution cosmological simulations of early structure formation that include radiative transfer, we show that a significant volume fraction of the IGM can be metal-polluted, as well as ionized, by massive Population III stars formed in small-mass (10^6-10^7 Msun) halos early on. If most of the early generation stars die as pair-instability supernovae with energies up to 10^{53} ergs, the volume-averaged mean metallicity will quickly reach Z ~ 10^{-4}Zsun by a redshift of 15-20, possibly causing a prompt transition to the formation of a stellar population that is dominated by low-mass stars. In this scenario, the early chemical enrichment history should closely trace the reionization history of the IGM, and the end of the Population III era is marked by the completion of reionization and pre-enrichment by z=15. We conclude that, while the pre-enrichment may partially account for the ``metallicity-floor'' in high-redshift Lyman-alpha clouds, it does not significantly affect the elemental abundance in the intracluster medium.

We examine the ability of the Sunyaev-Zel'dovich effect (SZE) statistics proposed by Lee (2002) using numerical simulations. The statistics describe the distribution of the peak heights in the noise-contaminated SZE sky maps, and provide an analysis technique to sort out the noise-contamination effect and estimate the number density of real clusters. The method is devised to be suitable for the interferometric SZE observations in drift-scanning mode like AMiBA experiment. We apply the proposed method to a set of realistic SZE sky maps constructed from large-scale cosmological simulations, and show that the method indeed allows us to estimate the number density of clusters efficiently. The efficiency of the method is demonstrated in two aspects: (i) it can reconstruct the number density of clusters even though the clusters contribute only a small fraction (10%) of the total number of peaks in each SZE map; (ii) it can count even those clusters whose amplitudes are low enough to be compatible with the noise level while conventional method of cluster identification using the threshold cutoff cannot count them properly within the same observation time. Thus, the proposed method may be useful for the study of cluster abundance exploiting future interferometric SZE surveys such as AMiBA.

We analyse the properties of substructures within dark matter halos (subhalos) using a set of high-resolution numerical simulations of the formation of structure in a Lambda-CDM Universe. Our simulation set includes 11 high-resolution simulations of massive clusters as well as a region of mean density, allowing us to study the spatial and mass distribution of substructures down to a mass resolution limit of 10^9 h^(-1)Mo. We also investigate how the properties of substructures vary as a function of the mass of the `parent' halo in which they are located. We find that the substructure mass function depends at most weakly on the mass of the parent halo and is well described by a power-law. The radial number density profiles of substructures are steeper in low mass halos than in high mass halos. More massive substructures tend to avoid the centres of halos and are preferentially located in the external regions of their parent halos. We also study the mass accretion and merging histories of substructures, which we find to be largely independent of environment. We find that a significant fraction of the substructures residing in clusters at the present day were accreted at redshifts z < 1. This implies that a significant fraction of present-day `passive' cluster galaxies should have been still outside the cluster progenitor and more active at z~1.

We investigate the supernova explosions that end the lives of massive Population III stars in low-mass minihalos (M~10^6 M_sun) at redshifts z~20. Employing the smoothed particle hydrodynamics method, we carry out numerical simulations in a cosmological set-up of pair-instability supernovae with explosion energies of E_SN=10^51 and 10^53 ergs. We find that the more energetic explosion leads to the complete disruption of the gas in the minihalo, whereas the lower explosion energy leaves much of the halo intact. The higher energy supernova expels > 90% of the stellar metals into a region ~1 kpc across over a timescale of 3-5 Myr. Due to this burst-like initial star formation episode, a large fraction of the universe could have been endowed with a metallicity floor, Z_min>10^-4 Z_sun, already at z>15.

The first-year Wilkinson Microwave Anisotropy Probe data suggest a high optical depth for Thomson scattering of 0.17 +/- 0.04, implying that the universe was reionized at an early epoch, z ~ 20. Such early reionization is likely to be caused by UV photons from first stars, but it appears that the observed high optical depth can be reconciled within the standard structure formation model only if star-formation in the early universe was extremely efficient. With normal star-formation efficiencies, cosmological models with non-Gaussian density fluctuations may circumvent this conflict as high density peaks collapse at an earlier epoch than in models with Gaussian fluctuations. We study cosmic reionization in non-Gaussian models and explore to what extent, within available constraints, non-Gaussianities affect the reionization history. For mild non-Gaussian fluctuations at redshifts of 30 to 50, the increase in optical depth remains at a level of a few percent and appears unlikely to aid significantly in explaining the measured high optical depth. On the other hand, within available observational constraints, increasing the non-Gaussian nature of density fluctuations can easily reproduce the optical depth and may remain viable in underlying models of non-Gaussianity with a scale-dependence.

(abridged) We study high redshift structure formation and reionization in a LCDM universe under the assumption that the spectral power index of primordial density fluctuations is a function of length scale. We adopt a particular formulation of the running spectral index (RSI) model as suggested by the recent WMAP data. While early structure forms hierarchically in the RSI model, the reduced power on small scales causes a considerable delay in the formation epoch of low mass (~ 10^6 Msun) ``mini-halos'' compared to the LCDM model. The extremely small number of gas clouds in the RSI model indicates that reionization is initiated later than z<15, generally resulting in a smaller total Thomson optical depth than in the LCDM model. By carrying out radiative transfer calculations, we also study reionization by stellar populations formed in galaxies. Even with a top-heavy intial mass function representing an early population of massive stars and/or an extraordinarily high photon emission rate from galaxies, the total optical depth can only be as large as tau ~ 0.1 for reasonable models of early star-formation. The RSI model is thus in conflict with the large Thomson optical depth inferred by the WMAP satellite.

We critically examine how the evolution of the matter density field in cosmological simulations is affected by details of setting up initial conditions. We show that it is non-trivial to realise an initial distribution of matter in N-body/hydrodynamic simulations so that the baryon and dark matter density fluctuations and their velocities evolve consistently as theoretically predicted. We perform a set of cosmological simulations and use them to distinguish and verify an appropriate method for generating initial conditions. We show that a straightforward way of applying the Zel'dovich approximation to each component using distinct transfer functions results in an incorrect growth of density fluctuations and that it is necessary to correct velocities at the initial epoch. The unperturbed uniform particle distribution must be also generated appropriately to avoid tight coupling of the baryonic and dark matter components. We recommend using independent ``glass'' particle distributions, using distinct transfer functions for baryons and dark matter, and taking into account the difference in the velocity fields at the initialisation epoch. The proposed method will be useful for studies of the evolution of the intergalactic medium and the formation of the first cosmological objects using numerical simulations.

We study first structure formation in Lambda-dominated universes using large cosmological N-body/SPH simulations. We consider a standard LCDM model and a LWDM model in which the mass of the dark matter particles is taken to be m_X=10 keV. The linear power spectrum for the LWDM model has a characteristic cut-off at a wavenumber k=200 /Mpc, suppressing the formation of low mass (< 10^6 Msun) nonlinear objects early on. The absence of low mass halos in the WDM model makes the formation of primordial gas clouds with molecular hydrogen very inefficient at high redshifts. The first star-forming gas clouds form at z~21 in the WDM model, considerably later than in the CDM counterpart, and the abundance of these gas clouds differs by an order of magnitude between the two models. We carry out radiative transfer calculations by embedding massive Population III stars in the gas clouds. We show that the volume fraction of ionized gas rises up close to 100% by z=18 in the CDM case, whereas that of the WDM model remains extremely small at a level of a few percent. Thus the WDM model with m_X=10 keV is strongly inconsistent with the observed high optical depth by the WMAP satellite.

(abridged) We use large cosmological simulations to study the origin of primordial star-forming clouds in a Lambda CDM universe, by following the formation of dark matter halos and the cooling of gas within them. To model the physics of chemically pristine gas, we employ a non-equilibrium treatment of the chemistry of 9 species and include cooling by molecular hydrogen. We explore the hierarchical growth of bound structures forming at redshifts z = 25 - 30 with total masses in the range 10^5 - 10^6 Msun. The complex interplay between the gravitational formation of dark halos and the thermodynamic and chemical evolution of the gas clouds compromises analytic estimates of the critical H2 fraction. Dynamical heating from mass accretion and mergers opposes relatively inefficient cooling by molecular hydrogen, delaying the production of star-forming clouds in rapidly growing halos. We also investigate the impact of photo-dissociating ultra-violet (UV) radiation on the formation of primordial gas clouds. We consider two extreme cases by first including a uniform radiation field in the optically thin limit and secondly by accounting for the maximum effect of gas self-shielding in virialized regions. In both the cases we consider, the overall impact can be described by computing an equilibrium H2 abundance for the radiation flux and defining an effective shielding factor. Based on our numerical results, we develop a semi-analytic model of the formation of the first stars, and demonstrate how it can be coupled with large N-body simulations to predict the star formation rate in the early universe.

The systematic magnification of background sources by the weak gravitational-lensing effects of foreground matter, also called cosmic magnification, is becoming an efficient tool both for measuring cosmological parameters and for exploring the distribution of galaxies relative to the dark matter. We extend here the formalism of magnification statistics by estimating the contribution of second-order terms in the Taylor expansion of the magnification and show that the effect of these terms was previously underestimated. We test our analytical predictions against numerical simulations and demonstrate that including second-order terms allows the accuracy of magnification-related statistics to be substantially improved. We also show, however, that both numerical and analytical estimates can provide only lower bounds to real correlation functions, even in the weak lensing regime. We propose to use count-in-cells estimators rather than correlation functions for measuring cosmic magnification since they can more easily be related to correlations measured in numerical simulations.

We use a simple analytic model to compute the angular correlation function of clusters identified in upcoming thermal SZ effect surveys. We then compute the expected fraction of close pairs of clusters on the sky that are also close along the line of sight. We show how the expected number of cluster pairs as a function of redshift is sensitive to the assumed biasing relation between the cluster and the mass distribution. We find that, in a LambdaCDM model, the fraction of physically associated pairs is 70% for angular separations smaller than 20 arcmin and clusters with specific flux difference larger than 200 mJy at 143 GHz. The agreement of our analytic results with the Hubble volume N-body simulations is satisfactory. These results quantify the feasibility of using SZ surveys to compile catalogues of superclusters at any redshifts.

We compare two techniques for following the cooling of gas and its condensation into galaxies within high resolution simulations of cosmologically representative regions. Both techniques treat the dark matter using N-body methods. One follows the gas using smoothed particle hydrodynamics (SPH) while the other uses simplified recipes from semi-analytic (SA) models. We compare the masses and locations predicted for dense knots of cold gas (the `galaxies') when the two techniques are applied to evolution from the same initial conditions and when the additional complications of star formation and feedback are ignored. We find that above the effective resolution limit of the two techniques, they give very similar results both for global quantities such as the total amount of cooled gas and for the properties of individual `galaxies'. The SA technique has systematic uncertainties arising from the simplified cooling model adopted, while details of the SPH implementation can produce substantial systematic variations in the galaxy masses it predicts. Nevertheless, for the best current SPH methods and the standard assumptions of the SA model, systematic differences between the two techniques are remarkably small. The SA technique gives adequate predictions for the condensation of gas into `galaxies' at less than one percent of the computational cost of obtaining similar results at comparable resolution using SPH.

We critically examine the reliability of the dark matter clustering in high-resolution cosmological $N$-body simulations on scales below the mean separation length of particles. The particle discreteness effect imposes the two fundamental limitations on those scales; the lack of the initial fluctuation power and the finite mass resolution. We address this problem applying the dark halo approach and are able to discuss separately how those two limitations affect the dark matter clustering in $N$-body simulations at early epochs. We find that limitations of the dark matter clustering are primarily determined by the mass of particles. By a detailed comparison with three major cosmological simulations, we also find that in order to reproduce a proper amplitude of the dark matter clustering on small scales, halos with a characteristic nonlinear mass, $M_{\rm NL}(z)$ defined by $\sigma_R(M_{\rm NL}; z)=1$, must be resolved in the simulation. This leads to a critical redshift $z_{\rm crit}$ determined by $M_{\rm NL}(z_{\rm crit})=n_{\rm halo} m_{\rm part}$ where $n_{\rm halo}$ is the number of particles necessary to resolve the typical nonlinear mass halo ($\sim 10$). We conclude that, at least as far as the two-point correlation functions are concerned, the dark matter clustering in high-resolution $N$-body simulations on scales below the mean particle separation is reliable down to the gravitational force resolution length only for $z<z_{\rm crit}$, while it is strongly affected by the finite mass resolution for $z>z_{\rm crit}$.

We use giga-particle N-body simulations to study galaxy cluster populations in Hubble Volumes of LCDM (Omega_m=0.3, Omega_Lambda=0.7) and tCDM (Omega_m=1) world models. Mapping past light-cones of locations in the computational space, we create mock sky surveys of dark matter structure to z~1.4 over 10,000 sq deg and to z~0.5 over two full spheres. Calibrating the Jenkins mass function at z=0 with samples of ~1.5 million clusters, we show that the fit describes the sky survey counts to <~20% acccuracy over all redshifts for systems larger than poor groups (M>5e13 Msun/h). Fitting the observed local temperature function determines the ratio beta of specific thermal energies in dark matter and intracluster gas. We derive a scaling with power spectrum normalization beta \propto sigma8^{5/3}, and measure a 4% error on sigma8 arising from cosmic variance in temperature-limited cluster samples. Considering distant clusters, the LCDM model matches EMSS and RDCS X-ray-selected survey observations under economical assumptions for intracluster gas evolution. Using transformations of mass-limited cluster samples that mimic sigma8 variation, we explore SZ search expectations for a 10 sq deg survey complete above 10^{14} Msun/h. Cluster counts are shown to be extremely sensitive to sigma8 uncertainty while redshift statistics, such as the sample median, are much more stable. For LCDM, the characteristic temperature at fixed sky surface density is a weak function of redshift, implying an abundance of hot clusters at z>1. Assuming constant beta, four kT>8 keV clusters lie at z>2 and 40 kT>5 keV clusters lie at z>3 on the whole sky. Detection of Coma-sized clusters at z>1 violate LCDM at 95% confidence if their surface density exceeds 0.003 per sq deg, or 120 on the whole sky.

We develop a phenomenological model to predict the clustering of dark matter halos on the light-cone by combining several existing theoretical models. Assuming that the velocity field of halos on large scales is approximated by linear theory, we propose an empirical prescription of a scale-, mass-, and time-dependence of halo biasing. We test our model against the Hubble Volume $N$-body simulation and examine its validity and limitations. We find a good agreement in two-point correlation functions of dark matter halos between the phenomenological model predictions and measurements from the simulation for $R>5h^{-1}$Mpc both in the real and redshift spaces. Although calibrated on the mass scale of groups and clusters and for redshifts up to $z\sim2$, the model is quite general and can be applied to a wider range of astrophysical objects, such as galaxies and quasars, if the relation between dark halos and visible objects is specified.

We use numerical simulations of a (480 Mpc/h)^3 volume to show that the distribution of peak heights in maps of the temperature fluctuations from the kinematic and thermal Sunyaev-Zeldovich effects will be highly non-Gaussian, and very different from the peak height distribution of a Gaussian random field. We then show that it is a good approximation to assume that each peak in either SZ effect is associated with one and only one dark matter halo. This allows us to use our knowledge of the properties of haloes to estimate the peak height distributions. At fixed optical depth, the distribution of peak heights due to the kinematic effect is Gaussian, with a width which is approximately proportional to optical depth; the non-Gaussianity comes from summing over a range of optical depths. The optical depth is an increasing function of halo mass, and the distribution of halo speeds is Gaussian, with a dispersion which is approximately independent of halo mass. This means that observations of the kinematic effect can be used to put constraints on how the abundance of massive clusters evolves, and on the evolution of cluster velocities. The non-Gaussianity of the thermal effect, on the other hand, comes primarily from the fact that, on average, the effect is larger in more massive haloes, and the distribution of halo masses is highly non-Gaussian. We also show that because haloes of the same mass may have a range of density and velocity dispersion profiles, the relation between halo mass and the amplitude of the thermal effect is not deterministic, but has some scatter.

We describe the newly written code GADGET which is suitable both for cosmological simulations of structure formation and for the simulation of interacting galaxies. GADGET evolves self-gravitating collisionless fluids with the traditional N-body approach, and a collisional gas by smoothed particle hydrodynamics. Along with the serial version of the code, we discuss a parallel version that has been designed to run on massively parallel supercomputers with distributed memory. While both versions use a tree algorithm to compute gravitational forces, the serial version of GADGET can optionally employ the special-purpose hardware GRAPE instead of the tree. Periodic boundary conditions are supported by means of an Ewald summation technique. The code uses individual and adaptive timesteps for all particles, and it combines this with a scheme for dynamic tree updates. Due to its Lagrangian nature, GADGET thus allows a very large dynamic range to be bridged, both in space and time. So far, GADGET has been successfully used to run simulations with up to 7.5e7 particles, including cosmological studies of large-scale structure formation, high-resolution simulations of the formation of clusters of galaxies, as well as workstation-sized problems of interacting galaxies. In this study, we detail the numerical algorithms employed, and show various tests of the code. We publically release both the serial and the massively parallel version of the code.

We create mock pencil-beam redshift surveys from very large cosmological $N$-body simulations of two Cold Dark Matter cosmogonies, an Einstein-de Sitter model ($\tau$CDM) and a flat model with $\Omega_0 =0.3$ and a cosmological constant ($\Lambda$CDM). We use these to assess the significance of the apparent periodicity discovered by Broadhurst et al. (1990). Simulation particles are tagged as `galaxies' so as to reproduce observed present-day correlations. They are then identified along the past light-cones of hypothetical observers to create mock catalogues with the geometry and the distance distribution of the Broadhurst et al. data. We produce 1936 (2625) quasi-independent catalogues from our $\tau$CDM ($\Lambda$CDM) simulation. A couple of large clumps in a catalogue can produce a high peak at low wavenumbers in the corresponding one-dimensional power spectrum, without any apparent large-scale periodicity in the original redshift histogram. Although the simulated redshift histograms frequently display regularly spaced clumps, the spacing of these clumps varies between catalogues and there is no `preferred' period over our many realisations. We find only a 0.72 (0.49) per cent chance that the highest peak in the power spectrum of a $\tau$CDM ($\Lambda$CDM) catalogue has a peak-to-noise ratio higher than that in the Broadhurst et al. data. None of the simulated catalogues with such high peaks shows coherently spaced clumps with a significance as high as that of the real data. We conclude that in CDM universes, the kind of regularity observed by Broadhurst et al. has a priori probability well below $10^{-3}$.

We carry out ray tracing through five high resolution simulations of a galaxy cluster to study how its ability to produce giant gravitationally lensed arcs is influenced by the collision cross-section of its dark matter. In three cases typical dark matter particles in the cluster core undergo between 1 and 100 collisions per Hubble time; two more explore the long (``collisionless'') and short (``fluid'') mean free path limits. We study the size and shape distributions of arcs and compute the cross-section for producing ``extreme'' arcs of various sizes. Even a few collisions per particle modify the core structure enough to destroy the cluster's ability to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to unrealistically large masses before it regains the ability to produce giant arcs. None of our models with self-interacting dark matter (except the ``fluid'' limit) is able to produce radial arcs; even the case with the smallest scattering cross-section must be scaled to the upper limit of observed cluster masses before it produces radial arcs. Apparently the elastic collision cross-section of dark matter in clusters must be very small, below 0.1 cm^2 g^-1, to be compatible with the observed ability of clusters to produce both radial arcs and giant arcs.

We combine data from a number of N-body simulations to predict the abundance of dark halos in Cold Dark Matter universes over more than 4 orders of magnitude in mass. A comparison of different simulations suggests that the dominant uncertainty in our results is systematic and is smaller than 10--30% at all masses, depending on the halo definition used. In particular, our ``Hubble Volume'' simulations of \tcdm and \lcdm cosmologies allow the abundance of massive clusters to be predicted with uncertainties well below those expected in all currently planned observational surveys. We show that for a range of CDM cosmologies and for a suitable halo definition, the simulated mass function is almost independent of epoch, of cosmological parameters, and of initial power spectrum when expressed in appropriate variables. This universality is of exactly the kind predicted by the familiar Press-Schechter model, although this model predicts a mass function shape which differs from our numerical results, overestimating the abundance of ``typical'' halos and underestimating that of massive systems.

We are making available on the WWW a selection of the archived data from N-body simulations carried out by the Virgo consortium and related groups. This currently includes: (i) time-slice, lightcone and cluster data from the two $10^9$-particle Hubble volume simulations described by Evrard 1998; (ii) time-slice data from simulations of 4 different cold dark matter cosmological models with $256^3$ particles analysed by Jenkins et al 1998; (iii) Dark halo catalogs, merger trees and galaxy catalogs from the GIF project described by Kauffmann et al 1999. Basic software is supplied to read the data. The data can be accessed from: http://www.mpa-garching.mpg.de/Virgo/data_download.html

We study the formation of dark halos in a $\Lambda$CDM universe under the assumption that Cold Dark Matter particles have a finite cross-section for elastic collisions. We compare evolution when CDM mean free paths are comparable to halo sizes with the collisionless and fluid limits. We show that a few collisions per particle per Hubble time at halo centre can substantially affect the central density profile. Cross-sections an order of magnitude larger produce sufficient relaxation for rich clusters to develop core radii in the range 100-200 $h^{-1}$kpc. The structural evolution of halos is a competition between collisional relaxation caused by individual particle interactions and violent relaxation resulting from the infall and merging processes by which clusters grow. Although our simulations concentrate on systems of cluster size, we can scale our results to address the halo structure expected for dwarf galaxies. We find that collision cross-sections sufficiently large to significantly modify the cores of such galaxies produce cluster cores which are too large and/or too round to be consistent with observation. Thus the simplest model for self-interacting dark matter is unable to improve fits to published dwarf galaxy rotation curves without violating other observational constraints.

We use very large cosmological N--body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ($\tau$CDM) and a flat low density model with a cosmological constant ($\Lambda$CDM). Our simulations each use $10^9$ particles to follow the mass distribution within cubes of side $2h^{-1}$Gpc ($\tau$CDM) and $3h^{-1}$Gpc ($\Lambda$CDM) with a force resolution better than $10^{-4}$ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from $r_0=12$ -- 13$h^{-1}$Mpc for samples with mean separation $d_{\rm c}=30h^{-1}$Mpc to $r_0=22$-- 27$h^{-1}$Mpc for samples with $d_{\rm c}=100h^{-1}$Mpc. The lower value here corresponds to $\tau$CDM and the upper to $\Lambda$CDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10% using the simple model of Sheth, Mo & Tormen (2000). This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the APM cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for $\Lambda$CDM. The observed power spectrum (for a cluster sample of mean separation $d_{\rm c}=31h^{-1}$Mpc) lies significantly above the predictions of both models.

We study how the internal structure of dark halos is affected if Cold Dark Matter particles are assumed to have a large cross-section for elastic collisions. We identify a cluster halo in a large cosmological N-body simulation and resimulate its formation with progressively increasing resolution. We compare the structure found in the two cases where dark matter is treated as collisionless or as a fluid. For the collisionless case our results agree with those of other workers. Collisional dark matter results in a cluster which is more nearly spherical and has a more singular central density profile. Substructure within the cluster is only weakly suppressed relative to the collisionless case. The observed structure of dwarf galaxies argues against self-interacting dark matter if, as seems likely, intermediate cross-sections produce structure lying between the extremes we have simulated.