The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is over-simplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds (PhEW), that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a sub-grid prescription to model galactic winds more robustly both physically and numerically.

We interpret simulations of secularly-evolving disc galaxies through orbit morphology. Using a new algorithm that measures the volume of orbits in real space using a tessellation, we rapidly isolate commensurate (resonant) orbits. We identify phase-space regions occupied by different orbital families. Compared to spectral methods, the tessellation algorithm can identify resonant orbits within a few dynamical periods, crucial for understanding an evolving galaxy model. The flexible methodology accepts arbitrary potentials, enabling detailed descriptions of the orbital families. We apply the machinery to four different potential models, including two barred models, and fully characterise the orbital membership. We identify key differences in the content of orbit families, emphasizing the presence of orbit families indicative of the bar evolutionary state and the shape of the dark matter halo. We use the characterisation of orbits to investigate the shortcomings of analytic and self-consistent studies, comparing our findings to the evolutionary epochs in self-consistent barred galaxy simulations. Using insight from our orbit analysis, we present a new observational metric that uses spatial and kinematic information from integral field spectrometers that may reveal signatures of commensurabilities and allow for a differentiation between models.

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$\alpha$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $\eta$ with velocity dispersion $\sigma$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $\eta \sim \sigma^5$ in low mass haloes with $\sigma < 106\ \mathrm{km\ s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc/h}$, $2\times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z\sim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique.

We track the angular momentum transfer in n-body simulations of barred galaxies by measuring torques to understand the dynamical mechanisms responsible for the evolution of the bar-disc-dark matter halo system. We find evidence for three distinct phases of barred galaxy evolution: assembly, secular growth, and steady-state equilibrium. Using a decomposition of the disc into orbital families, we track bar mass and angular momentum through time and correlate the quantities with the phases of evolution. We follow the angular momentum transfer between particles and identify the dominant torque channels. We find that the halo model mediates the assembly and growth of the bar for a high central density halo, and the outer disc mediates the assembly and growth of the bar in a low central density halo model. Both galaxies exhibit a steady-state equilibrium phase where the bar is neither lengthening nor slowing. The steady-state equilibrium results from the balance of torque between particles that are gaining and losing angular momentum. We propose observational metrics for barred galaxies that can be used to help determine the evolutionary phase of a barred galaxy, and discuss the implications of the phases for galaxy evolution as a whole.

We survey the highly ionized circumgalactic media (CGM) of 29 blindly selected galaxies at 0.49 < z_(gal) < 1.44 based on high-S/N ultraviolet spectra of z > 1 QSOs and the galaxy database from the COS Absorption Survey of Baryon Harbors (CASBaH). We detect the Ne VIII doublet in nine of the galaxies, and for gas with N(Ne VIII) > 10^13.3 cm^-2 (> 10^13.5 cm^-2), we derive a Ne VIII covering fraction f_c = 75 +15/-25% (44 +22/-20%) within impact parameter (rho) < 200 kpc of M_* = 10^(9.5-11.5) Msol galaxies and f_c = 70 +16/-22% (f_c = 42 +20/-17%) within rho < 1.5 virial radii. We estimate the mass in Ne VIII-traced gas to be M_gas(Ne VIII) > 10^9.5 Msol (Z/Zsol)^-1, or 6-20% of the expected baryonic mass if the Ne VIII absorbers have solar metallicity. Ionizing Ne VII to Ne VIII requires 207 eV, and photons with this energy are scarce in the CGM. However, for the median halo mass and redshift of our sample, the virial temperature is close to the peak temperature for the Ne VIII ion, and the Ne VIII-bearing gas is plausibly collisionally ionized near this temperature. Moreover, we find that photoionized Ne VIII requires cool and low-density clouds that would be highly underpressured (by approximately two orders of magnitude) relative to the putative, ambient virialized medium, complicating scenarios where such clouds could survive. Thus, more complex (e.g., non-equilibrium) models may be required; this first statistical sample of Ne VIII absorber/galaxy systems will provide stringent constraints for future CGM studies.

We study the mechanisms and evolutionary phases of bar formation in n-body simulations of a stellar disc and dark matter halo system using harmonic basis function expansion analysis to characterize the dynamical mechanisms in bar evolution. We correlate orbit families with phases of bar evolution by using empirical orthogonal functions that act as a spatial filter and form the gravitational potential basis. In both models we find evidence for three phases in evolution with unique harmonic signatures. We recover known analytic results, such as bar slowdown owing to angular momentum transfer. We also find new dynamical mechanisms for bar evolution: a steady-state equilibrium configuration and harmonic interaction resulting in harmonic mode locking, both of which may be observable. Additionally, we find that ellipse fitting may severely overestimate measurements of bar length by a factor of two relative to the measurements based on orbits that comprise the true backbone supporting the bar feature. The bias will lead to overestimates of both bar mass and bar pattern speed, affecting inferences about the evolution of bars in the real universe, such as the fraction of bars with fast pattern speeds. We propose a direct observational technique to compute the radial extent of trapped orbits and determine a dynamical length for the bar.

We explore the application of Bayesian image analysis to infer the properties of an SDSS early-type galaxy sample including AGN. We use GALPHAT (Yoon et al. 2010) with a Bayes-factor model comparison to photometrically infer an AGN population and verify this using spectroscopic signatures. Our combined posterior sample for the SDSS sample reveals distinct low and high concentration modes after the point-source flux is modeled. This suggests that ETG parameters are intrinsically bimodal. The bimodal signature was weak when analyzed by GALFIT (Peng et al. 2002, 2010). This led us to create several ensembles of synthetic images to investigate the bias of inferred structural parameters and compare with GALFIT. GALPHAT inferences are less biased, especially for high-concentration profiles: GALPHAT S\'ersic index $n$, $r_{e}$ and MAG deviate from the true values by $6\%$, $7.6\%$ and $-0.03 \,\mathrm{mag}$, respectively, while GALFIT deviates by $15\%$, $22\%$ and $-0.09$\, mag, respectively. In addition, we explore the reliability for the photometric detection of AGN using Bayes factors. For our SDSS sample with $r_{e}\ge 7.92\,$arcsec, we correctly identify central point sources with $\mathrm{Mag_{PS}}-\mathrm{Mag_{Sersic}}\le 5$ for $n\le6$ and $\mathrm{Mag_{PS}}-\mathrm{Mag_{Sersic}}\le 3$ for $n>6$. The magnitude range increases and classification error decreases with increasing resolution, suggesting that this approach will excel for upcoming high-resolution surveys. Future work will extend this to models that test hypotheses of galaxy evolution through the cosmic time.

We investigate the complex interactions between the stellar disc and the dark-matter halo during bar formation and evolution using N-body simulations with fine temporal resolution and optimally chosen spatial resolution. We find that the forming stellar bar traps dark matter in the vicinity of the stellar bar into bar-supporting orbits. We call this feature the shadow bar. The shadow bar modifies both the location and magnitude of the angular momentum transfer between the disc and dark matter halo and adds 10 per cent to the mass of the stellar bar over 4 Gyr. The shadow bar is potentially observable by its density and velocity signature in spheroid stars and by direct dark matter detection experiments. Numerical tests demonstrate that the shadow bar can diminish the rate of angular momentum transport from the bar to the dark matter halo by more than a factor of three over the rate predicted by dynamical friction with an untrapped dark halo, and thus provides a possible physical explanation for the observed prevalence of fast bars in nature.

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

Using Hubble Space Telescope Cosmic Origins Spectrograph observations of 89 QSO sightlines through the Sloan Digital Sky Survey footprint, we study the relationships between C IV absorption systems and the properties of nearby galaxies as well as large-scale environment. To maintain sensitivity to very faint galaxies, we restrict our sample to 0.0015 < z < 0.015, which defines a complete galaxy survey to L > 0.01 L* or stellar mass log M_* > 8 Msun. We report two principal findings. First, for galaxies with impact parameter rho < 1 rvir, C IV detection strongly depends on the luminosity/stellar mass of the nearby galaxy. C IV is preferentially associated with galaxies with log M_* > 9.5 Msun; lower mass galaxies rarely exhibit significant C IV absorption (covering fraction f = 9 +12-6% for 11 galaxies with log M_* < 9.5 Msun). Second, C IV detection within the log M_* > 9.5 Msun population depends on environment. Using a fixed-aperture environmental density metric for galaxies with rho < 160 kpc at z < 0.055, we find that 57+/-12% (8/14) of galaxies in low-density regions (regions with fewer than seven L > 0.15 L* galaxies within 1.5 Mpc) have affiliated C IV absorption; however, none (0/7) of the galaxies in denser regions show C IV. Similarly, the C IV detection rate is lower for galaxies residing in groups with dark-matter halo masses of log Mhalo > 12.5 Msun. In contrast to C IV, H I is pervasive in the CGM without regard to mass or environment. These results indicate that C IV absorbers with log N(C IV) > 13.5 cm^-2 trace the halos of log M_* > 9.5 Msun galaxies but also reflect larger scale environmental conditions.

Over a handful of rotation periods, dynamical processes in barred galaxies induce non-axisymmetric structure in dark matter halos. Using n-body simulations of a Milky Way-like barred galaxy, we identify both a trapped dark-matter component, a shadow bar, and a strong response wake in the dark-matter distribution that affects the predicted dark-matter detection rates for current experiments. The presence of a baryonic disk together with well-known dynamical processes (e.g. spiral structure and bar instabilities) increase the dark matter density in the disk plane. We find that the magnitude of the combined stellar and shadow bar evolution, when isolated from the effect of the axisymmetric gravitational potential of the disk, accounts for >30% of this overall increase in disk-plane density. This is significantly larger that of previously claimed deviations from the standard halo model. The dark-matter density and kinematic wakes driven by the Milky Way bar increase the detectability of dark matter overall, especially for the experiments with higher $v_{min}$. These astrophysical features increase the detection rate by more than a factor of two when compared to the standard halo model and by a factor of ten for experiments with high minimum recoil energy thresholds. These same features increase (decrease) the annual modulation for low (high) minimum recoil energy experiments. We present physical arguments for why these dynamics are generic for barred galaxies such as the Milky Way rather than contingent on a specific galaxy model.

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.

The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback.

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

To investigate the evolution of metal-enriched gas over recent cosmic epochs as well as to characterize the diffuse, ionized, metal-enriched circumgalactic medium (CGM), we have conducted a blind survey for C IV absorption systems in 89 QSO sightlines observed with the Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS). We have identified 42 absorbers at z < 0.16, comprising the largest uniform blind sample size to date in this redshift range. Our measurements indicate an increasing C IV absorber number density per comoving path length (dN/dX = 7.5 +/- 1.1) and modestly increasing mass density relative to the critical density of the Universe (Omega(C IV) = 10.0 +/- 1.5 x 10^-8 ) from z ~ 1.5 to the present epoch, consistent with predictions from cosmological hydrodynamical simulations. Furthermore, the data support a functional form for the column density distribution function that deviates from a single power-law, also consistent with independent theoretical predictions. As the data also probe heavy element ions in addition to C IV at the same redshifts, we identify, measure, and search for correlations between column densities of these species where components appear aligned in velocity. Among these ion-ion correlations, we find evidence for tight correlations between C II and Si II, C II and Si III, and C IV and Si IV, suggesting that these pairs of species arise in similar ionization conditions. However, the evidence for correlations decreases as the difference in ionization potential increases. Finally, when controlling for observational bias, we find only marginal evidence for a correlation (86.8% likelihood) between the Doppler line width b(C IV) and column density N(C IV).

We quantitatively examine the effects of accretion and environment on the HI content of galaxies within a cosmological hydrodynamic simulation that reproduces basic observed trends of HI in galaxies. We show that our model broadly reproduces the observed scatter in HI at a given stellar mass as quantified by the HI mass function in bins of stellar mass, as well as the HI richness versus local galaxy density. This shows that the predicted HI fluctuations and environmental effects are roughly consistent with data with few minor discrepancies. For satellite galaxies in >= 10^12M_* halos, the HI richness distribution is bimodal and drops towards the largest halo masses. The depletion rate of HI once a galaxy enters a more massive halo is more rapid at higher halo mass, in contrast to the specific star formation rate which shows much less variation in the attenuation rate versus halo mass. This suggests that, up to halo mass scales probed here (<= 10^14M_*), star formation is mainly attenuated by starvation, but HI is additionally removed by stripping once a hot gaseous halo is present. In low mass halos, the HI richness of satellites is independent of radius, while in high mass halos they become gas-poor towards the center. We show that the gas fraction of satellite and central galaxies decreases from z=5 to z=0,tracking each other until z~1 after which the satellites' HI content drops much more quickly, particularly for the highest halo masses. Mergers somewhat increase the HI richness and its scatter about the mean relation, tracking the metallicity in a way consistent with arising from inflow fluctuations, while star formation is significantly boosted relative to HI.

We investigate the time evolution of non-parametric morphological quantities and their relationship to major mergers between $4\geq z \geq 2$ in high-resolution cosmological zoom simulations of disk galaxies that implement kinetic wind feedback, $H_2$-based star formation, and minimal ISM pressurisation. We show that the resulting galaxies broadly match basic observed physical properties of $z\sim 2$ objects. We measure the galaxies' concentrations ($C$), asymmetries ($A$), and $Gini$ ($G$) and $M_{20}$ coefficients, and correlate these with major merger events identified from the mass growth history. We find that high values of asymmetry provide the best indicator for identifying major mergers of $>1:4$ mass ratio within our sample, with $Gini$-$M_{20}\,$ merger classification only as effective for face-on systems and much less effective for edge-on or randomly-oriented galaxies. The canonical asymmetry cut of $A\geq0.35$, however, is only able to correctly identify major mergers $\sim 10\%$ of the time, while a higher cut of $A\geq 0.8$ more efficiently picks out mergers at this epoch. We further examine the temporal correlation between morphological statistics and mergers, and show that for randomly-oriented galaxies, half the galaxies with $A\geq0.8$ undergo a merger within $\pm0.2\,{\rm Gyr}$, whereas $Gini$-$M_{20}\,$ identification only identifies about a third correctly. The fraction improves further using $A\geq 1.5$, but about the half the mergers are missed by this stringent cut.

The empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids.

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical 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 $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass.

We combine cosmological hydrodynamic simulations with analytic models to evaluate the role of galaxy-scale gravitational torques on the evolution of massive black holes at the centers of star-forming galaxies. We confirm and extend our earlier results to show that torque-limited growth yields black holes and host galaxies evolving on average along the Mbh-Mbulge relation from early times down to z = 0 and that convergence onto the scaling relation occurs independent of the initial conditions and with no need for mass averaging through mergers or additional self-regulation processes. Smooth accretion dominates the long-term evolution, with black hole mergers with mass ratios >1:5 representing typically a small fraction of the total growth. Winds from the accretion disk are required to eject significant mass to suppress black hole growth, but there is no need for coupling this wind to galactic-scale gas to regulate black holes in a non-linear feedback loop. Torque-limited growth yields a close-to-linear relation for the star formation rate and the black hole accretion rate averaged over galaxy evolution time scales. However, the SFR-AGN connection has significant scatter owing to strong variability of black hole accretion at all resolved time scales. Eddington ratios can be described by a broad lognormal distribution with median value evolving roughly as (1 + z)^1.9, suggesting a main sequence for black hole growth similar to the cosmic evolution of specific SFRs. Our results offer an attractive scenario consistent with available observations in which cosmological gas infall and transport of angular momentum in the galaxy by gravitational instabilities regulate the long-term co-evolution of black holes and star-forming galaxies.

We report new observations of circumgalactic gas from the COS-Dwarfs survey, a systematic investigation of the gaseous halos around 43 low-mass z $\leq$ 0.1 galaxies using background QSOs observed with the Cosmic Origins Spectrograph. From the projected 1D and 2D distribution of C IV absorption, we find that C IV absorption is detected out to ~ 0.5 R$_{vir}$ of the host galaxies. The C IV absorption strength falls off radially as a power law and beyond 0.5 R$_{vir}$, no C IV absorption is detected above our sensitivity limit of ~ 50-100 m$\AA$. We find a tentative correlation between detected C IV absorption strength and star formation, paralleling the strong correlation seen in highly ionized oxygen for L~L* galaxies by the COS-Halos survey. The data imply a large carbon reservoir in the CGM of these galaxies, corresponding to a minimum carbon mass of $\gtrsim$ 1.2$\times 10^6$ $M_\odot$ out to ~ 110 kpc. This mass is comparable to the carbon mass in the ISM and more than the carbon mass currently in stars of these galaxies. The C IV absorption seen around these sub-L* galaxies can account for almost two-thirds of all $W_r$> 100 m$\AA$ C IV absorption detected at low z. Comparing the C IV covering fraction with hydrodynamical simulations, we find that an energy-driven wind model is consistent with the observations whereas a wind model of constant velocity fails to reproduce the CGM or the galaxy properties.

We analyze the physical conditions of the cool, photoionized (T $\sim 10^4$ K) circumgalactic medium (CGM) using the COS-Halos suite of gas column density measurements for 44 gaseous halos within 160 kpc of $L \sim L^*$ galaxies at $z \sim 0.2$. These data are well described by simple photoionization models, with the gas highly ionized (n$_{\rm HII}$/n$_{\rm H} \gtrsim 99\%$) by the extragalactic ultraviolet background (EUVB). Scaling by estimates for the virial radius, R$_{\rm vir}$, we show that the ionization state (tracked by the dimensionless ionization parameter, U) increases with distance from the host galaxy. The ionization parameters imply a decreasing volume density profile n$_{\rm H}$ = (10$^{-4.2 \pm 0.25}$)(R/R$_{\rm vir})^{-0.8\pm0.3}$. Our derived gas volume densities are several orders of magnitude lower than predictions from standard two-phase models with a cool medium in pressure equilibrium with a hot, coronal medium expected in virialized halos at this mass scale. Applying the ionization corrections to the HI column densities, we estimate a lower limit to the cool gas mass M$_{\rm CGM}^{\rm cool} > 6.5 \times 10^{10}$ M$_{\odot}$ for the volume within R $<$ R$_{\rm vir}$. Allowing for an additional warm-hot, OVI-traced phase, the CGM accounts for at least half of the baryons purported to be missing from dark matter halos at the 10$^{12}$ M$_{\odot}$ scale.

We infer mechanisms of galaxy formation for a broad family of semi-analytic models (SAMs) constrained by the K-band luminosity function and HI mass function of local galaxies using tools of Bayesian analysis. Even with a broad search in parameter space the whole model family fails to match to constraining data. In the best fitting models, the star formation and feedback parameters in low-mass haloes are tightly constrained by the two data sets, and the analysis reveals several generic failures of models that similarly apply to other existing SAMs. First, based on the assumption that baryon accretion follows the dark matter accretion, large mass-loading factors are required for haloes with circular velocities lower than 200 km/s, and most of the wind mass must be expelled from the haloes. Second, assuming that the feedback is powered by Type-II supernovae with a Chabrier IMF, the outflow requires more than 25% of the available SN kinetic energy. Finally, the posterior predictive distributions for the star formation history are dramatically inconsistent with observations for masses similar to or smaller than the Milky-Way mass. The inferences suggest that the current model family is still missing some key physical processes that regulate the gas accretion and star formation in galaxies with masses below that of the Milky Way.

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

We present a budget and accounting of metals in and around star-forming galaxies at $z\sim 0$. We combine empirically derived star formation histories with updated supernova and AGB yields and rates to estimate the total mass of metals produced by galaxies with present-day stellar mass of $10^{9.3}$--$10^{11.6} M_{\odot}$. On the accounting side of the ledger, we show that a surprisingly constant 20--25% mass fraction of produced metals remain in galaxies' stars, interstellar gas and interstellar dust, with little dependence of this fraction on the galaxy stellar mass (omitting those metals immediately locked up in remnants). Thus, the bulk of metals are outside of galaxies, produced in the progenitors of today's $L^*$ galaxies. The COS-Halos survey is uniquely able to measure the mass of metals in the circumgalactic medium (to impact parameters of $< 150$ kpc) of low-redshift $\sim L^*$ galaxies. Using these data, we map the distribution of CGM metals as traced by both the highly ionized OVI ion and a suite of low-ionization species; combined with constraints on circumgalactic dust and hotter X-ray emitting gas out to similar impact parameters, we show that $\sim 40$% of metals produced by $M_{\star}\sim 10^{10}M_{\odot}$ galaxies can be easily accounted for out to 150 kpc. With the current data, we cannot rule out a constant mass of metals within this fixed physical radius. This census provides a crucial boundary condition for the eventual fate of metals in galaxy evolution models.

We develop an empirical approach to infer the star formation rate in dark matter halos from the galaxy stellar mass function (SMF) at different redshifts and the local cluster galaxy luminosity function (CGLF), which has a steeper faint end relative to the SMF of local galaxies. As satellites are typically old galaxies which have been accreted earlier, this feature can cast important constraint on the formation of low-mass galaxies at high-redshift. The evolution of the SMFs suggests the star formation in high mass halos ($>10^{12}M_{\odot}h^{-1}$) has to be boosted at high redshift beyond what is expected from a simple scaling of the dynamical time. The faint end of the CGLF implies a characteristic redshift $z_c\approx2$ above which the star formation rate in low mass halos with masses $< 10^{11}M_{\odot}h^{-1}$ must be enhanced relative to that at lower z. This is not directly expected from the standard stellar feedback models. Also, this enhancement leads to some interesting predictions, for instance, a significant old stellar population in present-day dwarf galaxies with $M_* < 10^8 M_{\odot}h^{-2}$ and steep slopes of high redshift stellar mass and star formation rate functions.

We study the high-ionization phase and kinematics of the circumgalactic medium around low-redshift galaxies using a sample of 23 Lyman Limit Systems (LLSs) at 0.08<z<0.93 observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. In Lehner et al. (2013), we recently showed that low-z LLSs have a bimodal metallicity distribution. Here we extend that analysis to search for differences between the high-ion and kinematic properties of the metal-poor and metal-rich branches. We find that metal-rich LLSs tend to show higher O VI columns and broader O VI profiles than metal-poor LLSs. The total H I line width (dv90 statistic) in LLSs is not correlated with metallicity, indicating that the H I kinematics alone cannot be used to distinguish inflow from outflow and gas recycling. Among the 17 LLSs with O VI detections, all but two show evidence of kinematic sub-structure, in the form of O VI-H I centroid offsets, multiple components, or both. Using various scenarios for how the metallicity in the high-ion and low-ion phases of each LLS compare, we constrain the ionized hydrogen column in the O VI phase to lie in the range log N(H II)~17.6-20. The O VI phase of LLSs is a substantial baryon reservoir, with M(high-ion)~10^{8.5-10.9}(r/150 kpc)^2 solar masses, similar to the mass in the low-ion phase. Accounting for the O VI phase approximately doubles the contribution of low-z LLSs to the cosmic baryon budget.

Recent theoretical work has implicated fast bar formation modes and subsequent evolution as the creation mechanism for the observed peanut-shaped bulges in some edge-on disk galaxies. We demonstrate an N-body simulation of a disk undergoing a contrasting slow mode of bar growth, unsubjected to a buckling instability, which nonetheless grows the 4:1 orbit family responsible for a peanut-shaped bulge. We also present a simulation with fast mode bar growth, which exhibits thickening similar to other work. A novel orbit classification method that finds dynamically distinct families is presented, allowing for a detailed analysis of angular momentum transfer channels within the disk.

We study the physical conditions of the circum-galactic medium (CGM) around z=0.25 galaxies as traced by HI and metal line absorption, using cosmological hydrodynamic simulations that include galactic outflows. Using lines of sight targeted at impact parameters from 10 kpc to 1 Mpc around galaxies with halo masses from 10^11-10^13 M_solar, we study the physical conditions and their variation with impact parameter b and line-of-sight velocity delta v in the CGM as traced by HI, MgII, SiIV, CIV, OVI, and NeVIII absorbers. All ions show a strong excess of absorption near galaxies compared to random lines of sight. The excess continues beyond 1 Mpc, reflecting the correlation of metal absorption with large-scale structure. Absorption is particularly enhanced within about v<300 km/sec and roughly 300 kpc of galaxies (with distances somewhat larger for the highest ion), approximately delineating the CGM; this range contains the majority of global metal absorption. Low ions like MgII and SiIV predominantly arise in denser gas closer to galaxies and drop more rapidly with b, while high ions OVI and NeVIII trace more diffusely distributed gas with a comparatively flat radial profile; CIV is intermediate. All ions predominantly trace T~10^4-4.5 K photo-ionised gas at all b, but when hot CGM gas is present (mostly in larger halos), we see strong collisionally-ionised OVI and NeVIII at b <= 100 kpc. Larger halo masses generally produce more absorption, though overall the trends are not as strong as that with impact parameter. These findings arise using our favoured outflow scalings as expected for momentum-driven winds; with no winds, the CGM gas remains mostly unenriched, while our outflow model with a constant velocity and mass loading factor produce hotter, more widely dispersed metals.

We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisation level of absorbers was a primary factor in determining the physical conditions of absorbing gas, we show here that it is also a governing factor for its dynamical state. Low-ionisation metal absorbers (e.g. MgII) tend to arise in gas that will fall onto galaxies within several Gyr, while high-ionisation metal absorbers (e.g. OVI) generally trace material that was deposited by outflows many Gyr ago. Inflowing gas is dominated by enriched material that was previously ejected in an outflow, hence accretion at low redshifts is typically substantially enriched. Recycling wind material is preferentially found closer to galaxies, and is more dominant in lower-mass halos since high-mass halos have more hot gas that is able to support itself against infall. Low-mass halos also tend to re-eject more of their accreted material, owing to our outflow prescription that employs higher mass loading factors for lower-mass galaxies. Typical HI absorbers trace unenriched ambient material that is not participating in the baryon cycle, but stronger HI absorbers arise in cool, enriched inflowing gas. Instantaneous radial velocity measures of absorbers are generally poor at distinguishing between inflowing and outflowing gas, except in the case of very recent outflows. These results suggest that probing halo gas using a range of absorbers can provide detailed information about the amount and physical conditions of material that is participating in the baryon cycle.

We present the design and methods of the COS-Halos survey, a systematic investigation of the gaseous halos of 44 z = 0.15-0.35 galaxies using background QSOs observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. This survey has yielded 39 spectra of z_em ~ 0.5 QSOs with S/N ~ 10-15 per resolution element. The QSO sightlines pass within 150 physical kpc of the galaxies, which span early and late types over stellar mass log M* / Msun= 9.5 - 11.5. We find that the CGM exhibits strong HI, averaging 1 Ang in Lya equivalent width out to 150 kpc, with 100% covering fraction for star-forming galaxies and 75% covering for passive galaxies. We find good agreement in column densities between this survey and previous studies over similar range of impact parameter. There is weak evidence for a difference between early- and late-type galaxies in the strength and distribution of HI. Kinematics indicate that the detected material is bound to the host galaxy, such that >~90% of the detected column density is confined within +/-200 km s^-1 of the galaxies. This material generally exists well below the halo virial temperatures at T<~ 10^5 K. We evaluate a number of possible origin scenarios for the detected material, and in the end favor a simple model in which the bulk of the detected HI arises in a bound, cool, low-density photoionized diffuse medium that is generic to all L* galaxies and may harbor a total gaseous mass comparable to galactic stellar masses.

We examine the global HI properties of galaxies in quarter-billion particle cosmological simulations using Gadget-2, focusing on how galactic outflows impact HI content. We consider four outflow models, including a new one (ezw) motivated by recent interstellar medium simulations in which the wind speed and mass loading factor scale as expected for momentum-driven outflows for larger galaxies and energy-driven outflows for dwarfs (sigma<75 km/s). To obtain predicted HI masses, we employ a simple but effective local correction for particle self-shielding, and an observationally-constrained transition from neutral to molecular hydrogen. Our ezw simulation produces an HI mass function whose faint-end slope of -1.3 agrees well with observations from the ALFALFA survey; other models agree less well. Satellite galaxies have a bimodal distribution in HI fraction versus halo mass, with smaller satellites and/or those in larger halos more often being HI-deficient. At a given stellar mass, HI content correlates with star formation rate and inversely correlates with metallicity, as expected if driven by stochasticity in the accretion rate. To higher redshifts, massive HI galaxies disappear and the mass function steepens. The global cosmic HI density conspires to remain fairly constant from z~5-0, but the relative contribution from smaller galaxies increases with redshift.

Weinberg (2012) described a constructive algorithm for computing the marginal likelihood, Z, from a Markov chain simulation of the posterior distribution. Its key point is: the choice of an integration subdomain that eliminates subvolumes with poor sampling owing to low tail-values of posterior probability. Conversely, this same idea may be used to choose the subdomain that optimizes the accuracy of Z. Here, we explore using the simulated distribution to define a small region of high posterior probability, followed by a numerical integration of the sample in the selected region using the volume tessellation algorithm described in Weinberg (2012). Even more promising is the resampling of this small region followed by a naive Monte Carlo integration. The new enhanced algorithm is computationally trivial and leads to a dramatic improvement in accuracy. For example, this application of the new algorithm to a four-component mixture with random locations in 16 dimensions yields accurate evaluation of Z with 5% errors. This enables Bayes-factor model selection for real-world problems that have been infeasible with previous methods.

We report new observations of circumgalactic gas in the halos of early type galaxies obtained by the COS-Halos Survey with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We find that detections of HI surrounding early type galaxies are typically as common and strong as around star-forming galaxies, implying that the total mass of circumgalactic material is comparable in the two populations. For early type galaxies, the covering fraction for HI absorption above 10^16 cm^2 is ~40-50% within ~150 kpc. Line widths and kinematics of the detected material show it to be cold (T ~< 10^5 K) in comparison to the virial temperature of the host halos. The implied masses of cool, photoionized CGM baryons may be up to 10^9 --- 10^11 Msun. Contrary to some theoretical expectations, strong halo HI absorbers do not disappear as part of the quenching of star-formation. Even passive galaxies retain significant reservoirs of halo baryons which could replenish the interstellar gas reservoir and eventually form stars. This halo gas may feed the diffuse and molecular gas that is frequently observed inside ETGs.

Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearly monotonic relation between stellar mass and halo mass provided that, for satellite galaxies, we use the mass of the subhalo at the epoch when it became a satellite. In each simulation, the median relation for central and satellite galaxies is nearly identical, though a somewhat larger fraction of satellites are outliers. SHAM-assigned masses (at z=0-2), luminosities (R-band at z=0), or star formation rates (at z=2) have a 68% scatter of 0.09-0.15 dex relative to the true simulation values. When we apply SHAM to the subhalo population of collisionless N-body simulation with the same initial conditions as the SPH runs, we find generally good agreement for the halo occupation distributions and halo radial profiles of galaxy samples defined by thresholds in stellar mass. However, because a small fraction of SPH galaxies suffer severe stellar mass loss after becoming satellites, SHAM slightly overpopulates high mass halos; this effect is more significant for the wind simulation, which produces galaxies that are less massive and more fragile. SHAM recovers the two-point correlation function of the SPH galaxies in the no-wind simulation to better than 10% at scales 0.1 < r < 10 Mpc/h. For the wind simulation, agreement is better than 15% at r > 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales.

We conduct Bayesian model inferences from the observed K-band luminosity function of galaxies in the local Universe, using the semi-analytic model (SAM) of galaxy formation introduced in Lu et al (2011). The prior distributions for the 14 free parameters include a large range of possible models. We find that some of the free parameters, e.g. the characteristic scales for quenching star formation in both high-mass and low-mass halos, are already tightly constrained by the single data set. The posterior distribution includes the model parameters adopted in other SAMs. By marginalising over the posterior distribution, we make predictions that include the full inferential uncertainties for the colour-magnitude relation, the Tully-Fisher relation, the conditional stellar mass function of galaxies in halos of different masses, the HI mass function, the redshift evolution of the stellar mass function of galaxies, and the global star formation history. Using posterior predictive checking with the available observational results, we find that the model family (i) predicts a Tully-Fisher relation that is curved; (ii) significantly over predicts the satellite fraction; (iii) vastly over predicts the HI mass function; (iv) predicts high-z stellar mass functions that have too many low mass galaxies and too few high mass ones. and (v) predicts a redshift evolution of the stellar mass density and the star formation history that are in moderate disagreement. These results suggest that some important processes are still missing in the current model family and we discuss a number of possible solutions to solve the discrepancies, such as interactions between galaxies and dark matter halos, tidal stripping, the bimodal accretion of gas, preheating, and a redshift-dependent initial mass function.

The circumgalactic medium (CGM) is fed by galaxy outflows and accretion of intergalactic gas, but its mass, heavy element enrichment, and relation to galaxy properties are poorly constrained by observations. In a survey of the outskirts of 42 galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, we detected ubiquitous, large (150 kiloparsec) halos of ionized oxygen surrounding star-forming galaxies, but we find much less ionized oxygen around galaxies with little or no star formation. This ionized CGM contains a substantial mass of heavy elements and gas, perhaps far exceeding the reservoirs of gas in the galaxies themselves. It is a basic component of nearly all star-forming galaxies that is removed or transformed during the quenching of star formation and the transition to passive evolution.

We investigate the metallicity evolution and content of the intergalactic medium (IGM) and galactic halo gas from z=2->0 using 110-million particle cosmological hydrodynamic simulations. We focus on the detectability and physical properties of UV resonance metal-line absorbers observable with Hubble's Cosmic Origins Spectrograph (COS). We confirm that galactic superwind outflows are required to enrich the IGM to observed levels down to z=0 using three wind prescriptions contrasted to a no-wind simulation. Our favoured momentum-conserved wind prescription deposits metals closer to galaxies owing to its moderate energy input, while the more energetic constant wind model enriches the warm-hot IGM 6.4x more. Despite these significant differences, all wind models produce metal-line statistics within a factor of two of existing observations. This is because OVI, CIV, SiIV, and NeVIII absorbers primarily arise from T<10^5 K, photo-ionised gas that is enriched to similar levels in the three feedback schemes. OVI absorbers trace the diffuse phase with delta<100, which is enriched to ~1/50 Zsol at z=0, although the absorbers themselves usually exceed 0.3 Zsol and arise from inhomogeneously distributed, un-mixed winds. CIV and SiIV absorbers trace primarily T~10^4 K gas inside haloes. We predict COS will observe a population of NeVIII photo-ionised absorbers tracing T<10^5 K, delta~10 gas. MgX and SiXII are rarely detected in COS S/N=30 simulated sight lines although detections trace T=10^(6-7) K halo gas. In general, the IGM is enriched in an outside-in manner, where wind-blown metals released at higher redshift reach lower overdensities, resulting in higher ionisation species tracing lower-density, older metals. At z=0, the 90% of baryons outside of galaxies are enriched to 0.096 Zsol, but the 65% of unbound baryons in the IGM have 0.018 Zsol and contain only 4% of all metals.

We believe that a wide range of physical processes conspire to shape the observed galaxy population but we remain unsure of their detailed interactions. The semi-analytic model (SAM) of galaxy formation uses multi-dimensional parameterisations of the physical processes of galaxy formation and provides a tool to constrain these underlying physical interactions. Because of the high dimensionality, the parametric problem of galaxy formation may be profitably tackled with a Bayesian-inference based approach, which allows one to constrain theory with data in a statistically rigorous way. In this paper we develop a SAM in the framework of Bayesian inference. We show that, with a parallel implementation of an advanced Markov-Chain Monte-Carlo algorithm, it is now possible to rigorously sample the posterior distribution of the high-dimensional parameter space of typical SAMs. As an example, we characterise galaxy formation in the current $\Lambda$CDM cosmology using the stellar mass function of galaxies as an observational constraint. We find that the posterior probability distribution is both topologically complex and degenerate in some important model parameters, suggesting that thorough explorations of the parameter space are needed to understand the models. We also demonstrate that because of the model degeneracy, adopting a narrow prior strongly restricts the model. Therefore, the inferences based on SAMs are conditional to the model adopted. Using synthetic data to mimic systematic errors in the stellar mass function, we demonstrate that an accurate observational error model is essential to meaningful inference.

We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations.

Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volume. Our primary model uses a simulation with galactic winds and assumes that dust traces the intergalactic metals. The predicted galaxy-dust correlation function is similar in form to the galaxy-mass correlation function, and reproducing the MSFR data requires a dust-to-metal mass ratio of 0.24, about half the value in the Galactic ISM. Roughly half of the reddening arises in dust that is more than 100Kpc/h from the nearest massive galaxy. We also examine a simulation with no galactic winds, which predicts a much smaller fraction of intergalactic metals (3% vs. 35%) and therefore requires an unphysical dust-to-metal ratio of 2.18 to reproduce the MSFR data. In both models, the signal is dominated by sightlines with E(g-i)=0.001-0.1. The no-wind simulation can be reconciled with the data if we also allow reddening to arise in galaxies up to several x 10^10 Msun. The wind model predicts a mean visual extinction of A_V ~0.0133 mag out to z=0.5, with a sightline-to-sightline dispersion similar to the mean, which could be significant for future supernova cosmology studies. Reproducing the MSFR results in these simulations requires that a large fraction of ISM dust survive its expulsion from galaxies and its residence in the intergalactic medium. Future observational studies that provide higher precision and measure the dependence on galaxy type and environment will allow detailed tests for models of enriched galactic outflows and the survival of IG dust.

We introduce a new galaxy image decomposition tool, GALPHAT (GALaxy PHotometric ATtributes), to provide full posterior probability distributions and reliable confidence intervals for all model parameters. GALPHAT is designed to yield a high speed and accurate likelihood computation, using grid interpolation and Fourier rotation. We benchmark this approach using an ensemble of simulated Sersic model galaxies over a wide range of observational conditions: the signal-to-noise ratio S/N, the ratio of galaxy size to the PSF and the image size, and errors in the assumed PSF; and a range of structural parameters: the half-light radius $r_e$ and the Sersic index $n$. We characterise the strength of parameter covariance in Sersic model, which increases with S/N and $n$, and the results strongly motivate the need for the full posterior probability distribution in galaxy morphology analyses and later inferences. The test results for simulated galaxies successfully demonstrate that, with a careful choice of Markov chain Monte Carlo algorithms and fast model image generation, GALPHAT is a powerful analysis tool for reliably inferring morphological parameters from a large ensemble of galaxies over a wide range of different observational conditions. (abridged)

The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies.

The intergalactic medium (IGM) is the dominant reservoir of baryons at all cosmic epochs. We investigate the evolution of the IGM from z=2-0 in 48 Mpc/h, 110-million particle cosmological hydrodynamic simulations using three prescriptions for galactic outflows. We focus on the evolution of IGM physical properties, and how such properties are traced by Ly-alpha absorption as detectable using HST/COS. Our results broadly confirm the canonical picture that most Ly-alpha absorbers arise from highly ionized gas tracing filamentary large-scale structure. Growth of structure causes gas to move from the diffuse photoionized IGM into other cosmic phases, namely stars, cold and hot gas within galaxy halos, and the unbound and shock-heated warm-hot intergalactic medium (WHIM). By today, baryons are roughly equally divided between bound phases (35%), the diffuse IGM (41%), and the WHIM (24%). Here we (re)define the WHIM as gas with overdensities lower than that in halos and temperatures >10^5 K, in order to more closely align it with "missing baryons". When we tune our photoionizing background to match the observed evolution of the Ly-alpha mean flux decrement, we obtain a line count evolution that broadly agrees with available data. We predict a column density distribution slope of -1.70 for our favored momentum-driven wind model, in agreement with recent observations, and it becomes shallower with redshift. With improved statistics, the frequency of strong lines can be a valuable diagnostic of outflows, and our favored wind model matches existing data best among our models. The relationship between column density and physical density is fairly tight from z=2-0, and evolves as rho N_HI^0.74 10^(-0.37z) for diffuse absorbers. Linewidths only loosely reflect the temperature of the absorbing gas, which will hamper attempts to quantify the WHIM using broad Ly-alpha absorbers. [Abridged]

We study the nature of rapidly star-forming galaxies at z=2 in cosmological hydrodynamic simulations, and compare their properties to observations of sub-millimetre galaxies (SMGs). We identify simulated SMGs as the most rapidly star-forming systems that match the observed number density of SMGs. In our models, SMGs are massive galaxies sitting at the centres of large potential wells, being fed by smooth infall and gas-rich satellites at rates comparable to their star formation rates (SFR). They are not typically undergoing major mergers that significantly boost their quiescent SFR, but they still often show complex gas morphologies and kinematics. Our simulated SMGs have stellar masses of log M*/Mo~11-11.7, SFRs of ~180-500 Mo/yr, a clustering length of 10 Mpc/h, and solar metallicities. The SFRs are lower than those inferred from far-IR data by a factor of 3, which we suggest may owe to one or more systematic effects in the SFR calibrations. SMGs at z=2 live in ~10^13 Mo halos, and by z=0 they mostly end up as brightest group galaxies in ~10^14 Mo halos. We predict that higher-M* SMGs should have on average lower specific SFRs, less disturbed morphologies, and higher clustering. We also predict that deeper far-IR surveys will smoothly join SMGs onto the massive end of the SFR-M* relationship defined by lower-mass z=2 galaxies. Overall, our simulated rapid star-formers provide as good a match to available SMG data as merger-based scenarios, offering an alternative scenario that emerges naturally from cosmological simulations.

We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the "hot" and "cold" modes emphasised in recent theoretical studies. Because of the interaction between outflows and gas in and around halos, the recycling timescale of wind material (t_rec) is shorter in higher-mass systems, which reside in denser gaseous environments. In these simulations, this differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF). If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed; galaxy masses are suppressed both by the direct removal of gas and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10^9 < M < 5x10^10 Msolar. At higher masses, wind recycling leads to excessive galaxy masses and excessive star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF.

We examine the impact of gas pressure on the transverse coherence of high-redshift (2 <= z <= 4) Lyman-alpha forest absorption along neighboring lines of sight that probe the gas Jeans scale (projected separation Delta r <= 500 kpc/h comoving; angular separation Delta theta <= 30"). We compare predictions from two smoothed particle hydrodynamics (SPH) simulations that have different photoionization heating rates and thus different temperature-density relations in the intergalactic medium (IGM). We also compare spectra computed from the gas distributions to those computed from the pressureless dark matter. The coherence along neighboring sightlines is markedly higher for the hotter, higher pressure simulation, and lower for the dark matter spectra. We quantify this coherence using the flux cross-correlation function and the conditional distribution of flux decrements as a function of transverse and line-of-sight (velocity) separation. Sightlines separated by Delta theta <= 15" are ideal for probing this transverse coherence. Higher pressure decreases the redshift-space anisotropy of the flux correlation function, while higher thermal broadening increases the anisotropy. In contrast to the longitudinal (line-of-sight) structure of the Lya forest, the transverse structure on these scales is dominated by pressure effects rather than thermal broadening. With the rapid recent growth in the number of known close quasar pairs, paired line-of-sight observations offer a promising new route to probe the IGM temperature-density relation and test the unexpectedly high temperatures that have been inferred from single sightline analyses.

In the low density intergalactic medium (IGM) that gives rise to the Lyman-alpha forest, gas temperature and density are tightly correlated. The velocity scale of thermal broadening and the Hubble flow across the gas Jeans scale are of similar magnitude (Hlambda_J ~ sigma_th). To separate the effects of gas pressure support and thermal broadening on the Lya forest, we compare spectra extracted from two smoothed particle hydrodynamics (SPH) simulations evolved with different photoionization heating rates (and thus different Jeans scales), imposing different temperature-density relations on the evolved particle distributions. The turnover scales in the flux power spectrum and flux autocorrelation function are determined mainly by thermal broadening rather than pressure. However, the insensitivity to pressure arises partly from a cancellation effect with a sloped temperature-density relation (T ~ rho^{0.6} in our simulations): the high density peaks in the colder, lower pressure simulation are less smoothed by pressure support than in the hotter simulation, and it is this higher density gas that experiences the strongest thermal broadening. Changes in thermal broadening and pressure support have comparably important effects on the flux probability distribution (PDF), which responds directly to the gas overdensity distribution rather than the scale on which it is smooth. Tests on a lower resolution simulation show that our statistical results are converged even at this lower resolution. While thermal broadening generally dominates the longitudinal structure in the Lya forest, we show in Paper II that pressure support determines the transverse coherence of the forest observed towards close quasar pairs. [ABRIDGED]

We investigate the physical mechanisms of tidal heating and satellite disruption in cold dark matter host haloes using N-body simulations based on cosmological initial conditions. We show the importance of resonant shocks and resonant torques with the host halo to satellite heating. A resonant shock (torque) couples the radial (tangential) motion of a satellite in its orbit to its phase space. For a satellite on a circular orbit, an ILR-like resonance dominates the heating and this heating results in continuous satellite mass loss. We estimate the requirements for simulations to achieve these dynamics using perturbation theory. Both resonant shocks and resonant torques affect satellites on eccentric orbits. We demonstrate that satellite mass loss is an outside-in process in energy space; a satellite's stars and gas are thus protected by their own halo against tidal stripping. We simulate the evolution of a halo similar to the Large Magellanic Cloud (LMC) in our Galactic dark matter halo and conclude that the LMC stars have not yet been stripped. Finally, we present a simple algorithm for estimating the evolution of satellite mass that includes both shock heating and resonant torques.

We present predictions for the fluorescent Lyman-alpha emission signature arising from photoionized, optically thick structures in Smoothed Particle Hydrodynamic (SPH) cosmological simulations of a Lambda-CDM universe using a Monte Carlo Lyman-alpha radiative transfer code. We calculate the expected Lyman-alpha image and 2-dimensional spectra for gas exposed to a uniform ultraviolet ionizing background as well as gas exposed additionally to the photoionizing radiation from a local quasar, after correcting for the self-shielding of hydrogen. As a test of our numerical methods and for application to current observations, we examine simplified analytic structures that are uniformly or anisotropically illuminated. We compare these results with recent observations. We discuss future observing campaigns on large telescopes and realistic strategies for detecting fluorescence owing to the ambient metagalactic ionization and in regions close to bright quasars. While it will take hundreds of hours on the current generation of telescopes to detect fluorescence caused by the ultraviolet background (UVB) alone, our calculations suggest that of order ten sources of quasar-induced fluorescent Lyman-alpha emission should be detectable after a 10 hour exposure in a 10 arcmin^2 field around a bright quasar. These observations will help probe the physical conditions in the densest regions of the intergalactic medium as well as the temporal light curves and isotropy of quasar radiation.

We compare the properties of galaxies that form in a cosmological simulation without strong feedback to observations at z=0. We confirm previous findings that models without strong feedback overproduce the observed galaxy baryonic mass function, especially at the low and high mass extremes. Through post-processing we investigate what kinds of feedback would be required to reproduce observed galaxy masses and star formation rates. To mimic an extreme form of "preventive" feedback (e.g., AGN radio mode) we remove all baryonic mass that was originally accreted via "hot mode" from shock-heated gas. This does not bring the high mass end of the galaxy mass function into agreement with observations because much of the stellar mass in these systems formed at high redshift from baryons that originally accreted via "cold mode" onto lower mass progenitors. An efficient "ejective" feedback mechanism, such as supernova driven winds, must reduce the masses of these progenitors. Feedback must also reduce the masses of lower mass z=0 galaxies, which assemble at lower redshifts and have much lower star formation rates. If we monotonically re-map galaxy masses to reproduce the observed mass function, but retain the simulation's predicted star formation rates, we obtain fairly good agreement with the observed sequence of star-forming galaxies but fail to recover the observed population of passive, low star formation rate galaxies. Suppressing all hot mode accretion improves agreement for high mass galaxies but worsens the agreement at intermediate masses. Reproducing these z=0 observations requires a feedback mechanism that dramatically suppresses star formation in a fraction of galaxies, increasing with mass, while leaving star formation rates of other galaxies essentially unchanged.

Using dark matter halos traced by galaxy groups selected from the Sloan Digital Sky Survey Data Release 4, we find that about 1/4 of the faint galaxies ($\rmag >-17.05$, hereafter dwarfs) that are the central galaxies in their own halo are not blue and star forming, as expected in standard models of galaxy formation, but are red. In contrast, this fraction is about 1/2 for dwarf satellite galaxies. Many red dwarf galaxies are physically associated with more massive halos. In total, about $\sim 45$% of red dwarf galaxies reside in massive halos as satellites, while another $\sim 25$% have a spatial distribution that is much more concentrated towards their nearest massive haloes than other dwarf galaxies. We use mock catalogs to show that the reddest population of non-satellite dwarf galaxies are distributed within about 3 times the virial radii of their nearest massive halos. We suggest that this population of dwarf galaxies are hosted by low-mass halos that have passed through their massive neighbors, and that the same environmental effects that cause satellite galaxies to become red are also responsible for the red colors of this population of galaxies. We do not find any significant radial dependence of the population of dwarf galaxies with the highest concentrations, suggesting that the mechanisms operating on these galaxies affect color more than structure. However, over 30% of dwarf galaxies are red and isolated and their origin remains unknown.

We study the formation of galaxies in a (50 Mpc/h)^3 cosmological simulation (2x288^3 particles), evolved using the entropy conserving SPH code Gadget-2. Most of the baryonic mass in galaxies of all masses is originally acquired through filamentary "cold mode" accretion of gas that was never shock heated to its halo virial temperature, confirming the key feature of our earlier results obtained with a different SPH code (Keres et al. 2005). Atmospheres of hot, virialized gas develop in halos above ~2.5e11 Msun, a transition mass that is nearly constant from z=3 to z=0. Cold accretion persists in halos above the transition mass, especially at z>=2. It dominates the growth of galaxies in low mass halos at all times, and it is the main driver of the cosmic star formation history. Satellite galaxies have accretion rates similar to central galaxies of the same baryonic mass at high redshifts, but they have less accretion than comparable central galaxies at low redshift. Relative to our earlier results, the Gadget-2 simulations predict much lower rates of "hot mode" accretion from the virialized gas component of massive halos. At z<=1, typical hot accretion rates in halos above 5e12 Msun are below 1 Msun/yr, even though our simulation does not include AGN heating or other forms of "preventive" feedback. The inner density profiles of hot gas in these halos are shallow, with long associated cooling times. The cooling recipes typically used in semi-analytic models can overestimate the accretion rates in these halos by orders of magnitude, so such models may overemphasize the role of preventive feedback in producing observed galaxy masses and colors. A fraction of the massive halos develop cuspy profiles and significant cooling rates between z=1 and z=0, a redshift trend similar to the observed trend in the frequency of cooling flow clusters.

We introduce a novel image decomposition package, GALPHAT, that provides robust estimates of galaxy surface brightness profiles using Bayesian Markov Chain Monte Carlo. The GALPHAT-determined posterior distribution of parameters enables us to assign rigorous statistical confidence intervals to maximum a posteriori estimates and to test complex galaxy formation and evolution hypotheses. We describe the GALPHAT algorithm, assess its performance using test image data, and demonstrate that it has sufficient speed for production analysis of a large galaxy sample. Finally we briefly introduce our ongoing science program to study the distribution of galaxy structural properties in the local universe using GALPHAT.

Using a representative sample of 911 central galaxies (CENs) from the SDSS DR4 group catalogue, we study how the structure of the most massive members in groups and clusters depend on (1) galaxy stellar mass (Mstar), (2) dark matter halo mass of the host group (Mhalo), and (3) their halo-centric position. We establish and thoroughly test a GALFIT-based pipeline to fit 2D Sersic models to SDSS data. We find that the fitting results are most sensitive to the background sky level determination and strongly recommend using the SDSS global value. We find that uncertainties in the background translate into a strong covariance between the total magnitude, half-light size (r50), and Sersic index (n), especially for bright/massive galaxies. We find that n depends strongly on Mstar for CENs, but only weakly or not at all on Mhalo. Less (more) massive CENs tend to be disk (spheroid)-like over the full Mhalo range. Likewise, there is a clear r50-Mstar relation for CENs, with separate slopes for disks and spheroids. When comparing CENs with satellite galaxies (SATs), we find that low mass (<10e10.75 Msun/h^2) SATs have larger median n than CENs of similar Mstar. Low mass, late-type SATs have moderately smaller r50 than late-type CENs of the same Mstar. However, we find no size differences between spheroid-like CENs and SATs, and no structural differences between CENs and SATs matched in both mass and colour. The similarity of massive SATs and CENs shows that this distinction has no significant impact on the structure of spheroids. We conclude that Mstar is the most fundamental property determining the basic structure of a galaxy. The lack of a clear n-Mhalo relation rules out a distinct group mass for producing spheroids, and the responsible morphological transformation processes must occur at the centres of groups spanning a wide range of masses. (abridged)

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7e9 M[Sun]. Most friends-of-friends halos in the simulation have a distinct central galaxy, typically two to five times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive halos form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Since z=1, 27% of central galaxies (above 3e10 M[Sun]) and 22% of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain "central" objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with the observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.

We investigate the incidence of major mergers creating >10e11 Msun galaxies in present-day groups and clusters more massive than 2.5e13 Msun. We identify 38 pairs of massive galaxies with mutual tidal interaction signatures selected from >5000 galaxies with >5e10 Msun that reside in 845 such groups. We fit the images of each galaxy pair as the line-of-sight projection of symmetric models and identify mergers by the presence of residual asymmetries around each progenitor, such as off-center isophotes, broad tidal tails, and dynamical friction wakes. At the resolution and sensitivity of the SDSS, such mergers are found in 16% of high-mass, galaxy-galaxy pairs with magnitude differences of <1.5 and <30 kpc projected separations. We find that 90% of these mergers have nearly equal-mass progenitors with red-sequence colors and centrally-concentrated morphologies, the hallmarks of dissipationless merger simulations. Mergers at group centers are more common than between 2 satellites, but both are morphologically indistinguishable and we tentatively conclude that the latter are likely located at the dynamical centers of recently accreted subhalos. The frequency of central and satellite merging diminishes with group mass consistent with dynamical friction expectations. Based on reasonable assumptions, the centers of these massive halos are growing in stellar mass by 1-9% per Gyr, on average. Compared to all LRG-LRG mergers, we find a 2-9 times higher rate for their merging when restricted to these dense environments. Our results imply that the massive end of the galaxy population continues to evolve hierarchically at a measurable level, and that the centers of massive groups are the preferred environment for merger-driven galaxy assembly. (abridged)

(Abbreviated) We investigate the dynamical mechanisms responsible for producing tidal tails from dwarf satellites using N-body simulations. We identify two important dynamical co-conspirators: 1) the points where the attractive force of the host halo and satellite are balanced do not occur at equal distances from the satellite centre or at the same equipotential value for massive satellites, breaking the morphological symmetry of the leading and trailing tails; and 2) the escaped ejecta in the leading (trailing) tail continues to be decelerated (accelerated) by the satellite's gravity leading to large offsets of the ejecta orbits from the satellite orbit. The effect of the satellite's self gravity decreases only weakly with a decreasing ratio of satellite mass to host halo mass, demonstrating the importance of these effects over a wide range of subhalo masses. Not only will the morphology of the leading and trailing tails for massive satellites be different, but the observed radial velocities of the tails will be displaced from that of the satellite orbit; both the displacement and the peak radial velocity is proportional to satellite mass. If the tails are assumed to follow the progenitor satellite orbits, the tails from satellites with masses greater than 0.0001 of the host halo virial mass in a spherical halo will appear to indicate a flattened halo. Therefore, a constraint on the Milky Way halo shape using tidal streams requires mass-dependent modelling. Similarly, we compute the the distribution of tail orbits both in E_{r}-r^{-2} space and in E-L_{z} space, advocated for identifying satellite stream relics. The acceleration of ejecta by a massive satellite during escape spreads the velocity distribution and obscures the signature of a well-defined ``moving group'' in phase space.

We present an analysis of star formation and feedback recipes appropriate for galactic smoothed particle hydrodynamics simulations. Using an isolated Milky Way-like galaxy, we constrain these recipes based on well-established observational results. Our star formation recipe is based on that of Katz (1992) with the additional inclusion of physically motivated supernova feedback recipes. We propose a new feedback recipe in which type II supernovae are modelled using an analytical treatment of blastwaves. With this feedback mechanism and a tuning of other star formation parameters, the star formation in our isolated Milky Way-like galaxy is constant and follows the slope and normalisation of the observed Schmidt law. In addition, we reproduce the low density cutoff and filamentary structure of star formation observed in disk galaxies. Our final recipe will enable better comparison of cosmological N-body simulations with observations.

Galaxy-galaxy lensing measures the mean excess surface density DS(r) around a sample of lensing galaxies. We develop a method for combining DS(r) with the galaxy correlation function xi_gg(r) to constrain Omega_m and sigma_8, going beyond the linear bias model to reach the level of accuracy demanded by current and future measurements. We adopt the halo occupation distribution (HOD) framework, and we test its applicability to this problem by examining the effects of replacing satellite galaxies in the halos of an SPH simulation with randomly selected dark matter particles from the same halos. The difference between dark matter and satellite galaxy radial profiles has a ~10% effect on DS(r) at r<1 Mpc/h. However, if radial profiles are matched, the remaining impact of individual subhalos around satellite galaxies and environmental dependence of the HOD at fixed halo mass is <5% in DS(r) for 0.1<r<15 Mpc/h. We develop an analytic approximation for DS(r) that incorporates halo exclusion and scale-dependent halo bias, and we demonstrate its accuracy with tests against a suite of populated N-body simulations. We use the analytic model to investigate the dependence of DS(r) and the galaxy-matter correlation function xi_gm(r) on Omega_m and sigma_8, once HOD parameters for a given cosmological model are pinned down by matching xi_gg(r). The linear bias prediction is accurate for r>2 Mpc/h, but it fails at the 30-50% level on smaller scales. The scaling of DS(r) ~ Omega_m^a(r) sigma_8^b(r) approaches the linear bias expectation a=b=1 at r>10 Mpc/h, but a(r) and b(r) vary from 0.8 to 1.6 at smaller r. We calculate a fiducial DS(r) and scaling indices a(r) and b(r) for two SDSS galaxy samples; galaxy-galaxy lensing measurements for these samples can be combined with our predictions to constrain Omega_m and sigma_8.

We have used the HST/ACS images to identify 4700 Lyman break galaxies (LBGs) in GOODS. We present the results from a parametric analysis of the 2-D surface brightness profiles, for 1333 LBGs at z > 2.5 with rest-frame UV(1600 Angstrom) AB magnitude < -20.5. Based on the Sersic index, n, which measures the profile shape, we find that about 40% of LBGs at z=3 have light profiles close to exponential, and only 30% have the high concentrations seen for spheroids. About 30% of LBGs appear to have multiple cores or disturbed morphologies suggestive of close pairs or on-going mergers. The fraction of spheroid-like (n > 2.5) LBGs decreases by about 15% from z = 5 to 3. A comparison of LBGs with the starburst galaxies at z = 1.2, shows that the fraction of spheroid-like profiles is about 20% higher among LBGs. The ellipticity distribution for LBGs exhibits a pronounced skew towards high ellipticities (> 0.5), which cannot be explained by morphologies similar to the local disks and spheroids viewed at random orientations. The peak of the distribution evolves toward lower ellipticities, from 0.7 at z = 4 to 0.5 at z = 3. At z = 1.2 the distribution is relatively flat as seen among the present-day galaxies. The dominance of elongated morphologies among LBGs suggests that in a significant fraction of them we may be witnessing star-formation in clumps along gas-rich filaments, or the earliest gas-rich bars that encompass essentially the entire visible galaxy. Similar features are found to be ubiquitous in hydrodynamical simulations in which galaxy formation at high redshifts occurs in filamentary inflows of dynamically cold gas within the dark matter halos, and involves gas- rich mergers.

We compare the galaxy population of an SPH simulation to those predicted by the GalICS semi-analytic model and a stripped down version without supernova and AGN feedback. The SPH simulation and the no-feedback GalICS model make similar predictions for the baryonic mass functions of galaxies and for the dependence of these mass functions on environment and redshift. The two methods also make similar predictions for the galaxy content of dark matter haloes as a function of halo mass and for the gas accretion history of galaxies. Both the SPH and no-feedback GalICS models predict a bimodal galaxy population at z=0. The "red'' sequence of gas poor, old galaxies is populated mainly by satellite systems while, contrary to observations, the central galaxies of massive haloes lie on the "blue'' star-forming sequence as a result of continuing hot gas accretion at late times. Furthermore, both models overpredict the observed baryonic mass function, especially at the high mass end. In the full GalICS model, supernova-driven outflows reduce the masses of low and intermediate mass galaxies by about a factor of two. AGN feedback suppresses gas cooling in large haloes, producing a sharp cut-off in the baryonic mass function and moving the central galaxies of these massive haloes to the red sequence. Our results imply that the observational failings of the SPH simulation and the no-feedback GalICS model are a consequence of missing input physics rather than computational inaccuracies, that truncating gas accretion by satellite galaxies automatically produces a bimodal galaxy distribution with a red sequence, but that explaining the red colours of the most massive galaxies requires a mechanism like AGN feedback that suppresses the accretion onto central galaxies in large haloes.

The total extragalactic background radiation can be an important test of the global star formation history (SFH). Using direct observational estimates of the SFH, along with standard assumptions about the initial mass function (IMF), we calculate the total extragalactic background radiation and the observed stellar density today. We show that plausible SFHs allow a significant range in each quantity, but that their ratio is very tightly constrained. Current estimates of the stellar mass and extragalactic background are difficult to reconcile, as long as the IMF is fixed to the Salpeter slope above 1 Msun. The joint confidence interval of these two quantities only agrees with that determined from the allowed range of SFH fits at the 3-sigma level, and for our best-fit values the discrepancy is about a factor of two. Alternative energy sources that contribute to the background, such as active galactic nuclei (AGN), Population III stars, or decaying particles, appear unlikely to resolve the discrepancy. However, changes to the IMF allow plausible solutions to the background problem. The simplest is an average IMF with an increased contribution from stars around 1.5--4 Msun. A ``paunchy'' IMF of this sort could emerge as a global average if low mass star formation is suppressed in galaxies experiencing rapid starbursts. Such an IMF is consistent with observations of star-forming regions, and would help to reconcile the fossil record of star formation with the directly observed SFH.

By comparing a collisionless cosmological N-body simulation (DM) to an SPH simulation with the same initial conditions, we investigate the correspondence between the dark matter subhalos produced by collisionless dynamics and the galaxies produced by dissipative gas dynamics in a dark matter background. When galaxies in the SPH simulation become satellites in larger groups, they retain local dark matter concentrations (SPH subhalos) whose mass is typically five times their baryonic mass. The more massive subhalos of the SPH simulation have corresponding subhalos of similar mass and position in the DM simulation; at lower masses, there is fairly good correspondence, but some DM subhalos are in different spatial positions and some suffer tidal stripping or disruption. The halo occupation statistics of DM subhalos -- the mean number of subhalos, pairs, and triples as a function of host halo mass -- are very similar to those of SPH subhalos and SPH galaxies. Gravity of the dissipative baryon component amplifies the density contrast of subhalos in the SPH simulation, making them more resistant to tidal disruption. Relative to SPH galaxies and SPH subhalos, the DM subhalo population is depleted in the densest regions of the most massive halos. The good agreement of halo occupation statistics between the DM subhalo and SPH galaxy populations leads to good agreement of their two-point correlation functions and higher order moments on large scales. The depletion of DM subhalos in dense regions depresses their clustering at R<1 Mpc/h. In these simulations, the "conversation" between dark matter and baryons is mostly one-way, with dark matter dynamics telling galaxies where to form and how to cluster, but the "back talk" of the baryons influences small scale clustering by enhancing the survival of substructure in the densest environments.

N-body simulations predict that CDM halo-assembly occurs in two phases: 1) a fast accretion phase with a rapidly deepening potential well; and 2) a slow accretion phase characterised by a gentle addition of mass to the outer halo with little change in the inner potential well. We demonstrate, using one-dimensional simulations, that this two-phase accretion leads to CDM halos of the NFW form and provides physical insight into the properties of the mass accretion history that influence the final profile. Assuming that the velocities of CDM particles are effectively isotropised by fluctuations in the gravitational potential during the fast accretion phase, we show that gravitational collapse in this phase leads to an inner profile rho(r) ~ r^{-1}. Slow accretion onto an established potential well leads to an outer profile with rho(r) ~ r^{-3}. The concentration of a halo is determined by the fraction of mass that is accreted during the fast accretion phase. Using an ensemble of realistic mass accretion histories, we show that the model predictions of the dependence of halo concentration on halo formation time, and hence the dependence of halo concentration on halo mass, and the distribution of halo concentrations all match those found in cosmological N-body simulations. Using a simple analytic model that captures much of the important physics we show that the inner r^{-1} profile of CDM halos is a natural result of hierarchical mass assembly with a initial phase of rapid accretion.

This paper explores resonance-driven secular evolution between a bar and dark-matter halo using N-body simulations. We make direct comparisons to our analytic theory (Weinberg & Katz 2005) to demonstrate the great difficulty that an N-body simulation has representing these dynamics for realistic astronomical interactions. In a dark-matter halo, the bar's angular momentum is coupled to the central density cusp (if present) by the Inner Lindblad Resonance. Owing to this angular momentum transfer and self-consistent re-equilibration, strong realistic bars WILL modify the cusp profile, lowering the central densities within about 30% of the bar radius in a few bar orbits. Past results to the contrary (Sellwood 2006, McMillan & Dehnen 2005) may be the result of weak bars or numerical artifacts. The magnitude depends on many factors and we illustrate the sensitivity of the response to the dark-matter profile, the bar shape and mass, and the galaxy's evolutionary history. For example, if the bar length is comparable to the size of a central dark-matter core, the bar may exchange angular momentum without changing its pattern speed significantly. We emphasise that this apparently simple example of secular evolution is remarkably subtle in detail and conclude that an N-body exploration of any astronomical scenario requires a deep investigation into the underlying dynamical mechanisms for that particular problem to set the necessary requirements for the simulation parameters and method (e.g. particle number and Poisson solver). Simply put, N-body simulations do not divinely reveal truth and hence their results are not infallible. They are unlikely to provide useful insight on their own, particularly for the study of even more complex secular processes such as the production of pseudo-bulges and disk heating.

We cross correlate the well-defined and very complete spectroscopic Main Galaxy Sample (MGS) of 156,000 bright (r<17.5 mag) galaxies from the SDSS with 2MASS sources to explore the nature and completeness of the 2MASS K-band selection of nearby galaxies. 2MASS detects 90% of the MGS brighter than r=17 mag. For r<16, 93.1% of the MGS is found in the 2MASS Extended Source Catalog (XSC). These detections span the representative range of optical and near-IR galaxy properties, but with a surface brightness-dependent bias to preferentially miss the most blue and low-concentration sources, consistent with the most morphologically late-type galaxy population. An XSC completeness of 97.5% is achievable at bright magnitudes, with blue LSBs being the only major source of incompleteness, if one follows our careful matching criteria and weeds out spurious SDSS sources. We conclude that the rapid drop in XSC completeness at r>16 reflects the sharp surface-brightness limit of the extended source detection algorithm in 2MASS. As a result, the r>16 galaxies found in the XSC are over-representative in red early types and under-representative in blue latetypes. At r>16 the XSC suffers an additional selection effect from the 2-3" spatial resolution limit of 2MASS. Therefore, 2MASS continues to detect 90% of of the MGS at 16<r<17, but with a growing fraction found in the Point Source Catalog (PSC) only. Overall, one third of the MGS is detected in the 2MASS PSC but not the XSC. A combined K<13.57 and r<16 selection provides the most representative inventory of galaxies in the local cosmos with near-IR and optical measurements, and 90.8% completeness. Using data from SDSS-DR2, this sample contains 19,156 galaxies with a median redshift of 0.052. (abridged)

We present an overview of some of the issues surrounding current models of galaxy formation, highlighting recent insights obtained from cosmological hydrodynamic simulations. Detailed examination of gas accretion processes show a hot mode of gas cooling from near the halo's virial temperature, and a previously underappreciated cold mode where gas flows in along filaments on dynamical timescales, emitting its energy in line radiation. Cold mode dominates in systems with halo masses slightly smaller than the Milky Way and below, and hence dominates the global accretion during the heydey of galaxy formation. This rapid accretion path enables prompt assembly of massive galaxies in the early universe, and results in $z\sim 4$ galaxy properties in broad agreement with observations, with the most massive galaxies being the most rapid star formers. Massive galaxies today are forming stars at a much reduced rate, a trend called downsizing. The trend of downsizing is naturally reproduced in simulations, owing to a transition from cold mode accretion in the early growth phase to slower hot mode accretion once their halos grow large. However, massive galaxies at the present epoch are still observed to have considerably redder colors than simulations suggest, suggesting that star formation is not sufficiently truncated in models by the transition to hot mode, and that another process not included in current simulations is required to suppress star formation.

We construct merger trees for galaxies identified in a cosmological hydrodynamical simulation and use them to characterize predicted merger rates as a function of redshift, galaxy mass, and merger mass ratio. At z=0.3, we find a mean rate of 0.054 mergers per galaxy per Gyr above a 1:2 mass ratio threshold for massive galaxies (baryonic mass above 6.4E10 \Msun), but only 0.018 / Gyr for lower mass galaxies. The mass ratio distribution is proportional to R_merg^-1.2 for the massive galaxy sample, so high mass mergers dominate the total merger growth rate. The predicted rates increase rapidly with increasing redshift, and they agree reasonably well with observational estimates. A substantial fraction of galaxies do not experience any resolved mergers during the course of the simulation, and even for the high mass sample only 50% of galaxies experience a greater than 1:4 merger since z=1. Typical galaxies thus have fairly quiescent merger histories. We assign bulge-to-disk ratios to simulated galaxies by assuming that mergers above a mass ratio threshold R_major convert stellar disks into spheroids. With R_major of 1:4, we obtain a fairly good match to the observed dependence of early-type fraction on galaxy mass. However, the predicted fraction of truly bulge-dominated systems (f_bulge} > 0.8) is small, and producing a substantial population of bulge-dominated galaxies may require a mechanism that shuts off gas accretion at late times and/or additional processes (besides major mergers) for producing bulges.

Using high resolution cosmological simulations, we study hydrogen and helium gravitational cooling radiation. We focus on the HeII cooling lines, which arise from gas with a different temperature history (T_max ~ 10^5K) than HI line emitting gas. We examine whether three major atomic cooling lines, HI 1216A, HeII 1640A and HeII 304A are observable, finding that HI 1216A and HeII 1640A cooling emission at z=2-3 are potentially detectable with deep narrow band (R>100) imaging and/or spectroscopy from the ground. While the expected strength of HI 1216A cooling emission depends strongly on the treatment of the self-shielded phase of the IGM in the simulations, our predictions for the HeII 1640A line are more robust because the HeII 1640A emissivity is negligible below T~10^4.5 K and less sensitive to the UV background. Although HeII 1640A cooling emission is fainter than HI 1216A by at least a factor of 10 and, unlike HI 1216A, might not be resolved spatially with current observational facilities, it is more suitable to study gas accretion in the galaxy formation process because it is optically thin and less contaminated by the recombination lines from star-forming galaxies. The HeII 1640A line can be used to distinguish among mechanisms for powering the so-called "Lyman alpha blobs" -- including gravitational cooling radiation, photoionization by stellar populations, and starburst-driven superwinds -- because (1) HeII 1640A emission is limited to very low metallicity (log(Z/Z_sun) < -5.3) and Population III stars, and (2) the blob's kinematics are probed unambiguously through the HeII 1640A line width, which, for cooling radiation, is narrower (sigma < 400 km/s) than typical wind speeds.

We use recent observations of the HI-mass function to constrain galaxy formation. The data conflicts with the standard model where most of the gas in a low-mass dark matter halo is assumed to settle into a disk of cold gas that is depleted by star formation and supernova-driven outflows until the disk becomes gravitationally stable. A consistent model can be found if low-mass haloes are embedded in a preheated medium, with a specific gas entropy ~ 10Kev cm^2. Such a model simultaneously matches the faint-end slope of the galaxy luminosity function. We propose a preheating model where the medium around low-mass haloes is preheated by gravitational pancaking. Since gravitational tidal fields suppress the formation of low-mass haloes while promoting that of pancakes, the formation of massive pancakes precedes that of the low-mass haloes within them. We demonstrate that the progenitors of present-day dark matter haloes with M<10^{12}h^{-1}\msun were embedded in pancakes of masses $~5x10^{12}h^{-1}\msun$ at z~2. The formation of such pancakes heats the gas to a temperature of 5x10^5K and compresses it to an overdensity of ~10. Such gas has a cooling time that exceeds the age of the Universe at z~2, and has a specific entropy of ~15Kev cm^2, almost exactly the amount required to explain the stellar and HI mass functions. (Abridged)

We analyze the flux power spectrum and its covariance using simulated Lyman alpha forests. We find that pseudo-hydro techniques are good approximations of hydrodynamical simulations at high redshift. However, the pseudo-hydro techniques fail at low redshift because they are insufficient for characterizing some components of the low-redshift intergalactic medium, notably the warm-hot intergalactic medium. Hence, to use the low-redshift Lyman alpha flux power spectrum to constrain cosmology, one would need realistic hydrodynamical simulations. By comparing one-dimensional mass statistics with flux statistics, we show that the nonlinear transform between density and flux quenches the fluctuations so that the flux power spectrum is much less sensitive to cosmological parameters than the one-dimensional mass power spectrum. The covariance of the flux power spectrum is nearly Gaussian. As such, the uncertainties of the underlying mass power spectrum could still be large, even though the flux power spectrum can be precisely determined from a small number of lines of sight.

Only through resonances can non-axisymmetric features such as spiral arms and bars exert torques over large scales and change the overall structure of a near-equilibrium galaxy. We describe the resonant interaction mechanism in detail and derive explicit criteria for the particle number required to simulate these dynamical processes accurately using N-body simulations and illustrate them with numerical experiments. To do this, we perform direct numerical solution of perturbation theory and make detailed comparisons with N-body simulations. The criteria include: sufficient particle coverage in phase space near the resonance and enough particles to minimize gravitational potential fluctuations that will change the dynamics of the resonant encounter. Some of our more surprising findings are as follows. First, the Inner-Lindblad-like resonance (ILR), responsible for coupling the bar to the central halo cusp, requires almost 10^9 equal mass particles within the virial radius for a Milky-Way-like bar in an NFW profile. Second, orbits that linger near the resonance receive more angular momentum than orbits that move through the resonance quickly. Small-scale fluctuations present in state-of-the-art particle-particle simulations can knock orbits out of resonance, preventing them from lingering and, thereby, decrease the torque. The required particle numbers are sufficiently high for scenarios of interest that apparent convergence in particle number is misleading: the convergence is in the noise-dominated regime. State-of-the-art simulations are not adequate to follow all aspects of secular evolution driven by the bar-halo interaction. We present a procedure to test the requirements for individual N-body codes for the actual problem of interest. [abridged]

We examine the temperature history of gas accreted by forming galaxies in SPH simulations. About half the gas shock heats to roughly the virial temperature of the galaxy potential well before cooling, condensing, and forming stars, but the other half radiates its acquired gravitational energy at much lower temperatures, typically T<10^5 K, and the histogram of maximum gas temperatures is clearly bimodal. The "cold mode" of gas accretion dominates for low mass galaxies (M_baryon < 10^{10.3}Msun or M_halo < 10^{11.4}Msun), while the conventional "hot mode" dominates the growth of high mass systems. Cold accretion is often directed along filaments, allowing galaxies to efficiently draw gas from large distances, while hot accretion is quasi-spherical. The galaxy and halo mass dependence leads to redshift and environment dependence of cold and hot accretion rates, with cold mode dominating at high redshift and in low density regions today, and hot mode dominating in group and cluster environments at low redshift. Star formation rates closely track accretion rates, and we discuss the physics behind the observed environment and redshift dependence of galactic scale star formation. If we allowed hot accretion to be suppressed by conduction or AGN feedback, then the simulation predictions would change in interesting ways, perhaps resolving conflicts with the colors of ellipticals and the cutoff of the galaxy luminosity function. The transition between cold and hot accretion at M_h ~ 10^{11.4}Msun is similar to that found by Birnboim & Dekel (2003) using 1-d simulations and analytic arguments. The corresponding baryonic mass is tantalizingly close to the scale at which Kauffmann et al. (2003) find a marked shift in galaxy properties. We speculate on connections between these theoretical and observational transitions.

The halo occupation distribution (HOD) describes the relation between galaxies and dark matter at the level of individual dark matter halos. The properties of galaxies residing at the centers of halos differ from those of satellite galaxies because of differences in their formation histories. Using a smoothed particle hydrodynamics (SPH) simulation and a semi-analytic (SA) galaxy formation model, we examine the separate contributions of central and satellite galaxies to the HOD, more specifically to the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class. In agreement with earlier results for dark matter subhalos, we find that the mean occupation function <N> for galaxies above a baryonic mass threshold can be approximated by a step function for central galaxies plus a power law for satellites, and that the distribution of satellite numbers is close to Poisson at fixed halo mass. For galaxy samples defined by different baryonic mass thresholds, there is a nearly linear relation between the minimum halo mass Mmin required to host a central galaxy and the mass M1 at which an average halo hosts one satellite, with M1 ~ 14 Mmin (SPH) or M1 ~ 18 Mmin (SA). The mean occupation number of young galaxies exhibits a local minimum at M ~ 10 Mmin where halos are too massive to host a young central galaxy but not massive enough to host satellites. We show that the conditional galaxy mass function at fixed halo mass cannot be described by a Schechter function because central galaxies produce a "bump" at high masses. We suggest parameterizations for the HOD and the conditional luminosity function that can be used to model observed galaxy clustering. Many of our predictions are in good agreement with recent results inferred from clustering in the Sloan Digital Sky Survey.

[Abridged] We test the reliability of a method to measure the mean halo mass of Damped Ly-alpha absorbers (DLAs). The method is based on measuring the ratio of the cross-correlation between DLAs and galaxies to the auto-correlation of the galaxies themselves ($w_{\rm dg}/w_{\rm gg}$), which is (in linear theory) the ratio of their bias factor. This is shown to be true irrespective of the galaxy redshift distribution, provided that one uses the same galaxies for the two correlation functions. The method is applicable to all redshifts. Here, we focus on z=3 DLAs and we demonstrate that the method robustly constrains the mean DLA halo mass using smoothed particle hydrodynamics (SPH) cosmological simulations. If we use the bias formalism of Mo & White with the DLA and galaxy mass distributions of these simulations, we predict a bias ratio of 0.771. Direct measurement from the simulations of $w_{\rm dg}/w_{\rm gg}$ st yields a ratio of 0.73+/-0.08, in excellent agreement with that prediction. Equivalently, inverting the measured correlation ratio to infer a mean DLA halo mass yields (log. averaging, in solar units) <log(M_DLA)> =11.13+/-013, in excellent agreement with the true value in the simulations: 11.16. The cross- correlation method thus appears to yield a robust estimate of the average host halo mass even though the DLAs and the galaxies occupy a broad mass spectrum of halos, and massive halos contain multiple galaxies with DLAs. We show that the inferred mean DLA halo mass is independent of the galaxy sub-sample used, i.e. the cross-correlation technique is also reliable. Our results imply that the cross-correlation length between DLAs and LBGs is predicted to be, at most, 2.85 Mpc. Future observations will soon distinguish models in which DLAs are in low mass halos from those in which DLAs are in massive halos.

We study the soft X-ray emission (0.1 - 1 keV) from the Warm-Hot Intergalactic Medium (WHIM) in a hydrodynamic simulation of a Cold Dark Matter universe. Our main goal is to investigate how such emission can be explored with a combination of imaging and spectroscopy, and to motivate future X-ray missions. We first present high resolution images of the X-ray emission in several energy bands, in which emission from different ion species dominates. We pick three different areas to study the high resolution spectra of X-rays from the warm-hot IGM: (a) a galaxy group; (b) a filament and (c) an underluminous region. By taking into account the background X-ray emission from AGNs and foreground emission from the Galaxy, we compute composite X-ray spectra of the selected regions. We briefly investigate angular clustering of the soft-X-ray emission, finding a strong signal. Most interestingly, the combination of high spectral resolution and angular information allows us to map the emission from the WHIM in 3 dimensions. We cross-correlate the positions of galaxies in the simulation with this redshift map of emission and detect the presence of six different ion species (Ne IX, Fe XVII, O VII, O VIII, N VII, C VI) in the large-scale structure traced by the galaxies. Finally we show how such emission can be detected and studied with future X-ray satellites, with particular attention to a proposed mission, the Missing Baryon Explorer, or MBE. We present simulated observations of the WHIM gas with MBE.

We calculate the angular correlation function of galaxies in the Two Micron All Sky Survey. We minimize the possible contamination by stars, dust, seeing and sky brightness by studying their cross correlation with galaxy density, and limiting the galaxy sample accordingly. We measure the correlation function at scales between 1-18 arcdegs using a half million galaxies. We find a best fit power law to the correlation function has a slope of 0.76 and an amplitude of 0.11. However, there are statistically significant oscillations around this power law. The largest oscillation occurs at about 0.8 degrees, corresponding to 2.8 h^{-1} Mpc at the median redshift of our survey, as expected in halo occupation distribution descriptions of galaxy clustering. We invert the angular correlation function using Singular Value Decomposition to measure the three-dimensional power spectrum and find that it too is in good agreement with previous measurements. A dip seen in the power spectrum at small wavenumber k is statistically consistent with CDM-type power spectra. A fit of CDM-type power spectra to k < 0.2 h Mpc^{-1} give constraints of \Gamma_{eff}=0.116 and \sigma_8=0.96. This suggest a K_s-band linear bias of 1.1+/-0.2. This \Gamma_{eff} is different from the WMAP CMB derived value. On small scales the power-law shape of our power spectrum is shallower than that derived for the SDSS. These facts together imply a biasing of these different galaxies that might be nonlinear, that might be either waveband or luminosity dependent, and that might have a nonlocal origin.

We study the relationship between galaxy luminosity, color, and environment in a cosmological simulation of galaxy formation. We compare the predicted relationship with that found for SDSS galaxies and find that the model successfully predicts most of the qualitative features seen in the data, but also shows some interesting differences. Specifically, the simulation predicts that the local density around bright red galaxies is a strong increasing function of luminosity, but does not depend much on color at fixed luminosity. Moreover, we show that this is due to central galaxies in dark matter halos whose baryonic masses correlate strongly with halo mass. The simulation also predicts that the local density around blue galaxies is a strong increasing function of color, but does not depend much on luminosity at fixed color. We show that this is due to satellite galaxies in halos whose stellar ages correlate with halo mass. Finally, the simulation fails to predict the luminosity dependence of environment observed around low luminosity red galaxies. However, we show that this is most likely due to the simulation's limited resolution. A study of a higher resolution, smaller volume simulation suggests that this dependence is caused by the fact that all low luminosity red galaxies are satellites in massive halos, whereas intermediate luminosity red galaxies are a mixture of satellites in massive halos and central galaxies in less massive halos.

We use a large sample of galaxies from the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) to calculate galaxy luminosity and stellar mass functions in the local Universe. We estimate k-corrections, evolution corrections, and stellar mass-to-light ratios (M/Ls) by fitting the galaxy colors with simple models. Our optical and near-infrared luminosity functions agree with most recent literature optical and near-infrared determinations within the uncertainties. We argue that 2MASS is biased against low surface brightness galaxies, and use SDSS plus our knowledge of stellar populations to estimate that the true K-band luminosity function has a steeper faint end slope than the direct estimate. Assuming a universally-applicable stellar initial mass function (IMF), we derive the stellar mass function of galaxies. The faint end slope slope for the stellar mass function is steeper than -1.1, reflecting the low stellar M/Ls characteristic of low-mass galaxies. We estimate an upper limit to the stellar mass density in the local Universe Omega* h = 0.0020+/-0.0006 by assuming an IMF as rich in low-mass stars as allowed by observations of galaxy dynamics. Finally, we find that the characteristic luminosity or mass of early-type galaxies is larger than for later types, and the faint end slope is steeper for later types than for earlier types. Accounting for typing uncertainties, we estimate that at least half, and perhaps as much as 3/4, of the stellar mass in the Universe is in early-type galaxies. We present also SDSS/2MASS color-M/L correlations, an updated discussion of near-infrared stellar M/L estimates, and the volume-corrected distribution of g and K-band stellar M/Ls as a function of stellar mass. [Abridged]

The evolution of a stellar bar transforms not only the galactic disk, but also the host dark matter halo. We present high resolution, fully self-consistent N-body simulations that clearly demonstrate that dark matter halo central density cusps flatten as the bar torques the halo. This effect is independent of the bar formation mode and occurs even for rather short bars. The halo and bar evolution is mediated by resonant interactions between orbits in the halo and the bar pattern speed, as predicted by linear Hamiltonian perturbation theory. The bar lengthens and slows as it loses angular momentum, a process that occurs even in rather warm disks. We demonstrate that the bar and halo response can be critically underestimated for experiments that are unable to resolve the relevant resonant dynamics; this occurs when the phase space in the resonant region is under sampled or plagued by noise.

We use hydrodynamic simulations to predict correlations between Lya forest absorption and galaxies at redshift z~3. The probability distribution function (PDF) of Lya flux decrements shifts systematically towards higher values in the vicinity of galaxies, reflecting the overdense environments in which these galaxies reside. The predicted signal remains strong in spectra smoothed over 50-200 km/s, allowing tests with moderate resolution quasar spectra. The strong bias of high redshift galaxies towards high density regions imprints a clear signature on the flux PDF, but the predictions are not sensitive to galaxy baryon mass or star formation rate, and they are similar for galaxies and for dark matter halos. The dependence of the flux PDF on galaxy proximity is sensitive to redshift determination errors, with rms errors of 150-300 km/s substantially weakening the predicted trends. On larger scales, the mean galaxy overdensity in a cube of 5 or 10 Mpc/h (comoving) is strongly correlated with the mean Lya flux decrement on a line of sight through the cube center. The slope of the correlation is ~3 times steeper for galaxies than for dark matter as a result of galaxy bias. The predicted large scale correlation is in qualitative agreement with recently reported observational results. However, observations also show a drop in absorption in the immediate vicinity of galaxies, which our models do not predict even if we allow the galaxies or AGNs within them to be ionizing sources. This decreased absorption could be a signature of galaxy feedback on the surrounding IGM, perhaps via galactic winds. Peculiar velocities often allow gas at comoving distances ~1.5 Mpc/h to produce saturated absorption at the galaxy redshift, so any feedback mechanism must suppress neutral hydrogen out to these radii to match the data. (Abridged)

Using a hydrodynamic simulation of a LCDM universe, we investigate the "X-ray forest" absorption imprinted on the spectra of background quasars by the intervening intergalactic medium (IGM). In agreement with previous studies, we find that OVII and OVIII produce the strongest absorption features. The strong oxygen absorbers that might be detectable with Chandra or XMM-Newton arise in gas with T ~ 10^6 K and overdensities delta >~ 100 that are characteristic of galaxy groups. Future X-ray missions could detect weaker oxygen absorption produced by gas with a wider range of temperatures and the lower densities of unvirialized structures; they could also detect X-ray forest absorption by C, N, Ne, Fe, and possibly Si. If the IGM metallicity is 0.1 solar, then the predicted number of systems strong enough for a ~5\sigma detection with Chandra or XMM-Newton is extremely low, though scatter in metallicity would increase the number of strong absorbers even if the mean metallicity remained the same. Our simulation reproduces the high observed incidence of OVI absorbers (in the UV), and the most promising strategy for finding the X-ray forest is to search at the redshifts of known OVI systems, thus reducing the signal-to-noise threshold required for a significant detection. However, while many OVI absorbers have associated OVII or OVIII absorption, the OVI systems trace only the low temperature phases of the X-ray forest, and a full accounting of the strong OVII and OVIII systems will require a mission with the anticipated capabilities of Constellation-X. The large effective area of the XEUS satellite would make it an extremely powerful instrument for studying the IGM, measuring X-ray forest absorption by a variety of elements and revealing the shock-heated filaments that may be an important reservoir of cosmic baryons.

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (a) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (b) the relative distributions of galaxies and dark matter within halos. We calculate predicted HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation and a semi-analytic (SA) galaxy formation model. Although the two methods predict different galaxy mass functions, their HOD predictions agree remarkably well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more steeply rising high occupancy regime. At low <N>, the mean pair and triple counts are well below Poisson expectations, with important consequences for small scale behavior of 2- and 3-point correlation functions. The HOD depends strongly on galaxy age, with high mass halos populated mainly by old galaxies and low mass halos by young galaxies. The SPH simulation supports several simplifying assumptions about HOD bias: the most massive galaxy in a halo usually lies close to the center and moves near the halo's mean velocity; satellite galaxies have the same radial profile and velocity dispersion as the dark matter; and the mean occupation at fixed halo mass is independent of the halo's larger scale environment. By applying the SPH and SA HODs to a large volume N-body simulation, we show that both methods predict slight, observable departures from a power-law galaxy correlation function. The predicted HODs are closely tied to the underlying galaxy formation physics, they offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. (Shortened)

The unprecedented sky coverage and photometric uniformity of the Two Micron All Sky Survey (2MASS) provides a rich resource for investigating the galaxies populating the local Universe. A full characterization of the large-scale clustering distribution is important for theoretical studies of structure formation. 2MASS offers an all-sky view of the local galaxy population at 2.15 micron, unbiased by young stellar light and minimally affected by dust. We use 2MASS to map the local distribution of galaxies, identifying the largest structures in the nearby universe. The inhomogeneity of these structures causes an acceleration on the Local Group of galaxies, which can be seen in the dipole of the Cosmic Microwave Background (CMB). We find that the direction of the 2MASS clustering dipole is 11 degrees from the CMB dipole, confirming that the local galaxy distribution accelerates the Local Group. From the magnitude of the dipole we find a value of the linear bias parameter b=1.37 +/- 0.3 in the K_s-band. The 2MASS clustering dipole is 19 degrees from the latest measurement of the dipole using galaxies detected by the Infrared Astronomical Satellite (IRAS) suggesting that bias may be non-linear in some wavebands.

A common feature of the X-ray bubbles observed in Chandra images of some "cooling flow" clusters is that they appear to be surrounded by bright, cool shells. Temperature maps of a few nearby luminous clusters reveal that the shells consist of the coolest gas in the clusters -- much cooler than the surrounding medium. Using simple models, we study the effects of this cool emission on the inferred cooling flow properties of clusters. We find that the introduction of bubbles into model clusters that do NOT have cooling flows results in temperature and surface brightness profiles that resemble those seen in nearby "cooling flow" clusters. They also approximately reproduce the recent XMM-Newton and Chandra observations of a high minimum temperature of 1-3 keV. Hence, bubbles, if present, must be taken into account when inferring the physical properties of the ICM. In the case of some clusters, bubbles may account entirely for these observed features, calling into question their designation as clusters with cooling flows. However, since not all nearby "cooling flow" clusters show bubble-like features, we suggest that there may be a diverse range of physical phenomena that give rise to the same observed features.

We estimate the baryonic (stellar+cold gas) mass function of galaxies in the local Universe by assigning a complete sample of Two Micron All Sky Survey and Sloan Digital Sky Survey galaxies a gas fraction based on a statistical sample of the entire population, under the assumption of a universally-applicable stellar initial mass function. The baryonic mass function is similar to the stellar mass function at the high mass end, and has a reasonably steep faint-end slope owing to the typically high cold gas fractions and low stellar mass-to-light ratios characteristic of low-mass galaxies. The Schechter Function fit parameters are phi* = 0.0108(6) per cubic Mpc per dex of mass, M* = 5.3(3)x10^10 solar masses, and alpha = -1.21(5), with formal error estimates given in parentheses (for a Hubble constant of 100 km/s per Mpc). We show that the neutral and molecular hydrogen mass functions derived using this indirect route are in agreement with direct estimates, validating our indirect method. Integrating under the baryonic mass function and incorporating all sources of uncertainty, we find that the baryonic (stellar+cold gas) mass density implied by this estimate is Omega_cold baryon = 0.0024+0.0007-0.0014, or 8% +4% -5% of the Big Bang nucleosynthesis expectation.

We review recent developments in the theory of the Lyman-alpha forest and their implications for the role of the forest as a test of cosmological models. Simulations predict a relatively tight correlation between the local Lya optical depth and the local gas or dark matter density. Statistical properties of the transmitted flux can constrain the amplitude and shape of the matter power spectrum at high redshift, test the assumption of Gaussian initial conditions, and probe the evolution of dark energy by measuring the Hubble parameter H(z). Simulations predict increased Lya absorption in the vicinity of galaxies, but observations show a Lya deficit within Delta_r ~ 0.5 Mpc/h (comoving). We investigate idealized models of "winds" and find that they must eliminate neutral hydrogen out to comoving radii ~1.5 Mpc/h to marginally explain the data. Winds of this magnitude suppress the flux power spectrum by \~0.1 dex but have little effect on the distribution function or threshold crossing frequency. In light of the stringent demands on winds, we consider the alternative possibility that extended Lya emission from target galaxies replaces absorbed flux, but we conclude that this explanation is unlikely. Taking full advantage of the data coming from large telescopes and from the Sloan Digital Sky Survey will require more complete understanding of the galaxy proximity effect, careful attention to continuum determination, and more accurate numerical predictions, with the goal of reaching 5-10% precision on key cosmological quantities.

Lyman alpha absorber correlations contain information about the underlying density distribution associated with a particular class of absorbers. As such, they provide an opportunity to independently measure the "bias" of the Lyman alpha forest, i.e. the relationship between HI column density and underlying dark matter density. In these proceedings we use hydrodynamic simulations to investigate whether the evolution of this bias is measurable from observable correlations. Unfortunately, the increasingly complex physics in the IGM at z<1 makes a direct measurement of the bias difficult. Nevertheless, current simulations do make predictions for HI absorber correlations that are in broad agreement with observations at both high and low redshift, thus reinforcing the bias evolution predictions given by these models.

We investigate galaxy clustering and galaxy-mass correlations in the LCDM cosmological model using a large volume SPH simulation. For the most part, the predicted biases between galaxies and dark matter lead to good agreement with current observations, including: (1) a nearly constant comoving correlation length from z=3 to z=0 for mass-selected galaxy samples of constant comoving space density; (2) an rms bias factor b~1 at z=0; (3) a scale-dependent bias on small scales that transforms the curved dark matter correlation function into a nearly power-law galaxy correlation function; (4) galaxy pairwise dispersion and hierarchical skewness ratio S_3 in good agreement with observed values, and lower than values for the dark matter by ~20%; (5) a ratio of galaxy-galaxy to galaxy-mass correlation functions consistent with recent measurements from the RCS survey; and (6) a mean excess mass dM(260/h kpc) approximately proportional to galaxy baryon mass M_b, in agreement with lensing estimates from the SDSS. These clustering properties vary with galaxy baryon mass and, more strongly, with stellar population age, in good qualitative agreement with the observed dependence of galaxy clustering and galaxy-mass correlations on galaxy type. The predicted ratio dM(260)/M_b is lower than the SDSS estimates for galaxies with M_b>~2e11 Msun, but numerical resolution effects account for most of this discrepancy, leaving little room for feedback or other astrophysical processes to reduce the stellar masses of luminous galaxies, at least given our adopted cosmological parameters. Our results show that the LCDM model and the galaxy formation physics incorporated in the SPH simulation give a good account of observed galaxy clustering, but anticipated improvements in clustering and weak lensing measurements will soon test this picture in much more detail.

Hydrodynamic cosmological simulations predict that the average opacity of the Ly-alpha forest should increase in the neighborhood of galaxies because galaxies form in dense environments. Recent observations (Adelberger et al. 2002) confirm this expectation at large scales, but they show a decrease of absorption at comoving separations Delta_r <~ 1 Mpc/h. We show that this discrepancy is statistically significant, especially for the innermost data point at Delta_r <= 0.5 Mpc/h, even though this data point rests on three galaxy-quasar pairs. Galaxy redshift errors of the expected magnitude are insufficient to resolve the conflict. Peculiar velocities allow gas at comoving distances >~ 1 Mpc/h to produce saturated absorption at the galaxy redshift, putting stringent requirements on any ``feedback'' solution. Local photoionization is insufficient, even if we allow for recurrent AGN activity that keeps the neutral hydrogen fraction below its equilibrium value. A simple ``wind'' model that eliminates all neutral hydrogen in spheres around the observed galaxies can marginally explain the data, but only if the winds extend to comoving radii ~1.5 Mpc/h.

Not the way one might have thought. In hydrodynamic simulations of galaxy formation, some gas follows the traditionally envisioned route, shock heating to the halo virial temperature before cooling to the much lower temperature of the neutral ISM. But most gas enters galaxies without ever heating close to the virial temperature, gaining thermal energy from weak shocks and adiabatic compression, and radiating it just as quickly. This ``cold mode'' accretion is channeled along filaments, while the conventional, ``hot mode'' accretion is quasi-spherical. Cold mode accretion dominates high redshift growth by a substantial factor, while at z<1 the overall accretion rate declines and hot mode accretion has greater relative importance. The decline of the cosmic star formation rate at low z is driven largely by geometry, as the typical cross section of filaments begins to exceed that of the galaxies at their intersections.

The unprecedented sky coverage and photometric uniformity of the Two Micron All Sky Survey (2MASS) provides a rich resource for obtaining a detailed understanding of the galaxies populating our local (z<0.1) Universe. A full characterization of the physical structure of nearby galaxies is essential for theoretical and observational studies of galaxy evolution and structure formation. We have begun a quantified description of the internal structure and morphology of 10,000 bright (10<Ks<11) 2MASS galaxies through multi-component model fits to the 2D surface brightness profiles using GIM2D. From our initial Monte Carlo tests on 77 galaxies drawn at random from the RC3, we find that the model derived structural parameter errors due to sky uncertainies are typically less than 10%.

We measure the filling factor, correlation function, and power spectrum of transmitted flux in a large sample of Lya forest spectra, comprised of 30 Keck HIRES spectra and 23 Keck LRIS spectra. We infer the linear matter power spectrum P(k) from the flux power spectrum P_F(k), using an improved version of the method of Croft et al. (1998) that accounts for the influence of z-space distortions, non- linearity, and thermal broadening on P_F(k). The evolution of the shape and amplitude of P(k) over the range z= 2-4 is consistent with gravitational instability, implying that non-gravitational fluctuations do not make a large contribution. Our fiducial measurement of P(k) comes from data with <z> = 2.72. It has amplitude Delta^2(k_p)=0.74^0.20_-0.16 at wavenumber k_p=0.03 (km/s)^-1 and is well described by a power-law of index -2.43 +/- 0.06 or by a CDM-like power spectrum with shape parameter Gamma'=1.3^+0.7_-0.5*10^-3 (km/s) at z=2.72. For Omega_m=0.4, Omega_Lam=0.6, the best-fit Gamma =0.16 (h^-1mpc)^-1, in good agreement with the 2dF Galaxy Redshift Survey, and the best-fit sigma_8=0.82 (Gamma/0.15)^-0.44. Matching the observed cluster mass function and our Delta^2(k_p) in spatially flat models requires Omega_m=0.38^+0.10_-0.08 + 2.2 (Gamma-0.15). Matching Delta^2(k_p) in COBE-normalized, flat CDM models with no tensor fluctuations requires Omega_m = (0.29 +/-0.04) n^-2.89 h_65^-1.9. The Lya forest complements other probes of P(k) by constraining a regime of redshift and lengthscale not accessible by other means, and the consistency of these inferred parameters with independent estimates provides further support for inflation, cold dark matter, and vacuum energy (abridged).

We examine the scalings of X-ray luminosity, temperature, and dark matter or galaxy velocity dispersion for galaxy groups in an LCDM cosmological simulation, which incorporates gravity, gas dynamics, radiative cooling, and star formation, but no substantial non-gravitational heating. In agreement with observations, the simulated L_X-sigma and L_X-T_X relations are steeper than those predicted by adiabatic simulations or self-similar models, with L_X\propto sigma^{4.4} and L_X\propto T_X^{2.6} for massive groups and significantly steeper relations below a break at sigma~180 km/s (T_X~0.7 keV), but T_X-sigma is fairly close to self-similar, with T_X\propto sigma^{1.75}. The entropy of hot gas in low mass groups is higher than predicted by self-similar scaling, and agrees with observations that suggest an "entropy floor". The steeper scalings of the luminosity relations are driven by radiative cooling, which reduces the hot gas fraction from 50% of the total baryons at sigma \approx 500 km/s to 20% at sigma \approx 100 km/s. A secondary effect is that hot gas in smaller systems is less clumpy, further driving down L_X. A smaller volume simulation with eight times higher mass resolution predicts nearly identical X-ray luminosities at a given group mass, demonstrating the insensitivity of the predicted scaling relations to numerical resolution. There remain some quantitative discrepancies: the predicted mass scale of the L_X-T_X and L_X-sigma breaks is somewhat too low, and the luminosity-weighted temperatures are too high at a given sigma. We conclude that radiative cooling has an important quantitative impact on group X-ray properties and can account for many of the observed trends that have been interpreted as evidence for non-gravitational heating. (abridged)

Simulations predict that the dark matter halos of galaxies should have central cusps, while those inferred from observed galaxies do not have cusps. We demonstrate, using both linear perturbation theory and n-body simulations, that a disk bar, which should be ubiquitous in forming galaxies, can produce cores in cuspy CDM dark matter profiles within five bar orbital times. Simulations of forming galaxies suggest that one of Milky Way size could have a 10 kpc primordial bar; this bar will remove the cusp out to approximately 5 kpc in approximately 1.5 gigayears, while the disk only loses approximately 8% of its original angular momentum. An inner Lindblad-like resonance couples the rotating bar to orbits at all radii through the cusp, transferring the bar pattern angular momentum to the dark matter cusp, rapidly flattening it. This resonance disappears for profiles with cores and is responsible for a qualitative difference in bar driven halo evolution with and without a cusp. This bar induced evolution will have a profound effect on the structure and evolution of almost all galaxies. Hence, both to understand galaxy formation and evolution and to make predictions from theory it is necessary to resolve these dynamical processes. Unfortunately, correctly resolving these important dynamical processes in ab initio calculations of galaxy formation is a daunting task, requiring at least 4,000,000 halo particles using our SCF code, and probably requiring many times more particles when using noisier tree, direct summation, or grid based techniques, the usual methods employed in such calculations.

We describe techniques for comparing spectra extracted from cosmological simulations and observational data, using the same methodology to link Lyman-alpha properties derived from the simulations with properties derived from observational data. The eventual goal is to measure the coherence or clustering properties of Lyman-alpha absorbers using observations of quasar pairs and groups. We quantify the systematic underestimate in opacity that is inherent in the continuum fitting process of observed spectra over a range of resolution and signal-to-noise ratio. We present an automated process for detecting and selecting absorption features over the range of resolution and signal-to-noise of typical observational data on the Lyman-alpha "forest". Using these techniques, we detect coherence over transverse scales out to 500 h^{-1}_{50} kpc in spectra extracted from a cosmological simulation at z = 2.

We investigate the X-ray vs. optical scaling relations of poor groups to small clusters (sigma~100-700 km/s) identified in a cosmological hydrodynamic simulation of a Lambda-CDM universe, with cooling and star formation but no pre-heating. We find that the scaling relations between X-ray luminosity, X-ray temperature, and velocity dispersion show significant departures from the relations predicted by simple hydrostatic equilibrium models or simulations without cooling, having steeper L_X-sigma and L_X-T_X slopes and a "break" at \~200 km/s (~0.3 keV). These departures arise because the hot (X-ray emitting) gas fraction varies substantially with halo mass in this regime. Our predictions roughly agree with observations. Thus radiative cooling is a critical physical process in modeling galaxy groups, and may present an alternative to ad hoc models such as pre-heating or entropy floors for explaining X-ray group scaling relations.

This paper discusses a semi-numerical method of investigating the enrichment of the intergalactic medium by galactic winds. We find that most galaxies at z >~ 3 should be driving winds, and that (if these winds are similar to those at low-z) these winds should escape to large distances. Our calculations -- which permit exploration of a large region of model parameter space -- indicate that the wind velocity, the mass of the wind-driving galaxies, the fraction of ambient material entrained, and the available time (between wind launch and the observed redshift) all affect wind propagation significantly; other physical effects can be important but are sub-dominant. We find that under reasonable assumptions, the enrichment by 3 <~ z <~ 6 galaxies could account for the quantity of metals seen in the Ly-alpha forest, though it is presently unclear whether this enrichment is compatible with the intergalactic medium's detailed metal distribution or relative quiescence.

Observations have established that the diffuse intergalactic medium (IGM) at z ~ 3 is enriched to ~0.1-1% solar metallicity and that the hot gas in large clusters of galaxies (ICM) is enriched to 1/3-1/2 solar metallicity at z=0. Metals in the IGM may have been removed from galaxies (in which they presumably form) during dynamical encounters between galaxies, by ram-pressure stripping, by supernova-driven winds, or as radiation-pressure driven dust efflux. This study develops a method of investigating the chemical enrichment of the IGM and of galaxies, using already completed cosmological simulations. To these simulations, we add dust and (gaseous) metals, distributing the dust and metals in the gas according to three simple parameterized prescriptions, one for each enrichment mechanism. These prescriptions are formulated to capture the basic ejection physics, and calibrated when possible with empirical data. Our results indicate that dynamical removal of metals from >~ 3*10^8 solar mass galaxies cannot account for the observed metallicity of low-column density Ly-alpha absorbers, and that dynamical removal from >~ 3*10^10 solar mass galaxies cannot account for the ICM metallicities. Dynamical removal also fails to produce a strong enough mass-metallicity relation in galaxies. In contrast, either wind or radiation-pressure ejection of metals from relatively large galaxies can plausibly account for all three sets of observations (though it is unclear whether metals can be distributed uniformly enough in the low-density regions without overly disturbing the IGM, and whether clusters can be enriched quite as much as observed). We investigate in detail how our results change with variations in our assumed parameters, and how results for the different ejection processes compare. (Abridged)

We analyze high-redshift structure in three hydrodynamic simulations that have identical initial conditions and cosmological parameters and differ only in the value of the baryon density parameter, Omega_b=0.02, 0.05, 0.125. Increasing Omega_b does not change the fraction of baryons in the diffuse (unshocked) phase of the intergalactic medium, but it increases cooling rates and therefore transfers some baryons from the shocked intergalactic phase to the condensed phase associated with galaxies. Predictions of Lyman-alpha forest absorption are almost unaffected by changes of Omega_b provided that the UV background intensity is adjusted so that the mean opacity of the forest matches the observed value. The required UV background intensity scales as Omega_b^1.7, and the higher photoionization rate increases the gas temperature in low density regions. Damped Lyman-alpha absorption and Lyman limit absorption both increase with increasing Omega_b, though the impact is stronger for damped absorption and is weaker at z=4 than at z=2-3. The mass of cold gas and stars in high-redshift galaxies increases faster than Omega_b but slower than Omega_b^2, and the global star formation rate scales approximately as Omega_b^1.5. In the higher Omega_b models, the fraction of baryonic material within the virial radius of dark matter halos is usually higher than the universal fraction, indicating that gas dynamics and cooling can lead to over-representation of baryons in virialized systems. On the whole, our results imply a fairly intuitive picture of the influence of Omega_b on high-redshift structure, and we provide scalings that can be used to estimate the impact of Omega_b uncertainties on the predictions of hydrodynamic simulations.

We study the predicted sub-mm emission from massive galaxies in a Lambda-CDM universe, using hydrodynamic cosmological simulations. Assuming that most of the emission from newly formed stars is absorbed and reradiated in the rest-frame far-IR, we calculate the number of galaxies that would be detected in sub-mm surveys conducted with SCUBA. The predicted number counts are strongly dependent on the assumed dust temperature and emissivity law. With plausible choices for SED parameters (e.g., T=35 K, beta=1.0), the simulation predictions reproduce the observed number counts above ~ 1 mJy. The sources have a broad redshift distribution with median z ~ 2, in reasonable agreement with observational constraints. However, the predicted count distribution may be too steep at the faint end, and the fraction of low redshift objects may be larger than observed. In this physical model of the sub-mm galaxy population, the objects detected in existing surveys consist mainly of massive galaxies (several M_*) forming stars fairly steadily over timescales ~ 10^8-10^9 years, at moderate rates ~100 Msun/yr. The typical descendants of these sub-mm sources are even more massive galaxies, with old stellar populations, found primarily in dense environments. While the resolution of our simulations is not sufficient to determine galaxy morphologies, these properties support the proposed identification of sub-mm sources with massive ellipticals in the process of formation. The most robust and distinctive prediction of this model, stemming directly from the long timescale and correspondingly moderate rate of star formation, is that the far-IR SEDs of SCUBA sources have a relative high 850 micron luminosity for a given bolometric luminosity. [Abridged]

We use hydrodynamic simulations to examine how the baryonic components of galaxies are assembled, focusing on the relative importance of mergers and smooth accretion in the formation of ~L_* systems. In our primary simulation, which models a (50\hmpc)^3 comoving volume of a Lambda-dominated cold dark matter universe, the space density of objects at our (64-particle) baryon mass resolution threshold, M_c=5.4e10 M_sun, corresponds to that of observed galaxies with L~L_*/4. Galaxies above this threshold gain most of their mass by accretion rather than by mergers. At the redshift of peak mass growth, z~2, accretion dominates over merging by about 4:1. The mean accretion rate per galaxy declines from ~40 M_sun/yr at z=2 to ~10 M_sun/yr at z=0, while the merging rate peaks later (z~1) and declines more slowly, so by z=0 the ratio is about 2:1. We cannot distinguish truly smooth accretion from merging with objects below our mass resolution threshold, but extrapolating our measured mass spectrum of merging objects, dP/dM ~ M^a with a ~ -1, implies that sub-resolution mergers would add relatively little mass. The global star formation history in these simulations tracks the mass accretion rate rather than the merger rate. At low redshift, destruction of galaxies by mergers is approximately balanced by the growth of new systems, so the comoving space density of resolved galaxies stays nearly constant despite significant mass evolution at the galaxy-by-galaxy level. The predicted merger rate at z<~1 agrees with recent estimates from close pairs in the CFRS and CNOC2 redshift surveys.

We combine COBE-DMR measurements of cosmic microwave background anisotropy with a recent measurement of the mass power spectrum at redshift z=2.5 from Lya forest data to derive constraints on cosmological parameters and test the inflation+CDM scenario of structure formation. By treating the inflationary spectral index n as a free parameter, we can find successful fits to the COBE and Lya forest constraints in Omega_m=1 models with and without massive neutrinos and in low-Omega_m models with and without a cosmological constant. Within each class of model, the combination of COBE and the Lya forest P(k) constrains a parameter combination of the form (Omega_m h^a n^b Omega_b^c), with different indices for each case. This new constraint breaks some of the degeneracies in cosmological parameter determinations from other measurements. The Lya forest P(k) provides the first measurement of the slope of the linear mass power spectrum on ~Mpc scales, and it confirms a basic prediction of the inflationary CDM scenario: a nearly scale-invariant spectrum of primeval fluctuations (n~1) that bends towards k^{n-4} on small scales. Considering additional observational data, we find that COBE-normalized, Omega_m=1 models that match the Lya forest P(k) do not match the observed masses of rich galaxy clusters and that a low-Omega_m model with a cosmological constant provides the best overall fit, even without the direct evidence for cosmic acceleration from supernovae. Modest improvements in the Lya forest P(k) measurement could greatly restrict the allowable region of parameter space for CDM models, constrain the contribution of tensor fluctuations to CMB anisotropy, and achieve a more stringent test of the current consensus model of structure formation.

Studies of quasar absorption lines reveal that the low density intergalactic medium at z ~ 3 is enriched to 0.001-0.01 solar metallicity. This enrichment may have occurred in an early generation of Population III stars at redshift z >~ 10, by protogalaxies at 6 <~ z <~ 10, or by larger galaxies at 3 <~ z <~ 6. This paper addresses the third possibility by calculating the enrichment of the IGM at z >~ 3 by galaxies of baryonic mass >~ 10^8.5 solar masses. We use already completed cosmological simulations to which we add a prescription for chemical evolution and metal ejection by winds, assuming that the winds have properties similar to those observed in local starbursts and Lyman-break galaxies. Results are given for a number of representative models, and we also examine the properties of the galaxies responsible for the enrichment as well as the physical effects responsible for wind escape and propagation. We find that winds of velocity >~ 200-300 km/s are capable of enriching the IGM to the mean level observed, though many low-density regions would remain metal free. Calibrated by observations of Lyman-break galaxies, our calculations suggest that most galaxies at z >~ 3 should drive winds that can escape and propagate to large radii. The primary effect limiting the enrichment of low-density IG gas in our scenario is then the travel time from high- to low-density regions, implying that the metallicity of low-density gas is a strong function of redshift.

The presence of metals in hot cluster gas and in Ly-alpha absorbers, as well as the mass-metallicity relation of observed galaxies, suggest that galaxies lose a significant fraction of their metals to the intergalactic medium (IGM). Theoretical studies of this process have concentrated on metal removal by dynamical processes or supernova-driven winds. Here, we investigate the enrichment of the IGM by the expulsion of dust grains from galaxies by radiation pressure. We use already completed cosmological simulations, to which we add dust assuming that most dust can reach the equilibrium point between radiation pressure and gravitational forces. We find that the expulsion of dust and its subsequent (partial) destruction in the IGM can plausibly account for the observed level of C and Si enrichment of the z=3 IGM. At low-z, dust ejection and destruction could explain a substantial fraction of the metals in clusters, but it cannot account for all of the chemical species observed. Dust expelled by radiation pressure could give clusters a visual opacity of up to 0.2-0.5 mag in their central regions even after destruction by the hot intracluster medium; this value is interestingly close to limits and claimed observations of cluster extinction. We also comment on the implications of our results for the opacity of the general IGM. Finally, we suggest a possible `hybrid' scenario in which winds expel gas and dust into galaxy halos but radiation pressure distributes the dust uniformly through the IGM.

We use hydrodynamic cosmological simulations to study damped Lyman-alpha (DLA) and Lyman limit (LL) absorption at redshifts z=2-4 in five variants of the cold dark matter scenario. Our standard simulations resolve the formation of dense concentrations of neutral gas in halos with circular velocity v_c roughly 140 km/s for Omega_m=1 and 90 km/s for Omega_m=0.4, at z=2; an additional LCDM simulation resolves halos down to v_c approximately 50 km/s at z=3. We find a clear relation between HI column density and projected distance to the center of the nearest galaxy, with DLA absorption usually confined to galactocentric radii less than 10-15 kpc and LL absorption arising out to projected separations of 30 kpc or more. Detailed examination provides evidence of non-equilibrium effects on absorption cross-section. If we consider only absorption in the halos resolved by our standard simulations, then all five models fall short of reproducing the observed abundance of DLA and LL systems at these redshifts. If we extrapolate to lower halo masses, we find all four models are consistent with the observed abundance of DLA systems if the the extrapolated behavior extends to circular velocities roughly 50-80 km/s, and they may produce too much absorption if the relation continues to 40 km/s. Our results suggest that LL absorption is closely akin to DLA absorption, arising in less massive halos or at larger galactocentric radii but not caused by processes acting on a radically different mass scale.

(Abridged) We use a hydrodynamic simulation of a LambdaCDM model to predict the extragalactic X-ray background (XRB), focussing on emission from the intergalactic medium (IGM). We also include X-rays from point sources associated with galaxies in the simulation, and make maps of the angular distribution of the emission. We find that filaments in the maps are not evident, being diluted by projection. In the soft (0.5-2 keV) band, the mean intensity of radiation from intergalactic and cluster gas is 2.3*10^-12 ergdeg^-2cm^-2s^-1, 35% of the total soft band emission. This is compatible at the ~1 sigma level with estimates of the unresolved soft background from ROSAT and {\it Chandra}. Only 4% of the hard (2-10 keV) emission is associated with the IGM. Relative to AGN flux, the IGM component peaks at a lower redshift (median z~0.45) so its clustering makes an important contribution to that of the total XRB. The angular correlations on 0.1-10 arcmin scales are significant, with an amplitude roughly consistent with an extrapolation of recent ROSAT results to small scales. A cross-correlation of the XRB against nearby galaxies taken from a simulated redshift survey also yields a strong signal from the IGM. Although some recent papers have argued that the expected soft band intensity from gas in galaxy, group, and cluster halos would exceed XRB limits unless much of the gas is expelled by supernova feedback, we obtain reasonable compatibility with current observations in a simulation that incorporates cooling, star formation, and only modest feedback. A prediction of our model is that the unresolved portion of the soft XRB will remain mostly unresolved.

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 examine the cooling radiation from forming galaxies in hydrodynamic simulations of the LCDM model (cold dark matter with a cosmological constant), focusing on the Ly-alpha line luminosities of high-redshift systems. Primordial composition gas condenses within dark matter potential wells, forming objects with masses and sizes comparable to the luminous regions of observed galaxies. As expected, the energy radiated in this process is comparable to the gravitational binding energy of the baryons, and the total cooling luminosity of the galaxy population peaks at z ~= 2. However, in contrast to the classical picture of gas cooling from the \sim 10^6 K virial temperature of a typical dark matter halo, we find that most of the cooling radiation is emitted by gas with T < 20,000 K. As a consequence, roughly 50% of this cooling radiation emerges in the Ly-alpha line. While a galaxy's cooling luminosity is usually smaller than the ionizing continuum luminosity of its young stars, the two are comparable in the most massive systems, and the cooling radiation is produced at larger radii, where the Ly-alpha photons are less likely to be extinguished by dust. We suggest, in particular, that cooling radiation could explain the two large (\sim 100 kpc), luminous (L_{Ly-alpha} \sim 10^{44} erg s^{-1}) ``blobs'' of Ly-alpha emission found in Steidel et al.'s (1999) narrow band survey of a z = 3 proto-cluster. Our simulations predict objects of the observed luminosity at about the right space density, and radiative transfer effects can account for the observed sizes and line widths. We discuss observable tests of this hypothesis for the nature of the Ly-alpha blobs, and we present predictions for the contribution of cooling radiation to the Ly-alpha luminosity function of galaxies as a function of redshift.

We use hydrodynamic cosmological simulations to predict the star formation properties of high-redshift galaxies (z=2-6) in five variants of the inflationary cold dark matter scenario, paying particular attention to z=3, the redshift of the largest "Lyman-break galaxy" (LBG) samples. Because we link the star formation timescale to the local gas density, the rate at which a galaxy forms stars is governed mainly by the rate at which it accretes cooled gas from the surrounding medium. At z=3, star formation in most of the simulated galaxies is steady on 200 Myr timescales, and the instantaneous star formation rate (SFR) is correlated with total stellar mass. However, there is enough scatter in this correlation that a sample selected above a given SFR threshold may contain galaxies with a fairly wide range of masses. The redshift history and global density of star formation in the simulations depend mainly on the amplitude of mass fluctuations in the underlying cosmological model. The three models whose mass fluctuation amplitudes agree with recent analyses of the Lyman-alpha forest also reproduce the observed luminosity function of LBGs reasonably well, though the dynamic range of the comparison is small and the theoretical and observational uncertainties are large. The models with higher and lower amplitudes appear to predict too much and too little star formation, respectively, though they are not clearly ruled out. The intermediate amplitude models predict SFR ~ 30-40 Msun/yr for galaxies with a surface density 1 per arcmin^2 per unit redshift at z=3. They predict much higher surface densities at lower SFR, and significant numbers of galaxies with SFR > 10 Msun/yr at z >= 5.

Using a high-resolution cosmological hydrodynamic simulation, we present a method to constrain extinction due to intergalactic grey dust based on the observed magnitudes of distant Type IA supernovae. We apply several simple prescriptions to relate the intergalactic dust density to the gas density in the simulation, thereby obtaining dust extinctions that may be directly compared to the observed distribution of supernova magnitudes. Our analysis is sensitive to the spatial distribution of grey dust, but is not dependent on its intrinsic properties such as its opacity or grain size. We present an application of our technique to the supernova data of Perlmutter et al., who find that their high redshift sample is ~0.2 magnitudes fainter than the expectation for a non-accelerating, low-density universe. We find that for grey dust to be responsible, it must be distributed quite smoothly, e.g., tracing intergalactic gas. More realistic dust distributions, such as dust tracing the metal density, are inconsistent with observations at the 1.5-2 sigma level. Upcoming observations and improved modelling of the dust distribution should lead to stronger constraints on intergalactic grey dust extinction.

We present the analysis of a suite of simulations of a Virgo mass galaxy cluster. Undertaken within the framework of standard cold dark matter cosmology, these simulations were performed at differing resolutions and with increasingly complex physical processes, with the goal of identifying the effects of each on the evolution of the cluster. We focus on the cluster at the present epoch and examine properties including the radial distributions of density, temperature, entropy and velocity. We also map `observable' projected properties such as the surface mass density, X-ray surface brightness and SZ signature. We identify significant differences between the simulations, which highlights the need for caution when comparing numerical simulations to observations of galaxy clusters. While resolution affects the inner density profile in dark matter simulations, the addition of a gaseous component, especially one that cools and forms stars, affects the entire cluster. We conclude that both resolution and included physical processes play an important role in simulating the formation and evolution of galaxy clusters. Therefore, physical inferences drawn from simulations that do not include a gaseous component that can cool and form stars present a poor representation of reality. (Abridged)

What type of objects are being detected as $z\sim 3$ "Lyman break galaxies"? Are they predominantly the most massive galaxies at that epoch, or are many of them smaller galaxies undergoing a short-lived burst of merger-induced star formation? We attempt to address this question using high-resolution cosmological hydrodynamic simulations including star formation and feedback. Our $\Lambda$CDM simulation, together with Bruzual-Charlot population synthesis models, reproduces the observed number density and luminosity function of Lyman break galaxies when dust is incorporated. The inclusion of dust is crucial for this agreement. In our simulation, these galaxies are predominantly the most massive objects at this epoch, and have a significant population of older stars. Nevertheless, it is possible that our simulations lack the resolution and requisite physics to produce starbursts, despite having a physical resolution of $\la 700$ pc at z=3. Thus we cannot rule out merger-induced starburst galaxies also contributing to the observed population of high-redshift objects.

We determine the evolution from $z=3\to 0$ of the galaxy and mass correlation functions and bias factor in a $50h^{-1}$Mpc $\Lambda$CDM hydrodynamic simulation with $10h^{-1}$kpc resolution. The mass correlation function grows with time, but the galaxy correlation function shows little evolution and is well described by a power law. At early times, galaxies are biased traces of mass, with bias being higher on smaller scales. By z=0, galaxies trace the mass, and the bias shows little scale dependence.

We review theoretical approaches to the study of galaxy formation, with emphasis on the role of hydrodynamic simulations in modeling the high redshift galaxy population. We present new predictions for the abundance of star-forming galaxies in the Lambda + cold dark matter model (Omega_m=0.4, Omega_L=0.6), combining results from several simulations to probe a wide range of redshift. At a threshold density of one object per arcmin^2 per unit z, these simulations predict galaxies with star formation rates of 2 msun/yr (z=10), 5 msun/yr (z=8), 20 msun/yr (z=6), 70-100 msun/yr (z=4-2), and 30 msun/yr (z=0.5). For galaxies selected at a fixed comoving space density n=0.003 h^3 Mpc^{-3], a (50 Mpc/h)^3 simulation predicts a galaxy correlation function (r/5 Mpc/h)^{-1.8} in comoving coordinates, essentially independent of redshift from z=4 to z=0.5. Different cosmological models predict global histories of star formation that reflect their overall histories of mass clustering, but robust numerical predictions of the comoving space density of star formation are difficult because the simulations miss the contribution from galaxies below their resolution limit. The LCDM model appears to predict a star formation history with roughly the shape inferred from observations, but it produces too many stars at low redshift, predicting Omega_* ~ 0.015 at z=0. We conclude with a brief discussion of this discrepancy and three others that suggest gaps in our current theory of galaxy formation: small disks, steep central halo profiles, and an excess of low mass dark halos. While these problems could fade as the simulations or observations improve, they could also guide us towards a new understanding of galactic scale star formation, the spectrum of primordial fluctuations, or the nature of dark matter.

Recent experimental results about the formation of molecular hydrogen on astrophysically relevant surfaces under conditions close to those encountered in the interstellar medium are analyzed using rate equations. The parameters of the rate equation model are fitted to temperature-programmed desorption curves obtained in the laboratory. These parameters are the activation energy barriers for atomic hydrogen diffusion and desorption, the barrier for molecular hydrogen desorption, and the probability of spontaneous desorption of a hydrogen molecule upon recombination. The model is a generalization of the Polanyi-Wigner equation and provides a description of both first and second order kinetic processes within a single model. Using the values of the parameters that fit best the experimental results, the efficiency of hydrogen recombination on olivine and amorphous carbon surfaces is obtained for a range of hydrogen flux and surface temperature pertinent to a wide range of interstellar conditions.

We analyse parallel N-body simulations of three Cold Dark Matter (CDM) universes to study the abundance and clustering of galaxy clusters. The simulations cover a volume comparable to the forthcoming SDSS. We are able to make robust measurements of cluster properties to a redshift larger than unity. We extract halos using two independent, public domain group finders (FOF & HOP) and find consistent results. The correlation function of clusters is in very good agreement with a simple analytic prescription based upon a Lagrangian biasing scheme developed by Mo & White (1996) and the Press-Schechter (PS) formalism for the mass function. The R_0--D_c relation for the open CDM model is in good agreement with the results from the APM Cluster Survey. The SCDM universe shows a robust deviation in the shape and evolution of the mass function when compared with that predicted by the PS formalism. Critical models with a low sigma_8 normalization or small shape parameter Gamma show an excess of massive clusters compared with the PS prediction. When cluster normalized, the SCDM universe at z =1 contains 10 times more clusters with temperatures greater than 7keV, compared with the PS prediction. The agreement between the analytic and N-body mass functions of SCDM can be improved if the value of the delta_c (the extrapolated linear theory threshold for collapse) is revised to be $ delta_c(z) = 1.685[(0.7/sigma_8)(1+z)]^{-0.125}. Our best estimate for the amplitude of fluctuations inferred from the local cluster abundance for SCDM is sigma_{8} = 0.5 \pm 0.04. However, the discrepancy between the temperature function predicted in a critical density universe and that observed at z=0.33 (Henry et al. 1998) remains. (abridged)

We outline the physical picture of the high-redshift Ly-alpha forest that has emerged from cosmological simulations, describe statistical characteristics of the forest that can be used to test theories of structure formation, present a preliminary comparison between simulation results and measurements from Keck HIRES spectra, and discuss a recent determination of the slope and amplitude of the linear mass power spectrum P(k) at z=2.5 from moderate resolution spectra. The physical picture is simple if each QSO spectrum is viewed as a continuous non-linear map of the line-of-sight density field rather than a collection of discrete absorption lines. The distribution of flux decrements depends mainly on the amplitude and PDF (Gaussian vs. non-Gaussian) of the primordial density fluctuations. The threshold crossing frequency, analogous to the 3-d genus curve, responds to the shape and amplitude of P(k) and to the values of Omega and Lambda. Open and Lambda-CDM models agree well with the measured flux decrement distribution at smoothing lengths of 25 km/s and 100 km/s and with the threshold crossing frequency at 100 km/s. Discrepancy with the observed threshold crossing frequency at 25 km/s may reflect the combined effects of noise in the data and limited mass resolution of the simulations. The slope of the measured P(k) agrees with the predictions of inflation+CDM models. Combining the amplitude with COBE normalization imposes a constraint on these models of the form Omega h^x n^y Omega_b^z = constant. Assuming Gaussian primordial fluctuations and a power spectrum shape parameter Gamma~0.2, consistency of the measured P(k) with the observed cluster mass function at z=0 requires Omega=0.46^{+0.12}_{-0.10} for an open universe and Omega=0.34^{+0.13}_{-0.09} for a flat universe (1-sigma errors). (Shortened)

We estimate the value of the matter density parameter \Omega_0 by combining constraints from the galaxy cluster mass function with Croft et al.'s recent measurement of the mass power spectrum, P(k), from \lya forest data. The key assumption of the method is that cosmic structure formed by gravitational instability from Gaussian primordial fluctuations. For a specified value of \Omega_0, matching the observed cluster mass function then fixes the value of \sigma_8, the rms amplitude of mass fluctuations in 8\hmpc spheres, and it thus determines the normalization of P(k) at z=0. The value of \Omega_0 also determines the ratio of P(k) at z=0 to P(k) at z=2.5, the central redshift of the \lya forest data; the ratio is different for an open universe (\Lambda=0) or a flat universe. Because the \lya forest measurement only reaches comoving scales 2\pi/k ~ 15-20\hmpc, the derived value of \Omega_0 depends on the value of the power spectrum shape parameter \Gamma, which determines the relative contribution of larger scale modes to \sigma_8. Adopting \Gamma=0.2, a value favored by galaxy clustering data, we find \Omega_0 = 0.46^{+0.12}_{-0.10} for an open universe and \Omega_0=0.34^{+0.13}_{-0.09} for a flat universe (1\sigma errors, not including the uncertainty in cluster normalization). Cluster-normalized models with \Omega_0=1 predict too low an amplitude for P(k) at z=2.5, while models with \Omega_0=0.1 predict too high an amplitude. The more general best fit parameter combination is approximately \Omega_0 + 0.2\Lambda_0 = 0.46 + 1.3(\Gamma-0.2). Analysis of larger, existing samples of QSO spectra could greatly improve the measurement of P(k) from the \lya forest, allowing a determination of \Omega_0 by this method with a precision of ~15%, limited mainly by uncertainty in the cluster mass function.

We measure the linear power spectrum of mass density fluctuations at redshift z=2.5 from the \lya forest absorption in a sample of 19 QSO spectra, using the method introduced by Croft et al. (1998). The P(k) measurement covers the range 2\pi/k ~ 450-2350 km/s (2-12 comoving \hmpc for \Omega=1). We examine a number of possible sources of systematic error and find none that are significant on these scales. In particular, we show that spatial variations in the UV background caused by the discreteness of the source population should have negligible effect on our P(k) measurement. We obtain consistent results from the high and low redshift halves of the data set and from an entirely independent sample of nine QSO spectra with mean redshift z=2.1. A power law fit to our measured P(k) yields a logarithmic slope n=-2.25 +/- 0.18 and an amplitude \Delta^2(k_p) = 0.57^{+0.26}_{-0.18}, where $\Delta^2$ is the contribution to the density variance from a unit interval of lnk and k_p=0.008 (km/s)^{-1}. Direct comparison of our mass P(k) to the measured clustering of Lyman Break Galaxies shows that they are a highly biased population, with a bias factor b~2-5. The slope of the linear P(k), never previously measured on these scales, is close to that predicted by models based on inflation and Cold Dark Matter (CDM). The P(k) amplitude is consistent with some scale-invariant, COBE-normalized CDM models (e.g., an open model with \Omega_0=0.4) and inconsistent with others (e.g., \Omega=1). Even with limited dynamic range and substantial statistical uncertainty, a measurement of P(k) that has no unknown ``bias factors'' offers many opportunities for testing theories of structure formation and constraining cosmological parameters. (Shortened)

We study the physical origin of the low-redshift Lyman alpha forest in hydrodynamic simulations of four CDM cosmologies. Our main conclusions are insensitive to the cosmological model but depend on our assumption that the UV background declines at low redshift. We find that the expansion of the universe drives rapid evolution of dN/dz (the number of absorbers per unit z) at z > 1.7, but that at lower redshift the fading of the UV background counters the influence of expansion, leading to slow evolution. At every redshift, weaker lines come primarily from moderate fluctuations of the diffuse, unshocked IGM, and stronger lines originate in shocked or radiatively cooled gas of higher overdensity. However, the neutral hydrogen column density associated with structures of fixed overdensity drops as the universe expands, so an absorber at z = 0 is dynamically analogous to an absorber with neutral hydrogen column density 10 to 50 times higher at z = 2-3. We find no clear distinction between lines arising in "galaxy halos" and lines arising in larger scale structures; however, galaxies tend to lie near the dense regions of the IGM that produce strong Lyman alpha lines. The simulations provide a unified physical picture that accounts for the most distinctive observed properties of the low redshift Lyman alpha forest: (1) a sharp transition in the evolution of dN/dz at z ~ 1.7, (2) stronger evolution for absorbers of higher equivalent width, (3) a correlation of increasing Lyman alpha equivalent width with decreasing galaxy impact parameter, and (4) a tendency for stronger lines to arise in close proximity to galaxies while weaker lines trace more diffuse large scale structure. (Abridged)

An analysis of the kinetics of H2 formation on interstellar dust grains is presented using rate equations. It is shown that semi-empirical expressions that appeared in the literature represent two different physical regimes. In particular, it is shown that the expression given by Hollenbach, Werner and Salpeter [ApJ, 163, 165 (1971)] applies when high flux, or high mobility, of H atoms on the surface of a grain, makes it very unlikely that H atoms evaporate before they meet each other and recombine. The expression of Pirronello et al.\ [ApJ, 483, L131 (1997)] -- deduced on the basis of accurate measurements on realistic dust analogue -- applies to the opposite regime (low coverage and low mobility). The implications of this analysis for the understanding of the processes dominating in the Interstellar Medium are discussed.

We investigate the clustering of high redshift galaxies in five variants of the cold dark matter (CDM) scenario, using hydrodynamic cosmological simulations that resolve the formation of systems with circular velocities $v_c \geq 100 \kms$ ($\Omega=1$) or $v_c \geq 70 \kms$ ($\Omega=0.4$). Although the five models differ in their cosmological parameters and in the shapes and amplitudes of their mass power spectra, they predict remarkably similar galaxy clustering at $z=2$, 3, and 4. The galaxy correlation functions show almost no evolution over this redshift range, even though the mass correlation functions grow steadily in time. Despite the fairly low circular velocity threshold of the simulations, the high redshift galaxies are usually highly biased tracers of the underlying mass distribution; the bias factor evolves with redshift and varies from model to model. Predicted correlation lengths for the resolved galaxy population are $2-3\hmpc$ (comoving) at $z=3$. More massive galaxies tend to be more strongly clustered. These CDM models have no difficulty in explaining the strong observed clustering of Lyman-break galaxies, and some may even predict excessive clustering. Because the effects of bias obscure differences in mass clustering, it appears that Lyman-break galaxy clustering will not be a good test of cosmological models but will instead provide a tool for constraining the physics of galaxy formation.

We analyze hierarchical structure formation based on scale-free initial conditions in an Einstein-de Sitter universe, including a baryonic component. We present three independent, smoothed particle hydrodynamics (SPH) simulations, performed with two different SPH codes (TreeSPH and P3MSPH) at two resolutions. Each simulation is based upon identical initial conditions, which consist of Gaussian distributed initial density fluctuations that have an n=-1 power spectrum. The baryonic material is modeled as an ideal gas subject only to shock heating and adiabatic heating and cooling. The evolution is expected to be self-similar in time, and under certain restrictions we identify the expected scalings for many properties of the distribution of collapsed objects in all three realizations. The distributions of dark matter masses, baryon masses, and mass and emission weighted temperatures scale quite reliably. However, the density estimates in the central regions of these structures are determined by the degree of numerical resolution. As a result, mean gas densities and luminosities obey the expected scalings only when calculated within a limited dynamic range in density contrast. The temperatures and luminosities of the groups show tight correlations with the baryon masses, which can be well-represented by power-laws. The Press-Schechter (PS) approximation predicts the distribution of group dark matter and baryon masses fairly well, though it tends to overestimate the baryon masses. Combining the PS mass distribution with the measured relations for T(M) and L(M) predicts the temperature and luminosity distributions reasonably, though there are some discrepancies at high temperatures/luminosities. The three simulations agree well for the properties of groups that are resolved by 32 or more particles.

We search for OVI absorption in a Keck HIRES spectrum of the z=3.62 quasar Q1422+231. Comparison of CIV measurements to cosmological simulations shows that \lya forest absorbers with N_HI > 10^{14.5} have [C/H]~=-2.5, for the UV background spectrum of Haardt & Madau (HM). Lower column density absorption arises in lower density gas, where OVI is the most sensitive metal tracer. Since OVI lines lie at wavelengths contaminated by Lyman series absorption, we interpret our Q1422 results by comparing to artificial spectra drawn from an SPH simulation of a Lambda-dominated CDM model. A search for deep, narrow features in Q1422 yields only a few candidate OVI lines, statistically consistent with the number in artificial spectra with no metals; spectra generated with the HM background and [O/H] >= -2.5 predict too many narrow lines. However, applying the optical depth ratio technique of Songaila (1998), we DO find significant OVI associated with CIV systems; matching Q1422 requires [O/C]~=+0.5, implying [O/H]~=-2.0. Taken together these results imply that (a) the metallicity in the low density IGM is at least a factor of three below that in the overdense regions where CIV absorption is detectable, and (b) oxygen is overabundant in these regions, consistent with the enrichment pattern of old halo stars. If the UV background is heavily truncated above 4 Ry, an implausibly high oxygen overabundance ([O/C]>+2) is required by the data; thus a majority of the volume of the universe must have undergone helium reionization by z=3.(Abridged)

A convincing physical picture for the Lyman-alpha forest has emerged from simulations and related semi-analytic studies of structure formation models. Observations can be be used in the context of this picture to study cosmology using the structure of the forest. With the availability of well motivated predictions, not only has it become possible to test models directly, but the physical processes involved appear to be simple enough that we can attempt to reconstruct aspects of the underlying cosmology from observations. We briefly summarise the method of Croft et al (1997) (astro-ph/9708018) for recovering the primordial mass power spectrum from Lyman-alpha forest data, emphasising the physical reasons that the derived P(k) is independent of unknown "bias factors". We present an illustrative application of the method to four quasar "spectra" reconstructed from published line lists.

We discuss simulations of gas at high redshift. We briefly review the methods used and the results for quasar absorption lines. We present gas mass functions and galaxy correlation functions for 5 different cosmological models. Galaxies should be detectable at redshifts greater than 2 by SKAI, and measurements of the gas mass functions and galaxy correlation functions could be used to discriminate between different cosmological models.

We develop a prescription for characterizing the strengths of metal lines associated with Lyman-alpha forest absorbers (LYFAs) of a given neutral hydrogen column density N_HI and metallicity [Fe/H]. This Line Observability Index (LOX) is line-specific and translates, for weak lines, into a measure of the equivalent width. It can be evaluated quickly for thousands of transitions within the framework of a given model of the Lyman-alpha forest, providing a ranking of the lines in terms of their strengths and enabling model builders to select the lines that should be detectable in observed spectra of a given resolution and signal-to-noise ratio. We compute the LOX for a large number of elements and transitions in two cosmological models of the Lyman-alpha forest at z=3 derived from a hydrodynamic simulation of structure formation, and we discuss how the LOX depends on redshift and on model parameters such as the mean baryonic density and radiation field. We find that the OVI (1032,1038) doublet is the best probe of the metallicity in low column density LYFAs N_{HI} \approx 10^{14.5} cm^{-2}). Metallicities down to [O/H] \sim -3 ([Fe/H] \sim -3.5 with the assumed [O/Fe] ratio) yield OVI absorption features that should be detectable in current high-quality spectra, provided that the expected position of the OVI feature is not contaminated by HI absorption. The strongest transitions in lower ionisation states of oxygen are OV(630), OIV(788), and OIII(833), and are likely to be detected with next generation UV instruments. Of the lines with rest wavelengths \lambda_r > 1216, which can potentially be observed redwards of the \lya forest, the CIV(1548,1551) doublet is expected to dominate in all LYFAs, regardless of the value of N_HI.

We have simulated disk galaxies undergoing continual bombardment by other galaxies in a rich cluster. "Galaxy harassment" leads to dramatic evolution of smaller disk galaxies and provides an extremely effective mechanism to fuel a central quasar. Within a few billion years after a small disk galaxy enters the cluster environment, up to 90% of its gas can be driven into the inner 500 pc. Up to half of the mass can be transferred in a burst lasting just 100-200 Myr. This transport of gas to the center of galaxy is far more efficient than any mechanism proposed before. Galaxy harassment was first proposed to explain the disturbed blue galaxies in clusters seen in clusters at ($z \gsim 0.3$), the "Butcher-Oemler effect". Quasars at the same reshifts lie in more clustered environments than those at lower redshift. Recent HST observations find that roughly half of all observed quasar host galaxiess are fainter than \l*, with many of these less luminous hosts occuring at redshifts $z \gsim 0.3$. We examine 5 quasars that are claimed to have low luminosity hosts and find that 3 are in rich clusters of galaxies, the fourth may be in a cluster but the evidence for this is marginal. The environment of the fifth has not been studied.

Galaxy morphologies in clusters have undergone a remarkable transition over the past several billion years. Distant clusters at $z \sim 0.4$ are filled with small spiral galaxies, many of which are disturbed and show evidence of multiple bursts of star-formation. This population is absent from nearby clusters where spheroidals comprise the faint end of the luminosity function. Our numerical simulations follow the evolution of disk galaxies in a rich cluster owing to encounters with brighter galaxies and the cluster's tidal field---galaxy harassment. After a bursting transient phase, they undergo a complete morphological transformation from "disks" to"spheroidals". We examine the remnants and find support for our theory in detailed comparisons of the photometry and kinematics of the spheroidal galaxies in clusters. Our model naturally accounts for the intermediate age stellar population seen in these spheroidals as well as the trend in dwarf to giant ratio with cluster richness. The final shapes are typically prolate and are flattened primarily by velocity anisotropy. Their mass to light ratios are in the range 3---8 in good agreement with observations.

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.

Leapfrog integration has been the method of choice in N-body simulations owing to its low computational cost for a symplectic integrator with second order accuracy. We introduce a new leapfrog integrator that allows for variable timesteps for each particle in large N-body simulations. Tests with single particles in fixed potentials show that it behaves as a symplectic integrator. We then examine the results of both standard leapfrog and our temporally adaptive leapfrog on full N-body integrations of clusters and large scale structure establishing accuracy criteria for both methods. The adaptive method shows significant speed-ups over single step integrations---but the integrator no longer appears to be symplectic or, in the case of large scale structure simulations, accurate. This loss of accuracy appears to be caused by the way that the timestep is chosen, not by the integrator itself. We present a related integration technique that does retain sufficient accuracy. Although it is not symplectic, it is apparently better than previous implementations and is our current integrator of choice for large astrophysical simulations. We also note that the standard leapfrog difference equations used in cosmological N-body integrations in comoving coordinates are not symplectic. We derive an implementation of leapfrog that is in comoving canonical coordinates to correct for this deficiency.

Lyman-alpha forest lines in QSO spectra have typical widths of 20-50 km/s. Low column density absorbers in cosmological simulations are large, diffuse structures, and the Hubble flow across the spatially extended absorber is usually the dominant contribution to the width of its associated absorption line. Thermal broadening is unimportant over most of the spectrum, and peculiar velocities tend to make absorption features narrower rather than broader. As a consequence of Hubble flow broadening, there is a close relation between local Lyman-alpha optical depth and local neutral hydrogen density, which is well approximated by the Gunn-Peterson formula. The physics that governs the unshocked intergalactic medium leads to a tight correlation between the neutral hydrogen density and the underlying gas and dark matter overdensity. For many purposes, it is simpler to regard a Lyman-alpha forest spectrum as a continuous, non-linear map of the density field rather than a collection of discrete lines. This continuous field view of the Lyman-alpha forest can be applied to measurement of the baryon density parameter, testing of cosmological models, and robust determination of the shape and amplitude of the primordial mass power spectrum.

The metal absorption lines found in association with \lya absorbers of moderate to high HI column density contain valuable information about the metallicity and ionization conditions within the absorbers and offer a stronger test of models of the intergalactic medium at \sim 3$ than HI absorption lines alone. We have developed a method to predict the strengths of metal absorption lines within the framework of cosmological models for the \lya forest. The method consists of evaluating a quantity, the Line Observability Index, for a database of hundreds of candidate metal lines, allowing a comprehensive identification of the lines the model predicts to be detectable associated with a \lya absorber of a given HI column density and metallicity. Applying this technique to a particular class of models at \sim 2-4$, we predict that the OVI(1032 \AA, 1038 \AA) doublet is the only practical probe of the metallicity of low column density absorbers (HI \simlt 10^{14.5} cm^{-2}$), that CIV (1548 \AA) is the strongest line with rest wavelength $\lambda_r > 1216$ \AA{} regardless of HI$, and that the strongest metal lines should be CIII(977 \AA) and SiIII(1206.5 \AA), which peak at HI \sim 10^{17} cm^{-2}$.

We describe cosmological simulation techniques and their application to studies of cosmic structure formation, with particular attention to recent hydrodynamic simulations of structure in the high redshift universe. Collisionless N-body simulations with Gaussian initial conditions produce a pattern of sheets, filaments, tunnels, and voids that resembles the observed large scale galaxy distribution. Simulations that incorporate gas dynamics and dissipation form dense clumps of cold gas with sizes and masses similar to the luminous parts of galaxies. Models based on inflation and cold dark matter predict a healthy population of high redshift galaxies, including systems with star formation rates of 20 M_{\sun}/year at z=6. At z~3, most of the baryons in these models reside in the low density intergalactic medium, which produces fluctuating Lyman-alpha absorption in the spectra of background quasars. The physical description of this ``Lyman-alpha forest'' is particularly simple if the absorption spectrum is viewed as a 1-dimensional map of a continuous medium instead of a collection of lines. The combination of superb observational data and robust numerical predictions makes the Lyman-alpha forest a promising tool for testing cosmological models.

We present a method to recover the shape and amplitude of the power spectrum of mass fluctuations, P(k), from observations of the high redshift \lya forest. The method is motivated by the physical picture of the \lya forest that has emerged from hydrodynamic cosmological simulations and related semi-analytic models, which predicts a tight correlation between the \lya optical depth and the underlying matter density. We monotonically map the QSO spectrum to a Gaussian density field, measure its 3-d P(k), and normalize by evolving cosmological simulations with this P(k) until they reproduce the observed power spectrum of the QSO flux. Imposing the observed mean \lya opacity as a constraint makes the derived P(k) normalization insensitive to the choice of cosmological parameters, ionizing background spectrum, or reionization history. Thus, in contrast to estimates of P(k) from galaxy clustering, there are no uncertain "bias parameters" in the recovery of the mass power spectrum. We test the full procedure on SPH simulations of 3 cosmological models and show that it recovers their true mass power spectra on comoving scales ~1-10/h Mpc, the upper scale being set by the size of the simulation boxes. The procedure works even for noisy (S/N ~ 10), moderate resolution (~40 km/s pixels) spectra. We present an illustrative application to Q1422+231; the recovered P(k) is consistent with an \Omega=1, \sigma_8=0.5 CDM model. Application to large QSO samples should yield the power spectrum of mass fluctuations on small scales at z ~ 2-4. (Compressed)

Observations from the HUT and the HST have recently detected HeII absorption along the lines of sight to two high redshift quasars. We use cosmological simulations with gas dynamics to investigate HeII absorption in the cold dark matter (CDM) theory of structure formation. We consider two Omega=1 CDM models with different normalizations and one Omega_0=0.4 CDM model, all incorporating the photoionizing UV background spectrum computed by Haardt & Madau (1996). The simulated gas distribution, combined with the H&M spectral shape, accounts for the relative observed values of taubar_HI and taubar_HeII, the effective mean optical depths for HI and HeII absorption. If the background intensity is as high as H&M predict, then matching the absolute values of taubar_HI and taubar_HeII requires a baryon abundance larger (by factors between 1.5 and 3 for the various CDM models) than our assumed value of Omega_b h^2=0.0125. The simulations reproduce the evolution of taubar_heII over the observed redshift range, 2.2 < z < 3.3, if the HeII photoionization rate remains roughly constant. HeII absorption in the CDM simulations is produced by a diffuse, fluctuating, intergalactic medium, which also gives rise to the HI ly-alpha forest. Much of the HeII opacity arises in underdense regions where the HI optical depth is very low. We compute statistical properties of the HeII and HI absorption that can be used to test the CDM models and distinguish them from an alternative scenario in which the HeII absorption is caused by discrete, compact clouds. The CDM scenario predicts that a substantial amount of baryonic material resides in underdense regions at high redshift. HeII absorption is the only sensitive probe of such extremely diffuse, intergalactic gas, so it can provide a vital test of this fundamental prediction.

We calculate the number of damped Lyman-alpha absorbers expected in various popular cosmological models as a function of redshift and compare our predictions with observed abundances. The Press-Schechter formalism is used to obtain the distribution of halos with circular velocity in different cosmologies, and we calibrate the relation between circular velocity and absorption cross-section using detailed gas dynamical simulations of a ``standard'' cold dark matter (CDM) model. Because of this calibration, our approach makes more realistic assumptions about the absorption properties of collapsed objects than previous, analytic calculations of the damped Lyman-alpha abundance. CDM models with Omega_0=1, H_0=50, baryon density Omega_b=0.05, and scale-invariant primeval fluctuations reproduce the observed incidence and redshift evolution of damped Lyman-alpha absorption to within observational uncertainty, for both COBE normalization (sigma_8=1.2) and a lower normalization (sigma_8=0.7) that better matches the observed cluster abundance at z=0. A tilted (n=0.8, sigma_8=0.7) CDM model tends to underproduce absorption, especially at z=4. With COBE normalization, a CDM model with Omega_0=0.4, Omega_{Lambda}=0.6 gives an acceptable fit to the observed absorption; an open CDM model is marginally acceptable if Omega_0 is at least 0.4 and strongly inconsistent with the z=4 data if Omega_0=0.3. Mixed dark matter models tend not to produce sufficient absorption, being roughly comparable to tilted CDM models if Omega_{nu} = 0.2 and failing drastically if Omega_{nu} = 0.3.

Any viable cosmological model must produce enough structure at early epochs to explain the amount of gas associated with high-redshift damped Ly$\alpha$ systems. We study the evolution of damped Ly$\alpha$ systems at redshifts $z\ge 2$ in cold dark matter (CDM) and cold+hot dark matter (CDM+HDM) models using both N-body and hydrodynamic simulations. Our approach incorporates the effects of gas dynamics, and we find that all earlier estimates which assumed that all the baryons in dark matter halos would contribute to damped Ly$\alpha$ absorption have overestimated the column density distribution $f(N)$ and the fraction of neutral dense gas $\Omega_g$ in damped Ly$\alpha$ systems. The differences are driven by ionization of hydrogen in the outskirts of galactic halos and by gaseous dissipation near the halo centers, and they tend to exacerbate the problem of late galaxy formation in CDM+HDM models. We only include systems up to the highest observed column density $N\sim 10^{21.8}$ cm$^{-2}$ in the estimation of $\Omega_g$ for a fair comparison with data. If the observed $f(N)$ and $\Omega_g$ inferred from a small number of confirmed and candidate absorbers are robust, the amount of gas in damped Ly$\alpha$ systems at high redshifts in the $\Omega_\nu=0.2$ CDM+HDM model falls well below the observations.

In gravitational instability models, \lya absorption arises from a continuous fluctuating medium, so that spectra provide a non-linear one-dimensional ``map'' of the underlying density field. We characterise this continuous absorption using statistical measures applied to the distribution of absorbed flux. We describe two simple members of a family of statistics which we apply to simulated spectra in order to show their sensitivity as probes of cosmological parameters (H$_{0}$, $\Omega$, the initial power spectrum of matter fluctuations) and the physical state of the IGM. We make use of SPH simulation results to test the flux statistics, as well as presenting a preliminary application to Keck HIRES data.

We study metal line absorption of CIV, CII, SiIV, and NV at redshifts z=3.5 to z=2 within the framework of a cosmological model for the Lyman alpha forest, comparing the results of numerical simulations to recent observations by Songaila & Cowie (1996, SC). In agreement with Rauch, Haehnelt & Steinmetz (1996), we find that the observed mean value of the CIV/HI ratio around z=3 is reproduced if a uniform metallicity of [C/H] = -2.5 is assumed in our model, but that the observed scatter in this ratio is larger than predicted, implying a scatter in the metallicity of the absorbing systems of roughly an order of magnitude. We further argue that absorbers with relatively low column densities (log(N_HI) < 15) likely have a mean metallicity [C/H] less than -2.5, a result which is basically independent of the model considered, and similar to what is predicted in the simulations by Gnedin & Ostriker (1997). Our model predicts no significant evolution in the mean values of metal line column densities between z=3.5 and z=2. Comparison of the predicted and observed numbers of SiIV and NV systems suggests that the photoionizing background radiation field around z=3 is somewhat softer than that proposed by Haardt & Madau (1996). Our model does not account for the increase in the SiIV/CIV ratio at z=3.2 found by SC. While SC suggested that the increase could be explained by a softening of the radiation field due to HeII absorption at z > 3, such a modification does not raise the mean value of SiIV/CIV in our simulation because it shifts numerous weak SiIV systems to just above the detection limit, thus keeping the mean column density of observable SiIV systems low.

We derive analytic lower bounds on the cosmic baryon density by requiring that the high-redshift IGM contain enough neutral hydrogen to produce the observed \lya absorption in quasar spectra. The key theoretical assumption is that absorbing structures are no more extended in redshift space than in real space. This assumption might not hold if \lya clouds are highly overdense and thermally broadened, but it is likely to hold in the gravitational instability picture for the \lya forest suggested by cosmological simulations, independently of the details of the cosmological model. The other ingredients in these bounds are an estimate of the UV background from quasars, a temperature $T~10^4K$ for the "warm" photoionized IGM that produces most of the \lya absorption, a value of the Hubble constant, and observational estimates of the mean \lya flux decrement $\Dbar$ or, for a more restrictive bound, the distribution function $P(\tau)$ of \lya optical depths. With plausible parameter values, the mean decrement bound implies a baryon density parameter $\Omb \ga 0.0125/h^2$. With conservative values, the bound weakens to $\Omb \ga 0.005/h^2$, but the required clustering of the IGM is then incompatible with other properties of quasar spectra. A recent observational determination of $P(\tau)$ implies $\Omb \ga 0.0125/h^2$ even for a conservative estimate of the UV background, and $\Omb \ga 0.018/h^2$ for a more reasonable estimate. These bounds are consistent with recent low estimates of the primordial deuterium-to-hydrogen ratio $\dtoh$, which imply $\Omb \approx 0.025/h^2$, but the $P(\tau)$ bound can only be reconciled with high $\dtoh$ estimates by abandoning standard big bang nucleosynthesis or the gravitational instability picture for the origin of the \lya forest. (Shortened abstract.)

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.

We use an automated Voigt-profile fitting procedure to extract statistical properties of the Ly$\alpha$ forest in a numerical simulation of an $\Omega=1$, cold dark matter (CDM) universe. Our analysis method is similar to that used in most observational studies of the forest, and we compare the simulations to recently published results derived from Keck HIRES spectra. With the Voigt-profile decomposition analysis, the simulation reproduces the large number of weak lines ($N_{\rm HI}\la 10^{13}\cdunits$) found in the HIRES spectra. The column density distribution evolves significantly between $z=3$ and $z=2$, with the number of lines at fixed column density dropping by a factor $\sim 1.6$ in the range where line blending is not severe. At $z=3$, the $b$-parameter distribution has a median of $35 \kms$ and a dispersion of $20 \kms$, in reasonable agreement with the observed values. The comparison between our new analysis and recent data strengthens earlier claims that the \lya forest arises naturally in hierarchical structure formation as photoionized gas falls into dark matter potential wells. However, there are two statistically signficant discrepancies between the simulated forest and the HIRES results: the model produces too many lines at $z=3$ by a factor $\sim 1.5-2$, and it produces more narrow lines ($b<20 \kms$) than are seen in the data. The first result is sensitive to our adopted normalization of the mean \lya optical depth, and the second is sensitive to our assumption that helium reionization has not significantly raised gas temperatures at $z=3$. It is therefore too early to say whether these discrepancies indicate a fundamental problem with the high-redshift structure of the $\Omega=1$ CDM model or reflect errors of detail in our modeling of the gas distribution or the observational procedure.

Lyman limit and damped Lyman-alpha absorption systems probe the distribution of collapsed, cold gas at high redshift. Numerical simulations that incorporate gravity and gas dynamics can predict the abundance of such absorbers in cosmological models. We develop a semi-analytical method to correct the numerical predictions for the contribution of unresolved low mass halos, and we apply this method to the Katz et al. (1996) simulation of the standard cold dark matter model ($\Omega=1$, $h=0.5$, $\Omega_b=0.05$, $\sigma_8=0.7$). Using this simulation and higher resolution simulations of individual low mass systems, we determine the relation between a halo's circular velocity $v_c$ and its cross section for producing Lyman limit or damped absorption. We combine this relation with the Press-Schechter formula for the abundance of halos to compute the number of absorbers per unit redshift. The resolution correction increases the predicted abundances by about a factor of two at z=2, 3, and 4, bringing the predicted number of damped absorbers into quite good agreement with observations. Roughly half of the systems reside in halos with circular velocities $v_c>100\kms$ and half in halos with $35\kms<v_c<100\kms$. Halos with $v_c>150\kms$ typically harbor two or more systems capable of producing damped absorption. Even with the resolution correction, the predicted abundance of Lyman limit systems is a factor of three below observational estimates, signifying either a failure of standard CDM or a failure of these simulations to resolve the systems responsible for most Lyman limit absorption. By comparing simulations with and without star formation, we find that depletion of the gas supply by star formation affects absorption line statistics at $z>=2$ only for column densities exceeding $N_{HI}=10^{22} cm^{-2}$.

We use an automated Voigt-profile fitting procedure to extract statistical properties of the Ly$\alpha$ forest in a numerical simulation of an $\Omega=1$, cold dark matter (CDM) universe. Our analysis method is similar to that used in most observational studies of the forest, and we compare the simulations to recently published results derived from Keck HIRES spectra. With the Voigt-profile decomposition analysis, the simulation reproduces the large number of weak lines ($N_{\rm HI}\la 10^{14} cm^{-2}$) found in the HIRES spectra. At $z=3$, the $b$-parameter distribution has a median of $35\kms$ and a dispersion of $20\kms$, in reasonable agreement with the observed values. The comparison between our new analysis and recent data strengthens earlier claims that the Ly$alpha$ forest arises naturally in hierarchical structure formation as photoionized gas falls into dark matter potential wells.

Using cosmological simulations that incorporate gas dynamics and gravitational forces, we investigate the influence of photoionization by a UV radiation background on the formation of galaxies. In our highest resolution simulations, we find that photoionization has essentially no effect on the baryonic mass function of galaxies at $z=2$, down to our resolution limit of $5e9 M_\sun$. We do, however, find a strong interplay between the mass resolution of a simulation and the microphysics included in the computation of heating and cooling rates. At low resolution, a photoionizing background can appear to suppress the formation of even relatively massive galaxies. However, when the same initial conditions are evolved with a factor of eight better mass resolution, this effect disappears. Our results demonstrate the need for care in interpreting the results of cosmological simulations that incorporate hydrodynamics and radiation physics. For example, we conclude that a simulation with limited resolution may yield more realistic results if it ignores some relevant physical processes, such as photoionization. At higher resolution, the simulated population of massive galaxies is insensitive to the treatment of photoionization and star formation, but it does depend significantly on the amplitude of the initial density fluctuations. By $z=2$, an $\Omega=1$ cold dark matter model normalized to produce the observed masses of present-day clusters has already formed galaxies with baryon masses exceeding $1e11 M_\sun$.

Cosmological simulations with gas dynamics suggest that the Lyman-alpha forest is produced mainly by "small scale structure" --- filaments and sheets that are the high redshift analog of today's galaxy superclusters. There is no sharp distinction between Lyman-alpha clouds and "Gunn-Peterson" absorption produced by the fluctuating IGM -- the Lyman-alpha forest {\it is} the Gunn-Peterson effect. Lyman limit and damped Lyman-alpha absorption arises in the radiatively cooled gas of forming galaxies. At $z~2-3$, most of the gas is in the photoionized, diffuse medium associated with the Lyman-alpha forest, but most of the {\it neutral} gas is in damped Lyman-alpha systems. We discuss generic evolution of cosmic gas in a hierarchical scenario of structure formation, with particular attention to the prospects for detecting 21cm emission from high redshift HI. A scaling argument based on the present-day cluster mass function suggests that objects with $M_{HI} >~ 5e11 h^{-1} \msun$ should be extremely rare at $z~3$, so detections with existing instruments will be difficult. An instrument like the proposed Square Kilometer Array could detect individual damped Lyman-alpha systems at high redshift, making it possible to map structure in the high redshift universe in much the same way that today's galaxy redshift surveys map the local large scale structure.

Disturbed spiral galaxies with high rates of star formation pervaded clusters of galaxies just a few billion years ago, but nearby clusters exclude spirals in favor of ellipticals. ``Galaxy harassment" (frequent high speed galaxy encounters) drives the morphological transformation of galaxies in clusters, provides fuel for quasars in subluminous hosts and leaves detectable debris arcs. Simulated images of harassed galaxies are strikingly similar to the distorted spirals in clusters at $z \sim 0.4$ observed by the Hubble Space Telescope.

It has been argued that a UV photoionizing background radiation field suppresses the formation of dwarf galaxies, and may even inhibit the formation of larger galaxies. In order to test this, we present gas-dynamical simulations of the formation of small objects in a CDM universe with and without a photoionizing background. The objects are selected from a collisionless simulation at a redshift of 2.4, and rerun at higher resolution including the effects of gas dynamics and using a hierarchical grid of particles. Five objects, each with a circular speed of 46 km/sec are simulated. The presence of the photoionizing background has only a small effect on the amount of gas that collapses in these objects, reducing the amount of cold collapsed gas by at most 30%. Analysis of the smaller objects found in the higher resolution simulation indicates that the photoionizing background only significantly affects the formation of objects with a virialized halo mass less than 10^9 soalr masses and circular speeds less than 23 km/sec. However, the ionization balance is greatly changed by the presence of the background radiation field. Typical lines of sight through the objects have 4 orders of magnitude less neutral hydrogen column density when the photoionizing background is included.

Cosmological simulations with gas provide a detailed description of the intergalactic medium, making possible predictions of neutral hydrogen absorption in the spectra of background QSOs. We present results from a high-resolution calculation of an $\Omega=1$ cold dark matter model. Our simulation reproduces many of the observed properties of the \lya forest surprisingly well. The distribution of HI column densities agrees with existing data to within a factor of $\sim$ two over most of the range from $10^{14}\cm^{-2}$ to $10^{22}\cm^{-2}$; i.e., from unsaturated \lya forest lines to damped \lya systems. The equivalent width distribution matches the observed exponential form with a characteristic width $W_* \approx 0.3$ angstroms. The distribution of b-parameters appears consistent with that derived from QSO spectra. Most of the low column density absorption arises in large, flattened structures of moderate or even relatively low overdensity, so there is no sharp distinction between the \lya forest and the ``Gunn-Peterson'' absorption produced by the smooth intergalactic medium. Our results demonstrate that a \lya forest like that observed develops naturally in a hierarchical clustering scenario with a photoionizing background. Comparison between simulations and high-resolution QSO spectra should open a new regime for testing theories of cosmic structure formation.

We study the formation of damped \lya and Lyman limit absorbers in a hierarchical clustering scenario using a gas dynamical simulation of an $\Omega = 1$, cold dark matter universe. In the simulation, these high column density systems are associated with forming galaxies. Damped \lya absorption, $N_{HI} \simgt 10^{20.2}\cm^{-2}$, arises along lines of sight that pass near the centers of relatively massive, dense protogalaxies. Lyman limit absorption, $10^{17}\cm^{-2} \simlt N_{HI} \simlt 10^{20.2}\cm^{-2}$, develops on lines of sight that pass through the outer parts of such objects or near the centers of smaller protogalaxies. The number of Lyman limit systems is less than observed, while the number of damped \lya systems is quite close to the observed abundance. Damped absorbers are typically $\sim 10$ kpc in radius, but the population has a large total cross section because the systems are much more numerous than present day $L_*$ galaxies. Our results demonstrate that high column density systems like those observed arise naturally in a hierarchical theory of galaxy formation and that it is now possible to study these absorbers directly from numerical simulations.

We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989) generalized to handle comoving coordinates and periodic boundary conditions, combines smoothed--particle hydrodynamics (SPH) with the hierarchical tree method for computing gravitational forces. The Lagrangian hydrodynamics approach and individual time steps for gas particles give the algorithm a large dynamic range, which is essential for studies of galaxy formation in a cosmological context. The code incorporates radiative cooling for an optically thin, primordial composition gas in ionization equilibrium with a user-specified ultraviolet background. We adopt a phenomenological prescription for star formation that gradually turns cold, dense, Jeans-unstable gas into collisionless stars, returning supernova feedback energy to the surrounding medium. In CDM simulations, some of the baryons that fall into dark matter potential wells dissipate their acquired thermal energy and condense into clumps with roughly galactic masses. The resulting galaxy population is insensitive to assumptions about star formation; we obtain similar baryonic mass functions and galaxy correlation functions from simulations with star formation and from simulations without star formation in which we identify galaxies directly from the cold, dense gas.

N-body simulations that follow only a collisionless dark matter component have failed to produce galaxy halos or substructure within dense environments. We investigate the `over-merging' problem analytically and with numerical simulations, by calculating dissolution timescales of halos due to physical and artificial dynamical effects. The numerical resolution that has recently been attained is such that mass-loss from two-body relaxation is negligible. We demonstrate that substructure is destroyed in present simulations as a result of large force softening combined with the heating sources of tides and encounters with dissolving substructure. In the limit of infinite numerical resolution, whether or not individual halos or substructure can survive depends sensitively on their inner density profiles. Singular isothermal halos will always survive at some level, however, if halos form with large core radii then the over-merging problem will always exist within dissipationless N-body simulations. In this latter case a dissipational component can increase the halos central density enabling galaxies to survive.

A small fraction of galaxies appear to reside in dense compact groups, whose inferred crossing times are much shorter than a Hubble time. These short crossing times have led to considerable disagreement about the dynamical state of these systems. We suggest that many of the observed groups are not physically bound but are chance projections of galaxies well-separated along the line of sight. Unlike earlier similar proposals, ours does not require that the galaxies in the compact group be members of a more diffuse, but physically bound entity. The probability of physically separated galaxies projecting into an apparent compact group is non-negligible if most galaxies are distributed, as observed, in thin filaments. We illustrate this general point with a specific example: a simulation of a cold dark matter universe, in which hydrodynamic effects are included to identify galaxies. The simulated galaxy distribution is filamentary, and end-on views of these filaments produce apparent galaxy associations that have sizes and velocity dispersions similar to those of observed compact groups. The frequency of such projections can explain the observed space-density of groups in the Hickson catalog. We discuss the implications of our proposal for the formation and evolution of groups and elliptical galaxies. The proposal can be tested by using redshift-independent distance estimators to measure the line-of-sight spatial extent of nearby compact groups.