The majority of baryons reside beyond the optical extent of a galaxy in the circumgalactic and intergalactic media (CGM/IGM). Gaseous halos are inextricably linked to the appearance of their host galaxies through a complex story of accretion, feedback, and continual recycling. The energetic processes, which define the state of gas in the CGM, are the same ones that 1) regulate stellar growth so that it is not over-efficient, and 2) create the diversity of today's galaxy colors, SFRs, and morphologies spanning Hubble's Tuning Fork Diagram. They work in concert to set the speed of growth on the star-forming Main Sequence, transform a galaxy across the Green Valley, and maintain a galaxy's quenched appearance on the Red Sequence. Most baryons in halos more massive than 10^12 Msolar along with their high-energy physics and dynamics remain invisible because that gas is heated above the UV ionization states. We argue that information on many of the essential drivers of galaxy evolution is primarily contained in this "missing" hot gas phase. Completing the picture of galaxy formation requires uncovering the physical mechanisms behind stellar and SMBH feedback driving mass, metals, and energy into the CGM. By opening galactic hot halos to new wavebands, we not only obtain fossil imprints of >13 Gyrs of evolution, but observe on-going hot-mode accretion, the deposition of superwind outflows into the CGM, and the re-arrangement of baryons by SMBH feedback. A description of the flows of mass, metals, and energy will only be complete by observing the thermodynamic states, chemical compositions, structure, and dynamics of T>=10^6 K halos. These measurements are uniquely possible with a next-generation X-ray observatory if it provides the sensitivity to detect faint CGM emission, spectroscopic power to measure absorption lines and gas motions, and high spatial resolution to resolve structures.

Massive galaxies today typically are not forming stars despite being surrounded by hot gaseous halos with short central cooling times. This likely owes to some form of "quenching feedback" such as merger-driven quasar activity or radio jets emerging from central black holes. Here we implement heuristic prescriptions for these phenomena on-the-fly within cosmological hydrodynamic simulations. We constrain them by comparing to observed luminosity functions and color-magnitude diagrams from SDSS. We find that quenching from mergers alone does not produce a realistic red sequence, because 1 - 2 Gyr after a merger the remnant accretes new fuel and star formation reignites. In contrast, quenching by continuously adding thermal energy to hot gaseous halos quantitatively matches the red galaxy luminosity function and produces a reasonable red sequence. Small discrepancies remain - a shallow red sequence slope suggests that our models underestimate metal production or retention in massive red galaxies, while a deficit of massive blue galaxies may reflect the fact that observed heating is intermittent rather than continuous. Overall, injection of energy into hot halo gas appears to be a necessary and sufficient condition to broadly produce red and dead massive galaxies as observed.

We conduct a detailed comparison of broad-band spectral energy distributions of six z >= 5.5 galaxies against galaxies drawn from cosmological hydrodynamic simulations. We employ a new tool called SPOC, which constrains the physical properties of observed galaxies through a Bayesian likelihood comparison with model galaxies. For five out of six observed z>=5.5 objects, our simulated galaxies match the observations at least as well as simple star formation histories such as tau-models, with similar favored values obtained for the intrinsic physical parameters such as stellar mass and star formation rate, but with substantially smaller uncertainties. Our results are broadly insensitive to simulation choices for galactic outflows and dust reddening. Hence the existence of early galaxies as observed is broadly consistent with current hierarchical structure formation models. However, one of the six objects has photometry that is best fit by a bursty SFH unlike anything produced in our simulations, driven primarily by a high K-band flux. These findings illustrate how SPOC provides a robust tool for optimally utilizing hydrodynamic simulations (or any model that predicts galaxy SFHs) to constrain the physical properties of individual galaxies having only photometric data, as well as identify objects that challenge current models. (abridged)

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.

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.

Using high S/N STIS echelle spectra (FWHM=7 km/s) of 3C 273, we constrain the metallicities of two Lya clouds in the vicinity of the Virgo cluster. We detect C II, Si II, and Si III absorption lines in the Lya absorber at z = 0.00530. Previous observations with FUSE have revealed Ly beta - Ly theta lines at this redshift, thereby accurately constraining N(H I). We model the ionization of the gas and derive [C/H] = -1.2^{+0.3}_{-0.2}, [Si/C] = 0.2+/-0.1, and log n_{H} = -2.8+/-0.3. The model implies a small absorber thickness, ~70 pc, and thermal pressure p/k ~ 40 cm^{-3} K. It is most likely that the absorber is pressure confined by an external medium because gravitational confinement would require a very high ratio of dark matter to baryonic matter. Based on Milky Way sight lines in which carbon and silicon abundances have been reliably measured in the same interstellar cloud (including new measurements presented herein), we argue that the overabundance of Si relative to C is not due to dust depletion. Instead, this probably indicates that the gas has been predominately enriched by Type II supernovae. Such enrichment is most plausibly provided by an unbound galactic wind, given the absence of galaxies within a projected distance of 100 kpc and the presence of galaxies capable of driving a wind at larger distances. We also constrain the metallicity and physical conditions of the Virgo absorber at z = 0.00337 based on detections of O VI and H I and an upper limit on C IV. If this absorber is collisionally ionized, the O VI/C IV limit requires T > 10^{5.3} K. For either collisional ionization or photoionization, we find that [O/H] > -2.0 at z = 0.00337.

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

The Lya forest in quasar spectra probes scales where massive neutrinos can strongly suppress the growth of mass fluctuations. Using hydrodynamic simulations with massive neutrinos, we successfully test techniques developed to measure the mass power spectrum from the forest. A recent observational measurement in conjunction with a conservative implementation of other cosmological constraints places upper limits on the neutrino mass: m_nu < 5.5 eV for all values of Omega_m, and m_nu < 2.4 (Omega_m/0.17 -1) eV, if 0.2 < Omega_m <0.5 as currently observationally favored (both 95 % C.L.).

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 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)

We describe PTreeSPH, a gravity treecode combined with an SPH hydrodynamics code designed for massively parallel supercomputers having distributed memory. Our computational algorithm is based on the popular TreeSPH code of Hernquist & Katz (1989). PTreeSPH utilizes a domain decomposition procedure and a synchronous hypercube communication paradigm to build self-contained subvolumes of the simulation on each processor at every timestep. Computations then proceed in a manner analogous to a serial code. We use the Message Passing Interface (MPI) communications package, making our code easily portable to a variety of parallel systems. PTreeSPH uses individual smoothing lengths and timesteps, with a communication algorithm designed to minimize exchange of information while still providing all information required to accurately perform SPH computations. We have additionally incorporated cosmology, periodic boundary conditions with forces calculated using a quadrupole Ewald summation method, and radiative cooling and heating from a parameterized ionizing background following Katz, Weinberg & Hernquist (1996). The addition of other physical processes, such as star formation, is straightforward. A cosmological simulation from z=49 to z=2 with 64^3 gas particles and 64^3 dark matter particles requires ~6000 node-hours on a Cray T3D, with a communications overhead of ~10% and is load balanced to a ~90% level. When used on the new Cray T3E, this code will be capable of performing cosmological hydrodynamical simulations down to z=0 with ~2x10^6 particles, or to z=2 with ~10^7 particles, in a reasonable amount of time. Even larger simulations will be practical in situations where the matter is not highly clustered or when periodic boundaries are not required.