A series of gravitational instabilities in a circumnuclear gas disk (CND) are required to trigger gas transport to a central supermassive black hole (SMBH) and ignite Active Galactic Nuclei (AGNs). A test of this scenario is to investigate whether an enhanced molecular gas mass surface density ($\Sigma_{\rm mol}$) is found in the CND-scale of quasars relative to a comparison sample of inactive galaxies. Here we performed sub-kpc resolution CO(2-1) observations with ALMA of four low-redshift ($z \sim 0.06$), luminous ($\sim 10^{45}$ erg s$^{-1}$) quasars with each matched to a different star-forming galaxy, having similar redshift, stellar mass, and star-formation rate. We detected CO(2-1) emission from all quasars, which show diverse morphologies. Contrary to expectations, $\Sigma_{\rm mol}$ of the quasar sample, computed from the CO(2-1) luminosity, tends to be smaller than the comparison sample at $r < 500$ pc; there is no systematic enhancement of $\Sigma_{\rm mol}$ in our quasars. We discuss four possible scenarios that would explain the lower molecular gas content (or CO(2-1) luminosity as an actual observable) at the CND-scale of quasars, i.e., AGN-driven outflows, gas-rich minor mergers, time-delay between the onsets of a starburst-phase and a quasar-phase, and X-ray-dominated region (XDR) effects on the gas chemical abundance and excitation. While not extensively discussed in the literature, XDR effects can have an impact on molecular mass measurements particularly in the vicinity of luminous quasar nuclei; therefore higher resolution molecular gas observations, which are now viable using ALMA, need to be considered.

Recent cosmological simulations have shown that turbulence should be generally prevailing in clusters because clusters are continuously growing through matter accretion. Using one-dimensional hydrodynamic simulations, we study the heating of cool-core clusters by the ubiquitous turbulence as well as feedback from the central active galactic nuclei (AGNs) for a wide range of cluster and turbulence parameters, focusing on the global stability of the core. We find that the AGN shows intermittent activities in the presence of moderate turbulence similar to the one observed with Hitomi. The cluster core maintains a quasi-equilibrium state for most of the time because the heating through turbulent diffusion is nearly balanced with radiative cooling. The balance is gradually lost because of slight dominance of the radiative cooling, and the AGN is ignited by increased gas inflow. Finally, when the AGN bursts, the core is heated almost instantaneously. Thanks to the pre-existing turbulence, the heated gas is distributed throughout the core without becoming globally unstable and causing catastrophic cooling, and the core recovers the quasi-equilibrium state. The AGN bursts can be stronger in lower-mass clusters. Predictions of our model can be easily checked with future X-ray missions like XRISM and Athena.

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation-hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation-hydrodynamic simulations to study the clumping factors needed to solve the analytical equations, and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve $<7\%$ difference from simulation results while the other models have $>20\%$ deviations. The Thomson optical depths from reionization calculated from our analytical models result in $<5\%$ deviation from simulations, while the previous analytical models have $>20\%$ difference in and could result in biased conclusions of the IGM reionization.

We present an HST Cosmic Origins Spectrograph (COS) spectrum of the QSO SDSS J095109.12+330745.8 ($z_{em} = 0.645$) whose sightline passes through the SMC-like dwarf galaxy UGC 5282 ($M_B = -16.0$, $cz = 1577$ km s$^{-1}$ ), 1.2 kpc in projection from the central HII region of the galaxy. Damped Ly$\alpha$ (DLA) absorption is detected at the redshift of UGC 5282 with $\log$[$N$(HI) cm$^{-2}\:$] = 20.89$^{+0.12}_{-0.21}\:\:$. Analysis of the accompanying SII, PII and OI metal lines yields a neutral gas metallicity, $Z_{HI}$, of [S/H] $\simeq$ [P/H] $= -0.80\pm0.24$. The metallicity of ionized gas from the central HII region, $Z_{HII}\:$, measured from its emission lines is [O/H]$= -0.37\pm0.10$, a difference of $+0.43\pm 0.26$ from $Z_{HI}\:$. This difference $\delta$ is consistent with that seen towards HII regions in other star-forming galaxies and supports the idea that ionized gas near star forming regions shows systematically higher metallicities than exist in the rest of a galaxy's neutral interstellar medium (ISM). The positive values of $\delta$ found in UGC 5282 (and the other star forming galaxies) is likely due to infalling low metallicity gas from the intergalactic medium that mixes with the galaxy's ISM on kpc scales. This model is also consistent with broad Ly$\alpha$ emission detected at the bottom of the DLA absorption, offset by $\sim 125$ km s$^{-1}$ from the absorption velocity. Models of galaxy evolution that attempt to replicate population characteristics, such as the mass-metallicity relation, may need to start with a galaxy metallicity represented by $Z_{HI}$ rather than that measured traditionally from $Z_{HII}\:$.

Lyman-$\alpha$ blobs (LABs) are spatially extended nebulae of emission in the Ly$\alpha$ line of hydrogen, seen at high redshifts$^{1,2}$, and most commonly found in the dense environment of star-forming galaxies$^{3,4}$. The origin of Ly$\alpha$ emission in the LABs is still unclear and under debate$^{5}$. Proposed powering sources generally fall into two categories: (1) photoionization, galactic super-winds/outflows, resonant scattering of Ly$\alpha$ photons from starbursts or active galactic nuclei (AGNs)$^{6,7,8,9,10}$ and (2) cooling radiation from cold streams of gas accreting onto galaxies$^{12}$. Here we analyze the gas kinematics within a LAB providing rare observational evidence for infalling gas. This is consistent with the release of gravitational accretion energy as cold streams radiate Ly$\alpha$ photons. It also provides direct evidence for possible cold streams feeding the central galaxies. The infalling gas is not important by mass but hints at more than one mechanism to explain the origin of the extended Ly$\alpha$ emission around young galaxies. It is also possible that the infalling gas may represent material falling back to the galaxy from where it originated, forming a galactic fountain.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

An analysis of the dynamics of a star formation event is performed. It is shown that galaxies able to drive leftover gas to sufficient altitudes in a few million years are characterized by two basic properties: small sizes (<1kpc) and high star formation rate surface densities (Sigma_SFR > 10 Msun/yr/kpc2). For the parameter space of relevance, the outflow is primarily driven by supernovae with radiation pressure being significant but subdominant. Our analysis provides the unifying physical origin for a diverse set of observed LyC leakers, including the green-peas galaxies, [SII]-weak galaxies, Lyman-alpha emitters, with these two characteristics as the common denominator. Among verifiable physical properties of LyC leakers, we predict that (1) the newly formed stellar masses are are typically in the range of 1e8-1e10 Msun, except perhaps ULIRGs, (2) the outflow velocities are typically in the range typically of 100-600km/s, but may exceed 1e3 km/s in ULIRGs, with a strong positive correlation between the stellar masses formed and the outflow velocities, (3) the overall escape fraction of galaxies is expected to increase with increasing redshift, given the cosmological trend that galaxies become denser and more compact with increasing redshift. In addition, two interesting by-product predictions are also borne out. First, ULIRGs appear to be in a parameter region where they should be prodigious LyC leakers, unless there is a large ram-pressure. Second, Lyman break galaxies (LBGs) are not supposed to be prodigious LyC leakers in our model, given their claimed effective radii exceeding 1kpc.

A starburst induced by a galaxy merger may create a relatively thin central stellar disk at radius $\le 100$pc. We calculate the rate of tidal disruption events (TDEs) by the inspiraling secondary supermassive black (SMBH) through the disk. With a small enough stellar velocity dispersion ($\sigma/v_c \le 0.1$) in the disk, it is shown that $10^5-10^6$ TDEs of solar-type main sequence stars per post-starburst galaxy (PSB) can be produced to explain their dominance in producing observed TDEs. Although the time it takes to bring the secondary SMBH to the disk apparently varies in the range of $\sim 0.1-1$Gyr since the starburst, depending on its landing location and subsequently due to dynamical friction with stars exterior to the central stellar disk in question, the vast majority of TDEs by the secondary SMBH in any individual PSB occurs within a space of time shorter than $\sim 30$Myr. Five unique testable predictions of this model are suggested.

The evolution of the galaxy size - stellar mass (Mstellar) relation has been a puzzle for over a decade. High redshift galaxies are significantly more compact than galaxies observed today, at an equivalent mass, but how much of this apparent growth is driven by progenitor bias, minor mergers, secular processes, or feedback from AGN is unclear. To help disentangle the physical mechanisms at work by addressing the latter, we study the galaxy size - Mstellar relation of 32 carefully-selected broad-line AGN hosts at 1.2 < z < 1.7 (7.5 < log M_BH < 8.5; L_bol/L_Edd > 0.1). Using HST with multi-band photometry and state-of-the-art modeling techniques, we measure half-light radii while accounting for uncertainties from subtracting bright central point sources. We find AGN hosts to have sizes ranging from 1 to 6 kpc at Mstellar ~ 0.3 - 1 x 10^11 Msun. Thus, many hosts have intermediate sizes as compared to equal-mass star-forming and quiescent galaxies. While inconsistent with the idea that AGN feedback may induce an increase in galaxy sizes, this finding is consistent with hypotheses in which AGNs preferentially occur in systems with prior concentrated gas reservoirs, or are involved in secular compaction processes perhaps responsible for simultaneously building bulges and shutting down star formation. If driven by minor mergers, which do not grow central black holes as fast as they do bulge-like stellar structures, such a process would explain both the galaxy size - mass relation observed here and the evolution in the black hole, bulge mass relation described in a companion paper.

The relationship between galaxies and the state/chemical enrichment of the warm-hot intergalactic medium (WHIM) expected to dominate the baryon budget at low-z provides sensitive constraints on structure formation and galaxy evolution models. We present a deep redshift survey in the field of 1ES1553+113, a blazar with a unique combination of UV+X-ray spectra for surveys of the circum-/intergalactic medium (CGM/IGM). Nicastro et al. 2018 reported the detection of two O VII WHIM absorbers at $z=0.4339$ and $0.3551$ in its spectrum, suggesting that the WHIM is metal-rich and sufficient to close the missing baryons problem. Our survey indicates that the blazar is a member of a $z=0.433$ group and that the higher-$z$ O VII candidate arises from its intragroup medium. The resulting bias precludes its use in baryon censuses. The $z=0.3551$ candidate occurs in an isolated environment 630 kpc from the nearest galaxy (with stellar mass $\log M_*/M_\odot \approx 9.7$) which we show is unexpected for the WHIM. Finally, we characterize the galactic environments of broad H I Ly$\alpha$ absorbers (Doppler widths of $b=40-80$ \kms; $T\lesssim4\times10^5$ K) which provide metallicity independent WHIM probes. On average, broad Ly$\alpha$, absorbers are ${\approx}2\times$ closer to the nearest luminous ($L>0.25 L_*$) galaxy (700 kpc) than narrow ($b<30$ \kms; $T\lesssim4\times10^5$ K) ones (1300 kpc) but ${\approx}2\times$ further than O\,VI absorbers (350 kpc). These observations suggest that gravitational collapse heats portions of the IGM to form the WHIM but with feedback that does not enrich the IGM far beyond galaxy/group halos to levels currently observable in UV/X-ray metal lines.

We report on the results of deep and wide-field (1.1 deg$^2$) narrow-band observations with Subaru/Hyper Suprime-Cam (HSC) of a field around a hyperluminous QSO (HLQSO), HS1549+1919, residing in a protocluster at $z=2.84$, to map the large-scale structure of Ly$\alpha$ emitters (LAEs). One HSC pointing enables us to detect 3490 LAEs and 76 extended Ly$\alpha$ blobs (LABs), probing diverse environments from voids to protoclusters. The HLQSO is found to be near the center of the protocluster, which corresponds to the intersection of $\sim$100 cMpc-scale structures of LAEs. LABs are basically distributed along the large-scale structure, with larger ones particularly clustered around the HLQSO, confirming a previously noted tendency of LABs to prefer denser environments. Moreover, the shapes of LABs near the HLQSO appear to be aligned with the large-scale structure. Finally, a deep Ly$\alpha$ image reveals a diffuse Ly$\alpha$ nebula along a filamentary structure with no luminous UV/sub-mm counterpart. We suggest that the diffuse nebula is due to a cold filament with high clumping factor illuminated by the QSO, with a required high clumpiness provided by unresolved residing halos of mass $\leq 10^{9-10}M_\odot$.

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations.

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we present new radiation-hydrodynamic simulations with updated high-redshift galaxy populations and varying radiation escape fractions. The simulations are designed to have fixed Thomson optical depth $\tau \approx 0.06$, consistent with recent Planck observations, and similar midpoints of reionization $7.5 \lesssim z \lesssim 8.0$, but with different ionization histories. The galaxy luminosity functions and ionizing photon production rates in our model are in good agreement with recent HST observations. Adopting a power-law form for the radiation escape fraction $f_{\text{esc}}(z) = f_8[(1+z)/9]^{a_{8}}$, we simulate the cases for $a_8 = 0$, 1, and 2 and find $a_8 \lesssim 2$ in order to end reionization in the range $5.5 \lesssim z \lesssim 6.5$ to be consistent with Lyman alpha forest observations. At fixed $\tau$ and as the power-law slope $a_8$ increases, the reionization process starts earlier but ends later with a longer duration $\Delta z$ and the decreased redshift asymmetry $Az$. We find a range of durations $3.9 \lesssim \Delta z \lesssim 4.6$ that is currently in tension with the upper limit $\Delta z < 2.8$ inferred from a recent joint analysis of Planck and South Pole Telescope observations.

Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$\alpha$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from clustering by cross-correlating the clusters with an external galaxy redshift catalog from the Sloan Digital Sky Survey III, the LOWZ sample. Our cross-correlations can constrain assembly bias up to a statistical uncertainty of 6%. Given our H$\alpha$ criteria for CC and NCC, we find no significant differences in their clustering amplitude. Interpreting this difference as the absence of halo assembly bias, our results rule out the possibility of having different large-scale (tens of Mpc) environments as the source of diversity observed in cluster cores. Combined with recent observations of the overall mild evolution of CC and NCC properties, such as central density and CC fraction, this would suggest that either the cooling properties of the cluster core are determined early on solely by the local (<200 kpc) gas properties at formation or that local merging leads to stochastic CC relaxation and disruption in a periodic way, preserving the average population properties over time. Studying the small-scale clustering in clusters at high redshift would help shed light on the exact scenario.

We study how the void environment affects galactic chemical evolution by comparing the oxygen and nitrogen abundances of dwarf galaxies in voids with dwarf galaxies in denser regions. Using spectroscopic observations from SDSS DR7, we estimate oxygen, nitrogen, and neon abundances of 889 void dwarf galaxies and 672 dwarf galaxies in denser regions. A substitute for the [OII] 3727 doublet is developed, permitting oxygen abundance estimates of SDSS dwarf galaxies at all redshifts with the Direct Te method. We find that void dwarf galaxies have about the same oxygen abundance and Ne/O ratio, slightly higher neon abundances, and slightly lower nitrogen abundance and N/O ratio than dwarf galaxies in denser environments. We conclude that the void environment has a slight influence on dwarf galaxy chemical evolution. Our mass-N/O relationship shows that the secondary production of nitrogen commences at a lower stellar mass in void dwarf galaxies than in dwarf galaxies in denser environments. Our dwarf galaxy sample demonstrates a strong anti-correlation between the sSFR and N/O ratio, providing evidence that oxygen is produced in higher mass stars than those which synthesize nitrogen. The lower N/O ratios and smaller stellar mass for secondary nitrogen production seen in void dwarf galaxies may indicate both delayed star formation and a dependence of cosmic downsizing on the large-scale environment. A shift toward slightly higher oxygen abundances in void dwarf galaxies could be evidence of larger ratios of dark matter halo mass to stellar mass in voids than in denser regions.

It has been known for decades that the observed number of baryons in the local universe falls about 30-40% short of the total number of baryons predicted by Big-Bang Nucleosynthesis, as inferred from density fluctuations of the Cosmic Microwave Background and seen during the first 2-3 billion years of the universe in the so called Lyman-alpha Forest. A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas - hydrogen - is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios. Indeed, despite the large observational efforts, only a few marginal claims of detection have been made so far. Here we report observations of two absorbers of highly ionized oxygen (OVII) in the high signal-to-noise-ratio X-ray spectrum of a quasar at redshift >0.4. These absorbers show no variability over a 2-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The OVII systems lie in regions characterized by large (x4 compared to average) galaxy over-densities and their number (down to the sensitivity threshold of our data), agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium (WHIM). We conclude that the missing baryons have been found.

We report an integral-field spectroscopic (IFS) observation of a gravitationally lensed spiral galaxy A1689B11 at redshift $z=2.54$. It is the most ancient spiral galaxy discovered to date and the second kinematically confirmed spiral at $z\gtrsim2$. Thanks to gravitational lensing, this is also by far the deepest IFS observation with the highest spatial resolution ($\sim$ 400 pc) on a spiral galaxy at a cosmic time when the Hubble sequence is about to emerge. After correcting for a lensing magnification of 7.2 $\pm$ 0.8, this primitive spiral disk has an intrinsic star formation rate of 22 $\pm$ 2 $M_{\odot}$ yr$^{-1}$, a stellar mass of 10$^{9.8 \pm 0.3}$$M_{\odot}$ and a half-light radius of $r_{1/2}=2.6 \pm 0.7$ kpc, typical of a main-sequence star-forming (SF) galaxy at $z\sim2$. However, the H\alpha\ kinematics show a surprisingly tranquil velocity field with an ordered rotation ($V_{\rm c}$ = 200 $\pm$ 12 km/s) and uniformly small velocity dispersions ($V_{\rm \sigma, mean}$ = 23 $\pm$ 4 km/s and $V_{\rm \sigma, outer-disk}$ = 15 $\pm$ 2 km/s). The low gas velocity dispersion is similar to local spiral galaxies and is consistent with the classic density wave theory where spiral arms form in dynamically cold and thin disks. We speculate that A1689B11 belongs to a population of rare spiral galaxies at $z\gtrsim2$ that mark the formation epoch of thin disks. Future observations with JWST will greatly increase the sample of these rare galaxies and unveil the earliest onset of spiral arms.

In La Plante et al. (2017), we presented a new suite of hydrodynamic simulations with the aim of accurately capturing the process of helium II reionization. In this paper, we discuss the observational signatures present in the He II Ly$\alpha$ forest. We show that the effective optical depth of the volume $\tau_\mathrm{eff}$ is not sufficient for capturing the ionization state of helium II, due to the large variance inherent in sightlines. However, the He II flux PDF can be used to determine the timing of helium II reionization. The amplitude of the one-dimensional flux power spectrum can also determine the ionization state of helium II. We show that even given the currently limited number of observations ($\sim$50 sightlines), measurements of the flux PDF can yield information about helium II reionization. Further, measurements using the one-dimensional power spectrum can provide clear indications of the timing of reionization, as well as the relative bias of sources of ionizing radiation.

We have run a new suite of simulations that solve hydrodynamics and radiative transfer simultaneously to study helium II reionization. Our suite of simulations employs various models for populating quasars inside of dark matter halos, which affect the He II reionization history. In particular, we are able to explore the impact that differences in the timing and duration of reionization have on observables. We examine the thermal signature that reionization leaves on the IGM, and measure the temperature-density relation. As previous studies have shown, we confirm that the photoheating feedback from helium II reionization raises the temperature of the IGM by several thousand kelvin. To compare against observations, we generate synthetic Ly$\alpha$ forest sightlines on-the-fly and match the observed effective optical depth $\tau_{\mathrm{eff}}(z)$ of hydrogen to recent observations. We show that when the simulations have been normalized to have the same values of $\tau_\mathrm{eff}$, the effect that helium II reionization has on observations of the hydrogen Ly$\alpha$ forest is minimal. Specifically, the flux PDF and the one-dimensional power spectrum are sensitive to the thermal state of the IGM, but do not show direct evidence for the ionization state of helium. We show that the peak temperature of the IGM typically corresponds to the time of 90%-95% helium ionization by volume, and is a relatively robust indicator of the timing of reionization. Future observations of helium reionization from the hydrogen Ly$\alpha$ forest should thus focus on measuring the temperature of the IGM, especially at mean density. Detecting the peak in the IGM temperature would provide valuable information about the timing of the end of helium II reionization.

We examine the statistics of the low-redshift Ly-alpha forest in an adaptive mesh refinement hydrodynamic cosmological simulation of sufficient volume to include distinct large-scale environments. We compare our HI column density distribution of absorbers both with recent work and between two highly-refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on scales of approximately 20 Mpc). We recover the average results presented in Kollmeier et al. (2014) using different simulation methods. We further break down these results as a function of environment to examine the detailed dependence of absorber statistics on large-scale density. We find that the slope of the HI column density distribution in the 10$^{12.5}$ $\le$ N$_{HI}$/cm$^{-2}$ $\le$ 10$^{14.5}$ range depends on environment such that the slope becomes steeper for higher environmental density, and this difference reflects distinct physical conditions of the intergalactic medium on these scales. We track this difference to the different temperature structures of filaments in varying environments. Specifically, filaments in the overdensity are hotter and, correspondingly, are composed of gas with lower HI fractions than those in underdense environments. Our results highlight that in order to understand the physics driving the HI CDD, we need not only improved accounting of the sources of ionizing UV photons, but also of the physical conditions of the IGM and how this may vary as a function of large-scale environment.

We systematically study, in the context of the standard cold dark matter model, star-formation suppression effects of two important known physical processes---photoheating due to reionization of the intergalactic medium and gravitational shock heating due to formation of massive halos and large-scale structure---on the global evolution of star formation rate (SFR) density and the so-called cosmic downsizing phenomenon in the redshift range z=0-6. We show that the steep decline of cosmic SFR density from z~2 to z=0 can be primarily explained by gravitational shock heating in two forms: massive halo self-quenching and hot environment. Simultaneously, we show a decreasing trend in the average SFR of star-forming galaxies from z=2 to z=0, reproducing the observed cosmic downsizing at z<2. Nevertheless, the average halo mass of star-forming galaxies is found to continue upsizing from z=2 to z=0. In stark contrast to z<2, both photoheating and gravitational shock heating effects are found to play a minor role in suppressing star formation. Additional negative feedback effects are required to reconcile our model with observations at z>2. Internal feedback from stellar evolution and supermassive black hole growth are the natural candidates for this role, as galaxies at z>2 are more moderate in mass but stronger in star formation and are thus more vulnerable. Our physical model can be used to treat star formation in cosmological N-body simulations.

There are well-observed differences between cool-core (CC) and non-cool-core (NCC) clusters, but the origin of this distinction is still largely unknown. Competing theories can be divided into internal (inside-out), in which internal physical processes transform or maintain the NCC phase, and external (outside-in), in which the cluster type is determined by its initial conditions, which in turn lead to different formation histories (i.e., assembly bias). We propose a new method that uses the relative assembly bias of CC to NCC clusters, as determined via the two-point cluster-galaxy cross-correlation function (CCF), to test whether formation history plays a role in determining their nature. We apply our method to 48 ACCEPT clusters, which have well resolved central entropies, and cross-correlate with the SDSS-III/BOSS LOWZ galaxy catalog. We find that the relative bias of NCC over CC clusters is $b = 1.42 \pm 0.35$ ($1.6\sigma$ different from unity). Our measurement is limited by the small number of clusters with core entropy information within the BOSS footprint, 14 CC and 34 NCC. Future compilations of X-ray cluster samples, combined with deep all-sky redshift surveys, will be able to better constrain the relative assembly bias of CC and NCC clusters and determine the origin of the bimodality.

An analysis of the physics-rich endgame of reionization at $z=5.7$ is performed, utilizing jointly the observations of the Ly$\alpha$ forest, the mean free path of ionizing photons, the luminosity function of galaxies and new physical insight. We find that an upper limit on ${\rm \tau_e}$ provides a constraint on the minimum mean free path (of ionizing photons) that is primarily due to dwarf galaxies, which in turn yields a new and yet the strongest constraint on the matter power spectrum on $10^6-10^9M_\odot$ scales. With the latest Planck measurements of ${\rm \tau_e = 0.055 \pm 0.009}$, we can place an upper limit of $(8.9\times 10^6, 3.8\times 10^7, 4.2\times 10^8)M_\odot$ on the lower cutoff mass of the halo mass function, or equivalent a lower limit on warm dark matter particle mass ${\rm m_x \ge (15.1, 9.8, 4.6)keV}$ or on sterile neutrino mass ${\rm m_s \ge (161, 90, 33)keV}$, at $(1, 1.4, 2.2)\sigma$ confidence level, respectively.

Utilizing high-resolution cosmological hydrodynamic simulations we investigate various ultra-violet absorption lines in the circumgalactic medium of star forming galaxies at low redshift, in hopes of checking and alleviating the claimed observational conundrum of the ratio of NV to OVI absorbers, among others. We find a satisfactory agreement between simulations and extant observational data with respect to the ratios of the following four line pairs examined, NV/OVI, SiIV/OVI, NIII/OVI and NII/OVI. For the pairs involving nitrogen lines, we examine two cases of nitrogen abundance, one with constant N/O ratio and the other with varying N/O ratio, with the latter motivated by theoretical considerations of two different synthetic sources of nitrogen that is empirically verified independently. Along a separate vector, for all line pairs, we examine two cases of radiation field, one with the Haardt-Madau background radiation field and the other with an additional local radiation field sourced by hot gas in the host galaxy. In all cases, two-sample Kolmogorov-Smirnov tests indicate excellent agreements. We find that the apparent agreements between simulations and observations will be strongly tested, if the bulk of current upper limits of various line ratios are turned into actual detections. We show that an increase in observational sensitivity by 0.2 dex will already start to significantly constrain the models.

Utilizing {\it ab initio} ultra-high resolution hydrodynamical simulations, we investigate the properties of the interstellar and circum-galactic medium of Ly$\alpha$ Blobs (LABs) at $z=3$, focusing on three important emission lines: Ly$\alpha$ 1216\AA, \heii 1640\AA\ and \civ 1449\AA. Their relative strengths provide a powerful probe of the thermodynamic properties of the gas when confronted with observations. By adjusting the dust attenuation effect using one parameter and matching the observed size-luminosity relation of LABs using another parameter, we show that our simulations can reproduce the observed \civ/\lya\ and \heii/\lya\ ratios adequately. This analysis provides the first successful physical model to account for simultaneously the LAB luminosity function, luminosity-size relation, and the \civ/Ly$\alpha$ and \heii/Ly$\alpha$ ratios, with only two parameters. The physical underpinning for this model is that, in addition to the stellar component for the \lya\ emission, the \lya\ and \civ\ emission lines due to shock heated gas are primarily collisional excitation driven and the \heii\ emission line collisional ionization driven. We find that the density, temperature and metallicity of the gas responsible for each emission line is significantly distinct, in a multi-phase interstellar and circumgalactic medium that is shock-heated primarily by supernovae and secondarily by gravitational accretion of gas.

While galaxies at $6 \lesssim z \lesssim 10$ are believed to dominate the epoch of cosmic reionization, the escape fraction of ionizing flux $f_\mathrm{esc}$ and the photon production rate $\dot n_\gamma$ from these galaxies must vary with redshift to simultaneously match CMB and low-redshift observations. We constrain $f_\mathrm{esc}(z)$ and $\dot n_\gamma(z)$ with Planck 2015 measurements of the Thomson optical depth $\tau$, recent low multipole E-mode polarization measurements from Planck 2016, SDSS BAO data, and $3 \lesssim z \lesssim 10$ galaxy observations. We compare different galaxy luminosity functions that are calibrated to HST observations, using both parametric and non-parametric statistical methods that marginalize over the effective clumping factor $C_\mathrm{HII}$, the LyC production efficiency $\xi_\mathrm{ion}$, and the time-evolution of the UV limiting magnitude $dM_\mathrm{SF}/dz$. Using a power-law model, we find $f_\mathrm{esc} \lesssim 0.5$ at $z=8$ with slope $\beta \gtrsim 2.0$ at $68\%$ confidence with little dependence on the galaxy luminosity function or data, although there is non-negligible probability for no redshift evolution $\beta \sim 0$ or small escape fraction $f_\mathrm{esc} \sim 10^{-2}$. A non-parametric form for $f_\mathrm{esc}(z)$ evolves significantly with redshift, yielding $f_\mathrm{esc} \sim 0.2, 0.3, 0.6$ at $z=6,9,12$, respectively. However, a model-independent reconstruction of $\dot n_\gamma(z)$ predicts a suppressed escaped photon production rate at $z=9$ for the latest Planck data compared to the other models, implying a quicker period of reionization. We find evidence for redshift evolution in the limiting magnitude of the galaxy luminosity function for empirical models of the galaxy luminosity function.

We perform a statistical analysis of strong gravitational lensing by quasar hosts of background galaxies, in the two competing models of dark matter halos of quasars, HOD and CS models. Utilizing the BolshoiP Simulation we demonstrate that strong gravitational lensing provides a potentially very powerful test of models of quasar hosting halos. For quasars at $z=0.5$, the lensing probability by quasars of background galaxies in the HOD model is higher than that of the CS model by two orders of magnitude or more for lensing image separations in the range of $\theta\sim 1.2-12~$arcsec. To observationally test this, we show that, as an example, at the depth of the CANDELS wide field survey and with a quasar sample of $1000$ at $z=0.5$, the two models can be differentiated at $3-4\sigma$ confidence level.

An analysis of negative radiative feedback from resident stars in minihalos is performed. It is found that the most effective mechanism to suppress star formation is provided by infrared photons from resident stars via photo-detachment of ${\rm H^-}$. It is shown that a stringent upper bound on (total stellar mass, metallicity) of ($\sim 1000{\rm M_\odot}$, $-3.3\pm 0.2$) in any newly minted atomic cooling halo can be placed, with the actual values possibly significantly lower. This has both important physical ramifications on formation of stars and supermassive black seeds in atomic cooling halos at high redshift, pertaining to processes of low temperature metal cooling, dust formation and fragmentation, and direct consequences on the faint end galaxy luminosity function at high redshift and cosmological reionization. The luminosity function of galaxies at the epoch of reionization may be substantially affected due to the combined effect of a diminished role of minihalos and an enhanced contribution from Pop III stars in atomic cooling halos. Upcoming results on reionization optical depth from Planck High-Frequency Instrument data may provide a significant constraint on and a unique probe of this star formation physical process in minihalos. As a numerical example, in the absence of significant contributions from minihalos with virial masses below $1.5\times 10^{8}{\rm M_\odot}$ the reionization optical depth is expected to be no greater than $0.065$, whereas allowing for minihalos of masses as low as ($10^7{\rm M_\odot}$, $10^{6.5}{\rm M_\odot}$) to form stars unconstrained by this self-regulation physical process, the reionization optical depth is expected to exceed $(0.075,0.085)$, respectively.

We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with $M_{\rm halo} \sim4\times10^7 M_\odot$ at $z>10$ using cosmological radiation-hydrodynamics simulations. We find that very compact ($\lesssim$ 1 pc) and massive ($\sim6\times10^5 M_\odot$) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Ly$\alpha$ emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short ($\ll 1\,{\rm Myr}$), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. The simulated clusters closely resemble the local GCs in mass and size but show a metallicity spread that is much wider than found in the local GCs. We discuss a role of pre-enrichment by Pop III and II stars as a potential solution to the latter issue. Although not without shortcomings, it is encouraging that a naive blind (not tuned) cosmological simulation presents a possible channel for the formation of at least some massive GCs.

Utilizing the high-resolution, large-scale LAOZI cosmological simulations we investigate the nature of the metal-poor (${\rm [Z/H]<-2}$) damped Lyman alpha systems (mpDLA) at $z=3$. The following physical picture of mpDLAs emerges. The majority of mpDLAs inhabit regions $\ge 20$~kpc from the host galaxy center on infalling cold gas streams originating from the intergalactic medium, with infall velocity of $\sim 100$ km/s and temperature of $\sim 10^{4}$ K. For each host galaxy, on average, about $1\%$ of the area within a radius $150$~kpc is covered by mpDLAs. The mpDLAs are relatively diffuse ($n_{\rm{gas}} \sim 10^{-2}$ cm$^{-3}$), Jeans quasi-stable, and have very low star formation rate ($\dot{\Sigma} \le 10^{-4} M_{\odot} \rm{\ yr}^{-1} \rm{\ kpc}^{-2}$). As mpDLAs migrate inward to the galaxy center, they mix with high metallicity gas and stellar outflows in the process, removing themselves from the metal-poor category and rendering the central ($\le 5$ kpc) regions of galaxies devoid of mpDLAs. Thus, the central regions of the host galaxies are populated by mostly metal-rich DLAs instead of mpDLAs. All observables of the simulated mpDLAs are in excellent agreement with observations, except the gas density, which is about a factor of ten lower than the value inferred observationally. However, the observationally inferred value is based on simplified assumptions that are not borne out in the simulations.

Three-dimensional hydrodynamic simulations, covering the spatial domain from hundreds of Schwarzschild radii to $2\ \mathrm{pc}$ around the central supermassive black hole of mass $10^8 M_\odot$, with detailed radiative cooling processes, are performed. Generically found is the existence of a significant amount of shock heated, high temperature ($\geq 10^8\ \mathrm{K}$) coronal gas in the inner ($\leq 10^4 r_\mathrm{sch}$) region. It is shown that the composite bremsstrahlung emission spectrum due to coronal gas of various temperatures are in reasonable agreement with the overall ensemble spectrum of AGNs and hard X-ray background. Taking into account inverse Compton processes, in the context of the simulation-produced coronal gas, our model can readily account for the wide variety of AGN spectral shape, which can now be understood physically. The distinguishing feature of our model is that X-ray coronal gas is, for the first time, an integral part of the inflow gas and its observable characteristics are physically coupled to the concomitant inflow gas. One natural prediction of our model is the anti-correlation between accretion disk luminosity and spectral hardness: as the luminosity of SMBH accretion disk decreases, the hard X-ray luminosity increases relative to the UV/optical luminosity.

SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $\dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $\approx 20\%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $\dot{M}$, but reaches a maximum value at a characteristic mass $\sim 2 \times 10^{11}\ M_\odot$ and a characteristic accretion rate $\sim 6 \times 10^2\ M_\odot/{\rm yr}$ at $z \approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $\phi_\star$ decreases, $M_\star$ is more positive (fainter), and $\alpha$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{\rm AB} \approx 31\ (32)$, it can observe $\gtrsim 11\ (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_\star$ at $z \lesssim 13\ (14)$.

The Lyman alpha forest power spectrum has been measured on large scales by the BOSS survey in SDSS-III at $z\sim 2.3$, has been shown to agree well with linear theory predictions, and has provided the first measurement of Baryon Acoustic Oscillations at this redshift. However, the power at small scales, affected by non-linearities, has not been well examined so far. We present results from a variety of hydrodynamic simulations to predict the redshift space non-linear power spectrum of the Lyman Alpha transmission for several models, testing the dependence on resolution and box size. A new fitting formula is introduced to facilitate the comparison of our simulation results with observations and other simulations. The non-linear power spectrum has a generic shape determined by a transition scale from linear to non-linear anisotropy, and a Jeans scale below which the power drops rapidly. In addition, we predict the two linear bias factors of the Lyman Alpha forest and provide a better physical interpretation of their values and redshift evolution. The dependence of these bias factors and the non-linear power on the amplitude and slope of the primordial fluctuations power spectrum, the temperature-density relation of the intergalactic medium, and the mean Lyman Alpha transmission, as well as the redshift evolution, is investigated and discussed in detail. A preliminary comparison to the observations shows that the predicted redshift distortion parameter is in good agreement with the recent determination of Blomqvist et al., but the density bias factor is lower than observed. We make all our results publicly available in the form of tables of the non-linear power spectrum that is directly obtained from all our simulations, and parameters of our fitting formula.

Supernovae (SN), the most energetic stellar feedback mechanism, are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing 3D hydrodynamical simulations at various SN rates, $S$, and ISM average densities, $\bar{n}$. The evolution of a SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By travelling faster and further in the low-density hot phase, the domain of a SNR increases by $>10^{2.5}$. Varying $\bar{n}$ and $S$, we find that a steady state can only be achieved when the hot gas volume fraction $f_{\rm{V,hot}}\lesssim 0.6 \pm 0.1 $. Above that level, overlapping SNRs render connecting topology of the hot gas, and the ISM is subjected to thermal runaway. Photoelectric heating (PEH) has a surprisingly strong impact on $f_{\rm{V,hot}}$. For $\bar{n}\gtrsim 3 \cm-3 $, a reasonable PEH rate is able to suppress the thermal runaway. Overall, we determine the critical SN rate for the onset of thermal runaway to be $S_{\rm{crit}} = 200 (\bar{n}/1\cm-3)^k (E_{\rm{SN}}/10^{51}\erg)^{-1} \kpc^{-3} \myr-1$, where $k = (1.2,2.7)$ for $\bar{n} \leq 1$ and $> 1\cm-3 $, respectively. We present a fitting formula of the ISM pressure $P(\bar{n}$, $S$), which can be used as an effective equation of state in cosmological simulations. Despite the 5 orders of magnitude span of $(\bar{n},S)$, the average Mach number varies little: $\mathcal{M} \approx \ 0.5\pm 0.2, \ 1.2\pm 0.3,\ 2.3\pm 0.9$ for the hot, warm and cold phases, respectively.

The connection between dark matter halos and galactic baryons is often not well-constrained nor well-resolved in cosmological hydrodynamical simulations. Thus, Halo Occupation Distribution (HOD) models that assign galaxies to halos based on halo mass are frequently used to interpret clustering observations, even though it is well-known that the assembly history of dark matter halos is related to their clustering. In this paper we use high-resolution hydrodynamical cosmological simulations to compare the halo and stellar mass growth of galaxies in a large-scale overdensity to those in a large-scale underdensity (on scales of about 20 Mpc). The simulation reproduces assembly bias, that halos have earlier formation times in overdense environments than in underdense regions. We find that the stellar mass to halo mass ratio is larger in overdense regions in central galaxies residing in halos with masses between 10$^{11}$-10$^{12.9}$ M$_{\odot}$. When we force the local density (within 2 Mpc) at z=0 to be the same for galaxies in the large-scale over- and underdensities, we find the same results. We posit that this difference can be explained by a combination of earlier formation times, more interactions at early times with neighbors, and more filaments feeding galaxies in overdense regions. This result puts the standard practice of assigning stellar mass to halos based only on their mass, rather than considering their larger environment, into question.

A statistical analysis of stacked Compton$-y$ maps of quasar hosts with a median redshift of $1.5$ using Millennium Simulation is performed to address two issues, one on the feedback energy from quasars and the other on testing dark matter halo models for quasar hosts. On the first, we find that, at the resolution of FWHM=$10$ arcmin obtained by Planck data, the observed thermal Sunyaev-Zeldovich (tSZ) effect can be entirely accounted for and explained by the thermal energy of halos sourced by gravitational collapse of halos, without a need to invoke additional, large energy sources, such as quasar or stellar feedback. Allowing for uncertainties of dust temperature in the calibration of observed Comton$-y$ maps, the maximum additional feedback energy is $\sim 25\%$ of that previously suggested. Second, we show that, with FWHM=$1$ arcmin beam, tSZ measurements will provide a potentially powerful test of quasar-hosting dark matter halo models, limited only by possible observational systematic uncertainties, not by statistical ones, even in the presence of possible quasar feedback.

Lyman alpha blobs (LABs) are spatially extended lyman alpha nebulae seen at high redshift. The origin of Lyman alpha emission in the LABs is still unclear and under debate. To study their heating mechanism(s), we present Australia Telescope Compact Array (ATCA) observations of the 20 cm radio emission and Herschel PACS and SPIRE measurements of the far-infrared (FIR) emission towards the four LABs in the protocluster J2143-4423 at z=2.38. Among the four LABs, B6 and B7 are detected in the radio with fluxes of 67+/-17 microJy and 77+/-16 microJy, respectively, and B5 is marginally detected at 3 sigma (51+/-16 microJy). For all detected sources, their radio positions are consistent with the central positions of the LABs. B6 and B7 are obviously also detected in the FIR. By fitting the data with different templates, we obtained redshifts of 2.20$^{+0.30}_{-0.35}$ for B6 and 2.20$^{+0.45}_{-0.30}$ for B7 which are consistent with the redshift of the lyman alpha emission within uncertainties, indicating that both FIR sources are likely associated with the LABs. The associated FIR emission in B6 and B7 and high star formation rates strongly favor star formation in galaxies as an important powering source for the lyman alpha emission in both LABs. However, the other two, B1 and B5, are predominantly driven by the active galactic nuclei or other sources of energy still to be specified, but not mainly by star formation. In general, the LABs are powered by quite diverse sources of energy.

Feedback from supernovae is essential to understanding the self-regulation of star formation in galaxies. However, the efficacy of the process in a cosmological context remains unclear due to excessive radiative losses during the shock propagation. To better understand the impact of SN explosions on the evolution of galaxies, we perform a suite of high-resolution (12 pc), zoom-in cosmological simulations of a Milky Way-like galaxy at z=3 with adaptive mesh refinement. We find that SN explosions can efficiently regulate star formation, leading to the stellar mass and metallicity consistent with the observed mass-metallicity relation and stellar mass-halo mass relation at z~3. This is achieved by making three important changes to the classical feedback scheme: i) the different phases of SN blast waves are modelled directly by injecting radial momentum expected at each stage, ii) the realistic time delay of SNe, commencing at as early as 3 Myr, is required to disperse very dense gas before a runaway collapse sets in at the galaxy centre via mergers of gas clumps, and iii) a non-uniform density distribution of the ISM is taken into account below the computational grid scale for the cell in which SN explodes. The last condition is motivated by the fact that our simulations still do not resolve the detailed structure of a turbulent ISM in which the fast outflows can propagate along low-density channels. The simulated galaxy with the SN feedback model shows strong outflows, which carry approximately ten times larger mass than star formation rate, as well as smoothly rising circular velocity. Other feedback models that do not meet the three conditions form too many stars, producing a peaked rotation curve. Our results suggest that understanding the structure of the turbulent ISM may be crucial to assess the role of SN and other feedback processes in galaxy formation theory.

We reason that, without physical fine-tuning, neither the supermassive black holes (SMBHs) nor the stellar bulges can self-regulate or inter-regulate by driving away already fallen cold gas to produce the observed correlation between them. We suggest an alternative scenario where the observed mass ratios of the SMBHs to bulges reflect the angular momentum distribution of infallen gas such that the mass reaching the stable accretion disc is a small fraction of that reaching the bulge region, averaged over the cosmological time scales. We test this scenario using high resolution, large-scale cosmological hydrodynamic simulations (without AGN feedback), assuming the angular momentum distribution of gas landing in the bulge region to yield a Mestel disc that is supported by independent simulations resolving the Bondi radii of SMBHs. A mass ratio of $0.1-0.3\%$ between the very low angular momentum gas that free-falls to the sub-parsec region to accrete to the SMBH and the overall star formation rate is found. This ratio is found to increase with increasing redshift to within a factor of $\sim 2$, suggesting that the SMBH to bulge ratio is nearly redshift independent, with a modest increase with redshift, a testable prediction. Furthermore, the duty cycle of active galactic nuclei (AGN) with high Eddington ratios is expected to increase significantly with redshift. Finally, while SMBHs and bulges are found to coevolve on $\sim 30-150$Myr time scales or longer, there is indication that, on shorer time scales, the SMBH accretion rate and star formation may be less correlated.

Ly$\alpha$ photons scattered by neutral hydrogen atoms in the circumgalactic media or produced in the halos of star-forming galaxies are expected to lead to extended Ly$\alpha$ emission around galaxies. Such low surface brightness Ly$\alpha$ halos (LAHs) have been detected by stacking Ly$\alpha$ images of high-redshift star-forming galaxies. We study the origin of LAHs by performing radiative transfer modeling of nine $z=3.1$ Lyman-Alpha Emitters (LAEs) in a high resolution hydrodynamic cosmological galaxy formation simulation. We develop a method of computing the mean Ly$\alpha$ surface brightness profile of each LAE by effectively integrating over many different observing directions. Without adjusting any parameters, our model yields an average Ly$\alpha$ surface brightness profile in remarkable agreement with observations. We find that observed LAHs cannot be accounted for solely by photons originating from the central LAE and scattered to large radii by hydrogen atoms in the circumgalactic gas. Instead, Ly$\alpha$ emission from regions in the outer halo is primarily responsible for producing the extended LAHs seen in observations, which potentially includes both star-forming and cooling radiation. With the limit on the star formation contribution set by the ultra-violet (UV) halo measurement, we find that cooling radiation can play an important role in forming the extended LAHs. We discuss the implications and caveats of such a picture.

Simulations have indicated that most of the escaped Lyman continuum photons escape through a minority of solid angles with near complete transparency, with the remaining majority of the solid angles largely opaque, resulting in a very broad and skewed probability distribution function (PDF) of the escape fraction when viewed at different angles. Thus, the escape fraction of Lyman continuum photons of a galaxy observed along a line of sight merely represents the properties of the interstellar medium along that line of sight, which may be an ill-representation of true escape fraction of the galaxy averaged over its full sky. Here we study how Lyman continuum photons escape from galaxies at $z=4-6$, utilizing high-resolution large-scale cosmological radiation-hydrodynamic simulations. We compute the PDF of the mean escape fraction ($\left<f_{\rm esc,1D}\right>$) averaged over mock observational samples, as a function of the sample size, compared to the true mean (had you an infinite sample size). We find that, when the sample size is small, the apparent mean skews to the low end. For example, for a true mean of 6.7%, an observational sample of (2,10,50) galaxies at $z=4$ would have have 2.5% probability of obtaining the sample mean lower than $\left<f_{\rm esc,1D}\right>=$(0.007%, 1.8%, 4.1%) and 2.5% probability of obtaining the sample mean being greater than (43%, 18%, 11%). Our simulations suggest that at least $\sim$ 100 galaxies should be stacked in order to constrain the true escape fraction within 20% uncertainty.

Morphology of the complex HI gas distribution can be quantified by statistics like the Minkowski functionals, and can provide a way to statistically study the large scale structure in the HI maps both at low redshifts, and during the epoch of reionization (EoR). At low redshifts, the 21cm emission traces the underlying matter distribution. Topology of the HI gas distribution, as measured by the genus, could be used as a "standard ruler". This enables the determination of distance-redshift relation and also the discrimination of various models of dark energy and of modified gravity. The topological analysis is also sensitive to certain primordial non-Gaussian features. Compared with two-point statistics, the topological statistics are more robust against the nonlinear gravitational evolution, bias, and redshift-space distortion. The HI intensity map observation naturally avoids the sparse sampling distortion, which is an important systematic in optical galaxy survey. The large cosmic volume accessible to SKA would provide unprecedented accuracy using such a measurement... [abridged]

Accretion is thought to primarily contribute to the mass accumulation history of supermassive black holes throughout cosmic time. While this may be true at high redshifts, at lower redshifts and for the most massive black holes mergers themselves might add significantly to the mass budget. We evolve SMBHs from $4 > z > 0$ using merger trees derived from hydrodynamical cosmological simulations of a cluster and void region, scaled to the observed value of the stellar mass fraction to account for overcooling. Mass gains from gas accretion proportional to bulge growth and BH-BH mergers are tracked, as are black holes that remain "orbiting" due to insufficient dynamical friction in a merger remnant, as well as those that are ejected due to gravitational recoil. We find that gas accretion remains the dominant source of mass accumulation in almost all SMBHs; mergers contribute $2.5\pm0.1\%$ for all SMBHs in the cluster and $1.0\pm0.1\%$ in the void since $z = 4$. However, mergers are significant for massive SMBHs. The fraction of mass accumulated from mergers for central BHs generally increases for larger values of the host bulge mass: in the void, the fraction is $2\%$ at $M_{*, bul} = 10^{10} M_{\odot}$, increasing to $4\%$ at $M_{*, bul} \gtrsim 10^{11} M_{\odot}$, and in the cluster it is $4\%$ at $M_{*, bul} = 10^{10} M_{\odot}$ and $23\%$ at $10^{12} M_{\odot}$. We find that $40\%$ of SMBHs and $\approx 8\%$ of the total SMBH mass is found orbiting in the cluster region at $z = 0$. The existence of orbiting and ejected SMBHs requires modification of the Soltan argument. We estimate this correction to the integrated accreted mass density of SMBHs to be in the range $6-21\%$, with a mean value of $11\pm3\%$. We also calculate the total energy output and strain from gravitational waves emitted by merging SMBHs, and obtain a signal potentially detectable by pulsar timing arrays.

A new model for quasar-hosting dark matter halos, meeting two physical conditions, is put forth. First, significant interactions are taken into consideration to trigger quasar activities. Second, satellites in very massive halos at low redshift are removed from consideration, due to their deficiency of cold gas. We analyze the {\em Millennium Simulation} to find halos that meet these two conditions and simultaneously match two-point auto-correlation functions of quasars and cross-correlation functions between quasars and galaxies at $z=0.5-3.2$. %The found halos have some distinct properties worth noting. The masses of found quasar hosts decrease with decreasing redshift, with the mass thresholds being $[(2-5)\times 10^{12}, (2-5)\times 10^{11}, (1-3)\times 10^{11}]\msun$ for median luminosities of $\sim[10^{46}, 10^{46}, 10^{45}]$erg/s at $z=(3.2, 1.4, 0.53)$, respectively, an order of magnitude lower than those inferred based on halo occupation distribution modeling. In this model quasar hosts are primarily massive central halos at $z\ge 2-3$ but increasingly dominated by lower mass satellite halos experiencing major interactions towards lower redshift. But below $z=1$ satellite halos in groups more massive than $\sim 2\times 10^{13}\msun$ do not host quasars. Whether for central or satellite halos, imposing the condition of significant interactions substantially boosts the clustering strength compared to the total population with the same mass cut. The inferred lifetimes of quasars at $z=0.5-3.2$ of $3-30$Myr are in agreement with observations. Quasars at $z\sim 2$ would be hosted by halos of mass $\sim 5\times 10^{11}\msun$ in this model, compared to $\sim 3\times 10^{12}\msun$ previously thought, which would help reconcile with the observed, otherwise puzzling high covering fractions for Lyman limit systems around quasars.

Analysis is performed on ultra-high resolution large-scale cosmological radiation-hydrodynamic simulations to, for the first time, quantify the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high density ($\ge 10^4$cm$^{-3}$) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density ($\sim 10^{-2.5}$cm$^{-3}$) high temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature and metallicity of intervening absorbers, vary strongly from sightline to sightline. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium of the host galaxy based on a very small number of (usually one) sightlines would be precarious. Utilizing high-z GRBs to probe interstellar medium and intergalactic medium should be undertaken properly taking into consideration the physical diversities of the interstellar medium.

We present predictions for the 21 cm brightness temperature power spectrum during the Epoch of Reionization (EoR). We discuss the implications of the "light cone" effect, which incorporates evolution of the neutral hydrogen fraction and 21 cm brightness temperature along the line of sight. Using a novel method calibrated against radiation-hydrodynamic simulations, we model the neutral hydrogen density field and 21 cm signal in large volumes ($L = 2$ Gpc/$h$). The inclusion of the light cone effect leads to a relative decrease of about 50% in the 21 cm power spectrum on all scales. We also find that the effect is more prominent at the midpoint of reionization and later. The light cone effect also can introduce an anisotropy along the line of sight. By decomposing the 3D power spectrum into components perpendicular to and along the line of sight, we find that in our fiducial reionization model, there is no significant anisotropy. However, parallel modes can contribute up to 40% more power for shorter reionization scenarios. The scales on which the light cone effect is relevant are comparable to scales where one measures the baryon acoustic oscillation. We argue that due to its large comoving scale and introduction of anisotropy, the light cone effect is important when considering redshift space distortions and future application to the Alcock-Paczynski test for the determination of cosmological parameters.

While the Sersic profile family provide adequate fits for the surface brightness profiles of observed galaxies, the physical origin is unknown. We show that, if the cosmological density field are seeded by random gaussian fluctuations, as in the standard cold dark matter model, galaxies with steep central profiles have simultaneously extended envelopes of shallow profiles in the outskirts, whereas galaxies with shallow central profiles are accompanied by steep density profiles in the outskirts. These properties are in accord with those of the Sersic profile family. Moreover, galaxies with steep central profiles form their central regions in smaller denser subunits that possibly merge subsequently, which naturally leads to formation of bulges. In contrast, galaxies with shallow central profiles form their central regions in a coherent fashion without significant substructure, a necessary condition for disk galaxy formation. Thus, the scenario is self-consistent with respect to the correlation between observed galaxy morphology and Sersic index. We predict further that clusters of galaxies should display a similar trend, which should be verifiable observationally.

Recent large surveys have found a reversal of the star formation rate (SFR)-density relation at z=1 from that at z=0 (e.g. Elbaz et al.; Cooper et al.), while the sign of the slope of the color-density relation remains unchanged (e.g. Cucciati et al.; Quadri et al.). We use state-of-the-art adaptive mesh refinement cosmological hydrodynamic simulations of a 21x24x20 (Mpc/h)$^3$ region centered on a cluster to examine the SFR-density and color-density relations of galaxies at z=0 and z=1. The local environmental density is defined by the dark matter mass in spheres of radius 1 Mpc/h, and we probe two decades of environmental densities. Our simulations produce a large increase of SFR with density at z=1, as in the observations of Elbaz et al. We also find a significant evolution to z=0, where the SFR-density relation is much flatter. The color-density relation in our simulations is consistent from z=1 to z=0, in agreement with observations. We find that the increase in the median SFR with local density at z=1 is due to a growing population of star-forming galaxies in higher-density environments. At z=0 and z=1 both the SFR and cold gas mass are tightly correlated with the galaxy halo mass, and therefore the correlation between median halo mass and local density is an important cause of the SFR-density relation at both redshifts. We also show that the local density on 1 Mpc/h scales affects galaxy SFRs as much as halo mass at z=0. Finally, we find indications that the role of the 1 Mpc/h scale environment reverses from z=0 to z=1: at z=0 high-density environments depress galaxy SFRs, while at z=1 high-density environments tend to increase SFRs.

We use high-resolution cosmological hydrodynamic galaxy formation simulations to gain insights on how galaxies lose their cold gas at low redshift as they migrate from the field to the high density regions of clusters of galaxies. We find that beyond three cluster virial radii, the fraction of gas-rich galaxies is constant, representing the field. Within three cluster-centric radii, the fraction of gas-rich galaxies declines steadily with decreasing radius, reaching $\mathbf{<10%}$ near the cluster center. Our results suggest that galaxies start to feel the impact of the cluster environment on their gas content well beyond the cluster virial radius. We show that almost all gas-rich galaxies at the cluster virial radius are falling in for the first time at nearly radial orbits. Furthermore, we find that almost no galaxy moving outward at the cluster virial radius is gas-rich (with gas to baryon ratio greater than $\mathbf{1%}$). These results suggest that galaxies that fall into clusters lose their cold gas within a single radial round-trip.

A new physical framework for the emergence of the Hubble sequence is outlined, based on novel analyses performed to quantify the evolution of cold streams of a large sample of galaxies from a state-of-the-art ultra-high resolution, large-scale adaptive mesh-refinement hydrodynamic simulation in a fully cosmological setting. It is found that the following three key physical variables of galactic cold inflows crossing the virial sphere substantially decrease with decreasing redshift: the number of streams N_{90} that make up 90% of concurrent inflow mass flux, average inflow rate per stream dot M_{90} and mean (mass flux weighted) gas density in the streams n_{gas}. Another key variable, the stream dimensionless angular momentum parameter lambda, instead is found to increase with decreasing redshift. Assimilating these trends and others leads naturally to a physically coherent scenario for the emergence of the Hubble sequence, including the following expectations: (1) the predominance of a mixture of disproportionately small irregular and complex disk galaxies at z>2 when most galaxies have multiple concurrent streams, (2) the beginning of the appearance of flocculent spirals at z~1-2 when the number of concurrent streams are about 2-3, (3) the grand-design spiral galaxies appear at z<1 when galaxies with only one major cold stream significantly emerge. These expected general trends are in good accord with observations. Early type galaxies are those that have entered a perennial state of zero cold gas stream, with their abundance increasing with decreasing redshift.

The fraction of hydrogen ionizing photons escaping from galaxies into the intergalactic medium is a critical ingredient in the theory of reionization. We use two zoomed-in, high-resolution (4 pc), cosmological radiation hydrodynamic simulations with adaptive mesh refinement to investigate the impact of two physical mechanisms (supernova feedback and runaway OB stars) on the escape fraction (f_esc) at the epoch of reionization (z>7). We implement a new, physically motivated supernova feedback model that can approximate the Sedov solutions at all (from the free expansion to snowplow) stages. We find that there is a significant time delay of about ten million years between the peak of star formation and that of escape fraction, due to the time required for the build-up and subsequent destruction of the star-forming cloud by supernova feedback. Consequently, the photon number-weighted mean escape fraction for dwarf galaxies in halos of mass 10^8-10^10.5 Msun is found to be <fesc>~11%, although instantaneous values of f_esc>20% are common when star formation is strongly modulated by the supernova explosions. We find that the inclusion of runaway OB stars increases the mean escape fraction by 22% to <fesc>~14%. As supernovae resulting from runaway OB stars tend to occur in less dense environments, the feedback effect is enhanced and star formation is further suppressed in halos with Mvir>10^9 Msun in the simulation with runaway OB stars compared with the model without them. While both our models produce enough ionizing photons to maintain a fully ionized universe at z>7 as observed, a still higher amount of ionizing photons at z>9 appears necessary to accommodate the high observed electron optical depth inferred from cosmic microwave background observations.

We study the evolution of angular momenta of ($M_*=10^{10}-10^{12}\msun$) galaxies utilizing large-scale ultra-high resolution cosmological hydrodynamic simulations and find that spin of the stellar component changes direction frequently, caused by major mergers, minor mergers, significant gas inflows and torques by nearby systems. The rate and nature of change of spin direction can not be accounted for by large-scale tidal torques, because the latter fall short in rates by orders of magnitude and because the apparent random swings of the spin direction are inconsistent with alignment by linear density field. The implications for galaxy formation as well as intrinsic alignment of galaxies are profound. Assuming the large-scale tidal field is the sole alignment agent, a new picture emerging is that intrinsic alignment of galaxies would be a balance between slow large-scale coherent torquing and fast spin reorientation by local interactions. What is still open is whether other processes, such as feeding galaxies with gas and stars along filaments or sheets, introduce coherence for spin directions of galaxies along the respective structures.

We report on the discovery of a relation between the number of star formation (SF) peaks per unit time, $\nu_{\rm peak}$, and the size of the temporal smoothing window function, $\Delta t$, used to define the peaks: $\nu_{\rm peak}\propto\Delta t^{1-\phi}$ ($\phi\sim 1.618$). This relation holds over the range of $\Delta t=10$ to $1000$Myr that can be reliably computed, using a large sample of galaxies obtained from a state-of-the-art cosmological hydrodynamic simulation. This means that the temporal distribution of SF peaks in galaxies as a population is fractal with a Hausdorff fractal dimension equal to $\phi-1$. This finding reveals, for the first time, that the superficially chaotic process of galaxy formation is underlined by a temporal self-organization up to at least one gigayear. It is tempting to suggest that, given the known existence of spatial fractals (such as the power-law two-point function of galaxies), there is a joint spatio-temporal self-organization in galaxy formation. From an observational perspective, it will be urgent to devise diagnostics to probe SF histories of galaxies with good temporal resolution to facilitate a test of this prediction. If confirmed, it would provide unambiguous evidence for a new picture of galaxy formation that is interaction driven, cooperative and coherent in and between time and space. Unravelling its origin may hold the key to understanding galaxy formation.

The outskirts of galaxy clusters are continuously disturbed by mergers and gas infall along filaments, which in turn induce turbulent flow motions and shock waves. We examine the properties of shocks that form within $r_{200}$ in sample galaxy clusters from structure formation simulations. While most of these shocks are weak and inefficient accelerators of cosmic rays (CRs), there are a number of strong, energetic shocks which can produce large amounts of CR protons via diffusive shock acceleration. We show that the energetic shocks reside mostly in the outskirts and a substantial fraction of them are induced by infall of the warm-hot intergalactic medium from filaments. As a result, the radial profile of the CR pressure in the intracluster medium is expected to be broad, dropping off more slowly than that of the gas pressure, and might be even temporarily inverted, peaking in the outskirts. The volume-integrated momentum spectrum of CR protons inside $r_{200}$ has the power-law slope of $4.25 - 4.5$, indicating that the average Mach number of the shocks of main CR production is in the range of $\left< M_s \right>_{\rm CR} \approx 3 - 4$. We suggest that some radio relics with relatively flat radio spectrum could be explained by primary electrons accelerated by energetic infall shocks with $M_s >~ 3$ induced in the cluster outskirts.

If a supermassive black hole has some material orbiting around it at close to its innermost stable circular orbit (ISCO), then, when it plunges into a second supermassive black hole, the orbiting material has a velocity dispersion of order of speed of light about the orbital velocity of its host black hole. It becomes plausible that some of the orbiting material will be "catapulted" to the negative-energy ergosphere orbits of the second black hole at the plunge. This may provide an astrophysically plausible way to extract energy from the black hole, originally suggested by Penrose.

An analysis of more than 3000 galaxies resolved at better than 114 pc/h at z=0.62 in a LAOZI cosmological adaptive mesh refinement hydrodynamic simulation is performed and insights gained on star formation quenching and color migration. The vast majority of red galaxies are found to be within three virial radii of a larger galaxy, at the onset of quenching when the specific star formation rate experiences the sharpest decline to fall below ~10^{-2}-10^{-1}/Gyr (depending on the redshift). We shall thus call this mechanism "environment quenching", which encompasses satellite quenching. Two physical processes are largely responsible: ram-pressure stripping first disconnects the galaxy from the cold gas supply on large scales, followed by a longer period of cold gas starvation taking place in high velocity dispersion environment, during the early part of which the existing dense cold gas in the central region (<10kpc) is consumed by in situ star formation. Quenching is found to be more efficient, but not faster, on average, in denser environment. Throughout this quenching period and the ensuing one in the red sequence galaxies follow nearly vertical tracks in the color-stellar-mass diagram. In contrast, individual galaxies of all masses grow most of their stellar masses in the blue cloud, prior to the onset of quenching, and progressively more massive blue galaxies with already relatively older mean stellar ages continue to enter the red sequence. Consequently, correlations among observables of red galaxies - such as the age-mass relation - are largely inherited from their blue progenitors at the onset of quenching. While the color makeup of the entire galaxy population strongly depends on environment, which is a direct result of environment quenching, physical properties of blue galaxies as a sub-population show little dependence on environment.

Three-dimensional panchromatic dust radiative transfer calculations are performed on a set of 198 galaxies of stellar masses in the range 5x10^8-3x10^10 Msun from a cosmological hydrodynamic simulation (resolved at 29pc/h) at z=7. In a companion paper (Kimm & Cen), the stellar mass and UV luminosity functions, and UV-optical and FUV-NUV colors are shown to be in good agreement with observations, if an SMC-type dust extinction curve is adopted. Here we make useful predictions, self-consistently, of the infrared properties of these z=7 simulated galaxies that can be confronted with upcoming ALMA data. Our findings are as follows. (1) The effective radius in the rest-frame MIPS 70 micron band is in the range of 80-400pc proper for z=7 galaxies with L_FIR=10^{11.3-12}Lsun. (2) The median of the peak wavelength of the far-infrared (FIR) spectral energy distribution is in the range of 45-60 micron, depending on the dust-to-metal ratio. (3) For star formation rate in the range 3-100 Msun/yr the median FIR to bolometric luminosity ratio is 60-90%. (4) The FIR luminosity function displays a power law in the high end with a slope of -3.1 +- 0.4, instead of the usual exponential decline.

Hubble Space Telescope observations with Wide Field Camera 3/IR reveal that galaxies at z~7 have very blue ultraviolet (UV) colors, consistent with these systems being dominated by young stellar populations with moderate or little attenuation by dust. We investigate UV and optical properties of the high-z galaxies in the standard cold dark matter model using a high-resolution adaptive mesh refinement cosmological hydrodynamic simulation. For this purpose, we perform panchromatic three-dimensional dust radiative transfer calculations on 198 galaxies of stellar mass 5x10^8-3x10^{10} Msun with three parameters, the dust-to-metal ratio, the extinction curve, and the fraction of directly escaped light from stars (\fesc). Our stellar mass function is found to be in broad agreement with Gonzalez et al., independent of these parameters. We find that our heavily dust-attenuated galaxies (A_V~1.8) can also reasonably match modest UV-optical colors, blue UV slopes, as well as UV luminosity functions, provided that a significant fraction (~10%) of light directly escapes from them. The observed UV slope and scatter are better explained with an Small Magellanic Cloud-type extinction curve, whereas Milky Way-type dust predicts too blue UV colors due to the 2175A bump. We expect that upcoming observations by ALMA will be able to test this heavily obscured model.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Halo gas in low-z (z<0.5) >0.1L* galaxies in high-resolution, large-scale cosmological hydrodynamic simulations is examined with respect to three components: (cold, warm, hot) with temperatures equal to (<10^5, 10^{5-6}, >10^6)K, respectively. The warm component is compared, utilizing O VI \lambda\lambda 1032, 1038 absorption lines, to observations and agreement is found with respect to the galaxy-O VI line correlation, the ratio of O VI line incidence rate in blue to red galaxies and the amount of O VI mass in star-forming galaxies. A detailed account of the sources of warm halo gas (stellar feedback heating, gravitational shock heating and accretion from the intergalactic medium), inflowing and outflowing warm halo gas metallicity disparities and their dependencies on galaxy types and environment is also presented. Having the warm component securely anchored, our simulations make the following additional predictions. First, cold gas is the primary component in inner regions, with its mass comprising 50% of all gas within galacto-centric radius r=(30,150)kpc in (red, blue) galaxies. Second, at r>(30,200)kpc in (red, blue) galaxies the hot component becomes the majority. Third, the warm component is a perpetual minority, with its contribution peaking at ~30% at r=100-300kpc in blue galaxies and never exceeding 5% in red galaxies. The significant amount of cold gas in low-z early-type galaxies found in simulations, in agreement with recent observations (Thom et al.), is intriguing, so is the dominance of hot gas at large radii in blue galaxies.

We present the first results from our pilot 500 ks Chandra-LETG Large Program observation of the soft X-ray brightest source in the z>=0.4 sky, the blazar 1ES 1553+113, aimed to secure the first uncontroversial detections of the missing baryons in the X-rays. We identify a total of 11 possible absorption lines, with single-line statistical significances between 2.2-4.1\sigma. Six of these lines are detected at high single-line statistical significance (3.6 <= sigma <= 4.1), while the remaining five are regarded as marginal detections in association with either other X-ray lines detected at higher significance and/or Far-Ultraviolet (FUV) signposts. In particular, five of these possible intervening absorption lines, are identified as CV and CVI K\alpha absorbers belonging to three WHIM systems at z_X = 0.312, z_X = 0.237 and <z_X> = 0.133, which also produce broad HI (and OVI for the z_X = 0.312 system) absorption in the FUV. For two of these systems (z_X = 0.312 and 0.237), the Chandra X-ray data led the a-posteriori discovery of physically consistent broad HI associations in the FUV, so confirming the power of the X-ray-FUV synergy for WHIM studies. The true statistical significances of these three X-ray absorption systems, after properly accounting for the number of redshift trials, are 5.8 sigma (z_X = 0.312; 6.3 sigma if the low-significance OV and CV K-beta associations are considered), 3.9 sigma (z_X = 0.237), and 3.8 sigma (\langle z_X \rangle = 0.133), respectively.

We examine the validity of the hydrostatic equilibrium (HSE) assumption for galaxy clusters using one of the highest-resolution cosmological hydrodynamical simulations. We define and evaluate several effective mass terms corresponding to the Euler equations of the gas dynamics, and quantify the degree of the validity of HSE in terms of the mass estimate. We find that the mass estimated under the HSE assumption (the HSE mass) deviates from the true mass by up to ~ 30 %. This level of departure from HSE is consistent with the previous claims, but our physical interpretation is rather different. We demonstrate that the inertial term in the Euler equations makes a negligible contribution to the total mass, and the overall gravity of the cluster is balanced by the thermal gas pressure gradient and the gas acceleration term. Indeed the deviation from the HSE mass is well explained by the acceleration term at almost all radii. We also clarify the confusion of previous work due to the inappropriate application of the Jeans equations in considering the validity of HSE from the gas dynamics extracted from cosmological hydrodynamical simulations.

We present a new method for modeling inhomogeneous cosmic reionization on large scales. Utilizing high-resolution radiation-hydrodynamic simulations with 2048^3 dark matter particles, 2048^3 gas cells, and 17 billion adaptive rays in a L = 100 Mpc/h box, we show that the density and reionization-redshift fields are highly correlated on large scales (>~ 1 Mpc/h). This correlation can be statistically represented by a scale-dependent linear bias. We construct a parametric function for the bias, which is then used to filter any large-scale density field to derive the corresponding spatially varying reionization-redshift field. The parametric model has three free parameters which can be reduced to one free parameter when we fit the two bias parameters to simulations results. We can differentiate degenerate combinations of the bias parameters by combining results for the global ionization histories and correlation length between ionized regions. Unlike previous semi-analytic models, the evolution of the reionization-redshift field in our model is directly compared cell by cell against simulations and preforms well in all tests. Our model maps the high-resolution, intermediate-volume radiation-hydrodynamic simulations onto lower-resolution, larger-volume N-body simulations (>~ 2 Gpc/h) in order to make mock observations and theoretical predictions.

We present new predictions for temperature (on small angular scales) and polarization (on large angular scales) CMB anisotropies induced during the epoch of reionization (EoR). Using a novel method calibrated from Radiation-Hydrodynamic simulations we model the EoR in large volumes (L >~ 2 Gpc/h) in the context of galactic reionization. We find that the EoR contribution to the kinetic Sunyaev- Zel'dovich power spectrum (patchy kSZ) ranges between ~0.6 - 2.8 muK^2 at ell = 3000, for the parameter space we explored. These patchy kSZ power spectra are calculated from large 15 Deg x 15 Deg maps that are found to be necessary. Decreasing the size of these maps biases the overall patchy kSZ power to higher values. We find that the amplitude of the patchy kSZ power spectrum at ell = 3000 follows simple scalings of D_ell=3000^kSZ propto <z> and D_ell=3000^kSZ propto Delz^0.47 for the mean redshift (<z>) of reionization and the duration (dz). Using the constraints on <z> from WMAP 7-year results and the lower limit on dz from EDGES we find a lower limit of ~ 0.4 muK^2 on the kSZ at ell = 3000. Planck will constrain the mean redshift and the Thomson optical depth from the low-ell polarization power spectrum. Future measurements of the high-ell CMB power spectrum from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) should detect the patchy kSZ signal if the cross correlation between the cosmic infrared background and the thermal Sunyaev Zel'dovich effect is constrained. We show that the combination of temperature and polarization measurements constrains both <z> and dz. The patchy kSZ maps, power spectra templates and the polarization power spectra will be publicly available.

We present a new model for the observed Lyman alpha blobs (LABs) within the context of the standard cold dark matter model. In this model, LABs are the most massive halos with the strongest clustering (proto-clusters) undergoing extreme starbursts in the high-z universe. Aided by calculations of detailed radiative transfer of Lya photons through ultra-high resolution (159pc) large-scale (>30Mpc) adaptive mesh-refinement cosmological hydrodynamic simulations with galaxy formation, this model is shown to be able to, for the first time, reproduce simultaneously the global Lya luminosity function and luminosity-size relation of the observed LABs. Physically, a combination of dust attenuation of Lya photons within galaxies, clustering of galaxies, and complex propagation of Lya photons through circumgalactic and intergalactic medium gives rise to the large sizes and frequently irregular isophotal shapes of LABs that are observed. A generic and unique prediction of this model is that there should be strong far-infrared (FIR) sources within each LAB, with the most luminous FIR source likely representing the gravitational center of the proto-cluster, not necessarily the apparent center of the Lya emission of the LAB or the most luminous optical source. Upcoming ALMA observations should unambiguously test this prediction. If verified, LABs will provide very valuable laboratories for studying formation of galaxies in the most overdense regions of the universe at a time when global star formation is most vigorous.

We adapt a modern scheme of smoothed particle hydrodynamics (SPH) to our tree N-body/SPH galactic chemodynamics code GCD+. The applied scheme includes imple- mentations of the artificial viscosity switch and artificial thermal conductivity pro- posed by Morris & Monaghan (1997), Rosswog & Price (2007) and Price (2008), to model discontinuities and Kelvin-Helmholtz instabilities more accurately. We first present hydrodynamics test simulations and contrast the results to runs undertaken without artificial viscosity switch or thermal conduction. In addition, we also explore the different levels of smoothing by adopting larger or smaller smoothing lengths, i.e. a larger or smaller number of neighbour particles, Nnb. We demonstrate that the new version of GCD+ is capable of modelling Kelvin-Helmholtz instabilities to a simi- lar level as the mesh code, Athena. From the Gresho vortex, point-like explosion and self-similar collapse tests, we conclude that setting the smoothing length to keep the number of neighbour particles as high as Nnb~58 is preferable to adopting smaller smoothing lengths. We present our optimised parameter sets from the hydrodynamics tests.

Radio observations of galaxy clusters show that there are $\mu$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($\sigma_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $\sigma_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

In a hierarchical universe, mergers may be an important mechanism not only in increasing the mass of galaxies but also in driving the color and morphological evolution of galaxies. We use a large sample of ~1000 simulated galaxies of stellar mass greater than 10^9.6 solar masses (for ~4800 observations at multiple redshifts) from a high-res (0.46 h^{-1} kpc) cosmological simulation to determine under what circumstances being a member of a pair influences galaxy properties at z <= 0.2. We identify gravitationally bound pairs, and find a relative fraction of blue-blue, red-red, and blue-red pairs that agrees with observations (Lin et al. 2010). Pairs tend to avoid the extreme environments of clusters and void centres. While pairs in groups can include galaxies that are both blue, both red, or one of each color, in the field it is rare for pair galaxies to both be red. We find that physically bound pairs closer than 250 h^{-1} kpc tend to have higher sSFRs than the galaxy population as a whole. However, the sSFR of a bound galaxy relative to galaxies in a comparable local density environment (determined by the distance to the fifth nearest neighbor, rho_5), differs depending on the local density. In regions of high rho_5 the bound population has a higher fraction of star-forming (bluer) galaxies, whereas there is very little difference between bound and unbound galaxies in low rho_5 regions. This effect on the star-forming fraction may be driven by the higher fraction of bound HI-rich galaxies compared to unbound galaxies, particularly at high local densities. It appears that being in a pair has an incremental, but not overwhelming, effect on the star formation rate of the paired galaxies, compared to the more pronounced trend where galaxies overall have low sSFR (are red) in clusters and higher sSFR (blue) at the centre of voids. This trend depends most strongly on rho_5.(abridged)

We examine galaxy formation in a cosmological AMR simulation, which includes two high resolution boxes, one centered on a 3 \times 10^14 M\odot cluster, and one centered on a void. We examine the evolution of 611 massive (M\ast > 10^10M\odot) galaxies. We find that the fraction of the final stellar mass which is accreted from other galaxies is between 15 and 40% and increases with stellar mass. The accreted fraction does not depend strongly on environment at a given stellar mass, but the galaxies in groups and cluster environments are older and underwent mergers earlier than galaxies in lower density environments. On average, the accreted stars are ~2.5 Gyrs older, and ~0.15 dex more metal poor than the stars formed in-situ. Accreted stellar material typically lies on the outskirts of galaxies; the average half-light radius of the accreted stars is 2.6 times larger than that of the in-situ stars. This leads to radial gradients in age and metallicity for massive galaxies, in qualitative agreement with observations. Massive galaxies grow by mergers at a rate of approximately 2.6% per Gyr. These mergers have a median (mass-weighted) mass ratio less than 0.26 \pm 0.21, with an absolute lower limit of 0.20, for galaxies with M\ast ~ 10^12 M\odot. This suggests that major mergers do not dominate in the accretion history of massive galaxies. All of these results agree qualitatively with results from SPH simulations by Oser et al. (2010, 2012).

A new model for coevolution of galaxies and supermassive black holes (SMBH) is presented that is physically based. The evolutionary track starts with an event that triggers a significant starburst in the central region of a galaxy. In this model, the main SMBH growth takes place in post-starburst phase fueled by recycled gas from inner bulge stars in a self-regulated fashion on a time scale that is substantially longer than 100Myrs and at a diminishing Eddington ratio with time. We argue that the SMBH cannot gorge itself during the starburst phase, despite the abundant supply of cold gas, because star formation is a preferred mode of gas consumption in such an environment than accretion to the central SMBH. We also show that feedback from star formation is at least as strong as that from AGN and thus, if star formation is in need of being quenched, AGN feedback generally does not play the primary role. The predicted relation between SMBH mass and bulge mass/velocity dispersion is consistent with observations. A clear prediction is that early-type galaxy hosts of high Eddingtion rate AGNs are expected to be light-blue to green in optical color, gradually evolving to the red sequences with decreasing AGN luminosity. A suite of falsifiable predictions and implications with respect to relationships between various types of galaxies and AGN, and others, are made. For those where comparisons to extant observations are possible, the model appears to be in good standing.

With high resolution (0.46kpc/h), adaptive mesh-refinement Eulerian cosmological hydrodynamic simulations we compute properties of O VI and O VII absorbers from the warm-hot intergalactic medium (WHIM). Our new simulations are in broad agreement with previous simulations, with ~40% of the intergalactic medium being in the WHIM at z=0. It is found (1) The amount of gas in the WHIM at temperature below and above 10^6K is about equal within uncertainties. (1) Our simulations are in excellent agreement with observed properties of O VI absorbers, with respect to the line incidence rate and Doppler width-column density relation. (2) Velocity structures within absorbing regions are a significant, and for large Doppler width clouds, a dominant contributor to the Doppler widths of both O VI and O VII absorbers. A non-negligible fraction (in number and mass) of O VI and O VII clouds can arise from gas of temperature lower than 10^5, until the Doppler width is well in excess of 100km/s. (3) Strong O VI absorbers are predominantly collisionally ionized. About (61%, 57%, 39%) of O VI absorbers in the column density ranges of log N(OVI) cm^2=(12.5-13,13-14,>14) have temperature lower than 10^5K. (4) Quantitative prediction is made for the presence of broad and shallow O VI lines, which current observations may have largely missed. Upcoming observations by COS may be able to provide a test. (5) The reported 3 sigma upper limit on the mean column density of coincidental O VII lines at the location of detected O VI lines by Yao et al is above the predicted value by a factor of 2.5-4. (6) The claimed observational detection of O VII lines by Nicastro et al, if true, is 2 sigma above what our simulations predict.

We present a new method for incorporating nonthermal pressure from bulk motions of gas into an analytic model of the intracluster medium in clusters of galaxies, which is based on a polytropic equation of state and hydrostatic equilibrium inside gravitational potential wells drawn from cosmological dark matter simulations. The pressure is allowed to have thermal and nonthermal components with different radial distributions; the overall level of nonthermal support is based on the dynamical state of the halo, such that it is lower in more relaxed clusters. This level is normalized by comparison to pressure profiles derived from X-ray observations, and to a high resolution hydrodynamical simulation. The nonthermal pressure fraction measured at r_500 is typically in the range 10-20%, increasing with cluster mass and with redshift. The resulting model cluster properties are in accord with Sunyaev-Zel'dovich (SZ) effect observations of clusters. Inclusion of nonthermal pressure reduces the expected angular power spectrum of SZ fluctuations in the microwave sky by 24%.

Utilizing a high-resolution (114 pc/h) adaptive mesh-refinement cosmological galaxy formation simulation of the standard cold dark matter model with a large (2000-3000 galaxies with stellar mass greater than 1e9 Msun) statistical sample, we examine the role of major mergers in driving star formation at z>1 in a cosmological setting, after validating that some of the key properties of simulated galaxies are in reasonable agreement with observations, including luminosity functions, SF history, effective sizes and damped Lyman alpha systems. We find that major mergers have a relatively modest effect on star formation, in marked contrast to previous idealized merger simulations of disk galaxies that show up to two orders of magnitude increase in star formation rate. At z=2.4-3.7, major mergers tend to increase the specific star formation rate by 10-25% for galaxies in the entire stellar mass range 10^9-10^12 Msun probed. Their effect appears to increase with decreasing redshift, but is capped at 60% at z=1.4-2.4. Two factors may account for this modest effect. First, SFR of galaxies not in major mergers are much higher at z>1 than local disk galaxy counterparts. Second, most galaxies at z>1 have small sizes and contain massive dense bulges, which suppress the merger induced structural effects and gas inflow enhancement. Various other predictions are also made that will provide verifiable tests of the model.

Utilizing state-of-the-art, adaptive mesh-refinement cosmological hydrodynamic simulations with ultra-high resolution (114h-1pc) and large sample size (>3300 galaxies of stellar mass >10^9Msun), we show how the stellar light of Lyman Break Galaxies at z=2 is distributed between optical/ultra-violet (UV) and far-infrared (FIR) bands. With a single scalar parameter for dust obscuration we can simultaneously reproduce the observed UV luminosity function for the entire range (3-100 Msun/yr) and extant FIR luminosity function at the bright end (>20Msun/yr). We quantify that galaxies more massive or having higher SFR tend to have larger amounts of dust obscuration mostly due to a trend in column density and in a minor part due to a mass (or SFR)-metallicity relation. It is predicted that the FIR luminosity function in the range SFR=1-100Msun/yr is a powerlaw with a slope about -1.7. We further predict that there is a "galaxy desert" at SFR(FIR) < 0.02 (SFR(UV)/10Msun/yr)^2.1 Msun/yr in the SFR(UV)-SFR(FIR) plane. Detailed distributions of SFR(FIR) at a fixed SFR(UV) are presented. Upcoming observations by ALMA should test this model. If confirmed, it validates the predictions of the standard cold dark matter model and has important implications on the intrinsic SFR function of galaxies at high redshift.

Lyman-alpha (Lya) photons that escape the interstellar medium of star-forming galaxies may be resonantly scattered by neutral hydrogen atoms in the circumgalactic and intergalactic media, thereby increasing the angular extent of the galaxy's Lya emission. We present predictions of this extended, low surface brightness Lya emission based on radiative transfer modeling in a cosmological reionization simulation. The extended emission can be detected from stacked narrowband images of Lya emitters (LAEs) or of Lyman break galaxies (LBGs). Its average surface brightness profile has a central cusp, then flattens to an approximate plateau beginning at an inner characteristic scale below ~0.2 Mpc (comoving), then steepens again beyond an outer characteristic scale of ~1 Mpc. The inner scale marks the transition from scattered light of the central source to emission from clustered sources, while the outer scale marks the spatial extent of scattered emission from these clustered sources. Both scales tend to increase with halo mass, UV luminosity, and observed Lya luminosity. The extended emission predicted by our simulation is already within reach of deep narrowband photometry using large ground-based telescopes. Such observations would test radiative transfer models of emission from LAEs and LBGs, and they would open a new window on the circumgalactic environment of high-redshift star-forming galaxies.

Utilizing high-resolution large-scale galaxy formation simulations of the standard cold dark matter model, we examine global trends in the evolution of galaxies due to gravitational shock heating by collapse of large halos and large-scale structure. We find two major global trends. (1) The mean specific star formation rate (sSFR) at a given galaxy mass is a monotonically increasing function with increasing redshift. (2) The mean sSFR at a given redshift is a monotonically increasing function of decreasing galaxy mass that steepens with decreasing redshift. The general dimming trend with time merely reflects the general decline of gas inflow rate with increasing time. The differential evolution of galaxies of different masses with redshift is a result of gravitational shock heating of gas due to formation of large halos (groups and clusters) and large-scale structure that move a progressively larger fraction of galaxies and their satellites into environments where gas has too high an entropy to cool to continue feeding resident galaxies. Overdense regions where larger halos are preferentially located begin to be heated earlier and have higher temperatures than lower density regions at any given time, causing sSFR of larger galaxies to fall below the general dimming trend at higher redshift than less massive galaxies and galaxies with high sSFR to gradually shift to lower density environments at lower redshift. We find that several noted cosmic downsizing phenomena are different manifestations of these general trends. We also find that the great migration of galaxies from blue cloud to red sequence as well as color-density relation, among others, may arise naturally in this picture.

Using hydrodynamic simulations we compute the metal enrichment history of the intergalactic medium (IGM). We show that galactic superwind (GSW) feedback can transport metals to the IGM and that the properties of simulated metal absorbers match observations. The distance of influence of GSW is typically limited to >0.5Mpc and within regions of overdensity >10. Most CIV and OVI absorbers are located within shocked regions of elevated temperature (T>2x10^4K), overdensity (>10), and metallicity ([-2.5,-0.5]). OVI absorbers have typically higher metallicity, lower density and higher temperature than CIV absorbers. For OVI absorbers collisional ionization dominates over the entire redshift range z=0-6, whereas for CIV absorbers the transition occurs at moderate redshift z~3 from collisionally dominated to photoionization dominated. We find that the observed column density distributions for CIV and OVI in the range log N cm^2=12-15 are reasonably reproduced by the simulations. The evolution of mass densities contained in CIV and OVI lines, Omega_CIV and Omega_OVI, is also in good agreement with observations, which shows a near constancy at low redshifts and an exponential drop beyond redshift z=3-4. For both CIV and OVI, most absorbers are transient and the amount of metals probed by CIV and OVI lines of column log N cm^2=12-15 is only ~2% of total metal density at any epoch. While gravitational shocks from large-scale structure formation dominate the energy budget (80-90%) for turning about 50% of IGM to the warm-hot intergalactic medium (WHIM) by z=0, GSW feedback shocks are energetically dominant over gravitational shocks at z > 1-2. Most of the so-called "missing metals" at z=2-3 are hidden in a warm-hot (T=10^{4.5-7}K) gaseous phase, heated up by GSW feedback shocks. Their mass distribution is broadly peaked at $\delta=1-10$ in the IGM, outside virialized halos.

We show that a quasi-spherical (QS) hot accretion flow is expected to operate in all SMBHs, with its rate being capped at dot m_QS,max~0.001 (M/1E8Msun) in units of the Eddington rate. It is then proposed that AGN jet power is proportional to the product of the hot accretion rate and a SMBH spin-dependent efficiency for energy extraction in the form of jets. Predictions from this model include (1) while radio jets should emerge from all SMBHs, the maximum jet power goes with SMBH mass approximately as 1E43.6(M/1E8Msun)^{2}erg/s, (2) even although there are two separate underlying populations of radio-quiet (RQ) and radio-loud (RL) AGNs, any bimodality in the observed AGN radio loudness distribution is likely due to a selection effect, such as some imposed optical magnitude limits, (3) the RL fraction of quasars is expected to decrease with redshift in the cold dark matter model, (4) host galaxies of RQ and RL quasars should be drawn from the same underlying elliptical galaxy population, although RL quasars may have SMBHs that are somewhat more massive than their RQ counterparts and RL quasars may reside predominantly in core elliptical galaxies, (5) RL low-luminosity AGNs and LINERs may represent the long, declining "trailing" phase following the initial, more luminous AGN phase, (6) a broad anti-correlation between radio- loudness and disk accretion rate is expected, (7) RL AGNs may be expected to be more abundant in Type Ia supernovae than RQ AGNs, (8) among RL AGNs a correlation between radio power and clustering strength is predicted, and (9) RQ and RL AGNs should have, on average, similar IR-optical-UV properties. ~

We analyze a (120 h^{-1} Mpc)^3 adaptive mesh refinement hydrodynamic simulation that contains a higher-resolution 31 x 31 x 35 h^{-3} Mpc subvolume centered on a ~30 Mpc diameter void. Our detailed ~1 kpc resolution allows us to identify 1300 galaxies within this void to a limiting halo mass of ~10^{10} M_sun. Nearly 1000 galaxies are found to be in underdense regions, with 300 galaxies residing in regions less than half the mean density of the simulation volume. We construct mock observations of the stellar and gas properties of these systems, and reproduce the range of colors and luminosities observed in the SDSS for nearby (z < 0.03) galaxies. We find no trends with density for the most luminous (M_r < -18) galaxies, however our dwarf void galaxies (M_r > -16), though they are less reliably resolved, typically appear bluer, with higher rates of star formation and specific star formation and lower mean stellar ages than galaxies in average density environments. We find a significant population of low luminosity (M_r ~ -14) dwarf galaxies that is preferentially located in low density regions and specifically in the void center. This population may help to reduce, but not remove, the discrepancy between the predicted and observed number of void galaxies.

Observations of anisotropies in the brightness temperature of the 21 cm line of neutral hydrogen from the period before reionization would shed light on the dawn of the first stars and galaxies. In this paper, we use large-scale semi-numerical simulations to analyse the imprint on the 21 cm signal of spatial fluctuations in the Lyman-alpha flux arising from the clustering of the first galaxies. We show that an experiment such as the Square Kilometer Array (SKA) can probe this signal at the onset of reionization, giving us important information about the UV emission spectra of the first stars and characterizing their host galaxies. SKA-pathfinders with ~ 10% of the full collecting area should be capable of making a statistical detection of the 21 cm power spectrum at redshifts z < 20 (corresponding to frequencies $\nu$ > 67 MHz). We then show that the SKA should be able to measure the three dimensional power spectrum as a function of the angle with the line of sight and discuss the use of the redshift space distortions as a way to separate out the different components of the 21 cm power spectrum. We demonstrate that, at least on large scales where the Lyman-alpha fluctuations are linear, they can be used as a model independent way to extract the power spectra due to these Lyman-alpha fluctuations.

We study the clustering properties of z~5.7 Lyman-alpha emitters (LAEs) in a cosmological reionization simulation with a full Lya radiative transfer calculation. Lya radiative transfer substantially modifies the intrinsic Lya emission properties, compared to observed ones, depending on the density and velocity structure environment around the Lya emitting galaxy. This environment-dependent Lya selection introduces new features in LAE clustering, suppressing (enhancing) the line-of-sight (transverse) density fluctuations and giving rise to scale-dependent galaxy bias. In real space, the contours of the three-dimensional two-point correlation function of LAEs appear to be prominently elongated along the line of sight on large scales, an effect that is opposite to and much stronger than the linear redshift-space distortion effect. The projected two-point correlation function is greatly enhanced in amplitude by a factor of up to a few, compared to the case without the environment dependent selection effect. The new features in LAE clustering can be understood with a simple, physically motivated model, where Lya selection depends on matter density, velocity, and their gradients. We discuss the implications and consequences of the effects on galaxy clustering from Lya selection in interpreting clustering measurements and in constraining cosmology and reionization from LAEs.

Using adaptive mesh-refinement cosmological hydrodynamic simulations with a physically motivated supernova feedback prescription we show that the standard cold dark matter model can account for extant observed properties of damped Lyman alpha systems (DLAs). We then examine the properties of DLA host galaxies. We find: (1) While DLA hosts roughly trace the overall population of galaxies at all redshifts, they are always gas rich. (2) The history of DLA evolution reflects primarily the evolution of the underlying cosmic density, galaxy size and galaxy interactions. With higher density and more interactions at high redshift DLAs are larger in both absolute terms and in relative terms with respect to virial radii of halos. (3) The variety of DLAs at high redshift is richer with a large contribution coming from galactic filaments, created through close galaxy interactions. The portion of gaseous disks of galaxies where most stars reside makes relatively small contribution to DLA incidence at z=3-4. (4) The vast majority of DLAs arise in halos of mass M_h=10^10-10^12 Msun at z=1.6-4. At z=3-4, 20-30% of DLA hosts are Lyman Break Galaxies (LBGs). (5) Galactic winds play an indispensable role in shaping the kinematic properties of DLAs. Specifically, the high velocity width DLAs are a mixture of those arising in high mass, high velocity dispersion halos and those arising in smaller mass systems where cold gas clouds are entrained to high velocities by galactic winds. (6) In agreement with observations, we see a weak but noticeable evolution in DLA metallicity. The metallicity distribution centers at [Z/H]=-1.5 to -1 at z=3-4, with the peak moving to [Z/H]=-0.75 at z=1.6 and [Z/H]=-0.5 by z=0. (7) The star formation rate of DLA hosts is concentrated in the range 0.3-30Msun/yr at z=3-4, gradually shifting lower to peak at ~0.5-1 Msun/yr by z=0.

We present the Smoothed Hessian Major Axis Filament Finder (SHMAFF), an algorithm that uses the eigenvectors of the Hessian matrix of the smoothed galaxy distribution to identify individual filamentary structures. Filaments are traced along the Hessian eigenvector corresponding to the largest eigenvalue, and are stopped when the axis orientation changes more rapidly than a preset threshold. In both N-body simulations and the Sloan Digital Sky Survey (SDSS) main galaxy redshift survey data, the resulting filament length distributions are approximately exponential. In the SDSS galaxy distribution, using smoothing lengths of 10 h^{-1} Mpc and 15 h^{-1} Mpc, we find filament lengths per unit volume of 1.9x10^{-3} h^2 Mpc^{-2} and 7.6x10^{-4} h^2 Mpc^{-2}, respectively. The filament width distributions, which are much more sensitive to non-linear growth, are also consistent between the real and mock galaxy distributions using a standard cosmology. In SDSS, we find mean filament widths of 5.5 h^{-1} Mpc and 8.4 h^{-1} Mpc on 10 h^{-1} Mpc and 15 h^{-1} Mpc smoothing scales, with standard deviations of 1.1 h^{-1} Mpc and 1.4 h^{-1} Mpc, respectively. Finally, the spatial distribution of filamentary structure in simulations is very similar between z=3 and z=0 on smoothing scales as large as 15 h^{-1} Mpc, suggesting that the outline of filamentary structure is already in place at high redshift.

In the context stellar reionization in the standard cold dark matter model, we analyze observations at z~6 and are able to draw three significant conclusions with respect to star formation and the state of the intergalactic medium (IGM) at z~6. (1) An initial stellar mass function (IMF) more efficient, by a factor of 10-20, in producing ionizing photons than the standard Salpeter IMF is required at z~6. This may be achieved by having either (A) a metal-enriched IMF with and a lower mass cutoff of >= 30Msun or (B) 2-4% of stellar mass being Population III massive metal-free stars at z~6. While there is no compelling physical reason or observational evidence to support (A), (B) could be fulfilled plausibly by continued existence of some pockets of uncontaminated, metal-free gas for star formation. (2) The volume-weighted neutral fraction of the IGM of <f_HI>_V~ 10^-4 at z=5.8 inferred from the SDSS observations of QSO absorption spectra provides enough information to ascertain that reionization is basically complete with at most ~0.1-1% of IGM that is un-ionized at z=5.8. (3) Barring some extreme evolution of the IMF, the neutral fraction of the IGM is expected to rise quickly toward high redshift from the point of HII bubble percolation, with the mean neutral fraction of the IGM expected to reach 6-12% at z=6.5, 13-27% at z=7.7 and 22-38% at z=8.8.

Although coherent large-scale structures such as filaments and walls are apparent to the eye in galaxy redshift surveys, they have so far proven difficult to characterize with computer algorithms. This paper presents a procedure that uses the eigenvalues and eigenvectors of the Hessian matrix of the galaxy density field to characterize the morphology of large-scale structure. By analysing the smoothed density field and its Hessian matrix, we can determine the types of structure - walls, filaments, or clumps - that dominate the large-scale distribution of galaxies as a function of scale. We have run the algorithm on mock galaxy distributions in a LCDM cosmological N-body simulation and the observed galaxy distributions in the Sloan Digital Sky Survey. The morphology of structure is similar between the two catalogues, both being filament-dominated on 10-20 h^{-1} Mpc smoothing scales and clump-dominated on 5 h^{-1} Mpc scales. There is evidence for walls in both distributions, but walls are not the dominant structures on scales smaller than ~25 h^{-1} Mpc. Analysis of the simulation suggests that, on a given comoving smoothing scale, structures evolve with time from walls to filaments to clumps, where those found on smaller smoothing scales are further in this progression at a given time.

We combine a cosmological reionization simulation with box size of 100Mpc/h on a side and a Monte Carlo Lyman-alpha (Lya) radiative transfer code to model Lyman Alpha Emitters (LAEs) at z~5.7. The model introduces Lya radiative transfer as the single factor for transforming the intrinsic Lya emission properties into the observed ones. Spatial diffusion of Lya photons from radiative transfer results in extended Lya emission and only the central part with high surface brightness can be observed. Because of radiative transfer, the appearance of LAEs depends on density and velocity structures in circumgalactic and intergalactic media as well as the viewing angle, which leads to a broad distribution of apparent (observed) Lya luminosity for a given intrinsic Lya luminosity. Radiative transfer also causes frequency diffusion of Lya photons. The resultant Lya line is asymmetric with a red tail. The peak of the Lya line shifts towards longer wavelength and the shift is anti-correlated with the apparent to intrinsic Lya luminosity ratio. The simple radiative transfer model provides a new framework for studying LAEs. It is able to explain an array of observed properties of z~5.7 LAEs in Ouchi et al. (2008), producing Lya spectra, morphology, and apparent Lya luminosity function (LF) similar to those seen in observation. The broad distribution of apparent Lya luminosity at fixed UV luminosity provides a natural explanation for the observed UV LF, especially the turnover towards the low luminosity end. The model also reproduces the observed distribution of Lya equivalent width (EW) and explains the deficit of UV bright, high EW sources. Because of the broad distribution of the apparent to intrinsic Lya luminosity ratio, the model predicts effective duty cycles and Lya escape fractions for LAEs.

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies.

We compare the predictions of four different algorithms for the distribution of ionized gas during the Epoch of Reionization. These algorithms are all used to run a 100 Mpc/h simulation of reionization with the same initial conditions. Two of the algorithms are state-of-the-art ray-tracing radiative transfer codes that use disparate methods to calculate the ionization history. The other two algorithms are fast but more approximate schemes based on iterative application of a smoothing filter to the underlying source and density fields. We compare these algorithms' resulting ionization and 21 cm fields using several different statistical measures. The two radiative transfer schemes are in excellent agreement with each other (with the cross-correlation coefficient of the ionization fields >0.8 for k < 10 h/Mpc and in good agreement with the analytic schemes (>0.6 for k < 1 h/Mpc). When used to predict the 21cm power spectrum at different times during reionization, all ionization algorithms agree with one another at the 10s of percent level. This agreement suggests that the different approximations involved in the ray tracing algorithms are sensible and that semi-numerical schemes provide a numerically-inexpensive, yet fairly accurate, description of the reionization process.

Using high redshift radio sources as background, the 21cm forest observations probe the neutral hydrogen absorption signatures of early structures along the lines of sight. Directly sensitive to the spin temperature, it complements the 21cm tomography observations, and provides information on the temperature as well as the ionization state of the intergalactic medium (IGM). We use a radiative transfer simulation to investigate the 21cm forest signals during the epoch of reionization. We first confirmed that the optical depth and equivalent width (EW) are good representations of the ionization and thermal state of the IGM. The features selected by their relative optical depth are excellent tracers of ionization fields, while those selected by absolute optical depth are very sensitive to the IGM temperature, so the IGM temperature information could potentially be extracted, thus breaking a degeneracy in 21cm tomographic observation. From the distributions of EWs and the number evolution of absorbers and leakers with different EWs, we see clearly the cosmological evolution of ionization state of the IGM. The number density of potentially observable features decreases rapidly with increasing gas temperature. The sensitivity of the proposed EW statistic to the IGM temperature makes it a unique and potentially powerful probe of reionization. Missing small-scale structures, such as small filaments and minihalos that are unresolved in our current simulation, and lack of an accurate calculation of the IGM temperature, however, likely have rendered the presented signals quantitatively inaccurate. Finally we discuss the requirements of the background radio sources for such observations, and find that signals with equivalent widths larger than 1kHz are hopeful to be detected.(Abridged)

The inhomogeneous cosmological reionization process leaves tangible imprints in the intergalactic medium down to z=4-5. The Lyman-alpha forest flux power spectrum provides a potentially powerful probe of the epoch of reionization. With the existing SDSS I/II quasar sample we show that two cosmological reionization scenarios, one completing reionization at z=6 and the other at z=9, can be distinguished at ~7 sigma level by utilizing Lyman alpha forest absorption spectra at z=4.5+-0.5, in the absence of other physical processes that may also affect the Lyman alpha flux power spectrum. The redshift range z=4-5 may provide the best window, because there is still enough transmitted flux and quasars to measure precise statistics of the flux fluctuations, and the IGM still retains a significant amount of memory of reionization.

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.

One of the next decade's most exciting prospects is to explore the cosmic "dark ages," during which the first stars in the Universe formed, with the 21 cm line of neutral hydrogen. At z>6, this light redshifts into the low-frequency radio regime and can be observed with new generations of wide-field arrays. These experiments have the potential to observe enormous volumes and to improve both cosmological and astrophysical constraints dramatically. In this white paper, we describe how the next decade will see the first steps toward unlocking this potential, allowing us to answer two key questions about the fundamental physics of our Universe and the intergalactic medium: Does the standard cosmological model describe the Universe during the "dark ages?" How did the intergalactic medium evolve during this important time, ending with the reionization of hydrogen?

The cosmic dark ages and the epoch of reionization, during which the first generations of stars and galaxies formed, are among the most compelling frontiers of extragalactic astrophysics and cosmology. Here we describe an exciting new probe of these eras: the 21 cm line of neutral hydrogen, which will allow us to map the neutral intergalactic medium throughout the era of first galaxy formation. In this white paper, we describe how these studies can answer two key questions about galaxy formation: What were the properties of high-redshift galaxies? How did they affect the Universe around them?

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

It has been argued that low-luminosity dwarf galaxies are the dominant source of ionizing radiation during cosmological reionization. The fraction of ionizing radiation that escapes into the intergalactic medium from dwarf galaxies with masses less than ~10^9.5 solar masses plays a critical role during this epoch. Using an extensive suite of very high resolution (0.1 pc), adaptive mesh refinement, radiation hydrodynamical simulations of idealized and cosmological dwarf galaxies, we characterize the behavior of the escape fraction in galaxies between 3 x 10^6 and 3 x 10^9 solar masses with different spin parameters, amounts of turbulence, and baryon mass fractions. For a given halo mass, escape fractions can vary up to a factor of two, depending on the initial setup of the idealized halo. In a cosmological setting, we find that the time-averaged photon escape fraction always exceeds 25% and reaches up to 80% in halos with masses above 10^8 solar masses with a top-heavy IMF. The instantaneous escape fraction can vary up to an order of magnitude in a few million years and tend to be positively correlated with star formation rate. We find that the mean of the star formation efficiency times ionizing photon escape fraction, averaged over all atomic cooling (T_vir > 8000 K) galaxies, ranges from 0.02 for a normal IMF to 0.03 for a top-heavy IMF, whereas smaller, molecular cooling galaxies in minihalos do not make a significant contribution to reionizing the universe due to a much lower star formation efficiency. These results provide the physical basis for cosmological reionization by stellar sources, predominately atomic cooling dwarf galaxies.

We study the impact of inhomogeneous hydrogen reionization on the thermal evolution of the intergalactic medium (IGM) using hydrodynamic + radiative transfer simulations where reionization is completed either early (z ~ 9) or late (z ~ 6). In general, we find that low-density gas near large-scale overdensities is ionized and heated earlier than gas in the large-scale, underdense voids. Furthermore, at a later time the IGM temperature is inversely related to the reionization redshift because gas that is heated earlier has more time to cool through adiabatic expansion and Compton scattering. Thus, at the end of reionization the median temperature-density relation is an inverted power-law with slope gamma-1 ~ -0.2, in both models. However, at fixed density, there is up to order unity scatter in the temperature due to the distribution of reionization redshifts. Because of the complex equation-of-state, the evolved IGM temperature-density relations for the redshift range 4 < z < 6 can still have significant curvature and scatter. These features must be taken into account when interpreting the Lyman alpha absorption in high redshift quasar spectra.

Using state-of-the-art adaptive mesh refinement cosmological hydrodynamic simulations with a spatial resolution of proper 0.21kpc/h in refined subregions embedded within a comoving cosmological volume (27.4Mpc/h)^3, we investigate the sizes of galaxies at z=3 in the standard cold dark matter model where reionization is assumed to complete at zri~6. Our simulated galaxies are found to be significantly smaller than the observed ones: while more than one half of the galaxies observed by HST and VLT ranging from rest-frame UV to optical bands with stellar masses larger than 2E10 Msun have half-light radii larger than 2kpc/h, none of the simulated massive galaxies in the same mass range have half-light radii larger than 2kpc/h, after taking into account dust extinction. Corroborative evidence is provided by the rotation curves of the simulated galaxies with total masses of 1E11-1E12Msun, which display values 300-1000km/s at small radii (0.5kpc/h) due to high stellar concentration in the central regions, larger than those of any well observed galaxies. Possible physical mechanisms to resolve this serious problem include: (1) an early reionization at zri>>6 to suppress gas condensation hence star formation, (2) a strong, internal energetic feedback from stars or central black holes to reduce the overall star formation efficiency, or (3) a substantial small-scale cutoff in the matter power spectrum.

We measure several properties of the reionization process and the corresponding low-frequency 21-cm signal associated with the neutral hydrogen distribution, using a large volume, high resolution simulation of cosmic reionization. The brightness temperature of the 21-cm signal is derived by post-processing this numerical simulation with a semi-analytical prescription. Our study extends to high redshifts (z ~ 25) where, in addition to collisional coupling, our post-processed simulations take into account the inhomogeneities in the heating of the neutral gas by X-rays and the effect of an inhomogeneous Lya radiation field. Unlike the well-studied case where spin temperature is assumed to be significantly greater than the temperature of the cosmic microwave background due to uniform heating of the gas by X-rays, spatial fluctuations in both the Lya radiation field and X-ray intensity impact predictions related to the brightness temperature at z > 10, during the early stages of reionization and gas heating. The statistics of the 21-cm signal from our simulation are then compared to existing analytical models in the literature and we find that these analytical models provide a reasonably accurate description of the 21-cm power spectrum at z < 10. Such an agreement is useful since analytical models are better suited to quickly explore the full astrophysical and cosmological parameter space relevant for future 21-cm surveys. We find, nevertheless, non-negligible differences that can be attributed to differences in the inhomogeneous X-ray heating and Lya coupling at z > 10 and, with upcoming interferometric data, these differences in return can provide a way to better understand the astrophysical processes during reionization.

We investigate the relationship between galaxies and metal-line absorption systems in a large-scale cosmological simulation with galaxy formation. Our detailed treatment of metal enrichment and non-equilibrium calculation of oxygen species allow us, for the first time, to carry out quantitative calculations of the cross-correlations between galaxies and O VI absorbers. We find the following: (1) The cross-correlation strength depends weakly on the absorption strength but strongly on the luminosity of the galaxy. (2) The correlation distance increases monotonically with luminosity from ~0.5-1h^-1 Mpc for 0.1L* galaxies to ~3-5h^-1 Mpc for L* galaxies. (3) The correlation distance has a complicated dependence on absorber strength, with a luminosity-dependent peak. (4) Only 15% of O VI absorbers lie near >=Lz* galaxies. The remaining 85%, then, must arise ``near'' lower-luminosity galaxies, though, the positions of those galaxies is not well-correlated with the absorbers. This may point to pollution of intergalactic gas predominantly by smaller galaxies. (5) There is a subtle trend that for >~0.5Lz* galaxies, there is a positive correlation between absorber strength and galaxy luminosity in the sense that stronger absorbers have a slightly higher probability of finding such a large galaxy at a given projection distance. For less luminous galaxies, there seems to be a negative correlation between luminosity and absorber strength.

We present general properties of ionized hydrogen (HII) bubbles and their growth based on a state-of-the-art large-scale (100 Mpc/h) cosmological radiative transfer simulation. The simulation resolves all halos with atomic cooling at the relevant redshifts and simultaneously performs radiative transfer and dynamical evolution of structure formation. Our major conclusions include: (1) for significant HII bubbles, the number distribution is peaked at a volume of $\sim 0.6 {\rm Mpc^{3}/h^{3}}$ at all redshifts. But, at $z\le 10$, one large, connected network of bubbles dominates the entire HII volume. (2) HII bubbles are highly non-spherical. (3) The HII regions are highly biased with respect to the underlying matter distribution with the bias decreasing with time. (4) The non-gaussianity of the HII region is small when the universe becomes 50% ionized. The non-gaussianity reaches its maximal near the end of the reionization epoch $z\sim 6$. But at all redshifts of interest there is a significant non-gaussianity in the HII field. (5) Population III galaxies may play a significant role in the reionization process. Small bubbles are initially largely produced by Pop III stars. At $z\ge 10$ even the largest HII bubbles have a balanced ionizing photon contribution from Pop II and Pop III stars, while at $z\le 8$ Pop II stars start to dominate the overall ionizing photon production for large bubbles, although Pop III stars continue to make a non-negligible contribution. (6) The relationship between halo number density and bubble size is complicated but a strong correlation is found between halo number density and bubble size for large bubbles.

Using the largest cosmological reionization simulation to-date (~24 billion particles), we use the genus curve to quantify the topology of neutral hydrogen distribution on scales > 1 Mpc as it evolves during cosmological reionization. We find that the reionization process proceeds primarily in an inside-out fashion, where higher density regions become ionized earlier than lower density regions. There are four distinct topological phases: (1) Pre-reionization at z ~ 15, when the genus curve is consistent with a Gaussian density distribution. (2) Pre-overlap at 10 < z < 15, during which the number of HII bubbles increases gradually with time, until percolation of HII bubbles starts to take effect, characterized by a very flat genus curve at high volume fractions. (3) Overlap at 8 < z < 10, when large HII bubbles rapidly merge, manifested by a precipitous drop in the amplitude of the genus curve. (4) Post-overlap at 6 < z < 8, when HII bubbles have mostly overlapped and the genus curve is consistent with a diminishing number of isolated neutral islands. After the end of reionization (z < 6), the genus of neutral hydrogen is consistent with Gaussian random phase, in agreement with observations.

Based on a simulation of galaxy formation in the standard cosmological model, we suggest that a consistent picture for Gamma-Ray Bursts and star formation may be found that is in broad agreement with observations: GRBs preferentially form in low metallicity environments and in galaxies substantially less luminous that L*. We find that the computed formation rate of stars with metallicity less than 0.1Zsun agrees remarkably well with the rate evolution of Gamma-Ray Bursts observed by Swift from z=0 to z=4, whereas the evolution of total star formation rate is weaker by a factor of about 4. Given this finding, we caution that any inference of star formation rate based on observed GRB rate may require a more involved exercise than a simple proportionality.

We present a new hybrid code for large volume, high resolution simulations of cosmic reionization, which utilizes a N-body algorithm for dark matter, physically motivated prescriptions for baryons and star formation, and an adaptive ray tracing algorithm for radiative transfer of ionizing photons. Two test simulations each with 3 billion particles and 400 million rays in a 50 Mpc/h box have been run to give initial results. Halos are resolved down to virial temperatures of 10^4 K for the redshift range of interest in order to robustly model star formation and clumping factors. This is essential to correctly account for ionization and recombination processes. We find that the halos and sources are strongly biased with respect to the underlying dark matter, re-enforcing the requirement of large simulation boxes to minimize cosmic variance and to obtain a qualitatively correct picture of reionization. We model the stellar initial mass function (IMF), by following the spatially dependent gas metallicity evolution, and distinguish between the first generation, Population III (PopIII) stars and the second generation, Population II (PopII) stars. The PopIII stars with a top-heavy IMF produce an order of magnitude more ionizing photons at high redshifts z>10, resulting in a more extended reionization. In our simulations, complete overlap of HII regions occurrs at z~6.5 and the computed mass and volume weighted residual HI fractions at 5<z<6.5 are both in good agreement with high redshift quasar absorption measurements from SDSS. The values for the Thomson optical depth are consistent within 1-sigma of the current best-fit value from third-year WMAP.

A significant fraction of nearby elliptical galaxies are known to have high density gas disks in their circumnuclear (CN) region (0.1 to a few kpc). Yet, ellipticals, especially luminous ones, show little signs of recent star formation (SF). To investigate the possible cause of the dearth of SF in these systems, we study the gravitational stability of CN gas disks embedded within the potentials of both the stellar bulge and the central massive black hole (BH) in ellipticals. We find that CN disks in higher mass galaxies are generally more stable than those in lower mass galaxies, because higher mass galaxies tend to have more massive BHs and more centrally concentrated stellar density profiles. We also consider the case in which the central stellar density profile has a core, which is often observed for ellipticals whose total stellar mass is higher than about 10^11 Msun. Such a cored stellar density profile leads to more unstable CN disks than the power-law density profile characteristic of less massive galaxies. However, the more massive BHs in high-mass galaxies act to stabilize the CN disk. Our results demonstrate that the gravitational potentials of both the central BH and the stellar component should be taken into account when studying the properties of CN disks, as their stability is sensitive to both the BH mass and the stellar density profile. Our results could explain the observed trend that less luminous ellipticals have a greater tendency to exhibit ongoing SF than giant ellipticals.

Extended 21cm absorption regions (dubbed ``21cm absorption halos'') around first galaxies at $z\sim 30$ are likely the first distinctive structures accessible to radio observations. Though the radio array capable of detecting and resolving them must have $\sim 200$ km$^2$ total collecting area, given the great impact of such detections to the understanding of the reionization process and cosmology, such radio survey would be extremely profitable. As an example, we point out a potentially useful byproduct of such survey. The resolved 21cm absorption ``halos'', likely close to spherical, can serve as (almost) ideal sources for measuring the {\it cosmic shear} and mapping the matter distribution to $z\sim 30$. We investigate the expected lensing signal and consider a variety of noise contributions on the shear measurement. We find that S/N $\sim 1$ can be achieved for individual ``halos''. Given millions of 21cm absorption ``halos'' across the sky, the total S/N will be comparable to traditional shear measurement of $\sim$$10^9$ galaxies at $z\sim 1$.

We differentiate between the metal enrichment of the gas in virialized minihalos and that of the intergalactic medium at high redshift, pertinent to cosmological reionization, with the initial expectation that gas in the high density regions within formed dark matter halos may be more robust thus resistant to mixing with lower density intergalactic medium. Using detailed hydrodynamic simulations of gas clouds in minihalos subject to destructive processes associated with the encompassing intergalactic shocks carrying metal-enriched gas, we find, as an example, that, for realistic shocks of velocities of 10-100km/s, more than (90%,65%) of the high density gas with rho>500 rhob inside a minihalo virialized at z=10 of mass (10^7,10^6)Msun remains at a metallicity lower than 3% of that of the intergalactic medium by redshift z=6. It may be expected that the high density gas in minihalos becomes fuel for subsequent star formation, when they are incorporated into larger halos where efficient atomic cooling can induce gas condensation hence star formation. Since minihalos virialize at high redshift when the universe is not expected to have been significantly reionized, the implication is that gas in virialized minihalos may provide an abundant reservoir of primordial gas to possibly allow for the formation of Population-III metal-free stars to extend to much lower redshift than otherwise expected based on the enrichment of intergalactic medium.

Cosmological shock waves result from supersonic flow motions induced by hierarchical clustering of nonlinear structures in the universe. These shocks govern the nature of cosmic plasma through thermalization of gas and acceleration of nonthermal, cosmic-ray (CR) particles. We study the statistics and energetics of shocks formed in cosmological simulations of a concordance $\Lambda$CDM universe, with a special emphasis on the effects of non-gravitational processes such as radiative cooling, photoionization/heating, and galactic superwind feedbacks. Adopting an improved model for gas thermalization and CR acceleration efficiencies based on nonlinear diffusive shock acceleration calculations, we then estimate the gas thermal energy and the CR energy dissipated at shocks through the history of the universe. Since shocks can serve as sites for generation of vorticity, we also examine the vorticity that should have been generated mostly at curved shocks in cosmological simulations. We find that the dynamics and energetics of shocks are governed primarily by the gravity of matter, so other non-gravitational processes do not affect significantly the global energy dissipation and vorticity generation at cosmological shocks. Our results reinforce scenarios in which the intracluster medium and warm-hot intergalactic medium contain energetically significant populations of nonthermal particles and turbulent flow motions.

In this era of precision cosmology, a detailed physical understanding on the evolution of cosmic baryons is required. Cosmic magnetic fields, though still poorly understood, may represent an important component in the global cosmic energy flow that affects the baryon dynamics. We have developed an Eulerian-based cosmological magnetohydrodynamics code (CosmoMHD) with modern shock capturing schemes to study the formation and evolution of cosmic structures in the presence of magnetic fields. The code solves the ideal MHD equations as well as the non-equilibrium rate equations for multiple species, the Vlasov equation for dynamics of collisionless particles, the Poisson's equation for the gravitational potential field and the equation for the evolution of the intergalactic ionizing radiation field. In addition, a detailed star formation prescription and feedback processes are implemented. Several methods for solving the MHD by high-resolution schemes with finite-volume and finite-difference methods are implemented. The divergence-free condition of the magnetic fields is preserved at a level of computer roundoff error via the constraint transport method. We have also implemented a high-resolution method via dual-equation formulations to track the thermal energy accurately in very high Mach number or high Alfven-Mach number regions. Several numerical tests have demonstrated the efficacy of the proposed schemes.

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

It is shown that, if gas accretion via a disk onto the central supermassive black hole is efficient only for surface density Sigma >= 10 g/cm^2, the black hole mass - galactic bulge velocity dispersion relation (Tremaine et al 2002) is naturally borne out, so is the modest dispersion in that relation, in the context of hierarchical structure formation theory. This relation is not expected to evolve with redshift in this model.

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

We study the chemical and kinematic properties of the first galaxies which formed at a high redshift, using high resolution cosmological numerical simulations, and compared them with the recent observational results for the Sculptor dwarf spheroidal galaxy by Tolstoy et al., who found two distinct stellar populations: the lower metallicity stars are more spatially extended and possess a higher velocity dispersion than the higher metallicity stars. Our calculations reproduce these observations as the result of a steep metallicity gradient, within a single populations, induced by dissipative collapse of the gas component. We also predict strong [N/O] enhancements in the lowest metallicity stars in dwarf spheroidals, due to the preferential retention of ejected gas from intermediate mass stars, compared to Type II supernovae.

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

Population-III (Pop-III) starformation (SF) is thought to be quenched when the metallicity of the star-forming gas reaches a critical level. At high z, when the general intergalactic medium (IGM) was enriched with metals, the fraction of primordial gas already collapsed in minihalos was significantly larger than the fraction of primordial gas that had already been involved in Pop-III SF. We argue that this minihalo gas remained largely in a metal-free state, until these minihalos merged into large systems and formed stars. As a result, the era of Pop-III SF was significantly prolonged, leading to an integrated Pop-III SF an order of magnitude larger than expected for an abrupt transition redshift. The contribution of Pop-III SF to the reionization of hydrogen could have been significant until z~10 and may have extended to z~6. Our modeling allows for gradual enrichment of the IGM, feedback from photo-ionization and screening of reionization by minihalos. Nevertheless, extended Pop-III SF can result in complex, multi-peaked reionization histories. The contribution of Pop-III stars to reionization will be tested by the three-year WMAP results: (1) if Pop-III stars do not contribute to reionization, tau_es<0.05-0.06 and a rapid reionization at z~6 is expected; (2) if the product of star formation efficiency and escape fraction for Pop-III stars is significantly larger than for Pop-II stars, then a maximum tau_es=0.21 is achievable; (3) in a scenario where the product of star formation efficiency and escape fraction for Pop-III stars is comparable to that for Pop-II stars, tau_es=0.09-0.12 would be observed, with reionization histories characterized by an extended ionization plateau from z=7-12. This result holds regardless of the redshift where the IGM becomes enriched with metals.

Numerical simulations of the intergalactic medium have shown that at the present epoch a significant fraction (40-50%) of the baryonic component should be found in the (T~10^6K) Warm-Hot Intergalactic Medium (WHIM) - with several recent observational lines of evidence indicating the validity of the prediction. We here recompute the evolution of the WHIM with the following major improvements: (1) galactic superwind feedback processes from galaxy/star formation are explicitly included; (2) major metal species (O V to O IX) are computed explicitly in a non-equilibrium way; (3) mass and spatial dynamic ranges are larger by a factor of 8 and 2, respectively, than in our previous simulations. We find: (1) non-equilibrium calculations produce significantly different results from ionization equilibrium calculations. (2) The abundance of O VI absorption lines based on non-equilibrium simulations with galactic superwinds is in remarkably good agreement with latest observations, implying the validity of our model, while the predicted abundances for O VII and O VIII absorption lines appear to be lower than observed but the observational errorbars are currently very large. The expected abundances for O VI (as well as Lyman alpha), O VII and O VIII absorption systems are in the range 50-100 per unit redshift at EW=1km/s decreasing to 10-20 per unit redshift at EW=10km/s. The number of O VI absorption lines with EW>100km/s is very small, while there are about 1-3 lines per unit redshift for O VII and O VIII absorption lines at EW=100km/s. (3) Emission lines, primarily O VI and \lya in the UV and O VII and O VIII in the soft X-rays are potentially observable by future missions. The number of emission lines per unit redshift that may be detectable by planned UV and soft X-ray missions are in the order of 0.1-1.

Numerical simulations of the intergalactic medium have shown that at the present epoch a significant fraction (40-50%) of the baryonic component should be found in the (T~10^6K) Warm-Hot Intergalactic Medium (WHIM) - with several recent observational lines of evidence indicating the validity of the prediction. We here recompute the evolution of the WHIM with the following major improvements: (1) galactic superwind feedback processes from galaxy/star formation are explicitly included; (2) major metal species (O V to O IX) are computed explicitly in a non-equilibrium way; (3) mass and spatial dynamic ranges are larger by a factor of 8 and 2, respectively, than in our previous simulations. Here are the major findings: (1) galactic superwinds have dramatic effects, increasing the WHIM mass fraction by about 20%, primarily through heating up warm gas near galaxies with density 10^{1.5}-10^4 times the mean density. (2) the fraction of baryons in WHIM is increased modestly from the earlier work but is ~40-50%. (3) the gas density of the WHIM is broadly peaked at a density 10-20 times the mean density, ranging from underdense regions to regions that are overdense by 10^3-10^4. (4) the median metallicity of the WHIM is 0.18 Zsun for oxygen with 50% and 90% intervals being (0.040,0.38) and (0.0017,0.83).

First galaxies formed within halos of mass M=E7.5-E9 Msun at z=30-40 in the standard cold dark matter (CDM) universe may each display an extended hydrogen 21-cm absorption halo against the cosmic microwave background with a brightness temperature decrement of del T=-(100-150)mK at a radius 0.3 < r < 3.0 comoving Mpc, corresponding to an angular size of 10-100 arcseconds. A 21-cm tomographic survey in the redshift shell z=30-40 (at 35-45MHz), which could be carried out by the next generation of radio telescopes, is expected to be able to detect millions of first galaxies and may prove exceedingly profitable in enabling (at least) four fundamental applications for cosmology and galaxy formation. First, it may yield direct information on star formation physics in first galaxies. Second, it could provide a unique and sensitive probe of small-scale power in the standard cosmological model hence physics of dark matter and inflation. Third, it would allow for an independent, perhaps "cleaner" characterization of interesting features on large scales in the power spectrum such as the baryonic oscillations. Finally, possibly the most secure, each 21-cm absorption halo is expected to be highly spherical and faithfully follow the Hubble flow. By applying the Alcock-Paczynski test to a significant sample of first galaxies, one may be able to determine the dark energy equation of state with an accuracy likely only limited by the accuracy with which the matter density can be determined independently.

We analyze the SDSS Ly-alpha forest P_F(k,z) measurement to determine the linear theory power spectrum. Our analysis is based on fully hydrodynamic simulations, extended using hydro-PM simulations. We account for the effect of absorbers with damping wings, which leads to an increase in the slope of the linear power spectrum. We break the degeneracy between the mean level of absorption and the linear power spectrum without significant use of external constraints. We infer linear theory power spectrum amplitude Delta^2_L(k_p=0.009s/km,z_p=3.0)=0.452_{-0.057-0.116}^{+0.069+0.141} and slope n_eff=-2.321_{-0.047-0.102}^{+0.055+0.131} (possible systematic errors are included through nuisance parameters in the fit - a factor >~5 smaller errors would be obtained on both parameters if we ignored modeling uncertainties). The errors are correlated and not perfectly Gaussian, so we provide a chi^2 table to accurately describe the results. The result corresponds to sigma_8=0.85, n=0.94, for a LCDM model with Omega_m=0.3, Omega_b=0.04, and h=0.7, but is most useful in a combined fit with the CMB. The inferred curvature of the linear power spectrum and the evolution of its amplitude and slope with redshift are consistent with expectations for LCDM models, with the evolution of the slope, in particular, being tightly constrained. We use this information to constrain systematic contamination, e.g., fluctuations in the UV background. This paper should serve as a starting point for more work to refine the analysis, including technical improvements such as increasing the size and number of the hydrodynamic simulations, and improvements in the treatment of the various forms of feedback from galaxies and quasars.

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

The physical basis for the belief that the abundant old dwarf galaxies seen in present-day galaxy clusters formed before reionization of the universe is very compelling, because (1) the observed faint end slope of the galaxy luminosity function ascends from a shallower slope at brighter magnitudes and tracks that of dark matter halo mass function, and (2) that steep power-law slope is observed to extend all the way to galaxies of inferred velocity dispersion of ~10km/s. We then show that the number of ionizing photons emitted by these dwarf galaxies should be able to produce an H II region at least as large as each proto-cluster region, occupying at least 20% of the entire IGM volume at z >= 6. A more likely scenario is that our estimate of the ionizing photon production rate based on these dwarfs under-estimates the true rate from proto-clusters by a factor of ~5 for a variety of physically plausible reasons, in which case a fully self- consistent picture emerges. Specific implications include (1) Strong clustering of sources would collectively produce very large individual H II regions of sizes up to ~100Mpc before overlapping with others, having important implications for upcoming radio and CMB observations. There should be large spatial fluctuations of neutral fraction in the intergalactic medium up to cluster scales of ~100Mpc. We give an H II region size function. (2) The ionizing sources at z~6 are expected to be highly biased with a bias factor >=15. (3) Large H II regions would enable largely unattenuated transmission of flux of a significant number of sources prior to reionization. When more sensitive observations become available, one should expect to see clusters of fainter sources in the vicinity of known galaxies at z>6.

The observed luminosity function (LF) of satellite galaxies shows several interesting features that require a better understanding of gas-thermodynamic processes and feedback effects related to reionization and galaxy formation. In galaxy clusters, the abundance of dwarf galaxies is in good agreement with the expectation based on the subhalo mass function, whereas in galaxy groups, the relatively small abundance of dwarfs conflicts with theoretical expectations. In all halo systems, there is a dip in the abundance of galaxies with luminosities in the range ~ 2x10^8 L_sun to 10^10 L_sun, corresponding to subhalo mass scales between ~ 5x10^10 M_sun to few times 10^11 M_sun. Photoionization from reionization has been used to explain statistics of the dwarf population, with larger systems forming prior to, and smaller systems forming subsequent to, reionization. The observed dip in the LF is an imprint of small dwarf galaxies (<2x10^8 L_sun) that formed prior to reionization. The galactic winds powered by supernovae in these dwarf galaxies propagate energy and metals to large distances such that the intergalactic medium is uniformly enriched to a level of 10^-3 zsun, as observed in the low-redshift Ly-alpha forest. The associated energy related to this metallicity raises the intergalactic medium temperature and the Jeans mass to a range 10^10-10^11 M_sun at z ~ 3.4-6.0. Because the epoch of nonlinearity for halos in this mass range is at z > 3.4-4.4, their gas content, hence star formation, is greatly suppressed on average and leads to a dip in the observed LF at $z=0$. Larger halos (M > 10^11 M_sun), becoming nonlinear at z < 3.4-4.4, have masses topping the Jeans mass, where subhalo mass function based LF is again in agreement with observations.

We present a simple method for tracing the spatial distribution and predicting the physical properties of the Warm-Hot Intergalactic Medium (WHIM), from the map of galaxy light in the local universe. Under the assumption that biasing is local and monotonic we map the ~ 2 Mpc/h smoothed density field of galaxy light into the mass density field from which we infer the spatial distribution of the WHIM in the local supercluster. Taking into account the scatter in the WHIM density-temperature and density-metallicity relation, extracted from the z=0 outputs of high-resolution and large box size hydro-dynamical cosmological simulations, we are able to quantify the probability of detecting WHIM signatures in the form of absorption features in the X-ray spectra, along arbitrary directions in the sky. To illustrate the usefulness of this semi-analytical method we focus on the WHIM properties in the Virgo Cluster region.

We present predictions of Ly-alpha forest-galaxy correlations at z=3 from Eulerian simulations that include the effects of galactic winds, driven primarily by supernova explosions. Galactic winds produce expanding bubbles of shock-heated gas within ~0.5 comoving Mpc/h of luminous galaxies in the simulation, which have space density similar to that of observed Lyman break galaxies (LBGs). However, most of the low-density intergalactic gas that determines the observed properties of the Ly-alpha forest is unaffected by winds. The impact of winds on the Ly-alpha optical depth near galaxies is less dramatic than their impact on gas temperature because winds heat only a small fraction of the gas present in the turnaround regions surrounding galaxies. Hence, Ly-alpha absorption from gas outside the wind bubbles is spread out over the same velocity range occupied by the wind-heated gas. In general, Ly-alpha absorption is expected to be stronger than average near galaxies because of the high gas density. Winds result in a modest reduction of this expected increase of Ly-alpha absorption. Our predictions can be compared to future observations to detect the wind effects and infer their strength, although with the caveat that the results are still dependent on the correspondence of simulated galaxies and observed LBGs. We find that wind effects in our simulations are not strong enough to reproduce the high Ly-alpha transmission within 0.5 Mpc/h comoving of galaxies that has been suggested by recent observations; powerful galactic explosions or ejecta with hyper-escape velocities would be required, but these are unlikely to be produced by ordinary star formation and supernovae alone

The transition from Population III to Population II stars is determined by the presence of a sufficient amount of metals, in particular, oxygen and carbon. The vastly different yields of these relevant metals between different initial stellar mass functions would then cause such a transition to occur at different times. We show that the transition from Pop III to Pop II stars is likely to occur before the universe can be reionized, if the IMF is entirely very massive stars (M > 140 M_sun). A factor of about 10 more ionizing photons would be produced in the case with normal top-heavy IMF (e.g., M ~ 10-100 M_sun), when such a transition occurs. Thus, a high Thomson optical depth (tau_e >= 0.11-0.14) may be indication that the Population III stars possess a more conventional top-heavy IMF.

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

At redshifts beyond z>6, as the mean fraction of neutral hydrogen x_HI in the intergalactic medium (IGM) increases, the line flux of Lyman alpha (Lya) emitters can be significantly suppressed, which can result in a decrease in the observed number of emitters at a given Lya flux. However, cosmological HII regions surrounding the Lya emitting galaxies alleviate these effects. We use simple models of the Lya line suppression that incorporate the presence of HII regions to predict the overall effect of the Lya absorption on the Lya luminosity function. We find, in agreement with other recent studies, that unless ionizing sources are unusually strongly clustered, a fully neutral IGM may be inconsistent with the large abundance of confirmed z=6.5 Lya emitters. However, the presence of local HII regions prohibits placing a tight constraint on the mean neutral fraction. We find x_HI < 0.25; the presence of strong winds and/or the clustering of ionizing sources would further weaken this constraint. We conclude that the evolution of the Lya LF is consistent with reionization occurring near this redshift. Finally, we suggest that a measurement of observed Lya line width as a function of the Lya luminosity, in a future, larger sample of Lya emitters, may serve as a robust diagnostic of the neutral fraction in the IGM.

Cosmological shock waves play a major role in heating baryonic gas in the large scale structure of the universe. In this contribution we study the shock-heated gas in connection with shocks themselves. The distributions of shock speed and temperature of shock-hated gas both should reflect the depth of gravitational potential of associated nonlinear structures, so their morphology. From their distributions we find, in addition to hot gas in and around clusters/groups and warm-hot intergalactic medium (WHIM) with T=10^5-10^7 K mostly in filaments, there is a significant amount of low temperature WHIM with T < 10^5 K distributed mostly as sheet-like structures. Not only the WHIM with T=10^5-10^7 K but also the WHIM with T < 10^5 K make up a significant fraction of gas mass, implying the low temperature WHIM could be important in resolving the missing baryon problem. The shock-heated gas in filaments and sheets are manifested best through emissions and absorptions in soft X-ray and far UV. We confirm that the WHIM with T=10^5-10^7 K makes significant contributions to the soft X-ray background, absorptions of highly ionized species such as OVII and OVIII in AGN spectra, and line emissions from OVII and OVIII ions, as pointed by previous studies. However, the WHIM with T < 10^5 K is the major contributor to absorptions of lower ionized species such as OV and OVI, because these photo-ionized ions are most abundant in sheets of low density and temperature. (abridged)

The confirmation of the cluster soft excess (CSE) by XMM-Newton has rekindled interest as to its origin. The recent detections of CSE emission at large cluster radii together with reports of OVII line emission associated with the CSE has led many authors to conjecture that the CSE is, in fact, a signature of the warm-hot intergalactic medium (WHIM). In this paper we test the scenario by comparing the observed properties of the CSE with predictions based on models of the WHIM. We find that emission from the WHIM in current models is 3 to 4 orders of magnitude too faint to explain the CSE. We discuss different possibilities for this discrepancy including issues of simulation resolution and scale, and the role of small density enhancements or galaxy groups. Our final conclusion is that the WHIM alone is unlikely to be able to accout for the observed flux of the CSE.

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

We explore several physical effects on the power spectrum of the Lyman-alpha forest transmitted flux. The effects we investigate here are usually not part of hydrodynamic simulations and so need to be estimated separately. The most important effect is that of high column density absorbers with damping wings, which add power on large scales. We compute their effect using the observational constraints on their abundance as a function of column density. Ignoring their effect leads to an underestimation of the slope of the linear theory power spectrum. The second effect we investigate is that of fluctuations in the ionizing radiation field. For this purpose we use a very large high resolution N-body simulation, which allows us to simulate both the fluctuations in the ionizing radiation and the small scale LyaF within the same simulation. We find an enhancement of power on large scales for quasars and a suppression for galaxies. The strength of the effect rapidly increases with increasing redshift, allowing it to be uniquely identified in cases where it is significant. We develop templates which can be used to search for this effect as a function of quasar lifetime, quasar luminosity function, and attenuation length. Finally, we explore the effects of galactic winds using hydrodynamic simulations. We find the wind effects on the LyaF power spectrum to be be degenerate with parameters related to the temperature of the gas that are already marginalized over in cosmological fits. While more work is needed to conclusively exclude all possible systematic errors, our results suggest that, in the context of data analysis procedures where parameters of the LyaF model are properly marginalized over, the flux power spectrum is a reliable tracer of cosmological information.

We combine the constraints from the recent Ly-alpha forest and bias analysis of the SDSS with previous constraints from SDSS galaxy clustering, the latest supernovae, and WMAP . Combining WMAP and the lya we find for the primordial slope n_s=0.98\pm 0.02. We see no evidence of running, dn/d\ln k=-0.003\pm 0.010, a factor of 3 improvement over previous constraints. We also find no evidence of tensors, r<0.36 (95% c.l.). A positive correlation between tensors and primordial slope disfavors chaotic inflation type models with steep slopes: V \propto \phi^4 is outside the 3-sigma contour. For the amplitude we find sigma_8=0.90\pm 0.03 from the lyaf and WMAP alone. We find no evidence of neutrino mass: for the case of 3 massive neutrino families with an inflationary prior, \sum m_{\nu}<0.42eV and the mass of lightest neutrino is m_1<0.13eV at 95% c.l. For the 3 massless + 1 massive neutrino case we find m_{\nu}<0.79eV for the massive neutrino, excluding at 95% c.l. all neutrino mass solutions compatible with the LSND results. We explore dark energy constraints in models with a fairly general time dependence of dark energy equation of state, finding Omega_lambda=0.72\pm 0.02, w(z=0.3)=-0.98^{+0.10}_{-0.12}, the latter changing to w(z=0.3)=-0.92^{+0.09}_{-0.10} if tensors are allowed. We find no evidence for variation of the equation of state with redshift, w(z=1)=-1.03^{+0.21}_{-0.28}. These results rely on the current understanding of the lyaf and other probes, which need to be explored further both observationally and theoretically, but extensive tests reveal no evidence of inconsistency among different data sets used here.

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

We quantify and discuss the footprints of neutral hydrogen in the intergalactic medium (IGM) on the spectra of high-redshift (z ~ 6) sources, using mock spectra generated from hydrodynamical simulations of the IGM. We show that it should be possible to extract relevant parameters, including the mean neutral fraction in the IGM, and the radius of the local cosmological Stromgren region, from the flux distribution in the observed spectra of distant sources. We focus on quasars, but a similar analysis is applicable to galaxies and gamma ray burst (GRB) afterglows. We explicitly include uncertainties in the spectral shape of the assumed source template near the Lyman alpha line. Our results suggest that a mean neutral hydrogen fraction, x(HI) of unity can be statistically distinguished from x(HI)<0.01, by combining the spectra of tens of bright (M = -27) quasars. Alternatively, the same distinction can be achieved using the spectra of several hundred sources that are ~100 times fainter. Furthermore, if the radius of the Stromgren sphere can be independently constrained to within ~10 percent, this distinction can be achieved using a single source. The information derived from such spectra will help in settling the current debate as to what extent the universe was reionized at redshifts near z=6.

A new model for cooling flows in X-ray clusters, capable of naturally explaining salient features observed, is proposed. The only requirement is that a significant relativistic component, in the form of cosmic rays (CR), be present in the intra-cluster medium and significantly frozen to the thermal gas. Such an addition qualitatively alters the conventional isobaric thermal instability criterion such that a fluid parcel becomes thermally stable when its thermal pressure drops below a threshold fraction of its CR pressure. Consequently, the lowest possible temperature at any radius is about one third of the ambient temperature {\it at that radius}, exactly as observed, In addition, we suggest that dissipation of internal gravity waves, excited by radial oscillatory motions of inward drifting cooling clouds about their radial equilibrium positions, may be responsible for heating up cooling gas. With the ultimate energy source for powering the cooling X-ray luminosity and heating up cooling gas being gravitational due to inward drifting cooling clouds as well as the general inward flow, heating is spatially distributed and energetically matched with cooling. One desirable property of this heating mechanism is that heating energy is strongly centrally concentrated, providing the required heating for emission-line nebulae.

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

The recently discovered z=10 galaxy (Pello et al 2004) has a strong Lyman alpha emission line that is consistent with being surprisingly symmetric, even given the relatively poor quality of its spectrum. The blue wing of a Lyman alpha line originating at high redshift should be strongly suppressed by resonant hydrogen absorption along the line of sight, an expectation borne out by the observed asymmetric shapes of the existing sample of Lyman alpha emitting sources at lower redshifts (3< z < 6.7). Absorption on the blue side of the line of the Pello et al. source could be reduced if the intergalactic medium (IGM) in the vicinity of the galaxy is highly ionized, but we show that this requires an unrealistically high ionizing emissivity. We suggest instead that the Lyman alpha emitting gas be receding relative to the surrounding gas with a velocity of >= 35km/s, a large velocity that is plausible only if the galaxy is part of a larger system (group of galaxies) with a velocity dispersion >= 35km/s. We find that with this velocity shift, the observed strength and shape of the line is still consistent with the galaxy being surrounded by its own Stromgren sphere embedded in a fully neutral IGM. More generally, we predict that at any given redshift, the bright Lyman alpha emitters with broader lines would exhibit stronger asymmetry than fainter ones. Bright galaxies with symmetric Lyman alpha lines may be signposts for groups and clusters of galaxies, within which they can acquire random velocities comparable to or larger than their linewidths.

We compute the properties of dark matter halos with mass $10^{6.5}-10^9\msun$ at redshift $z=6-11$ in the standard cold dark matter cosmological model, utilizing a very high resolution N-body simulation. We find that dark matter halos in these mass and redshift ranges are significantly biased over matter with a bias factor in the range 2-6. The dark matter halo mass function displays a slope of $2.05\pm 0.15$ at the small mass end. We do not find a universal dark matter density profile. Instead, we find a significant dependence of the central density profile of dark matter halos on halo mass and epoch with $\alpha_0=0.4-1.0$; the high-mass ($M\ge 10^8\msun$) low-redshift ($z\sim 6$) halos occupy the high end of the range and low-mass ($M\sim 10^{7}\msun$) high-redshift ($z\sim 11$) halos occupy the low end. Additionally, for fixed mass and epoch there is a significant dispersion in $\alpha_0$ due to the stochastic assembly of halos. Our results fit a relationship of the form $\alpha_0=0.75((1+z)/7.0)^{-1.25}(M/10^7\msun)^{0.11(1+z)/7.0}$ with a dispersion about this fit of $\pm 0.5$ and no systematic dependence of variance correlated with environment. The median spin parameter of dark matter halos is $0.03-0.04$ but with a large lognormal dispersion of $\sim 0.4$. Various quantities are tabulated or fitted with empirical formulae.

We show that large, high-redshift (z>10) galaxies with virial temperature in excess of 10^4K may be mostly comprised of cold atomic clouds which were formerly minihalos. These clouds move at a speed of 15-30km/s and collide with one another on a time scale of 10^7yr. The supersonic collisions result in spatially distributed star formation with high efficiency. We then show that most, subsequent star formation in cold atomic clouds may be triggered by shockwaves launched from the first stars formed in collisions. Those shockwave-compressed clouds are even more widespread spatially because of large imparted velocities and some can possibly escape into the intergalactic medium. More importantly, with respect to cosmological reionization, widespread star formation would allow a much higher ionizing photon escape fraction. These favorable conditions may form the physical basis to enable the standard cosmological model to produce a reasonably high Thomson optical depth \tau_e=0.10-0.14. A chain reaction of star formation in minihalos in the intergalactic space may be triggered by explosions in the intergalactic medium, if minihalos are strongly clustered. In this case, a still higher \tau_e would be achievable.

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

A scenario is put forth for the formation of supermassive black holes at the centers of galaxies. It depends upon the formation of a Mestel disk with a flat rotation curve, where the mass is proportional to r and thickness to 1/r. Such disks should form from the collapse of uniformly rotating, isolated, gaseous clouds, either proto-galactic, galaxy-mass damped Lyman alpha clouds or the gas that survives galaxy mergers. We propose that in any case the disk will be unstable to the Rossby vortex instability (RVI). This instability grows from any large, steep pressure gradient in an optically thick disk. Such pressure gradients either occur adjacent to compact objects or could be triggered by individual supernovae in and around the disk. Upon excitation, the RVI transports angular momentum far more efficiently than the turbulence of a viscosity-based Shakura-Sunyaev disk. The critical column density condition in a L* Mestel disk predicts that the black hole mass is ~ 10^8 Msolar and is also proportional to the 1/4 power of the velocity dispersion, in accordance with observations.

The intrinsic Lyman-alpha emission lines of normal galaxies before reionization are much absorbed by the damping wing of the Gunn-Peterson trough, rendering their direct detection nearly impossible, if their intrinsic line widths are less than ~100km/s. High redshift luminous quasars prior to the completion of cosmological reionization at z~6, on the other hand, are capable of producing large HII regions around them (Stromgren spheres) to allow their intrinsic Lyman-alpha emission lines to be transmitted without overwhelming absorption (Cen & Haiman 2000). We suggest that targeted observations at the Stromgren spheres of known luminous quasars at z >= 6 would be able to detect Lyman-alpha emission lines of galaxies inside the Stromgren spheres largely unattenuated. A tunable, very narrowband filter of \Delta\lambda\over \lambda ~ 0.1% or a narrowband filter of \Delta\lambda\over \lambda ~1% with follow-up spectroscopic identifications will be required. Such observations could directly observe the sources of cosmological reionization including possibly the Pop III galaxies at z=6-20 by JWST. Possible applications include determinations of the ionization state of the intergalactic medium, the sizes of the Stromgren spheres, the ages of the quasars, the luminosity function of high redshift galaxies and its evolution, the spatial distribution of galaxies and its evolution, the biased distribution of galaxies around quasars and the anisotropy of quasar emission. Observations using Keck-class telescopes may already be made to enable a differentiation between a fully neutral and a 10% neutral intergalactic medium at z>6.

Utilizing a new, high mass resolution hydrodynamic simulation we compute the metallicity evolution of damped Lyman alpha systems (DLAs) and find a reasonable agreement with observations. In particular, the observed slow evolution of the DLA metallicity occurs naturally in the simulation due to the combined effects of physical and observational selection. The slow metallicity evolution is caused by the steady transformation, with increasing time, of the highest metallicity systems to "galaxies", thus depleting this category, while all the lower metallicity systems show, individually, an increase in metallicity. Although the trend of DLA metallicity with redshift is in good agreement with observations, it appears that the average metallicity of simulated DLAs is higher than observed by 0.3-0.5 dex in the probed redshift range (z=0-5). Our study indicates that this difference may be attributed to observational selection effects due to dust obscuration. If we allow for a dust obscuration effect, our model reproduces the observed metallicity evolution in both amplitude and slope. We find that DLAs are not a simple population but probe a range of different systems and the mix changes with redshift. About 50% of all metals in the gaseous phase is in DLAs at all times from z=5 to z=1, making a rapid downturn at z\le 1 to ~20% by z=0, as metals are swept into the hotter components of the IGM as well as locked up in stars. While not the primary focus of this study, we find that the model provides good matches to observations with respect to column density distribution and evolution of neutral gas content, if the same dust obscuration is taken into account. We find Omega_{DLA,comp}=(1-3)E-3, depending on the effect of dust obscuration.

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

We investigate the promise of the Ly-alpha forest for high precision cosmology in the era of the Sloan Digital Sky Survey using low order N-point statistics. We show that with the existing data one can determine the amplitude, slope and curvature of the slope of the matter power spectrum with a few percent precision. Higher order statistics such as the bispectrum provide independent information that can confirm and improve upon the statistical precision from the power spectrum alone. The achievable precision is comparable to that from the cosmic microwave background with upcoming satellites, and complements it by measuring the power spectrum amplitude and shape at smaller scales. Since the data cover the redshift range 2<z<4, one can also extract the evolution of the growth factor and Hubble parameter over this range, and provide useful constraints on the presence of dark energy at z>2.

In an earlier paper (Cen 2003) we pointed out the strong likelihood for the universal reionization to occur twice, giving rise to a larger Thomson optical depth. Here we perform a more focused analysis of the Thomson optical depth in light of the WMAP observations. While the current uncertainties on the observed Thomson optical depth are still relatively large, with tau_e=0.17+-0.04 (68%) (Kogut et al 2003), important implications on Pop III star formation processes at high redshift can already be inferred. We are able to draw four conclusions: (1) in the absence of a top-heavy initial stellar mass function (IMF) for Pop III metal-free stars and without a dramatic upturn in the star formation efficiency and ionizing photon escape fraction at high redshift (z>6), we find tau_e =< 0.09; (2) with a top-heavy IMF for the Pop III metal-free stars and plausible star formation efficiency and ionizing photon escape fraction, it is expected that tau_e =< 0.12; (3) it is possible to reach tau_e = 0.15, if the metal enrichment efficiency of the intergalactic medium by Pop III stars is very low thus Pop III era is prolonged; (4) to reach tau_e >= 0.17 requires either of the following two conditions: the cosmological model power index n is positively tilted to n >= 1.03, Pop III star formation in minihalos with molecular hydrogen cooling has an efficiency c_*(H_2,III)>0.01 (with ionizing photon escape fraction greater than 30%). If the current observed value of Thomson optical depth withstands future data, we will have strong observational evidence that Pop III stars are massive and their formation efficiency may be much higher than current theoretical works suggest. Alternatively, there may be unknown, non-stellar ionizing sources at very high redshift.

We propose a possible way to detect baryons at low redshifts from the analysis of X-ray absorption spectra of bright AGN pairs. A simple semi-analytical model to simulate the spectra is presented. We model the diffuse warm-hot intergalactic medium (WHIM) component, responsible for the X-ray absorption, using inputs from high-resolution hydro-dynamical simulations and analytical prescriptions. We show that the number of OVII absorbers per unit redshift with column density larger than $10^{13.5}$ cm$^{-2}$ - corresponding to an equivalent width of $\sim$ 1 km/s - which will be possibly detectable by {\it XEUS}, is $\magcir 30$ per unit redshift. {\it Constellation-X} will detect $\sim 6$ OVII absorptions per unit redshift with an equivalent width of 10 km/s. Our results show that, in a $\Lambda$CDM Universe, the characteristic size of these absorbers at $z\sim 0.1$ is $\sim 1$ $h^{-1}$ Mpc. The filamentary structure of WHIM can be probed by finding coincident absorption lines in the spectra of background AGN pairs. We estimate that at least 20 AGN pairs at separation $\mincir 20$ arcmin are needed to detect this filamentary structure at a 3$\sigma$ level. Assuming observations of distant sources using {\it XEUS} for exposure times of 500 ksec, we find that the minimum source flux to probe the filamentary structure is $\sim 2\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, in the 0.1-2.4 keV energy band. Thus, most pairs of these extragalactic X-ray bright sources have already been identified in the {\it ROSAT} All-Sky Survey. Re-observation of these objects by future missions could be a powerful way to search for baryons in the low redshift Universe.

We show the universe was reionized twice, first at z~15-16 and second at z~6. Such an outcome appears inevitable, when normalizing to two well determined observational measurements, namely, the epoch of the final cosmological reionization at z~6 and the density fluctuations at z~6, which in turn are tight ly constrained by Lyman alpha forest observations at z~3. These two observations most importantly fix the product of star formation efficiency and ionizing photon escape fraction from galaxies at high redshift. To the extent that the relative star formation efficiencies in gaseous minihalos with H2 cooling and large halos with atomic cooling at high redshift are still unknown, the primary source for the first reionization could be Pop III stars either in minihalos or in large halos. We show that gas in minihalos can be cooled efficiently by H2 molecules and star formation can continue to take place largely unimpeded throughout the first reionization period, thanks to two new mechanisms for generating a high X-ray background during the Pop III era, put forth here. Moreover, an important process for producing a large number of H2 molecules in relic HII regions of Pop III galaxies, first pointed out by Ricotti, Gnedin, & Shull, is quantified here. It is shown that the Lyman-Werner background may never build up during the Pop III era. The long cosmological reionization and reheating history is complex. We discuss a wide range of implications and possible tests for this new reionization picture. In particular, Thomson scattering optical depth is increased to 0.10 +- 0.03, compared to 0.027 for the case of only one rapid reionization at z=6. Upcoming Microwave Anisotropy Probe observation of the polarization of the cosmic microwave background should be able to distinguish between these two scenarios.

The recent discovery of bright quasars around redshift z=6 suggests that black holes (BHs) with masses in excess of 10^9 Msun have already assembled at a very early stage in the evolution of the universe. An alternative interpretation is that these quasars are powered by less massive BHs, but their fluxes are strongly magnified through gravitational lensing by intervening galaxies. Here we analyze the flux distribution of the Ly alpha emission of the quasar with the highest known redshift, SDSS 1030+0524, at z=6.28. We show that this object could not have been magnified by lensing by more than a factor of five. The constraint arises from the large observed size, 30 (comoving) Mpc, of the ionized region around this quasar, and relies crucially only on the assumption that the quasar is embedded in a largely neutral IGM. Based on the line/continuum ratio of SDSS 1030+0524, we argue further that this quasar also cannot be beamed by a significant factor. We conclude that the minimum mass for its resident BH is 4 x 10^8 Msun (for magnification by a factor of five); if the mass is this low, then the quasars had to switch on prior to redshift z=9. From the size of the ionized region, we are also able to place an absolute lower bound on the age of this quasar at t > 2 x 10^7 years.

We present predictions of the correlation between the Lyman-alpha forest absorption in quasar spectra and the mass within \sim 5 Mpc/h (comoving) of the line of sight, using fully hydrodynamic and hydro-PM numerical simulations of the cold dark matter model supported by present observations. The observed correlation based on galaxies and the Lya forest can be directly compared to our theoretical results, assuming that galaxies are linearly biased on large scales. Specifically, we predict the average value of the mass fluctuation, <delta_m>, conditioned to a fixed value of the Lya forest transmitted flux delta_F, after they have been smoothed over a 10 Mpc/h cube and line of sight interval, respectively. We find that <delta_m>/sigma_m as a function of delta_F/sigma_F has a slope of 0.6 at this smoothing scale, where sigma_m and sigma_F are the rms dispersions (this slope should decrease with the smoothing scale). We show that this value is largely insensitive to the cosmological model and other Lya forest parameters. Comparison of our predictions to observations should provide a fundamental test of our ideas on the nature of the Lya forest and the distribution of galaxies, and can yield a measurement of the bias factor of any type of galaxies that are observed in the vicinity of Lya forest lines of sight.

We have developed an algorithm for transferring radiation in three-dimensional space. The algorithm computes radiation source and sink terms using the Fast Fourier Transform (FFT) method, based on a formulation in which the integral of any quantity (such as emissivity or opacity) over any volume may be written in the classic convolution form. The algorithm is fast with the computational time scaling as N (log N)^2, where N is the number of grid points of a simulation box, independent of the number of radiation sources. Furthermore, in this formulation one can naturally account for both local radiation sources and diffuse background as well as any extra external sources, all in a self-consistent fashion. Finally, the algorithm is completely stable and robust. While the algorithm is generally applicable, we test it on a set of problems that encompass a wide range of situations in cosmological applications, demonstrating that the algorithm is accurate. These tests show that the algorithm produces results that are in excellent agreement with analytic expectations in all cases. In particular, radiation flux is guaranteed to propagate in the right direction, with the ionization fronts traveling at the correct speed with an error no larger than one cell for all the cases tested. The total number of photons is conserved in the worst case at 10% level and typically at 1-5% level over hundreds of time steps. As an added advantage, the accuracy of the results depends weakly on the size of the time step, with a typical cosmological hydrodynamic time step being sufficient.

The temporal evolution of the ionizing UV background radiation field at high redshift provides a probe of the evolution of the early star formation rate. By comparing the observed levels of absorption in the highest redshift quasar spectra to the predictions of a hydrodynamic simulation, we determine the evolution of the photoionization rate Gamma for neutral hydrogen in the intergalactic medium, over the redshift range 4.9 < z < 6.1. After accounting for sampling variance, we infer a sharp increase in Gamma from z~6.1 to z~5.8, probably implying reionization at this redshift. We find a decrease in Gamma from z~5.6 to 5.2, at 3sigma significance. This feature is a generic signature in the aftermath of reionization, entirely consistent with the cosmological reionization process being completed at z~6.1. This generic feature is a result of a significant change in the star formation rate subsequent to the cosmological reionization. There is an abrupt rise of the temperature of the intergalactic medium due to photo-heating, when it is reionized. This translates to a correspondingly sudden jump in the Jeans mass and a dramatic suppression of gas accretion onto the most abundant (sub-galactic) halos at the epochs of interest. The star formation rate suffers a temporary setback in the aftermath of reionization, resulting in a temporary decrease in the amplitude of the ionizing radiation field.

Combining conservation of energy throughout nearly-spherical collapse of galaxy clusters with the virial theorem, we derive the mass-temperature relation for X-ray clusters of galaxies $T=CM^{2/3}$. The normalization factor $C$ and the scatter of the relation are determined from first principles with the additional assumption of initial Gaussian random field. We are also able to reproduce the recently observed break in the M-T relation at $T \sim 3 \keV$, based on the scatter in the underlying density field for a low density $\Lambda$CDM cosmology. Finally, by combining observational data of high redshift clusters with our theoretical formalism, we find a semi-empirical temperature-mass relation which is expected to hold at redshifts up to unity with less than 20% error.

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

Hydrodynamic simulations of growth of cosmic structure suggest that 30-50% of the total baryons at z=0 may be in a warm-hot intergalactic medium (WHIM) with temperatures ~10^5-10^7K. The O VI \lambda \lambda 1032, 1038 absorption line doublet in the FUV portion of QSO spectra provides an important probe of this gas. Utilizing recent hydrodynamic simulations, it is found that there should be ~5 O VI absorption lines per unit redshift with equivalent widths >= 35 mA, decreasing rapidly to ~0.5 per unit redshift at >= 350 mA. About 10% of the total baryonic matter or 20-30% of the WHIM is expected to be in the O VI absorption line systems with equivalent width >= 20 mA; the remaining WHIM gas may be too hot or have too low metallicity to be detected in O VI. We find that the simulation results agree well with observations with regard to the line abundance and total mass contained in these systems. Some of the O VI systems are collisionally ionized and some are photoionized, but most of the mass is in the collisionally ionized systems. We show that the gas that produces the O VI absorption lines does not reside in virialized regions such as galaxies, groups, or clusters of galaxies, but rather has an overdensity of 10-40 times the average density. These regions form a somewhat connected network of filaments. The typical metallicity of these regions is 0.1-0.3Zsun.

Gas rich sub-galactic halos with mass Mt <= 10^7.5 Msun, while incapable of forming stars due to lack of adequate coolants, contain a large fraction of baryonic mass at cosmological reionization. We show that the reionization of the universe at z=10-20 has an interesting physical effect on these halos. The external radiation field causes a synchronous inward propagation of an ionization front towards each halo, resulting in an inward, convergent shock. The resident gas of mass Mb~10^4-10^7 Msun in low spin (initial dimensionless spin parameter lambda <= 0.01) halos with a velocity dispersion sigmav <= 11km/s would be compressed by a factor of ~100 in radius and form self -gravitating baryonic systems. Under the assumption that such compressed gaseous systems fragment to form stars, the final stellar systems will have a size 2-40pc, velocity dispersion 1-10km/s and a total stellar mass of M* 10^3-10^6 Msun. The characteristics of these proposed systems seem to match the observed properties of halo globular clusters. The expected number density is consistent with the observed number density of halo globular clusters. The observed mass function of slope ~-2 at the high mass end is predicted by the model. Strong correlation between velocity dispersion and luminosity (or surface brightness) and lack of correlation between velocity dispersion and size, in agreement with observations, are expected. Metallicity is, on average, expected to be low and should not correlate with any other quantities of globular clusters, except that a larger dispersion of metallicity among globular clusters is expected for larger galaxies. The observed trend of specific frequency with galaxy type may be produced in the model. We suggest that these stellar systems are seen as halo globular clusters today.

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

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

There appears to exist a conflict between the standard structure formation theory and the existence of the faint dwarf spheroidal galaxies in the Local Group. Theoretical considerations suggest that a cold dark matter universe has been a very hostile place for the formation of small galaxies. In particular, gas would not have been able to cool and condense to form stars inside dark matter halos with a velocity dispersion <10km/s since before the cosmological reionization epoch. In other words, one should not expect to see any major star formation activities in dwarf galaxies with a velocity dispersion <10km/s during the past 12Gyr, according to the current structure formation theory. The real universe, on the other hand, shows that all dwarf spheroidal galaxies in the Local Group have velocity dispersions <10km/s and extended and recent star formation activities are quite common in them. This apparent conflict between cold dark matter models and local observations can be resolved, if one allows the cold dark matter particles to decay relatively recently. The resolution comes about in that the dwarf spheroidal galaxies with a velocity dispersion of \~10km/s seen today were born, in this picture, in a small fraction of halos with a velocity dispersion of ~20km/s, where gas can be retained and is able to cool and contract to form stars. The presently observed lower velocity dispersion ~10km/s of these dwarf spheroidal galaxies is a consequence of the decay of one half of the dark matter in and subsequent expansion of the halos since redshift z~2.

At low redshifts, a census of the baryons in all known reservoirs falls a factor of two to four below the total baryon density predicted from Big Bang nucleosynthesis arguments and observed light element ratios. Recent cosmological hydrodynamic simulations suggest that a significant fraction of these missing baryons could be in the form of warm/hot gas in the filaments and halos within which most field galaxies are embedded. With the release of source count results from Chandra and recent detections of this gas in O VI quasar absorption lines, it becomes interesting to examine the predictions and limits placed on this component of the X-ray background (XRB). We have used new hydrodynamical simulations to predict the total X-ray spectrum from the gas in the 100 eV to 10 keV range. We find that, when uncertainties in the normalization of the observed XRB and the value of Omega_b are taken into account, our results are consistent with current observational limits placed on the contribution of emission from gas to the XRB. In the 0.5-2 keV range, we expect the contribution from this component to be 0.63 10^{-12} erg s^-1 cm^-2 deg^-2 or between 6% and 18% of the extragalactic surface brightness. The peak fraction occurs in the 0.5-1 keV range where the predicted line emission mirrors a spectral bump seen in the latest ASCA/ROSAT XRB data.

Hydrodynamic simulations of Lyman alpha clouds based on ab inito cosmological models have produced results that are in broad agreement with observations. However, further analyses have revealed that, with progressively higher numerical resolution, the median or cutoff line width of the simulated Lyman alpha clouds (i.e. the Doppler parameter) appears to converge to a value significantly below what is observed at z~3 (by about a factor of 1.5). These convergence test simulations do not include feedback from star formation. Given the observed metallicity in the Lyman alpha clouds we suggest that supernovae, which presumably polluted the IGM with metals, may have deposited a sufficient amount of energy in the IGM to reconcile the theory with observations. Simple arguments immediately narrow the redshift range of pollution down to 4<z_{dep}<9. It seems quite certain that dwarf and sub-dwarf galaxies with total masses in the range 10^{6.5-9.0}Msun have to be largely responsible for the pollution. Furthermore, it is implied that either star formation is very efficient or metal yield is very high for these early dwarf galaxies, if the mean metallicity in the universe at z=3 is as high as <Z>=0.01Zsun. Finally, assuming the specific supernova heating energy is proportional to the metallicity of a gas, we note that the picture proposed here would be consistent with supernovae being the apparently needed heating source for the intra-cluster gas, if the required heating of the intra-cluster gas is no greater than 1 keV per particle.

The canonical cosmological constant dominated cold dark matter model (LCDM) may possess too much power on small scales at z=0, manifested as central over-concentration of dark matter and over-abundance of dwarf galaxies. We suggest an alternative model, LDCDM, where one half of the cold dark matter particles decay into relativistic particles by z=0. The model successfully lowers the concentration of dark matter in dwarf galaxies as well as in large galaxies like our own at low redshift, while simultaneously retaining the virtues of the LCDM model. The model solves the problem of over-production of small dwarf galaxies in the LCDM not by removing them but by identifying them with failed, "dark" galaxies, where star-formation is quenched due to dark matter evaporation and consequent halo expansion. A dramatic difference between the LDCDM model and other proposed variants of the LCDM model is that the small-scale power at high redshift (z>2) in the LDCDM model is enhanced compared to the LCDM model. A COBE-and-cluster normalized LDCDM model can be constructed with the following parameters: H0=60km/sec/Mpc, lambda0=0.60, Omega_{0,CDM}=0.234, Omega_{0,b}=0.044, n=1.0, and sigma8=1.06. A clean test of this model can be made by measuring the evolution of gas fraction in clusters. The prediction is that the gas fraction should decrease with redshift and is smaller by 31% at z=1 than at z=0. X-ray and Sunyaev-Zel'dovich effect observations should provide such a test.

Ionizing sources embedded in the neutral intergalactic medium (IGM) before cosmological reionization generate discrete HII regions. We show that a sufficiently bright quasar (for example, one tenth as luminous as that recently discovered by Fan et al.) can ionize a large volume, allowing the transmission of a substantial fraction of the flux of its intrinsic Lyman alpha emission line on both sides of the Lyman alpha wavelength. The observed line profile is richly informative. We show that a sufficiently accurate, high spectral resolution ($R\approx 10^4$) measurement of the line profile of a bright quasar is feasible using the Next Generation Space Telescope (NGST) as well as large ground-based telescopes. Such a measurement has two potentially important applications. First, the red side of the Lyman alpha emission line provides a way to measure the quasar lifetime. Second, the blue side provides a direct measure of the density fluctuations of the intergalactic medium at the quasar redshift. The estimate of the absorption of the red side is limited by the accuracy to which the intrinsic profile of the Lyman alpha emission line is known. The blue side, however, does not sensitively depend on the intrinsic profile, because the former is much narrower than the latter.

Approximately 30-40% of all baryons in the present day universe reside in a warm-hot intergalactic medium (WHIM), with temperatures between 10^5<T<10^7 K. This is a generic prediction from six hydrodynamic simulations of currently favored structure formation models having a wide variety of numerical methods, input physics, volumes, and spatial resolutions. Most of these warm-hot baryons reside in diffuse large-scale structures with a median overdensity around 10-30, not in virialized objects such as galaxy groups or galactic halos. The evolution of the WHIM is primarily driven by shock heating from gravitational perturbations breaking on mildly nonlinear, non-equilibrium structures such as filaments. Supernova feedback energy and radiative cooling play lesser roles in its evolution. WHIM gas is consistent with observations of the 0.25 keV X-ray background without being significantly heated by non-gravitational processes because the emitting gas is very diffuse. Our results confirm and extend previous work by Cen & Ostriker and Dave' et al.

We use large-scale cosmological simulations to estimate the mass-to-light ratio of galaxy systems as a function of scale, and compare the results with observations of galaxies, groups, clusters, and superclusters of galaxies. We find remarkably good agreement between observations and simulations. Specifically, we find that the simulated mass-to-light ratio increases with scale on small scales and flattens to a constant value on large scales, as suggested by observations. We find that while mass typically follows light on large scales, high overdensity regions --- such as rich clusters and superclusters of galaxies --- exhibit higher M/L_B values than average, while low density regions exhibit lower M/L_B values; high density regions are thus antibiased in M/L_B, with mass more strongly concentrated than blue light. This is true despite the fact that the galaxy mass density is unbiased or positively biased relative to the total mass density in these regions. The M/L_B antibias is likely due to the relatively old age of the high density regions, where light has declined significantly since their early formation time, especially in the blue band which traces recent star formation. Comparing the simulated results with observations, we place a powerful constraint on the mass density of the universe; using, for the first time, the entire observed mass-to-light function, from galaxies to superclusters, we find Omega =0.16+/-0.05.

The evolution of the temperature in the intergalactic medium is related to the reionization of hydrogen and helium, and has important consequences for our understanding of the Lya forest and of galaxy formation in gravitational models of large-scale structure. We measure the temperature-density relation of intergalactic gas from Lya forest observations of eight quasar spectra with high resolution and signal-to-noise ratio, using a new line fitting technique to obtain a lower cutoff of the distribution of line widths from which the temperature is derived. We carefully test the accuracy of this technique to recover the gas temperature with a hydrodynamic simulation. The temperature at redshift z=(3.9, 3.0, 2.4) is best determined at densities slightly above the mean: T_star=(20200\pm2700, 20200\pm1300, 22600\pm1900)K (statistical error bars) for gas density (in units of the mean density) Delta_star=(1.42\pm0.08, 1.37\pm0.11, 1.66\pm0.11). The power-law index of the temperature-density relation, defined by T=T_star(Delta/Delta_star)^{gamma-1}, is gamma-1= (0.43\pm0.45, 0.29\pm0.30, 0.52\pm0.14) for the same three redshifts. The temperature at the fixed over-density Delta=1.4 is T_1.4=(20100\pm2800, 20300\pm1400, 20700\pm1900)K. These temperatures are higher than expected for photoionized gas in ionization equilibrium with a cosmic background, and can be explained by a gradual additional heating due to on-going HeII reionization. The measurement of the temperature reduces one source of uncertainty in the lower limit to the baryon density implied by the observed mean flux decrement. We find that the temperature cannot be reliably measured for under-dense gas, because the velocities due to expansion always dominate the widths of the corresponding weak lines.

We have examined the properties of shock waves in simulations of large scale structure formation for two cosmological scenarios (a SCDM and a LCDM with Omega =1). Large-scale shocks result from accretion onto sheets, filaments and Galaxy Clusters (GCs) on a scale of circa 5 Mpc/h in both cases. Energetic motions, both residual of past accretion history and due to current asymmetric inflow along filaments, generate additional, common shocks on a scale of about 1 Mpc/h, which penetrate deep inside GCs. Also collisions between substructures inside GCs form merger shocks. Consequently, the topology of the shocks is very complex and highly connected. During cosmic evolution the comoving shock surface density decreases, reflecting the ongoing structure merger process in both scenarios. Accretion shocks have very high Mach numbers (10-10^3), when photo-heating of the pre-shock gas is not included. The typical shock speed is of order v_{sh}(z) =H(z)lambda_{NL}(z), with lambda_{NL}(z) the wavelength scale of the nonlinear perturbation at the given epoch. However, the Mach number for shocks occuring within clusters is usually smaller (3-10), due to the fact that the intracluster gas is already hot. Statistical fits of shock speed around GCs as a function of GCs temperature give power-law's in accord with 1-D predictions. However, a very different result is obtained for fits of the shock radius, reflecting the very complex shock structures forming in 3-D simulations. The in-flowing kinetic energy across such shocks, giving the power available for cosmic-ray acceleration, is comparable to the cluster X-ray luminosity emitted from a central region of radius 0.5 Mpc/h. Considering their large size and long lifetimes, those shocks are potentially interesting sites for cosmic-ray acceleration, if modest magnetic fields exist within them.

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

The formation rate of luminous galaxies seems to be roughly constant from z~2 to z~4 from the recent observations of Lyman break galaxies (LBGs) (Steidel et al 1999). The abundance of luminous quasars, on the other hand, appears to drop off by a factor of more than twenty from z~2 to z~5 (Warren, Hewett, & Osmer 1994; Schmidt, Schneider, & Gunn 1995). The difference in evolution between these two classes of objects in the overlapping, observed redshift range, z=2-4, can be explained naturally, if we assume that quasar activity is triggered by mergers of luminous LBGs and one quasar lifetime is ~10^{7-8} yrs. If this merger scenario holds at higher redshift, for the evolutions of these two classes of objects to be consistent at z>4, the formation rate of luminous LBGs is expected to drop off at least as rapidly as exp(-(z-4)^{6/5}) at z>4.

We analyze a new large-scale ($100h^{-1}$Mpc) numerical hydrodynamic simulation of the popular $\Lambda$CDM cosmological model, including in our treatment dark matter, gas and star-formation, on the basis of standard physical processes. The method, applied with a numerical resolution of $<200h^{-1}$kpc (which is still quite coarse for following individual galaxies, especially in dense regions), attempts to estimate where and when galaxies form. We then compare the smoothed galaxy distribution with the smoothed mass distribution to determine the "bias" defined as $b\equiv (\delta M/M)_{gal}/(\delta M/M)_{total}$ on scales large compared with the code numerical resolution (on the basis of resolution tests given in the appendix of this paper). We find that (holding all variables constant except the quoted one) bias increases with decreasing scale, with increasing galactic age or metallicity and with increasing redshift of observations. At the $8h^{-1}$Mpc fiducial comoving scale bias (for bright regions) is 1.35 at $z=0$ reaching to 3.6 at $z=3$, both numbers being consistent with extant observations. We also find that $(10-20)h^{-1}$Mpc voids in the distribution of luminous objects are as observed (i.e., observed voids are not an argument against CDM-like models) and finally that the younger systems should show a colder Hubble flow than do the early type galaxies (a testable proposition). Surprisingly, little evolution is found in the amplitude of the smoothed galaxy-galaxy correlation function (as a function of {\it comoving} separation). Testing this prediction vs observations will allow a comparison between this work and that of Kauffmann et al which is based on a different physical modelingmethod.

A sample of eight quasars observed at high resolution and signal-to-noise is used to determine the probability distribution function (PDF), the power spectrum, and the correlation function of the transmitted flux in the \lya forest, in three redshift bins centered at z=2.41, 3.00, and 3.89. All the results are presented in tabular form, with full error covariance matrices to allow for comparisons with any numerical simulations and with other data sets. The observations are compared with a numerical simulation of the \lya forest of a Lambda-CDM model with Omega=0.4, known to agree with other large-scale structure observational constraints. There is excellent agreement for the PDF, if the mean transmitted flux is adjusted to match the observations. A small difference between the observed and predicted PDF is found at high fluxes and low redshift, which may be due to the uncertain effects of fitting the spectral continuum. Using the numerical simulation, we show how the flux power spectrum can be used to recover the initial power spectrum of density fluctuations. From our sample of eight quasars, we measure the amplitude of the mass power spectrum to correspond to a linear variance per unit ln(k) of $\Delta^2_\rho(k)=0.72\pm0.09$ at k=0.04(km/s)^{-1} and z=3, and the slope of the power spectrum near the same k to be $n_p=-2.55\pm0.10$ (statistical error bars). The results are statistically consistent with Croft et. al. (1999), although our value for the rms fluctuation is lower by a factor 0.75. For the Lambda-CDM model we use, the implied primordial slope is $n=0.93\pm0.10$, and the normalization is $\sigma_8=0.68+1.16(0.95-n)\pm0.04$.

During hierarchical clustering, smaller masses generally collapse earlier than larger masses and so are denser on the average. The core of a small mass halo could be dense enough to resist disruption and survive undigested, when it is incorporated into a bigger object. We explore the possibility that a nested sequence of undigested cores in the center of the halo, which have survived the hierarchical, inhomogeneous collapse to form larger and larger objects, determines the halo structure in the inner regions. For a flat universe with $P(k) \propto k^n$, scaling arguments then suggest that the core density profile is, $\rho \propto r^{-\alpha}$ with $\alpha = (9+3n)/(5+n)$. But whether such behaviour obtains depends on detailed dynamics. We first examine the dynamics using a fluid approach to the self-similar collapse solutions for the dark matter phase space density, including the effect of velocity dispersions. We highlight the importance of tangential velocity dispersions to obtain density profiles shallower than $1/r^2$ in the core regions. If tangential velocity dispersions in the core are constrained to be less than the radial dispersion, a cuspy core density profile shallower than 1/r cannot obtain, in self-similar collapse. We then briefly look at the profiles of the outer halos in low density cosmological models where the total halo mass is convergent. Finally, we analyze a suite of dark halo density and velocity dispersion profiles obtained in cosmological N-body simulations of models with n= 0, -1 and -2. We find that the core-density profiles of dark halos, show considerable scatter in their properties, but nevertheless do appear to reflect a memory of the initial power spectrum, with steeper initial spectra producing flatter core profiles. (Abridged)

There was a bright, transient companion spot to SN1987A with a projected distance of about 17 light-days, observed by optical speckle interferometry one to two months after explosion. It is shown here that the bright spot may be due to a receding ultra-relativistic jet traveling at 53 degress to the observer-to-SN1987A vector, through a circumstellar medium of density profile rho(r) ~ r^-2. If it had approached us along the line of sight, a very bright gamma-ray burst would have been seen with an apparent isotropic energy of about 10^54 erg and an openning angle of a few degrees. The model provides an adequate explanation for the evolution of the spot, although there are still problems in explaining its observed color. This model implies that at least some GRBs would be seen as going through a medium with density rho(r) ~ r^-2 rather than a uniform medium, which is frequently adopted in GRB calculations. Improved analysis of the speckle data has revealed another and fainter spot on the opposite side.

It is shown that differential Doppler shift of different patches of the blastwave front in the fireball model at varying angles to the line of sight could provide an intrinsic smoothing mechanism for the spectra of gamma-ray bursts (GRBs) and the associated afterglows at lower energy bands. For the model parameters of interest, it is illustrated that a monochromatic spectrum at $\nu$ in blastwave comoving frame is smoothed and observed to have a half-width-half-maximum (HWHM) of $\sim (0.6-2.0)\nu$. Some other implications of this smoothing process are discussed. In particular, if the circumburster medium is uniform and electron and magnetic energies are fixed fractions of the total post-shock energy with time, the observed GRB or afterglow spectra cannot be steeper than $\nu^{-\alpha}$ with $\alpha=0.75-1.25$, regardless of the intrinsic spectra in the comoving shock frame; i.e., a very steep electron distribution function (with $p>3$) could produced an observed spectrum with $\alpha\sim 1$. In addition, a generic fast-rise-slow-decay type of GRB temporal profile is expected.

I present an overview of some of the recent advances in our understanding of the distribution and evolution of the ordinary, baryonic matter in the universe. Two observations that strongly suggest that most of the baryons seen at high redshift ($z\ge 2$) have turned into some forms yet undetected at $z=0$ are highlighted. With the aid of large-scale cosmological hydrodynamic simulations, it is shown that most of the baryons today are in a gaseous form with a temperature of $10^5-10^7~$Kelvin -- the "warm/hot gas", shock heated during the gravitational collapse and formation of the large-scale structure at low redshift. Primarily line emissions from this warm/hot gas may account for a large fraction of the residual (after removal of identifiable discrete sources) soft X-ray background at $h\nu < 1.0$keV. How this warm/hot gas may be detected by the next generation of EUV and soft X-ray instruments is indicated. Detection or non-detection of this warm/hot gas will have profound implications for cosmology.

Using large-scale hydrodynamic simulations with heuristic criteria for galaxy formation, we investigate how the galaxy field is related to the mass density and the gas temperature. We find that the bias b(R) = sigma_g(R)/sigma(R), where sigma_g(R) is the variance of galaxy counts in spheres of radius R and sigma(R) is the same for mass, varies from 2.6 at 1 h^{-1} Mpc to 1.2 at 30 h^{-1} Mpc. Including the dependence of the galaxy density on local gas temperature as well as on local mass density can fully account for this scale dependence. Galaxy density depends on temperature because gas which is too hot cannot cool to form galaxies; this causes scale dependence of b(R) because local gas temperature is related to the gravitational potential, and thus contains information about the large scale density field. We show that temperature dependence generally causes b(R) to vary on quasilinear and nonlinear scales, indicating that scale dependence of bias may be a generic effect in realistic galaxy formation scenarios. We find that the relationship between the galaxy and mass density fields is also a function of galaxy age. On large scales, the older galaxies are highly biased (b ~ 1.7) and highly correlated (r ~ 1.0) with the mass density field; younger galaxies are not biased (b ~ 0.8) and are poorly correlated (r ~ 0.5) with the mass. We argue that linear bias is inadequate to describe the relationship between galaxies and mass. (Abridged)

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

Numerical simulations of standard cosmological scenarios have now reached the degree of sophistication required to provide tentative answers to the fundamental question: Where and when were the heavy elements formed? Averaging globally, these simulations give a metallicity that increases from 1% of the solar value at $z=3$ to 20% at present. This conclusion is, in fact, misleading, as it masks the very strong dependency of metallicity on local density. At every epoch higher density regions have much higher metallicity than lower density regions. Moreover, the highest density regions quickly approach near solar metallicity and then saturate, while more typical regions slowly catch up. These results are much more consistent with observational data than the simpler picture (adopted by many) of gradual, quasi-uniform increase of metallicity with time.

The clustering of galaxies relative to the mass distribution declines with time because: first, nonlinear peaks become less rare events; second, the densest regions stop forming new galaxies because gas there becomes too hot to cool and collapse; third, after galaxies form, they are gravitationally ``debiased'' because their velocity field is the same as the dark matter. To show these effects, we perform a hydrodynamic cosmological simulation and examine the density field of recently formed galaxies as a function of redshift. We find the bias b_* of recently formed galaxies (the ratio of the rms fluctuations of these galaxies and mass), evolves from 4.5 at z=3 to around 1 at z=0, on 8 h^{-1} Mpc comoving scales. The correlation coefficient r_* between recently formed galaxies and mass evolves from 0.9 at z=3 to 0.25 at z=0. As gas in the universe heats up and prevents star formation, star-forming galaxies become poorer tracers of the mass density field. After galaxies form, the linear continuity equation is a good approximation to the gravitational debiasing, even on nonlinear scales. The most interesting observational consequence of the simulations is that the linear regression of the star-formation density field on the galaxy density field evolves from about 0.9 at z=1 to 0.35 at z=0. These effects also provide a possible explanation for the Butcher-Oemler effect, the excess of blue galaxies in clusters at redshift z ~ 0.5. Finally, we examine cluster mass-to-light ratio estimates of Omega, finding that while Omega(z) increases with z, one's estimate Omega_est(z) decreases. (Abridged)

A new method to determine the bias parameter of galaxies relative to matter is suggested. The method is based on the assumption that gravity is the dominating force which determines the formation of the structure in the Universe. Due to gravitational instability the galaxy formation is a threshold process: in low-density environments galaxies do not form and matter remains in primordial form. We investigate the influence of the presence of void and clustered populations to the power spectrum of matter and galaxies. The power spectrum of galaxies is similar to the power spectrum of matter; the fraction of total matter in the clustered population determines the difference between amplitudes of fluctuations of matter and galaxies, i.e. the bias factor. To determine the fraction of matter in voids and clustered population we perform numerical simulations. The fraction of matter in galaxies at the present epoch is found using a calibration through the sigma_8 parameter.

Observed power spectrum of matter found in Papers I and II is compared with analytical power spectra. Spatially flat cold and mixed dark matter models with cosmological constant and open models are considered. The primordial power spectrum of matter is determined using the power spectrum of matter and the transfer functions of analytical models. The primordial power spectrum has a break in amplitude. We conclude that a scale-free primordial power spectrum is excluded if presently available data on the distribution of clusters and galaxies represent the true mass distribution of the Universe.

We calculate the mean power spectrum of galaxies using published power spectra of galaxies and clusters of galaxies. On small scales we use the power spectrum derived from the 2-dimensional distribution of APM galaxies, on large scales we use power spectra derived from 3-dimensional data for galaxy and cluster samples. Spectra are reduced to real space and to the amplitude of the power spectrum of APM galaxies. Available data indicate the presence of two different populations in the nearby Universe. Clusters of galaxies sample a relatively large region in the Universe where rich, medium and poor superclusters are well represented. Their mean power spectrum has a spike on scale 120 h^{-1}Mpc, followed by an approximate power-law spectrum of index n = -1.9 towards small scales. The power spectrum found from LCRS and IRAS 1.2 Jy surveys is flatter around the maximum, which may represent regions of the Universe with medium-rich and poor superclusters.

We perform a variety of tests to determine the numerical resolution of the cosmological TVD eulerian code developed by Ryu et al (1993). Tests include 512^3 and 256^3 simulations of a Pk=k^{-1} spectrum to check for self-similarity and comparison of results with those from higher resolution SPH and grid-based calculations (Frenk et al 1998). We conclude that in regions where density gradients are not produced by shocks the code degrades resolution with a Gaussian smoothing (radius) length of 1.7 cells. At shock caused gradients (for which the code was designed) the smoothing length is 1.1 cells. Finally, for \beta model fit clusters, we can approximately correct numerical resolution by the transformation R^2_{core}\to R^2_{core}-(C\Delta l)^2, where \Delta l is the cell size and C=1.1-1.7. When we use these corrections on our previously published computations for the SCDM and \Lambda CDM models we find luminosity weighted, zero redshift, X-ray cluster core radii of (210\pm 86, 280\pm 67)h^{-1}kpc, respectively, which are marginally consistent with observed (Jones & Forman 1992) values of 50-200h^{-1}kpc. Using the corrected core radii, the COBE normalized SCDM model predicts the number of bright L_x>10^{43}erg/s clusters too high by a factor of \sim 20 and the \Lambda CDM model is consistent with observations.

New, high resolution, large-scale, cosmological hydrodynamic galaxy formation simulations of a standard cold dark matter model (with a cosmological constant) are utilized to predict the distribution of baryons at the present and at moderate redshift. It is found that the average temperature of baryons is an increasing function of time, with most of the baryons at the present time having a temperature in the range 10^{5-7} K. Thus, not only is the universe dominated by dark matter, but more than one half of the normal matter is yet to be detected. Detection of this warm/hot gas poses an observational challenge, requiring sensitive EUV and X-ray satellites. Signatures include a soft, cosmic X-ray background, apparent warm components in hot clusters due to both intrinsic warm intra-cluster gas and warm inter-cluster gas projected onto clusters along the line of sight, absorption lines in X-ray and UV quasar spectra [e.g., O VI (1032,1038)A lines, OVII 574 eV line], strong emission lines (e.g., O VIII 653 eV line) and low redshift, broad, low column density $\lya$ absorption lines. We estimate that approximately 1/4 of the extragalactic soft X-ray background (SXRB) (at 0.7 keV) arises from the warm/hot gas, half of it coming from $z<0.65$ and three-quarters from $z<1.00$, so the source regions should be identifiable on deep optical images.

A scenario is proposed that explains both the observed high pulsar velocities and extragalactic gamma-ray bursts (GRBs). The model involves an ultra- relativistic jet from a supernova (SN), that produces a GRB and its afterglow, whose characteristics are similar to an isotropic fireball GRB perhaps with some differences at late times in the afterglow once some significant transverse diffusion has occurred. The time scales and many other properties of GRBs and their afterglows in this model are consistent with observations. GRBs in this model have special intrinsic properties, that can either falsify or prove this model unambiguously by observations. The most direct proof is the detection of a SN about the same time as the luminous GRB event. Most GRBs and SNe are expected occur at moderate redshift ($z\sim 1-3$), if they follow the observed universal star formation history, as implied in this model. Searching for GRB/SN associations is a challenge, because majority of the SNe will be faint. Some additional, dramatic observable consequences are predicted, which can also be utilized to test the model.

The COBE microwave background temperature fluctuations and the abundance of local rich clusters of galaxies provide the two most powerful constraints on cosmological models. When all variants of the standard cold dark matter (CDM) model are subject to the combined constraint, the power spectrum of any model is fixed to $\sim 10%$ accuracy in both the shape and overall amplitude. These constrained models are not expected to differ dramatically in their local large-scale structure properties. However, their evolutionary histories differ, resulting in dramatic differences towards high redshift. We examine in detail six standardized, COBE and cluster normalized CDM models with respect to a large set of independent observations. The observations include correlation function of rich clusters of galaxies, galaxy power spectrum, evolution of rich cluster abundance, gravitational lensing by moderate -to-high redshift clusters, $\lya$ forest, damped $\lya$ systems, high redshift galaxies, reionization of the universe and future CMB experiments. It seems that each of the independent observations examined is or potentially is capable of distinguishing between at least some of the models. The combined power of several or all of these observations is tremendous. Thus, we appear to be on the verge of being able to make dramatic tests of all models in the near future using a rapidly growing set of observations, mostly at moderate to high redshift. Consistency or inconsistency between different observed phenomena on different scales and/or at different epochs with respect to the models will have profound implications for theory of growth of cosmic structure.

(shortened) We develop and test a method to compute mass and auto-correlation functions of rich clusters of galaxies from linear density fluctuations, based on the formalism of Gaussian peaks (Bardeen et al 1986). The essential, new ingredient in the current approach is a simultaneous and unique fixture of the size of the smoothing window for the density field, $r_f$, and the critical height for collapse of a density peak, $\delta_c$, for a given cluster mass (enclosed within the sphere of a given radius rather than the virial radius, which is hard to measure observationally). Of these two parameters, $r_f$ depends on both the mass of the cluster in question and $\Omega$, whereas $\delta_c$ is a function of only $\Omega$ and $\Lambda$. These two parameters have formerly been treated as adjustable and approximate parameters. Thus, for the first time, the Gaussian Peak Method (GPM) becomes unambiguous, and more importantly, accurate, as is shown here. We apply this method to constrain all variants of the Gaussian cold dark matter (CDM) cosmological model using the observed abundance of local rich clusters of galaxies and the microwave background temperature fluctuations observed by COBE. The combined constraint fixes the power spectrum of any model to $\sim 10%$ accuracy in both the shape and overall amplitude. To set the context for analyzing CDM models, we choose six representative models of current interest, including an $\Omega_0=1$ tilted cold dark matter model, a mixed hot and cold dark matter model with 20% mass in neutrinos, two lower density open models with $\Omega_0=0.25$ and $\Omega_0=0.40$, and two lower density flat models with $(\Omega_0=0.25, \Lambda_0=0.75)$ and $(\Omega_0=0.40,\Lambda_0=0.60)$. This suite of CDM models

This study shows one important effect of preexistent cosmic microwave background temperature fluctuations on the determination of the Hubble constant through Sunyaev-Zel'dovich effect of clusters of galaxies, especially when coupled with the gravitational lensing effect by the same clusters. The effect results in a broad distribution of the apparent Hubble constant. The combination of this effect with other systematic effects such as the Loeb-Refregier Effect seems to provide an explanation for the observationally derived values of the Hubble constant currently available based on the Sunyaev-Zel'dovich effect, if the true value of the Hubble constant is 60-80 km/s/Mpc. It thus becomes possible that the values of the Hubble constant measured by other techniques which generally give a value around 60-80 km/s/Mpc be reconciled with the SZ effect determined values of the Hubble constant, where are systematically lower than others and have a broad distribution.

We propose a new method of analysis for the \lya forest, namely to measure the 1-point and 2-point joint probability distribution of the transmitted flux. The results for a sample of seven observed quasars and from two simulations of structure formation are shown and compared. Statistically significant differences in the 2-point function between the results of the numerical simulations and the observations are easily found. The analysis we suggest is very simple to apply to observed data sets, and we discuss its superiority over the traditional Voigt-profile fitting algorithms for accurate comparison to the predictions of theoretical models.

We present a method for determining the rms mass fluctuations on 8 h^-1 Mpc scale, sigma8. The method utilizes the rate of evolution of the abundance of rich clusters of galaxies. Using the Press-Schechter approximation, we show that the cluster abundance evolution is a strong function of sigma8: d log n/dz ~ -1/sigma8^2; low sigma8 models evolve exponentially faster than high sigma8 models, for a given mass cluster. For example, the number density of Coma-like clusters decreases by a factor of ~10^3$ from z = 0 to z ~ 0.5 for sigma8=0.5 models, while the decrease is only a factor of ~5 for sigma8 ~ 1. The strong exponential dependence on sigma8 arises because clusters represent rarer density peaks in low sigma8 models. We show that the evolution rate at z < 1 is insensitive to the density parameter $\Omega$ or to the exact shape of the power spectrum. Cluster evolution therefore provides a powerful constraint on sigma8. Using available cluster data to z ~ 0.8, we find sigma8 = 0.83 +/- 0.15. This amplitude implies a bias parameter b ~ 1/sigma8 = 1.2 +/- 0.2, i.e., a nearly unbiased universe with mass approximately tracing light on large scales.

We examine the correlational property of $\lya$ clouds in detail and compare it to that of mass and galaxies. We show that the correlation strength of $\lya$ clouds is somewhat weaker than that of the underlying matter, which in turn is weaker than that of galaxies (biased galaxy formation). On the scales probed, $10-300$km/s, for $\lya$ clouds, we find that higher density, higher optical depth, higher column density regions are more strongly clustered than lower density, lower optical depth, lower column density regions. Thus, a consistent picture seems to emerge: the correlation strength for a given set of objects is positively correlated with their characteristic global density and the differences among the correlations of galaxies, $\lya$ clouds and mass reflect the differences in density that each trace. Significant positive correlations with a strength of $0.1-1.0$ are found for $\lya$ clouds in the velocity range $50-300$km/s. Among the correlational measures examined, an optical depth correlation function (Equation 5) proposed here may serve as the best correlational measure. It reasonably faithfully represents the true correlation of the underlying matter, enabling a better indication of both matter correlation and the relationship between galaxies and $\lya$ clouds. Furthermore, it appears to be an alternative to the conventional line-line correlation function with the virtue that it does not require ambiguous post-observation fitting procedures such as those commonly employed in the conventional line-finding methods. Neither does it depend sensitively on the observational resolution (e.g., FWHM), insofaras the clouds are resolved Conveniently, it can be easily measured with the current observational sensitivity without being contaminated significantly by the presence of noise.

Gaussian cosmological models, typified by the inflationary cold dark matter models, and non-Gaussian topological defect based cosmological models, such as the texture seeded model, differ in the origin of large-scale cosmic structures. In the former it is believed that peaks at appropriate scales in the initial high density field are the sites onto which matter accretes and collapses to form the present galaxies and clusters of galaxies, whereas in the latter these structures can form around the density perturbation seeds (which are textures in the texture model). Textures initially are randomly distributed on scales larger than their size, in sharp contrast to the initial high density peaks in the Gaussian models which are already strongly clustered before any gravitational evolution has occured. One thus expects that the resultant correlation of large cosmic objects such as clusters of galaxies in the texture model should be significantly weaker than its Gaussian counterpart. We show that an $\Omega_0=1$ biased $b=2$ (as required by cluster abundance observations) texture model (or any random seed model) predicts a two-point correlation length of $\le 6.0h^{-1}$Mpc for rich clusters, independent of richness. On the other hand, the observed correlation length for rich clusters is $\ge 10.0h^{-1}$Mpc at an approximately $2\sigma$ confidence level. It thus appears that the global texture cosmological model or any random seed cosmological models of this sort are ruled out at a very high confidence ($>3\sigma$).

The standard method of measuring the galaxy pairwise velocity dispersion on small scales is heavily weighted by the densest regions in a way that is difficult to calibrate. We propose a new statistic which measures the small-scale velocity dispersion as an explicit function of density. Computing this statistic for a volume-limited subsample of the Optical Redshift Survey, we find that the small-scale velocity dispersion rises from 220 to 760 km/s as density increases. We calculate this statistic for a series of mock catalogs drawn from a hydrodynamic simulation of an \Omega h = 0.5 Cold Dark Matter universe (standard CDM), and find that the observed velocity distribution lies ~1 \sigma below the simulations in each of eight density bins, formally ruling out this model at the 7.4 \sigma level, quantifying the well-known problem that this model produces too high a velocity dispersion. This comparison is insensitive to the normalization of the power spectrum, although it is quite sensitive to the density and velocity bias of galaxies relative to dark matter on small scales. (abridged)

We show that the evolution of the number density of rich clusters of galaxies breaks the degeneracy between Omega (the mass density ratio of the universe) and sigma_{8} (the normalization of the power spectrum), sigma_{8}Omega^{0.5} \simeq 0.5, that follows from the observed present-day abundance of rich clusters. The evolution of high-mass (Coma-like) clusters is strong in Omega=1, low-sigma_{8} models (such as the standard biased CDM model with sigma_{8} \simeq 0.5), where the number density of clusters decreases by a factor of \sim 10^{3} from z = 0 to z \simeq 0.5; the same clusters show only mild evolution in low-Omega, high-sigma_{8} models, where the decrease is a factor of \sim 10. This diagnostic provides a most powerful constraint on Omega. Using observations of clusters to z \simeq 0.5-1, we find only mild evolution in the observed cluster abundance. We find Omega = 0.3 \pm 0.1 and sigma_{8} = 0.85 \pm 0.15 (for Lambda = 0 models; for Omega + Lambda = 1 models, Omega = 0.34 \pm 0.13). These results imply, if confirmed by future surveys, that we live in a low-den sity, low-bias universe.

We study the projection effects on various observables of clusters of galaxies at redshift near zero, including cluster richness, velocity dispersion, X-ray luminosity, three total mass estimates (velocity-based, temperature-based and gravitational lensing derived), gas fraction and substructure. Unlike previous studies focusing on the Abell clusters, we conservatively assume that both optical and X-ray observations can determine the source (galaxy or hot X-ray gas) positions along the line of sight as well as in the sky plane accurately; hence we only include sources inside the velocity space defined by the cluster galaxies (filtered through the pessimistic 3 sigma clipping algorithm) as possible contamination sources. Projection effects are found to be important for some quantities but insignificant for others. We show that, on average, the gas to total mass ratio in clusters appears to be 30-40% higher than its corresponding global ratio. Independent of its mean value, the broadness of the observed distribution of gas to total mass ratio is adequately accounted for by projection effects, alleviating the need to invoke (though not preventing) other non gravitational physical processes. While the moderate boost in the ratio narrows the gap, it is still not quite sufficient to reconcile the standard nucleosynthesis value of $\Omega_b=0.0125(H_0/100)^{-2}$ and $\Omega=1$ with the observed gas to mass ratio value in clusters of galaxies, $0.05(H_0/100)^{-3/2}$, for any plausible value of $H_0$. However, it is worth noting that real observations of X-ray clusters, especially X-ray imaging observations, may be subject to more projection contaminations than we allow for in our analysis. In contrast, the X-ray luminosity of a cluster within a radius $\leq 1.0h^{-1}$Mpc is hardly

This study analyzes the sizes, shapes and correlations of $\lya$ clouds produced by a hydrodynamic simulation of a spatially flat CDM universe with a non-zero cosmological constant ($\Omega_0=0.4$, $\Lambda_0=0.6$, $\sigma_8 =0.79$), over the redshift range $2\le z \le 4$. The $\lya$ clouds range in size from several kiloparsecs to about a hundred kiloparsecs in proper units, and they range in shape from roundish, high column density regions with $\nhi\ge 10^{15} cm^{-2}$ to low column density sheet-like structures with $\nhi \le 10^{13} cm^{-2}$ at z=3. The most common shape found in the simulation resembles that of a flattened cigar. The physical size of a typical cloud grows with time roughly as $(1+z)^{-3/2}$ while its shape hardly evolves (except for the most dense regions $\rho_{cut}>30$). Our result indicates that any simple model with a population of spheres (or other shapes) of a uniform size is oversimplified; if such a model agrees with observational evidence, it is probably only by coincidence. We also illustrate why the use of double quasar sightlines to set lower limits on cloud sizes is useful only when the perpendicular sightline separation is small ($\Delta r \le 50h^{-1}$ kpc). Finally, we conjecture that high column density $\lya$ clouds ($\nhi\ge 10^{15} cm^{-2}$) may be the progenitors of the lower redshift faint blue galaxies. This seems plausible because their correlation length, number density (extrapolated to lower redshift) and their masses are in fair agreement with those observed.

We pay particular attention to the high end of the $\lya$ optical depth distribution of a quasar spectrum. Based on the flux distribution (Miralda-Escud\'e et al 1996), a simple yet seemingly cosmological model -differentiating statistic, $\Delta_{\tau_0}$ -- the cumulative probability of a quasar spectrum with $\lya$ optical depth greater than a high value $\tau_0$ -- is emphasized. It is shown that two different models -- the cold dark matter model with a cosmological constant and the mixed hot and cold dark matter model, both normalized to COBE and local galaxy cluster abundance -- yield quite different values of $\Delta_{\tau_0}$: 0.13 of the former versus 0.058 of the latter for $\tau_0=3.0$ at $z=3$. Moreover, it is argued that $\Delta_{\tau_0}$ may be fairly robust to compute theoretically because it does not seem to depend sensitively on small variations of simulations parameters such as radiation field, cooling, feedback process, radiative transfer, resolution and simulation volume within the plausible ranges of the concerned quantities. Furthermore, it is illustrated that $\Delta_{\tau_0}$ can be obtained sufficiently accurately from currently available observed quasar spectra for $\tau_0\sim 3.0-4.0$, when observational noise is properly taken into account. We anticipate that analyses of observations of quasar $\lya$ absorption spectra over a range of redshift may be able to constrain the redshift evolution of the amplitude of the density fluctuations on small-to-intermediate scales, therefore providing an independent constraint on $\Omega_0$, $\Omega_{0,HDM}$ and $\Lambda_0$.

We have measured the distribution function of the flux decrement caused by Lyman alpha forest absorption in a new sample of high resolution QSO spectra. The observations are compared to the results from two simulations of the Lya forest: an Eulerian Lambda-CDM model, and an SPH standard CDM model. Good agreement between the shapes of simulated and observed distributions is achieved by globally scaling the optical depth to match the mean flux decrements. This procedure amounts to a measurement of the parameter Omega_b^2 h^3 / Gamma (where Omega_b is the baryonic matter density and Gamma is the HI ionization rate). Estimating a lower limit Gamma > 7 10^{-13} s^{-1} from the abundance of known QSOs, we derive a lower limit to the baryon density, Omega_b h^2>0.021(0.017) for the Lambda-CDM (SCDM) model. In both cases the large values are inconsistent with some recent D/H determinations (Rugers & Hogan 1996a,b), favoring a low deuterium abundance as reported by Tytler, Fan & Burles (1996). Adopting a fixed Omega_b, we can determine the evolution of the ionizing radiation field. Our models predict the intensity to be approximately constant with redshift, consistent with the assumption that the ionizing background is produced by known quasars for z < 3. However, additional sources of ionizing photons are required at higher redshift.

The dynamics of the Local Group and its environment provide a unique challenge to cosmological models. The velocity field within 5h-1 Mpc of the Local Group (LG) is extremely ``cold''. The deviation from a pure Hubble flow, characterized by the observed radial peculiar velocity dispersion, is measured to be about 60km/s. We compare the local velocity field with similarly defined regions extracted from N-body simulations of Universes dominated by cold dark matter (CDM). This test is able to strongly discriminate between models that have different mean mass densities. We find that neither the Omega=1 (SCDM) nor Omega=0.3 (OCDM) cold dark matter models can produce a single candidate Local Group that is embedded in a region with such small peculiar velocities. For these models, we measure velocity dispersions between 500-700km/s and 150-300km/s respectively, more than twice the observed value. Although both CDM models fail to produce environments similar to those of our Local Group on a scale of a few Mpc, they can give rise to many binary systems that have similar orbital properties as the Milky Way--Andromeda system. The local, gravitationally induced bias of halos in the CDM ``Local Group'' environment, if defined within a sphere of 10 Mpc around each Local Group is about 1.5, independent of Omega. No biasing scheme could reconcile the measured velocity dispersions around Local Groups with the observed one. Identification of binary systems using a halo finder (named Skid (http://www-hpcc.astro.washington.edu/tools/DENMAX for a public version)) based on local density maxima instead of a simple linking algorithm, gives a much more complete sample. We show that a standard ``friend of friends'' algorithm would miss 40% of the LG candidates present in the simulations.

Gravitational lensing directly measures mass density fluctuations along the lines of sight to very distant objects. No assumptions need to be made concerning bias, the ratio of fluctuations in galaxy density to mass density. Hence, lensing is a very useful tool to study the universe at low to moderate redshifts. We describe in detail a new method to trace light rays from redshift zero through a three dimensional mass distribution to high redshift. As an example, this method is applied here to a standard cold dark matter universe. We obtain a variety of results, some of them statistical in nature, others from rather detailed case studies of individual ``lines of sight". Among the former are the frequency of multiply imaged quasars, the distribution of separation of the multiple quasars, and the redshift distribution of lenses: all that as a function of quasar redshift. We find effects from very weak lensing to highly magnified multiple images of high redshift objects, which, for extended background sources, (i.e. galaxies), range from slight deformations of the shapes through tangentially aligned arclets up to giant luminous arcs. Different cosmological models differ, increasingly with redshift, in their predictions for the mass (thus gravitational potential) distributions. Our ultimate goal is to apply this method to a number of cosmogonic models and to eliminate some models whose gravitational lensing properties are inconsistent with those observed.

It is demonstrated that strong magnetic fields are produced from a zero initial magnetic field during the pregalactic era, when galaxies are first forming. Their development proceeds in three phases. In the first phase, weak magnetic fields are created by the Biermann battery mechanism, acting in shocked parts of the intergalactic medium where caustics form and intersect. In the second phase, these weak magnetic fields are amplified to strong magnetic fields by the Kolmogoroff turbulence endemic to gravitational structure formation of galaxies. During this second phase, the magnetic fields reach saturation with the turbulent power, but they are coherent only on the scale of the smallest eddy. In the third phase, the magnetic field strength increases to equipartition with the turbulent energy, and the coherence length of the magnetic fields increases to the scale of the largest turbulent eddy, comparable to the scale of the entire galaxy. The resulting magnetic field represents a galactic magnetic field of primordial origin. No further dynamo action is necessary, after the galaxy forms, to explain the origin of magnetic fields. However, the magnetic field may be altered by dynamo action once the galaxy and the galactic disk have formed. It is first shown by direct numerical simulations, that thermoelectric currentsassociated with the Biermann battery, build the field up from zero to $10^{-21}$ G in the regions about to collapse into galaxies, by $z\sim3$. For weak fields, in the absence of dissipation, the cyclotron frequency ${\bf \omega_{cyc}}=e{\bf B } /m_H c $ and $ {\bf \omega}/(1+ \chi )$, where ${\bf \omega = \nabla \times v }$ is the vorticity and $\chi$ is the degree of ionization, satisfy the same equations, and initial conditions ${\bf \omega_{cyc}=\omega}=0$, so that, ${\bf \omega_{cyc}}({\bf

Weak gravitational lensing by large-scale structure affects the determination of the cosmological deceleration parameter $q_0$. We find that the lensing induced dispersions on truly standard candles are $0.04$ and $0.02$ mag at redshift $z=1$ and $z=0.5$, respectively, in a COBE-normalized cold dark matter universe with $\Omega_0=0.40$, $\Lambda_0=0.6$, $H=65$km/s/Mpc and $\sigma_8=0.79$. It is shown that one would observe $q_0=-0.44^{+0.17}_{-0.05}$ and $q_0=-0.45^{+0.10}_{-0.03}$ (the errorbars are $2\sigma$ limits) with standard candles with zero intrinsic dispersion at redshift $z=1$ and $z=0.5$, respectively, compared to the truth of $q_0=-0.40$ in this case, i.e., a 10\% error in $q_0$ will be made. A standard COBE normalized $\Omega_0=1$ CDM model would produce three times as much variance and a mixed (hot and cold) dark matter model would lead to an intermediate result. One unique signature of this dispersion effect is its non Gaussianity. Although the lensing induced dispersion at lower redshift is still significantly smaller than the currently best observed (total) dispersion of $0.12$ mag in a sample of type Ia supernovae, selected with the multicolor light curve shape method, it becomes significant at higher redshift. We show that there is an optimal redshift, in the range $z\sim 0.5-2.0$ depending on the amplitude of the intrinsic dispersion of the standard candles, at which $q_0$ can be most accurately determined.

The nonlinear clustering of dark matter particles in an expanding universe is usually studied by N-body simulations. One can gain some insight into this complex problem if simple relations between physical quantities in the linear and nonlinear regimes can be extracted from the results of N-body simulations. Hamilton et al. (1991) and Nityananda and Padmanabhan (1994) have made an attempt in this direction by relating the mean relative pair velocities to the mean correlation function in a useful manner. We investigate this relation and other closely related issues in detail for the case of six different power spectra: power laws with spectral indexes $n=-2,-1$, cold dark matter (CDM), and hot dark matter models with density parameter $\Omega=1$; CDM including a cosmological constant ($\Lambda$) with $\Omega_{CDM}=0.4$, $\Omega_{\Lambda}=0.6$; and $n=-1$ model with $\Omega=0.1$. We find that: (i) Power law spectra lead to self-similar evolution in an $\Omega=1$ universe. (ii) Stable clustering does not hold in an $\Omega=1$ universe to the extent our simulations can ascertain. (iii) Stable clustering is a better approximation in the case of $\Omega<1$ universe in which structure formation freezes out at some low redshift. (iv) The relation between dimensionless pair velocity and the mean correlation function, $\bar\xi$, is only approximately independent of the shape of the power spectrum. At the nonlinear end, the asymptotic value of the dimensionless pair velocity decreases with increasing small scale power, because the stable clustering assumption is not universally true. (v) The relation between the evolved $\bar\xi$ and the linear regime $\bar\xi$ is also not universal but shows a weak spectrum dependence. We present simple theoretical arguments for these conclusions.

We compute, including a current state-of-the-art treatment of hydrodynamical processes, heating and cooling, a variety of cosmological models into the extreme nonlinear phase to enable comparisons with observations. First, we note the common, model independent results. All have a mean ($z=0$) temperature of $10^{4.5}-10^{5.5}$\Kel, set essentially by photoheating processes. Most gas is in one of two components: either at the photoheating floor $10^{4.5}$\Kel and primarily in low density regions or else shock heated to $10^5-10^6$\Kel and in regions of moderate overdensity (in caustics and near groups and clusters). It presents a major observational challenge to observationally detect this second, abundant component as it is neither an efficient radiator nor absorber. About 2% to 10% of the baryons cool and collapse into galaxies forming on caustics and migrating to clusters. About 1%-2% of baryons are in the very hot X-ray emitting gas near cluster cores, in good agreement with observations. These correspondances between the simulations and the real world imply that there is some significant truth to the underlying standard scenarios for the growth of structure. The differences among model predictions may help us find the path to the correct model. For COBE normalized models the most relevant differences concern epoch of structure formation. In the open variants having $\Omega=0.3$, with or without a cosmological constant, structure formation on galactic scales is well advanced at redshift z=5, and reionization occurs early. But if observations require models for which most galaxy formation occurs more recently than $z=2$, then the flat $\Omega=1$ models are to be preferred. The velocity dispersion on the $1h^{-1}$Mpc scale also provides a strong discriminant with, as expected, the $\Omega=1$ models giving a much higher

We model the $\lya$ forest in a spatially flat, CDM model with $\Omega=0.4$, with an Eulerian hydrodynamic simulation, and find that the intergalactic, photoionized gas collapses into sheet-like and filamentary structures which produce absorption lines like those observed in the $\lya$ forest. A typical filament is $\sim 1 h^{-1} \mpc$ long with thickness $\sim 50-100 h^{-1} \kpc$ (in proper units), and baryonic mass $\sim 10^{10} h^{-1} \msun$. The gas temperature increases with time as structures with larger velocities collapse gravitationally. The predicted distributions of column densities and b-parameters of the absorption lines agree with the observed ones, and their evolution is consistent with our model if the ionizing background has an approximately constant intensity between $z=2$ and $z=4$. We also predict the distribution of transmitted flux and its correlation along a spectrum and on parallel spectra, and the \heii flux decrement as a function of redshift. We predict a correlation length of $\sim 80 h^{-1}\kpc$ perpendicular to the line of sight for features in the Lyman alpha forest. In order to reproduce the observed number of lines and average flux transmission, the baryon content of the clouds may need to be significantly higher than in previous models because of the low densities we predict. If the background intensity $\jhi$ is at least that predicted from the observed quasars, $\Omega_b$ needs to be $\sim 0.025 h^{-2}$, higher than expected by light element nucleosynthesis; the model also predicts that most of the baryons at $z>2$ are in $\lya$ clouds, and that the rate at which the baryons move to more overdense regions is slow.

The observed baryon fraction and velocity--temperature relation in clusters of galaxies are compared with hydrodynamic simulations in two cosmological models : standard (Omega = 1) and a low-density flat (Omega=0.45 and \lambda=0.55) CDM models, normalized to the COBE background fluctuations. The observed properties of clusters include the velocity dispersion versus temperature relation, the gas mass versus total mass relation, and the gas mass fraction versus velocity dispersion relation. We find that, while both cosmological models reproduce well the shape of these observed functions, only low-density CDM can reproduce the observed amplitudes. The cluster gas mass fraction reflects approximately the baryon fraction in the models, with a slight anti-bias. Therefore, due to the low baryon density given by nucleosynthesis, Omega = 1 models produce too few baryons in clusters compared with observations. Scaling our results as a function of Omega, we find that a low-density CDM model, with Omega approximately 0.3 - 0.4, best reproduces the observed mean baryon fraction in clusters. The observed beta parameter of clusters, beta = sigma^2/(k T/mu m_p) = 0.94 \pm 0.08 discriminates less well between the models; it is consistent with that produced by low-density CDM (1.10 \pm 0.22), while it is slightly larger than expected but still consistent with Omega = 1 (0.70 \pm 0.14).

The topology of large scale structure is studied as a function of galaxy type using the genus statistic. In hydrodynamical cosmological CDM simulations, galaxies form on caustic surfaces (Zeldovich pancakes) then slowly drain onto filaments and clusters. The earliest forming galaxies in the simulations (defined as ``ellipticals") are thus seen at the present epoch preferentially in clusters (tending toward a meatball topology), while the latest forming galaxies (defined as ``spirals") are seen currently in a spongelike topology. The topology is measured by the genus (= number of ``donut" holes - number of isolated regions) of the smoothed density-contour surfaces. The measured genus curve for all galaxies as a function of density obeys approximately the theoretical curve expected for random-phase initial conditions, but the early forming elliptical galaxies show a shift toward a meatball topology relative to the late forming spirals. Simulations using standard biasing schemes fail to show such an effect. Large observational samples separated by galaxy type could be used to test for this effect.

We present a new treatment of two popular models for the growth of structure, examining the X-ray emission from hot gas with allowance for spectral line emission from various atomic species, primarily ``metals". The X-ray emission from the bright cluster sources is not significantly changed from prior work and, as noted earlier, shows the CDM$+\Lambda$ model (LCDM) to be consistent but the standard, COBE normalized model (SCDM) to be inconsistent with existing observations --- after allowance for still the considerable numerical modelling uncertainties. But we find one important new result. Radiation in the softer band 0.5-1.0keV is predominantly emitted by gas far from cluster centers (hence ``background"). This background emission dominates over the cluster emission below 1keV and observations of it should show clear spectral signatures indicating its origin. In particular the ``iron blend" should be seen prominantly in this spectral bin from cosmic background hot gas at high galactic latitudes and should show shadowing against the SMC indicating its extragalactic origin. Certain OVII lines also provide a signature of this gas which emits a spectrum characteristic of $10^{6.6\pm 0.6}$K gas. Recent ASCA observations of the X-ray background tentatively indicate the presence of component with exactly the spectral features we predict here.

The distributions of peculiar velocities of rich clusters and of groups of galaxies are investigated for different cosmological models and are compared with observations. Four cosmological models are studied: standard ($\Omega=1$) CDM, low-density ($\Omega=0.3$) CDM, HDM ($\Omega=1$), and PBI ($\Omega=0.3$), all COBE normalized. We find that rich clusters of galaxies exhibit a Maxwellian distribution of peculiar velocities in all models,as expected from a Gaussian initial density field. The cluster 3-D velocity distribution typically peaks at $v \sim 600$ km s$^{-1}$, and extends to high cluster velocities of $v \sim 2000$ km s$^{-1}$. The low-density CDM model exhibits somewhat lower velocities: it peaks at $\sim 400$ km s$^{-1}$, and extends to $\sim 1200$ km s$^{-1}$. Approximately 10\% ($\sim$1\% for low- density CDM) of all model rich clusters move with high peculiar velocities of $v \ge 10^3$ km s$^{-1}$. The highest velocity clusters frequently originate in dense superclusters. The model velocity distributions of groups and clusters of galaxies are compared with observations. The data exhibit a larger high-velocity tail, to $v_r \geq 2000$ km s$^{-1}$, than seen in the model simulations (except HDM). Due to the large observational uncertainties, however, the data are consistent at a $\sim 1$ to $3\sigma$ level with the model predictions and with a Gaussian initial density field. Accurate observations of cluster peculiar velocities, especially at the high-velocity tail, should provide powerful constraints on the cosmological models.

Gravitational lensing provides a strict test of cosmogonic models. Detailed numerical propagation of light rays through a universe having a distribution of inhomogeneities derived from the standard CDM (cold dark matter) scenario with the aid of massive, fully nonlinear computer simulations is used to test the model. It predicts that we should have seen far more widely split quasar images than have been found. These and other inconsistencies rule out the COBE normalized CDM model with $\Omega=1$ and $H_0=50$km/s/Mpc; but variants of this model might be constructed, which could pass the stringent tests provided by strong gravitational lensing.

We determine the velocity correlation function, pairwise peculiar velocity difference, and root-mean-square pairwise peculiar velocity dispersion of rich clusters of galaxies, as a function of pair separation, for three cosmological models: Omega=1 and Omega=0.3 CDM, and Omega=0.3 PBI models (all flat and COBE-normalized). We find that close cluster pairs, with separation r<10Mpc/h, exhibit strong attractive peculiar velocities in all models; the cluster pairwise velocities depend sensitively on the model. The mean pairwise attractive velocity of clusters on 5Mpc/h scale ranges from 1700 km/s for Omega=1 CDM, to 1000 km/s for PBI, to 700 km/s for Omega=0.3 CDM. The small-scale pairwise velocities depend also on cluster mass: richer, more massive clusters exhibit stronger attractive velocities than less massive clusters. On large scales, from 20 to 200Mpc/h, the cluster peculiar velocities are increasingly dominated by bulk and random motions; they are independent of cluster mass. The cluster velocity correlation function, which reflects the bulk motion minus the relative motion of pairs, is negative on small scales for Omega=1 and Omega=0.3 CDM, and positive for PBI; this indicates stronger pairwise motion than bulk motion on small scales for CDM, and relatively larger bulk motions for PBI. The cluster velocity correlation function is positive on very large scales, from 10 to 200Mpc/h, for all models. These positive correlations, which decrease monotonically with scale, indicate significant bulk motions of clusters up to 200Mpc/h. The strong dependence of the cluster velocity functions on models, especially at small separations, makes them useful tools in constraining cosmological models when compared with observations.

If gravitational clustering is a hierarchical process, the present large-scale structure of the galaxy distribution implies that structures on smaller scales must have formed at high redshift. We simulate the formation of small-scale structure (average cell mass: $\Delta \bar m_b=10^{4.2}$M$_\odot$) and the evolution of photoionized gas, in the specific case of a CDM model with a cosmological constant. The photoionized gas has a natural minimal scale of collapse, the Jeans scale ($m_{b,J}\simeq 10^{9}$M$_\odot$). We find that low column density ($\nhi \le 10^{14}\cm^{-2}$) lines originate in regions resembling Zel'dovich pancakes, where gas with overdensities in the range $3 - 30$ is enclosed by two shocks but is typically re-expanding at approximately the Hubble velocity. However, higher column density ($\nhi \ge 10^{15}\cm^{-2}$) lines stem from more overdense regions where the shocked gas is cooling. We show that this model can probably account for the observed number of lines, their distribution in column density and b-parameters, as well as the cloud physical sizes as observed in gravitationally lensed quasars. We find a redshift evolution that is

Lauer \& Postman (LP) observe that all Abell clusters with redshifts less than 15,000\kms\ appear to be participating in a bulk flow of 689 km s$^{-1}$ with respect to the Cosmic Microwave Background. We find this result difficult to reconcile with all popular models for large-scale structure formation that assume Gaussian initial conditions. This conclusion is based on Monte-Carlo realizations of the LP data, drawn from large Particle-Mesh $N$-body simulations. We have taken special care to treat properly the longest-wavelength components of the power spectra. Bulk flows with amplitude as large as that reported by LP are not uncommon in the Monte-Carlo datasets. However, the $\chi^2$ of the observed bulk flow, taking into account the anisotropy of the error ellipsoid, is much more difficult to match in the simulations. The models examined are ruled out at confidence levels between 94\% and 98\%. Any model that has {\it intrinsic} flows of less than 480\kms\ on the scales probed by LP scales can be ruled out at a similar level.

What is the role of clusters of galaxies in probing the large-scale velocity field of the universe? We investigate the distribution of peculiar velocities of clusters of galaxies in the popular low-density ($\Omega=0.3$) flat Cold-Dark-Matter (CDM) cosmological model, which best fits many large-scale structure observations. An $\Omega=1$ CDM model is also studied for comparison. We find that clusters of galaxies are efficient tracers of the large-scale velocity field. The clusters exhibit a Maxwellian distribution of peculiar velocities, as expected from Gaussian initial density fluctuations. The cluster 3-D velocity distribution for the $\Omega=0.3$ model peaks at $v \sim 400$ km s$^{-1}$, and extends to high velocities of $v \sim 1200$ km s$^{-1}$. The rms peculiar velocity of the clusters is $440$ km s$^{-1}$. Approximately 10\% of all model clusters move with high peculiar velocities of $v \ge 700$ km s$^{-1}$. The observed velocity distribution of clusters of galaxies is compared with the predictions from cosmological models. The observed data exhibit a larger velocity tail than seen in the model simulations; however, due to the large observational uncertainties, the data are consistent at a $\sim 3\sigma$ level with the model predictions, and with a Gaussian initial density field. The large peculiar velocities reported for some clusters of galaxies ($v \geq 3000$ km s$^{-1}$) are likely to be overestimated, if the current model is viable.

The mass and velocity distributions in the outskirts ($0.5-3.0h^{-1}$Mpc) of clusters of galaxies are examined for a suite of cosmogonic models utilizing large-scale Particle-Mesh (PM) simulations ($500^3$ cells, $250^3$ particles and box size of $100h^{-1}$Mpc, giving a nominal resolution of $0.2h^{-1}$Mpc with the true resolution about $0.5h^{-1}$Mpc). Through a series of model computations, designed to isolate the different effects, we find that both $\Omega_0$ and $P_k$ ($\lambda\le 16h^{-1}$Mpc) are important to the mass and velocity distributions in clusters of galaxies. More details in the full abstract.

The distributions of peculiar velocities of rich clusters and of groups of galaxies are investigated for different cosmological models and are compared with observations. Four cosmological models are studied: standard ($\Omega=1$) CDM, low-density CDM, HDM ($\Omega=1$), and PBI. We find that rich clusters of galaxies exhibit a Maxwellian distribution of peculiar velocities in all models, as expected from a Gaussian initial density fluctuation field. The cluster 3-D velocity distribution is generally similar in the models: it peaks at $v \sim 500$ km s$^{-1}$, and extends to high cluster velocities of $v \sim 1500$ km s$^{-1}$. Approximately 10\% of all model rich clusters move with high peculiar velocities of $v \ge 10^3$ km s$^{-1}$. The highest velocity clusters frequently originate in dense superclusters. The group velocity distribution is, in general, similar to the velocity distribution of the rich clusters. In all but the low-density CDM model, the mass exhibits a longer tail of high velocities than do the clusters. This high-velocity tail originates mostly from the high velocities that exist within rich clusters. The model velocity distributions of groups and clusters of galaxies are compared with observations. The data are generally consistent with the models, but exhibit a somewhat larger high-velocity tail, to $v_r \sim 3000$ km s$^{-1}$. While this high-velocity tail is similar to the HDM model predictions, the data are consistent with the other models studied, including the low-density CDM model, which best fits most other large-scale structure observations. The observed velocity

We present a detailed comparison of the simulation results of various cosmological hydrodynamic codes. Starting with identical initial conditions based on the Cold Dark Matter scenario for the growth of structure, we integrate from redshift $z=20$ to $z=0$ to determine the physical state within a representative volume of size $L^3$ where $L=64 h^{-1} {\rm Mpc}$. Five independent codes are compared: three of them Eulerian mesh based and two variants of the Smooth Particle Hydrodynamics "SPH" Lagrangian approach. The Eulerian codes were run at $N^3=(32^3,~64^3,~128^3,~{\rm and},~256^3)$ cells, the SPH codes at $N^3= 32^3$ and $64^3$ particles. Results were then rebinned to a $16^3$ grid with the expectation that the rebinned data should converge, by all techniques, to a common and correct result as $N \rightarrow \infty$. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainties in even primitive quantities such as $\langle T \rangle$, $\langle \rho^2\rangle^{1/2}$ persists at the 3\%-17\% level after completion of very large simulations. The two SPH codes and the two shock capturing Eulerian codes achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high density, high temperature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by $\rho^2$) is only 4.5\%. Overall the comparison allows us to better estimate errors, it points to ways of improving this current generation of hydrodynamic codes and of suiting their use to problems which exploit their individually best features.

We exmaine the distribution of hot gas in a standard CDM model of the universe using high resolution hydrodynamic simulations. Adopting standard parameters determined from COBE and light element nucleosynthesis, $\sigma_8=1.05$, $\Omega_b=0.06$ and assuming $h=0.5$, we find the X-ray emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution and estimated sizes as well as the evolution of these quantities with redshift. This standard CDM model, normalized to COBE, produces approximately 5 times too much emission from clusters having $L_x>10^{43}$erg/s, a not unexpected result. If all other parameters were unchanged, we would expect adequate agreement for $\sigma_8=0.6$. This provides a new and independent argument for lower small scale power than standard CDM at the $8h^{-1}$Mpc scale. The background radiation field at 1keV due to clusters in this model is approximately $1/3$ of the observed background which, after correction for numerical effects, again indicates approximately 5 times too much emission and the appropriateness of $\sigma_8=0.6$. If we had used the observed ratio of gas to total mass in clusters, rather than basing the mean density on light element nucleosynthesis, then the computed luminosity of each cluster would have increased still further, by a factor of approximately ten. Examining the ratio of gas to total mass in the clusters (which we find to be anti-biased by a factor of approximately 0.6), normalized to $\Omega_b h^2=0.015$, and comparing to observations, we conclude, in agreement with S. White, that the cluster observations argue for an open universe.

A new, three-dimensional, shock capturing, TVD hydrodynamic code is utilized to determine the distribution of hot gas in a CDM$+\Lambda$ model universe. This open model succeeds in matching local observations of clusters in contrast to the standard $\Omega=1$, CDM model, which fails. It predicts an order of magnitude decline in the number density of bright ($h\nu = 2-10$keV) clusters from $z=0$ to $z=2$ in contrast to a slight increase in the number density for standard $\Omega=1$, CDM model. This COBE-normalized CDM$+\Lambda$ model produces approximately the same number of X-ray clusters having $L_x>10^{43}$erg/s as observed. The number density of bright clusters increases to $z\sim 0.2-0.5$ and declines, but the luminosity per typical cluster decreases monotonically with redshift, with the result that the number density of bright clusters shows a broad peak near $z=0.5$, and then a rapid decline as $z\rightarrow 3$. The most interesting point which we find is that the temperatures of clusters in this model freeze out at later times ($z\le 0.3$), while previously we found in the CDM model that there was a steep increase during the same interval of redshift. The trend should be detectable even with a relatively ``soft" X-ray instrument such as ROSAT, providing a powerful discriminant between $\Omega=1$ and $\Omega<1$ models. Examining the ratio of gas-to-total mass in the clusters, we find a slight antibias ($b=0.9$ or $({\Omega_{gas}\over\Omega_{tot}})_{cl} =0.083\pm 0.007$), which is consistent with observations [$({\Omega_{gas}\over\Omega_{tot}})_{obs}=0.097\pm 0.019$ for the Coma cluster for the given value of $h$, \cf, White 1991].

We compute the evolution of spatially flat, mixed cold and hot dark matter (``MDM") models containing both baryonic matter and two kinds of dark matter. The mean final Zeldovich-Sunyaev $y$ parameter is estimated to be $\bar y=(5.4\pm 2.7)\times 10^{-7}$, with a rms fluctuation of approximately $\bar{\delta y}=(6.0\pm 3.0)\times 10^{-7}$ on arc minute scales. The rate of galaxy formation peaks at an even later epoch ($z\sim 0.3$) than in the standard ($\Omega=1$, $\sigma_8=0.67$) CDM model ($z\sim 0.5$) and, at a redshift of z=4 is nearly a factor of a hundred lower than for the CDM model with the same value of $\sigma_8$. The model yields a cluster mass function which is about a factor of 4 higher than observed but a cluster-cluster correlation length lower by a factor of 2 than what is observed but both are closer to observations than in the COBE normalized CDM model. The one dimensional pairwise velocity dispersion is $605\pm 8$km/s at $1h^{-1}$ separation, lower than that of the CDM model normalized to COBE, but still significantly higher than observations (Davis \& Peebles 1983). A plausible velocity bias $b_v=0.8\pm 0.1$ on this scale will reduce but not remove the discrepancy. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of $\sigma_8$, but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter significantly improves the ability of COBE normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same $\sigma_8$ and is still probably unsatisfactory with regard to several critical tests.

We show that the spatial correlation function of a flux-limited sample of X-ray selected clusters of galaxies will exhibit a correlation scale that is smaller than the correlation scale of a volume-limited, richness-limited sample of comparable apparent spatial density. The flux-limited sample contains clusters of different richnesses at different distances: poor groups are found nearby and rich clusters at greater distances. Since the cluster correlation strength is known to increase with richness, the flux-limited sample averages over the correlations of poor and rich clusters. On the other hand, a volume-limited, richness-limited sample has a minimum richness threshold, and a constant mixture of richnesses with redshift. Using the observed correlation scale for rich ($R\ge 1$) clusters, $r_o (R\ge 1) = 21\pm 2 h^{-1}$Mpc that was determined from previous volume-limited studies, we derive for the ROSAT flux-limited X-ray cluster sample $r_o({\rm flux}{\rm -}{\rm limited}) \approx 14h^{-1}$Mpc, in agreement with the recently observed value of $13.7\pm 2.3 h^{-1}$Mpc.

The large-scale bulk flow implied by the measurements of distances to brightest cluster members within 15,000 km/s by Lauer and Postman has the potential to put strong constraints on, or to even rule out, various models of large-scale structure. Using PM simulations 400$h^{-1}$ Mpc on a side of five cosmological models (Standard CDM, Tilted CDM, HDM, $\Omega = 0.3$ CDM, and PBI), we present detailed Monte-Carlo simulations of the peculiar velocity data, which we use to assess the probability that any given model can explain the data. These models are ruled out at confidence levels ranging from 94\% (HDM) to 97\% (PBI). However, these results may be affected by the finite size of the simulation volume. If the observed bulk flow remains unchanged when calculated using a refined distance indicator with smaller scatter, all these models will be ruled out at high confidence levels.