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.

The circumgalactic medium (CGM) of nearby star-forming galaxies shows clear indications of OVI absorption accompanied by little to no detectable NV absorption. This unusual spectral signature, accompanied by highly non-uniform absorption from lower ionization state species, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. Motivated by previous isotropic turbulent simulations, we carry out chemodynamical simulations of stratified media in a Navarro-Frenk-White (NFW) gravitational potential with a total mass of $10^{12}$ solar masses and turbulence that decreases radially. The simulations assume a metallicity of 0.3 solar, a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling using the MAIHEM package. We compare a suite of ionic column densities with the COS-Halos sample of low-redshift star-forming galaxies. Turbulence with an average one-dimensional velocity dispersion approximately 40 km/s, corresponding to an energy injection rate of approximately $10^{49}$ erg/yr, produces a CGM that matches many of the observed ionic column densities and ratios. In this simulation, the NVI to OVI ratio is suppressed from its equilibrium value due to a combination of radiative cooling and cooling from turbulent mixing. This level of turbulence is consistent with expectations from observations of better constrained, higher-mass systems, and could be sustained by energy input from supernovae, gas inflows, and dynamical friction from dark matter subhalos. We also conduct a higher resolution run which yields smaller-scale structures, but remains in agreement with observations.

Motivated by observations of outflowing galaxies, we investigate the combined impact of magnetic fields and radiative cooling on the evolution of cold clouds embedded in a hot wind. We perform a collection of three-dimensional adaptive mesh refinement, magnetohydrodynamical simulations that span two resolutions, and include fields that are aligned and transverse to the oncoming, super-Alfv\'enic material. Aligned fields have little impact on the overall lifetime of the clouds over the non-magnetized case, although they do increase the mixing between the wind and cloud material by a factor of $\approx 3.$ Transverse fields lead to magnetic draping, which isolates the clouds, but they also squeeze material in the direction perpendicular to the field lines, which leads to rapid mass loss. A resolution study suggests that the magnetized simulations have somewhat better convergence properties than non-magnetized simulations, and that a resolution of 64 zones per cloud radius is sufficient to accurately describe these interactions. We conclude that the combined effects of radiative cooling and magnetic fields are dependent on field orientation, but are unlikely to enhance cloud lifetimes beyond the effect of radiative cooling alone.

We study a turbulent helical dynamo in a periodic domain by solving the ideal magnetohydrodynamic (MHD) equations with the FLASH code using the divergence-cleaning eight-wave method and compare our results with direct numerical simulations (DNS) using the Pencil Code. At low resolution, FLASH reproduces the DNS results qualitatively by developing the large-scale magnetic field expected from DNS, but at higher resolution, no large-scale magnetic field is obtained. In all those cases in which a large-scale magnetic field is generated, the ideal MHD results yield too little power at small scales. As a consequence, the small-scale current helicity is too small compared with that of the DNS. The resulting net current helicity has then always the wrong sign, and its statistical average also does not approach zero at late times, as expected from the DNS. Our results have implications for astrophysical dynamo simulations of stellar and galactic magnetism using ideal MHD codes.

In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring unique capabilities for probing astrophysical processes at high redshifts and out to the low-density regions in the outskirts of galaxy clusters. Advances in SZ observations have largely been driven by developments in centimetre-, millimetre-, and submillimetre-wave instrumentation on ground-based facilities, with notable exceptions including results from the Planck satellite. Here we review the utility of the thermal, kinematic, relativistic, non-thermal, and polarised SZ effects for studies of galaxy clusters and other large scale structures, incorporating the many advances over the past two decades that have impacted SZ theory, simulations, and observations. We also discuss observational results, techniques, and challenges, and aim to give an overview and perspective on emerging opportunities, with the goal of highlighting some of the exciting new directions in this field.

Single-phase photoionization equilibrium (PIE) models are often used to infer the underlying physical properties of galaxy halos probed in absorption with ions at different ionization potentials. To incorporate the effects of turbulence, we use the MAIHEM code to model an isotropic turbulent medium exposed to a redshift zero metagalactic UV background, while tracking the ionizations, recombinations, and species-by-species radiative cooling for a wide range of ions. By comparing observations and simulations over a wide range of turbulent velocities, densities, and metallicity with a Markov chain Monte Carlo technique, we find that MAIHEM models provide an equally good fit to the observed low-ionization species compared to PIE models, while reproducing at the same time high-ionization species such as \ion{Si}{4} and \ion{O}{6}. By including multiple phases, MAIHEM models favor a higher metallicity ($Z/Z_\odot \approx 40\%$) for the circumgalactic medium compared to PIE models. Furthermore, all of the solutions require some amount of turbulence ($\sigma_{\rm 3D} \geqslant 26\ {\rm km}\ {\rm s}^{-1}$). Correlations between turbulence, metallicity, column density, and impact parameter are discussed alongside mechanisms that drive turbulence within the halo.

Outflows are a pervasive feature of mechanical feedback from super star clusters (SSC) in starburst galaxies, playing a fundamental role in galaxy evolution. Observations are now starting to confirm that outflows can undergo catastrophic cooling, suppressing adiabatic superwinds. Here we present a suite of one-dimensional, hydrodynamic simulations that study the ionization structure of these outflows and the resulting line emission generated by the cooling gas. We use the non-equilibrium atomic chemistry package within MAIHEM, our modified version of FLASH, which evolves the ionization state of the gas and computes the total cooling rate on an ion-by-ion basis. We find that catastrophically cooling models produce strong nebular line emission compared to adiabatic outflows. We also show that such models exhibit non-equilibrium conditions, thereby generating more highly ionized states than equivalent equilibrium models. When including photoionization from the parent SSC, catastrophically cooling models show strong C IV {\lambda}1549 and O VI {\lambda}1037 emission. For density bounded photoionization, He II {\lambda}1640, {\lambda}4686, C III] {\lambda}1908, Si IV {\lambda}1206, and Si III {\lambda}1400 are also strongly enhanced. These lines are seen in extreme starbursts where catastrophic cooling is likely to occur, suggesting that they may serve as diagnostics of such conditions. The higher ionization generated by these flows may help to explain line emission that cannot be attributed to SSC photoionization alone.

Cosmic rays play a pivotal role in launching galactic winds, particularly in quiescently star-forming galaxies where the hot gas alone is not sufficient to drive a wind. Except for the Milky Way, not much is known about the transport of cosmic rays in galaxies. In this Letter, we present low-frequency observations of the nearby edge-on spiral galaxy NGC 4565 using the LOw-Frequency ARray (LOFAR). With our deep 144-MHz observations, we obtain a clean estimate of the emission originating from old cosmic-ray electrons (CRe), which is almost free from contamination by thermal emission. We measured vertical profiles of the non-thermal radio continuum emission that we fitted with Gaussian and exponential functions. The different profile shapes correspond to 1D cosmic-ray transport models of pure diffusion and advection, respectively. We detect a warp in the radio continuum that is reminiscent of the previously known HI warp. Because the warp is not seen at GHz-frequencies in the radio continuum, its minimum age must be about 100 Myr. The warp also explains the slight flaring of the thick radio disc that can otherwise be well described by a Gaussian profile with an FWHM of 65 arcsec (3.7 kpc). The diffusive radio halo together with the extra-planar X-ray emission may be remnants of enhanced star-forming activity in the past where the galaxy had a galactic wind, as GHz-observations indicate only a weak outflow in the last 40 Myr. NGC 4565 could be in transition from an outflow- to an inflow-dominated phase.

We investigate the ability of current and third-generation gravitational wave (GW) detectors to determine the delay time distribution (DTD) of binary neutron stars (BNS) through a direct measurement of the BNS merger rate as a function of redshift. We assume that the DTD follows a power law distribution with a slope $\Gamma$ and a minimum merger time $t_{\rm min}$, and also allow the overall BNS formation efficiency per unit stellar mass to vary. By convolving the DTD and mass efficiency with the cosmic star formation history, and then with the GW detector capabilities, we explore two relevant regimes. First, for the current generation of GW detectors, which are only sensitive to the local universe, but can lead to precise redshift determinations via the identification of electromagnetic counterparts and host galaxies, we show that the DTD parameters are strongly degenerate with the unknown mass efficiency and therefore cannot be determined uniquely. Second, for third-generation detectors such as Einstein Telescope (ET) and Cosmic Explorer (CE), which will detect BNS mergers at cosmological distances, but with a redshift uncertainty inherent to GW-only detections ($\delta(z)/z\approx 0.1z$), we show that the DTD and mass efficiency can be well-constrained to better than 10\% with a year of observations. This long-term approach to determining the DTD through a direct mapping of the BNS merger redshift distribution will be supplemented by more near term studies of the DTD through the properties of BNS merger host galaxies at $z\approx 0$ (Safarzadeh & Berger 2019).

We model the history of Galactic r-process enrichment using high-redshift, high-resolution zoom cosmological simulations of a Milky Way (MW) type halo. We assume that all r-process sources are neutron star mergers (NSMs) with a power law delay time distribution. We model the time to mix pollutants at subgrid scales, which allows us to to better compute the properties of metal poor (MP) and carbon enhanced metal poor (CEMP) stars, along with statistics of their r-process enhanced subclasses. Our simulations underpredict the cumulative ratios of r-process enhanced MP and CEMP stars (MP-r, CEMP-r) over MP and CEMP stars by about one order of magnitude, even when the minimum coalescence time of the double neutron stars ($t_{\rm min}$) is set to 1 Myr. No r-process enhanced stars form if $t_{\rm min}=100$ Myr. Our results show that even when we adopt the r-process yield estimates observed in GW170817, NSMs by themselves can only explain the observed frequency of r-process enhanced stars if either the birth rate of double neutron stars per unit mass of stars is boosted to $\approx10^{-4} M_\odot^{-1}$ or the Europium yield of each NSM event is boosted to $\approx 10^{-4} M_{\odot}$.

We present a suite of one-dimensional spherically-symmetric hydrodynamic simulations that study the atomic ionization structure of galactic outflows. We track the ionization state of the outflowing gas with a non-equilibrium atomic chemistry network that includes photoionization, photo-heating, and ion-by-ion cooling. Each simulation describes a steady-state outflow that is defined by its mass and energy input rates, sonic radius, metallicity, and UV flux from both the host galaxy and meta-galactic background. We find that for a large range of parameter choices, the ionization state of the material departs strongly from what it would be in photo-ionization equilibrium, in conflict with what is commonly assumed in the analysis of observations. In addition, nearly all the models reproduce the low N V to O VI column density ratios and the relatively high O VI column densities that are observed.

The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UFD) galaxies such as Reticulum II and Tucana III. Assuming that such dwarf systems have lost a significant fraction of their stellar mass through tidal stripping, we conclude that contrary to most current models, it is the DNSs with rather large natal kicks but very short merging timescales that can enrich UFD-type galaxies. These binaries are either on highly eccentric orbits, or form with very short separations due to an additional mass-transfer between the first-born neutron star and a naked helium star, progenitor of the second-born neutron star. These DNSs are born with a frequency that agrees with the statistics of the r-process UFDs, and merge well within the virial radius of their host halos, therefore contributing significantly to their r-process enrichment.

We study the observational consequences of several unknown properties of Population III (Pop III) stars using large-scale cosmological simulations that include a subgrid model to track the unresolved mixing of pollutants. Varying the value of the critical metallicity that marks the boundary between Pop III and Population II (Pop II) star formation across 2 dex has a negligible effect on the fraction of Pop III stars formed and the subsequent fraction of Pop III flux from high-redshift galaxies. However, adopting a log normal initial mass function (IMF) for Pop III stars, in place of a baseline Salpeter IMF, results in a Pop III star formation rate density (SFRD) that is 1/4 of the baseline rate. The flux from high-redshift galaxies modeled with this IMF is highly bimodal, resulting in a tiny fraction of $z \leq 8$ galaxies with more than 75\% of their flux coming from Pop III stars. However, at $z=9$, right before reionization in our simulations, $\approx$ 20\% of galaxies are Pop III-bright with $m_{\rm UV} \le 31.4$ mag and at least 75\% of their flux generated by Pop III stars . Additionally, the log normal Pop III IMF results in a population of carbon enhanced, metal poor stars in reasonable agreement with MW halo observations. Our analysis supports the conclusion that the Pop III IMF was dominated by stars in the 20-120$M_{\odot}$ range that generate SN with carbon-enhanced ejecta.

Most galaxies comparable to or larger than the mass of the Milky Way host hot, X-ray emitting atmospheres, and many such galaxies are radio sources. Hot atmospheres and radio jets and lobes are the ingredients of radio-mechanical active galactic nucleus (AGN) feedback. While a consensus has emerged that such feedback suppresses cooling of hot cluster atmospheres, less attention has been paid to massive galaxies where similar mechanisms are at play. Observation indicates that the atmospheres of elliptical and S0 galaxies were accreted externally during the process of galaxy assembly and augmented significantly by stellar mass loss. Their atmospheres have entropy and cooling time profiles that are remarkably similar to those of central cluster galaxies. About half display filamentary or disky nebulae of cool and cold gas, much of which has likely cooled from the hot atmospheres. We review the observational and theoretical perspectives on thermal instabilities in galactic atmospheres and the evidence that AGN heating is able to roughly balance the atmospheric cooling. Such heating and cooling may be regulating star formation in all massive spheroids at late times.

Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($\Sigma_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC 3184, 4736, and 5055) from the LOFAR Two-metre Sky Survey (LoTSS), and archival LOw-Frequency ARray (LOFAR) data of NGC 5194. Maps were created with the facet calibration technique and converted to radio $\Sigma_{\rm SFR}$ maps using the Condon relation. We compared these maps with hybrid $\Sigma_{\rm SFR}$ maps from a combination of GALEX far-ultraviolet and Spitzer 24 $\mu\rm m$ data using plots tracing the relation at $1.2\times 1.2$-kpc$^2$ resolution. The RC emission is smoothed with respect to the hybrid $\Sigma_{\rm SFR}$ owing to the transport of cosmic-ray electrons (CREs). This results in a sublinear relation $(\Sigma_{\rm SFR})_{\rm RC} \propto [(\Sigma_{\rm SFR})_{\rm hyb}]^{a}$, where $a=0.59\pm 0.13$ (140 MHz) and $a=0.75\pm 0.10$ (1365 MHz). Both relations have a scatter of $\sigma = 0.3~\rm dex$. If we restrict ourselves to areas of young CREs ($\alpha > -0.65$; $I_\nu \propto \nu^\alpha$), the relation becomes almost linear at both frequencies with $a\approx 0.9$ and a reduced scatter of $\sigma = 0.2~\rm dex$. We then simulate the effect of CRE transport by convolving the hybrid $\Sigma_{\rm SFR}$ maps with a Gaussian kernel until the RC-SFR relation is linearised; CRE transport lengths are $l=1$-5 kpc. Solving the CRE diffusion equation, we find diffusion coefficients of $D=(0.13$-$1.5) \times 10^{28} \rm cm^2\,s^{-1}$ at 1 GeV. A RC-SFR relation at $1.4$ GHz can be exploited to measure SFRs at redshift $z \approx 10$ using $140$ MHz observations.

Molecular clouds exhibit lognormal probability density functions (PDF) of mass densities, which are thought to arise as a consequence of isothermal, supersonic turbulence. Star formation is then widely assumed to occur in perturbations in which gravitational collapse is faster than the rate of change due to turbulent motions. Here we use direct numerical simulations to measure this rate as a function of density for a range of turbulent Mach numbers, and show it is faster at high densities than at low densities. Furthermore, we show that both the density PDF and rate of change arise naturally in a simple model of turbulence as a continuous Markov process. The one-dimensional Langevin equation that describes this evolution depends on only two parameters, yet it captures the full evolution seen in direct three-dimensional simulations. If it is modified to include gravity, the Langevin equation also reproduces the rate of material collapsing to high densities seen in turbulent simulations including self-gravity. When generalized to include both temperature and density, similar analyses are likely applicable throughout astrophysics.

The circumgalactic medium (CGM) of nearby star-forming galaxies show clear indications of \ion{O}{6} absorption accompanied by little to no \ion{N}{5} absorption. This unusual spectral signature, accompanied by absorption from lower ionization state species whose columns vary by orders of magnitude along \st{difference} \textbf{different} sightlines, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. To explore this possibility, we carry out a series of chemodynamical simulations of a isotropic turbulent media, using the MAIHEM package. The simulations assume a metallicity of $0.3\ Z_{\odot}$ and a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling for a wide range of elements. We find that turbulence with a one-dimensional velocity dispersion of $\sigma_{1D} \approx 60$ km/s replicates many of the observed features within the CGM, such as clumping of low ionization-state ions and the existence of \ion{O}{6} at moderate ionization parameters. However, unlike observations, \ion{N}{5} often arises in our simulations with derived column densities of a similar magnitude to those of \ion{O}{6}. While higher values of $\sigma_{1D}$ lead to a thermal runaway in our isotropic simulations, this would not be the case in stratified media, and thus we speculate that more complex models of the turbulence may well match the absence of \ion{N}{5} in the CGM of star-forming galaxies.

Energetic feedback from active galactic nuclei (AGNs) is often used in simulations to resolve several outstanding issues in galaxy formation, but its impact is still not fully understood. Here we derive new constraints on AGN feedback by comparing observations and simulations of the thermal Sunyaev-Zel'dovich (tSZ) effect. We draw on observational results presented in Spacek et al. (2016, 2017) who used data from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) to measure the tSZ signal from >= 10^11 M_Sun and >= 1 Gyr galaxies at z=0.5-1.0 (low-z) and z=1.0-1.5 (high-z). Using the large-scale cosmological hydrodynamical simulations Horizon-AGN and Horizon-NoAGN, which include and omit AGN feedback, we extract simulated tSZ measurements around galaxies equivalent to the observational work. We find that the Horizon-AGN results only differ from the SPT measurements at levels of 0.4 sigma (low-z) and 0.6 sigma (high-z), but differ from the ACT measurements by 3.4 sigma (low-z) and 2.3 sigma (high-z). The Horizon-NoAGN results provide a slightly better fit to the SPT measurements by differing by 0.2 sigma (low-z) and 0.4 sigma (high-z), but a significantly better match to the ACT measurements by differing by only 0.5 sigma (low-z) and 1.4 sigma (high-z). We conclude that, while the lower-mass (<~ 5 x 10^11 M_Sun) SPT results allow for the presence AGN feedback energy, the higher-mass (>~ 5 x 10^11 M_Sun) ACT results show significantly less energy than predicted in the simulation including AGN feedback, while more closely matching the simulation without AGN feedback, indicating that AGN feedback may be milder than often predicted in simulations.

Absorption line studies are essential to understanding the origin, nature, and impact of starburst-driven galactic outflows. Such studies have revealed a multiphase medium with a number of poorly-understood features leading to a need to study the ionization mechanism of this gas. To better interpret these observations, we make use of a suite of adaptive mesh refinement hydrodynamic simulations of cold, atomic clouds driven by hot, supersonic outflows, including the effect of radiative cooling, thermal conduction, and an ionizing background characteristic of a starbursting galaxy. Using a new analysis tool, Trident, we estimate the equilibrium column density distributions for ten different ions: H I, Mg II, C II, C III, C IV, Si III, Si IV, N V, O VI, and Ne VIII. These are fit to model profiles with two parameters describing the maximum column density and coverage, and for each ion we provide a table of these fit parameters, along with average velocities and line widths. Our results are most sensitive to Mach number and conduction efficiency, with higher Mach numbers and more efficient conduction leading to more compact, high column density clouds. We use our results to interpret down-the-barrel observations of outflows and find that the adopted ionization equilibrium model overpredicts column densities of ions such as Si IV and does not adequately capture the observed trends for N V and O VI, implying the presence of strong non equilibrium ionization effects.

The recent EDGES collaboration detection of an absorption signal at a central frequency of $\nu = 78 \pm 1$ MHz points to the presence of a significant Lyman-$\alpha$ background by a redshift of $z=18$. The timing of this signal constrains the dark matter particle mass ($m_\chi$) in the warm dark matter (WDM) cosmological model. WDM delays the formation of small-scale structures, and therefore a stringent lower limit can be placed on $m_\chi,$ based on the presence of a sufficiently strong Ly-$\alpha$ background due to star formation at $z=18$. Our results show that the coupling the spin temperature to the gas through Ly-$\alpha$ pumping requires a minimum mass of $m_\chi>3$ keV if atomic cooling halos dominate the star formation rate at $z=18,$ and $m_\chi>2$ keV if ${\rm\,H_2}$ cooling halos also form stars efficiently at this redshift. These limits match or exceed the most stringent limits cited to date in the literature, even in the face of the many uncertainties regarding star-formation at high redshift.

The smallest satellites of the Milky Way ceased forming stars during the epoch of reionization and thus provide archaeological access to galaxy formation at $z>6$. Numerical studies of these ultra-faint dwarf galaxies (UFDs) require expensive cosmological simulations with high mass resolution that are carried out down to $z=0$. However, if we are able to statistically identify UFD host progenitors at high redshifts \emph{with relatively high probabilities}, we can avoid this high computational cost. To find such candidates, we analyze the merger trees of Milky Way type halos from the high-resolution ${\it Caterpillar}$ suite of dark matter only simulations. Satellite UFD hosts at $z=0 $ are identified based on four different abundance matching (AM) techniques. All the halos at high redshifts are traced forward in time in order to compute the probability of surviving as satellite UFDs today. Our results show that selecting potential UFD progenitors based solely on their mass at z=12 (8) results in a 10\% (20\%) chance of obtaining a surviving UFD at $z=0$ in three of the AM techniques we adopted. We find that the progenitors of surviving satellite UFDs have lower virial ratios ($\eta$), and are preferentially located at large distances from the main MW progenitor, while they show no correlation with concentration parameter. Halos with favorable locations and virial ratios are $\approx 3$ times more likely to survive as satellite UFD candidates at $z=0.$

Direct observational searches for Population III (Pop III) stars at high redshift are faced with the question of how to select the most promising targets for spectroscopic follow-up. To help answer this, we use a large-scale cosmological simulation, augmented with a new subgrid model that tracks the fraction of pristine gas, to follow the evolution of high-redshift galaxies and the Pop III stars they contain. We generate rest-frame ultraviolet (UV) luminosity functions for our galaxies and find that they are consistent with current $z \ge 7 $ observations. Throughout the redshift range $7 \le z \le 15$ we identify "Pop III-bright" galaxies as those with at least 75% of their flux coming from Pop III stars. While less than 1% of galaxies brighter than $m_{\rm UV, AB} = 31.4$ mag are Pop III--bright in the range $7\leq z \leq8$, roughly 17% of such galaxies are Pop III--bright at $z=9$, immediately before reionization occurs in our simulation. Moving to $z=10$, $m_{\rm UV, AB} = 31.4$ mag corresponds to larger, more luminous galaxies and the Pop III-bright fraction falls off to 5%. Finally, at the highest redshifts, a large fraction (29% at $z=14$ and 41% at $z=15)$ of all galaxies are Pop III-bright regardless of magnitude. While $m_{\rm UV, AB} = 31.4$ mag galaxies are extremely rare during this epoch, we find that 13% of galaxies at $z = 14$ are Pop III-bright with $m_{\rm UV, AB} \le 33$ mag, an intrisic magnitude within reach of the James Webb Space Telescope using lensing. Thus, we predict that the best redshift to search for luminous Pop III--bright galaxies is just before reionization, while lensing surveys for fainter galaxies should push to the highest redshifts possible.

We investigate the stellar mass loss of gas rich galaxies falling into clusters due to the change in the gravitational potential caused by the ram pressure stripping of their gas. We model the satellites with exponential stellar and gas disk profiles, assume rapid ram pressure stripping, and follow the stellar orbits in the shocked potential. Due to the change of the potential, the stars move from circular orbits to elliptical orbits with apocenters that are often outside the tidal radius, causing those stars to be stripped. We explore the impact of the redshift of infall, gas fraction, satellite halo mass and cluster mass on this process. The puffing of the satellites makes them appear as ultra diffuse galaxies, and the stripped stars contribute to the intracluster light. Our results show these effects are most significant for less-massive satellites, which have larger gas fractions when they are accreted into clusters. The preferential destruction of low mass systems causes the red fraction of cluster galaxies to be smaller at lower masses, an observation that is otherwise difficult to explain.

Motivated by the observed differences in the nebular emission of nearby and high-redshift galaxies, we carry out a set of direct numerical simulations of turbulent astrophysical media exposed to a UV background. The simulations assume a metallicity of $Z/Z_{\odot}$=0.5 and explicitly track ionization, recombination, charge transfer, and ion-by-ion radiative cooling for several astrophysically important elements. Each model is run to a global steady state that depends on the ionization parameter $U$, and the one-dimensional turbulent velocity dispersion, $\sigma_{\rm 1D}$, and the turbulent driving scale. We carry out a suite of models with a T=42,000K blackbody spectrum, $n_e$ = 100 cm$^{-3}$ and $\sigma_{\rm 1D}$ ranging between 0.7 to 42 km s$^{-1},$ corresponding to turbulent Mach numbers varying between 0.05 and 2.6. We report our results as several nebular diagnostic diagrams and compare them to observations of star-forming galaxies at a redshift of $z\approx$2.5, whose higher surface densities may also lead to more turbulent interstellar media. We find that subsonic, transsonic turbulence, and turbulence driven on scales of 1 parsec or greater, have little or no effect on the line ratios. Supersonic, small-scale turbulence, on the other hand, generally increases the computed line emission. In fact with a driving scale $\approx 0.1$ pc, a moderate amount of turbulence, $\sigma_{\rm 1D}$=21-28 km s$^{-1},$ can reproduce many of the differences between high and low redshift observations without resorting to harder spectral shapes.

To explain the high observed abundances of r-process elements in local ultra- faint dwarf (UFD) galaxies, we perform cosmological zoom simulations that include r-process production from neutron star mergers (NSMs). We model star-formation stochastically and simulate two different halos with total masses $\approx 10^8 M_{\odot}$ at z = 6. We find that the final distribution of [Eu/H] vs. [Fe/H] is relatively insensitive to the energy by which the r-process material is ejected into the interstellar medium, but strongly sensitive to the environment in which the NSM event occurs. In one halo the NSM event takes place at the center of the stellar distribution, leading to high-levels of r-process enrichment such as seen in a local UFD, Reticulum II (Ret II). In a second halo, the NSM event takes place outside of the densest part of the galaxy, leading to a more extended r-process distribution. The subsequent star formation occurs in an interstellar medium with shallow levels of r-process enrichment which results in stars with low levels of [Eu/H] compared to Ret II stars even when the maximum possible r-process mass is assumed to be ejected. This suggests that the natal kicks of neutron stars may also play an important role in determining the r-process abundances in UFD galaxies, a topic that warrants further theoretical investigation.

I present a suite of three-dimensional simulations of the evolution of initially-hot material ejected by starburst-driven galaxy outflows. The simulations are conducted in a comoving frame that moves with the material, tracking atomic/ionic cooling, Compton cooling, and dust cooling and destruction. Compton cooling is most efficient of these processes, while the main role of atomic/ionic cooling is to enhance density inhomogeneities. Dust, on the other hand, has little effect on the outflow evolution, and is rapidly destroyed in all the simulations except the case with the smallest mass flux. I use the results to construct a simple steady-state model of the observed UV/optical emission from each outflow. The velocity profiles in this case are dominated by geometric effects, and the overall luminosities are extremely strong functions of the properties of the host system, as observed in ultra-luminous infrared galaxies (ULIRGs). Furthermore the luminosities and maximum velocities in several models are consistent with emission-line observations of ULIRGs, although the velocities are significantly greater than observed in absorption-line studies. It may be that absorption line observations of galaxy outflows probe entrained cold material at small radii, while emission-line observations probe cold material condensing from the initially hot medium at larger distances.

Recent Hubble Space Telescope (HST) observations of the extragalactic radio source Centaurus A (Cen A) display a young stellar population around the southwest tip of the inner filament 8.5 kpc from the Cen A galactic center, with ages in the range of 1-3 Myr. Crockett et al. (2012) argue that the transverse bow shock of the Cen A jet triggered this star formation as it impacted dense molecular cores of clouds in the filament. To test this hypothesis, we perform three-dimensional numerical simulations of induced star formation by the jet bow shock in the inner filament of Cen A, using a positivity preserving WENO method to solve the equations of gas dynamics with radiative cooling. We find that star clusters form inside a bow-shocked molecular cloud when the maximum initial density of the cloud is > 40 H2 molecules/cm^3. In a typical molecular cloud of mass 10^6 M_sun and diameter 200 pc, approximately 20 star clusters of mass 10^4 M_sun are formed, matching the HST images.

We directly measure the thermal energy of the gas surrounding galaxies through the thermal Sunyaev-Zel'dovich (tSZ) effect. We perform a stacking analysis of microwave background images from the Atacama Cosmology Telescope, around 1179 massive quiescent elliptical galaxies at 0.5 <= z <= 1.0 ('low-z') and 3274 galaxies at 1.0 <= z <= 1.5 ('high-z'), selected using data from the Wide-Field Infrared Survey Explorer All-Sky Survey and the Sloan Digital Sky Survey (SDSS) within the SDSS Stripe-82 field. The gas surrounding these galaxies is expected to contain energy from past episodes of active galactic nucleus (AGN) feedback, and after using modeling to subtract undetected contaminants, we detect a tSZ signal at a significance of 0.9-sigma for our low-z galaxies and 1.8-sigma for our high-z galaxies. We then include data from the high-frequency Planck bands for a subset of 227 low-z galaxies and 529 high-z galaxies and find low-z and high-z tSZ detections of 1.0-sigma and 1.5-sigma, respectively. These results indicate an average thermal heating around these galaxies of 5.6(+5.9/-5.6) x 10^60 erg for our low-z galaxies and 7.0(+4.7/-4.4) x 10^60 erg for our high-z galaxies. Based on simple heating models, these results are consistent with gravitational heating without additional heating due to AGN feedback.

We make use of new subgrid model of turbulent mixing to accurately follow the cosmological evolution of the first stars, the mixing of their supernova ejecta, and the impact on the chemical composition of the Galactic Halo. Using the cosmological adaptive mesh refinement code RAMSES, we implement a model for the pollution of pristine gas as described in Pan et al. Tracking the metallicity of Pop III stars with metallicities below a critical value allows us to account for the fraction of Z < Zcrit stars formed even in regions in which the gas' average metallicity is well above Zcrit. We demonstrate that such partially-mixed regions account for 0.5 to 0.7 of all Pop III stars formed up to z = 5. Additionally, we track the creation and transport of "primordial metals" (PM) generated by Pop III supernovae (SNe). These neutron-capture deficient metals are taken up by second-generation stars and likely lead to unique abundance signatures characteristic of carbon-enhanced, metal-poor (CEMP-no) stars. As an illustrative example, we associate primordial metals with abundance ratios used by Keller et al. to explain the source of metals in the star SMSS J031300.36-670839.3, finding good agreement with the observed [Fe/H], [C/H], [O/H], and [Mg/Ca] ratios in CEMP-no Milky Way halo stars. Similar future simulations will aid in further constraining the properties of Pop III stars using CEMP observations, as well as improve predictions of the spatial distribution of Pop III stars, as will be explored by the next generation of ground- and space-based telescopes.

We investigate the relation between the turbulent Mach number (\mach) and the escape fraction of Lyman continuum photons ($f_{\rm esc}$) in high-redshift galaxies. Approximating the turbulence as isothermal and isotropic, we show that the increase in the variance in column densities from $\mathcal{M}=1$ to $\mathcal{M}=10$ causes $f_{\rm esc}$ to increase by $\approx 25$\%, and the increase from $\mathcal{M}=1$ to $\mathcal{M}=20$ causes $f_{\rm esc}$ to increases by $\approx 50$\% for a medium with opacity $\tau\approx1$. At a fixed Mach number, the correction factor for escape fraction relative to a constant column density case scales exponentially with the opacity in the cell, which has a large impact for simulated star forming regions. Furthermore, in simulations of isotropic turbulence with full atomic/ionic cooling and chemistry, the fraction of HI drops by a factor of $\approx 2.5$ at $\mathcal{M}\approx10$ even when the mean temperature is $\approx5\times10^3 K$. If turbulence is unresolved, these effects together enhance $f_{\rm esc}$ by a factor $>3$ at Mach numbers above 10. Such Mach numbers are common at high-redshifts where vigorous turbulence is driven by supernovae, gravitational instabilities, and merger activity, as shown both by numerical simulations and observations. These results, if implemented in the current hydrodynamical cosmological simulations to account for unresolved turbulence, can boost the theoretical predictions of the Lyman Continuum photon escape fraction and further constrain the sources of reionization.

We study the origin of the cold molecular clumps in quasar outflows, recently detected in CO and HCN emission. We first describe the physical properties of such radiation-driven outflows and show that a transition from a momentum- to an energy-driven flow must occur at a radial distance of R ~ 0.25 kpc. During this transition, the shell of swept up material fragments due to Rayleigh-Taylor instabilities, but these clumps contain little mass and are likely to be rapidly ablated by the hot gas in which they are immersed. We then explore an alternative scenario in which clumps form from thermal instabilities at R >~ 1 kpc, possibly containing enough dust to catalyze molecule formation. We investigate this processes with 3D two-fluid (gas+dust) numerical simulations of a kpc^3 patch of the outflow, including atomic and dust cooling, thermal conduction, dust sputtering, and photoionization from the QSO radiation field. In all cases, dust grains are rapidly destroyed in ~10,000 years; and while some cold clumps form at later times, they are present only as transient features, which disappear as cooling becomes more widespread. In fact, we only find a stable two-phase medium with dense clumps if we artificially enhance the QSO radiation field by a factor 100. This result, together with the complete destruction of dust grains, renders the interpretation of molecular outflows a very challenging problem.

Energetic feedback due to active galactic nuclei (AGN) is likely to play an important role in the observed anti-hierarchical trend in the evolution of galaxies, and yet the energy injected into the circumgalactic medium by this process is largely unknown. One promising approach to constrain this feedback is through measurements of CMB spectral distortions due to the thermal Sunyaev-Zel'dovich (tSZ) effect, whose magnitude is directly proportional to the energy input by AGN. Here we co-add South Pole Telescope SZ (SPT-SZ) survey data around a large set of massive quiescent elliptical galaxies at z >= 0.5. We use data from the Blanco Cosmology Survey and VISTA Hemisphere Survey to create a large catalog of galaxies split up into two redshift bins, with 3394 galaxies at 0.5 <= z <= 1.0 and 924 galaxies at 1.0 <= z <= 1.5, with typical stellar masses of 1.5 x 10^11 M_Sun. We then co-add the emission around these galaxies, resulting in a measured tSZ signal at 2.2 sigma significance for the lower redshift bin and a contaminating signal at 1.1 sigma for the higher redshift bin. To remove contamination due to dust emission, we use SPT-SZ source counts to model a contaminant source population in both the SPT-SZ bands and Planck high-frequency bands for a subset of 937 low-redshift galaxies and 240 high-redshift galaxies. This increases our detection to 3.6 sigma for low redshifts and 0.9 sigma for high redshifts. We find the mean angularly-integrated Compton-y values to be 2.2 (-0.7+0.9) x 10^-7 Mpc^2 for low redshifts and 1.7 (-1.8+2.2) x 10^-7 Mpc^2 for high redshifts, corresponding to total thermal energies of 7.6 (-2.3+3.0) x 10^60 ergs and 6.0 (-6.3+7.7) x 10^60 ergs, respectively. These numbers are higher than expected from simple theoretical models that do not include AGN feedback, and serve as constraints that can be applied to current simulations of massive galaxy formation. (abridged)

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

We explore the impact of electron thermal conduction on the evolution of radiatively-cooled cold clouds embedded in flows of hot and fast material, as occur in outflowing galaxies. Performing a parameter study of three-dimensional adaptive mesh refinement hydrodynamical simulations, we show that electron thermal conduction causes cold clouds to evaporate, but it can also extend their lifetimes by compressing them into dense filaments. We distinguish between low column-density clouds, which are disrupted on very short times, and high-column density clouds with much-longer disruption times that are set by a balance between impinging thermal energy and evaporation. We provide fits to the cloud lifetimes and velocities that can be used in galaxy-scale simulations of outflows, in which the evolution of individual clouds cannot be modeled with the required resolution. Moreover, we show that the clouds are only accelerated to a small fraction of the ambient velocity because compression by evaporation causes the clouds to present a small cross-section to the ambient flow. This means that either magnetic fields must suppress thermal conduction, or that the cold clouds observed in galaxy outflows are not formed of cold material carried out from the galaxy.

High surface density, rapidly star-forming galaxies are observed to have $\approx 50-100\,{\rm km\,s^{-1}}$ line-of-sight velocity dispersions, which are much higher than expected from supernova driving alone, but may arise from large-scale gravitational instabilities. Using three-dimensional simulations of local regions of the interstellar medium, we explore the impact of high velocity dispersions that arise from these disk instabilities. Parametrizing disks by their surface densities and epicyclic frequencies, we conduct a series of simulations that probe a broad range of conditions. Turbulence is driven purely horizontally and on large scales, neglecting any energy input from supernovae. We find that such motions lead to strong global outflows in the highly-compact disks that were common at high redshifts, but weak or negligible mass loss in the more diffuse disks that are prevalent today. Substantial outflows are generated if the one-dimensional horizontal velocity dispersion exceeds $\approx 35\,{\rm km\,s^{-1}},$ as occurs in the dense disks that have star formation rate densities above $\approx 0.1\,{\rm M}_\odot\,{\rm yr}^{-1}\,{\rm kpc}^{-2}.$ These outflows are triggered by a thermal runaway, arising from the inefficient cooling of hot material coupled with successive heating from turbulent driving. Thus, even in the absence of stellar feedback, a critical value of the star-formation rate density for outflow generation can arise due to a turbulent heating instability. This suggests that in strongly self-gravitating disks, outflows may be enhanced by, but need not caused by, energy input from supernovae.

We carry out direct numerical simulations of turbulent astrophysical media exposed to the redshift zero metagalactic background. The simulations assume solar composition and explicitly track ionizations, recombinations, and ion-by-ion radiative cooling for hydrogen, helium, carbon, nitrogen, oxygen, neon, sodium, magnesium, silicon, sulfur, calcium, and iron. Each run reaches a global steady state that not only depends on the ionization parameter, $U,$ and mass-weighted average temperature, $T_{\rm MW},$ but also on the the one-dimensional turbulent velocity dispersion, \soned. We carry out runs that span a grid of models with $U$ ranging from 0 to 10$^{-1}$ and \soned\ ranging from 3.5 to 58 km s$^{-1}$, and we vary the product of the mean density and the driving scale of the turbulence, $nL,$ which determines the average temperature of the medium, from $nL =10^{16}$ to $nL =10^{20}$ cm$^{-2}$. The turbulent Mach numbers of our simulations vary from $M \approx 0.5$ for the lowest velocity dispersions cases to $M \approx 20$ for the largest velocity dispersion cases. When $M \lesssim1,$ turbulent effects are minimal, and the species abundances are reasonably described as those of a uniform photoionized medium at a fixed temperature. On the other hand, when $M \gtrsim 1,$ dynamical simulations such as the ones carried out here are required to accurately predict the species abundances. We gather our results into a set of tables, to allow future redshift zero studies of the intergalactic medium to account for turbulent effects.

We investigate the high-redshift evolution of the restframe UV-luminosity function (LF) of galaxies via hydrodynamical cosmological simulations, coupled with an emulated observational astronomy pipeline that provides a direct comparison with observations. We do this by creating mock images and synthetic galaxy catalogs of approximately 100 square arcminute fields from the numerical model at redshifts ~ 4.5 to 10.4. We include the effects of dust extinction and the point spread function (PSF) for the Hubble WFC3 camera for comparison with space observations. We also include the expected zodiacal background to predict its effect on space observations, including future missions such as the James Webb Space Telescope (JWST). When our model catalogs are fitted to Schechter function parameters, we predict that the faint-end slope alpha of the LF evolves as alpha = -1.16 - 0.12 z over the redshift range z ~ 4.5 to 7.7, in excellent agreement with observations from e.g., Hathi et al. (2010). However, for redshifts z ~ 6 to 10.4, alpha(z) appears to display a shallower evolution, alpha = -1.79 - 0.03 z. Augmenting the simulations with more detailed physics - specifically stellar winds and supernovae (SN) - produces similar results. The model shows an overproduction of galaxies, especially at faint magnitudes, compared with the observations, although the discrepancy is reduced when dust extinction is taken into account.

To better understand the nature of the multiphase material found in outflowing galaxies, we study the evolution of cold clouds embedded in flows of hot and fast material. Using a suite of adaptive-mesh refinement simulations that include radiative cooling, we investigate both cloud mass loss and cloud acceleration under the full range of conditions observed in galaxy outflows. The simulations are designed to track the cloud center of mass, enabling us to study the cloud evolution at long disruption times. For supersonic flows, a Mach cone forms around the cloud, which damps the Kelvin-Helmholtz instability but also establishes a streamwise pressure gradient that stretches the cloud apart. If time is expressed in units of the cloud crushing time, both the cloud lifetime and the cloud acceleration rate are independent of cloud radius, and we find simple scalings for these quantities as a function of the Mach number of the external medium. A resolution study suggests that our simulations have sufficient resolution to accurately describe the evolution of cold clouds in the absence of thermal conduction and magnetic fields, physical processes whose roles will be studied in forthcoming papers.

We carry out direct numerical simulations of turbulent astrophysical media that explicitly track ionizations, recombinations, and species-by-species radiative cooling. The simulations assume solar composition and follows the evolution of hydrogen, helium, carbon, oxygen, sodium, and magnesium, but they do not include the presence of an ionizing background. In this case, the medium reaches a global steady state that is purely a function of the one-dimensional turbulent velocity dispersion, $\sigma_{\rm 1D},$ and the product of the mean density and the driving scale of turbulence, $n L.$ Our simulations span a grid of models with $\sigma_{\rm 1D}$ ranging from 6 to 58 km s$^{-1}$ and $n L$ ranging from 10$^{16}$ to 10$^{20}$ cm$^{-2},$ which correspond to turbulent Mach numbers from $M=0.2$ to 10.6. The species abundances are well described by single-temperature estimates whenever $M$ is small, but local equilibrium models can not accurately predict the global equilibrium abundances when $M \gtrsim 1.$ To allow future studies to account for nonequilibrium effects in turbulent media, we gather our results into a series of tables, which we will extend in the future to encompass a wider range of elements, compositions, and ionizing processes.

We conduct simulations of turbulent mixing in the presence of a magnetic field, grown by the small-scale dynamo. We show that the scalar gradient field, $\nabla C$, which must be large for diffusion to operate, is strongly biased perpendicular to the magnetic field, ${\mathbf B}$. This is true both early-on, when the magnetic field is negligible, and at late times, when the field is strong enough to back react on the flow. This occurs because $\nabla C$ increases within the plane of a compressive motion, but ${\mathbf B}$ increases perpendicular to it. At late times the magnetic field resists compression, making it harder for scalar gradients to grow and likely slowing mixing.

Observations have revealed cold gas with large velocity dispersions (~300 km/s) within the hot outflows of ultra-luminous infrared galaxies (ULIRGs). This gas may trace its origin to the Rayleigh-Taylor (RT) fragmentation of a super-bubble or may arise on smaller scales. We model a ULIRG outflow at two scales to recreate this gas in three-dimensional hydrodynamic simulations using FLASH. Although resolution is limited, these models successfully produce cold gas in outflows with large velocity dispersions. Our small-scale models produce this cold gas through RT fragmentation of the super-bubble wall, but the large-scale models produce the cold gas after hot bubbles fragment the disc's gas into cold clouds which are then accelerated by thermal pressure, and supplemented by cooling within the outflow. We produce simple mock spectra to compare these simulations to observed absorption spectra and find line-widths of ~250 km/s, agreeing with the lower end of observations.

Turbulent motions are essential to the mixing of entrained fluids and are also capable of amplifying weak initial magnetic fields by small-scale dynamo action. Here we perform a systematic study of turbulent mixing in magnetized media, using three-dimensional magnetohydrodynamic simulations that include a scalar concentration field. We focus on how mixing depends on the magnetic Prandtl number, Pm, from 1 to 4 and the Mach number, M}, from 0.3 to 2.4. For all subsonic flows, we find that the velocity power spectrum has a k^-5/3 slope in the early, kinematic phase, but steepens due to magnetic back reactions as the field saturates. The scalar power spectrum, on the other hand, flattens compared to k^-5/3 at late times, consistent with the Obukohov-Corrsin picture of mixing as a cascade process. At higher Mach numbers, the velocity power spectrum also steepens due to the presence of shocks, and the scalar power spectrum again flattens accordingly. Scalar structures are more intermittent than velocity structures in subsonic turbulence while for supersonic turbulence, velocity structures appear more intermittent than the scalars only in the kinematic phase. Independent of the Mach number of the flow, scalar structures are arranged in sheets in both the kinematic and saturated phases of the magnetic field evolution. For subsonic turbulence, scalar dissipation is hindered in the strong magnetic field regions, probably due to Lorentz forces suppressing the buildup of scalar gradients, while for supersonic turbulence, scalar dissipation increases monotonically with increasing magnetic field strength. At all Mach numbers, mixing is significantly slowed by the presence of dynamically-important small-scale magnetic fields, implying that mixing in the interstellar medium and in galaxy clusters is less efficient than modeled in hydrodynamic simulations.

The early Universe hosted a large population of small dark matter `minihalos' that were too small to cool and form stars on their own. These existed as static objects around larger galaxies until acted upon by some outside influence. Outflows, which have been observed around a variety of galaxies, can provide this influence in such a way as to collapse, rather than disperse the minihalo gas. Gray & Scannapieco performed an investigation in which idealized spherically-symmetric minihalos were struck by enriched outflows. Here we perform high-resolution cosmological simulations that form realistic minihalos, which we then extract to perform a large suite of simulations of outflow-minihalo interactions including non-equilibrium chemical reactions. In all models, the shocked minihalo forms molecules through non-equilibrium reactions, and then cools to form dense chemically homogenous clumps of star-forming gas. The formation of these high-redshift clusters will be observable with the next generation of telescopes, and the largest of them should survive to the present day, having properties similar to halo globular clusters.

In the early universe, substantial relative "stream" velocities between the gas and dark matter arise due to radiation pressure and persist after recombination. To asses the impact of these velocities on high-redshift structure formation, we carry out a suite of high-resolution Adaptive Mesh Refinement (AMR) cosmological simulations, which use Smoothed Particle Hydrodynamic datasets as initial conditions, converted using a new tool developed for this work. These simulations resolve structures with masses as small as a few 100 M$_\odot$, and we focus on the $10^6$ M$_\odot$ "mini-halos" in which the first stars formed. At $z \approx 17,$ the presence of stream velocities has only a minor effect on the number density of halos below $10^6$ M$_\odot$, but it greatly suppresses gas accretion onto all halos and the dark matter structures around them. Stream velocities lead to significantly lower halo gas fractions, especially for $\approx 10^5$ M$_\odot$ objects, an effect that is likely to depend on the orientation of a halo's accretion lanes. This reduction in gas density leads to colder, more compact radial profiles, and it substantially delays the redshift of collapse of the largest halos, leading to delayed star formation and possibly delayed reionization. These many differences suggest that future simulations of early cosmological structure formation should include stream velocities to properly predict gas evolution, star-formation, and the epoch of reionization.

We conduct a comprehensive theoretical and numerical investigation of the pollution of pristine gas in turbulent flows, designed to provide new tools for modeling the evolution of the first generation of stars. The properties of such Population III (Pop III) stars are thought to be very different than later generations, because cooling is dramatically different in gas with a metallicity below a critical value Z_c, which lies between ~10^-6 and 10^-3 solar value. Z_c is much smaller than the typical average metallicity, <Z>, and thus the mixing efficiency of the pristine gas in the interstellar medium plays a crucial role in the transition from Pop III to normal star formation. The small critical value, Z_c, corresponds to the far left tail of the probability distribution function (PDF) of the metallicity. Based on closure models for the PDF formulation of turbulent mixing, we derive equations for the fraction of gas, P, lying below Z_c, in compressible turbulence. Our simulation data shows that the evolution of the fraction P can be well approximated by a generalized self-convolution model, which predicts dP/dt = -n/tau_con P (1-P^(1/n)), where n is a measure of the locality of the PDF convolution and the timescale tau_con is determined by the rate at which turbulence stretches the pollutants. Using a suite of simulations with Mach numbers ranging from M = 0.9 to 6.2, we provide accurate fits to n and tau_con as a function of M, Z_c/<Z>, and the scale, L_p, at which pollutants are added to the flow. For P>0.9, mixing occurs only in the regions surrounding the pollutants, such that n=1. For smaller P, n is larger as mixing becomes more global. We show how the results can be used to construct one-zone models for the evolution of Pop III stars in a single high-redshift galaxy, as well as subgrid models for tracking the evolution of the first stars in large cosmological simulations.

Intergalactic filaments form the foundation of the cosmic web that connect galaxies together, and provide an important reservoir of gas for galaxy growth and accretion. Here we present very high resolution two-dimensional simulations of the thermal and chemical evolution of such filaments, making use of a 32 species chemistry network that tracks the evolution of key molecules formed from hydrogen, oxygen, and carbon. We study the evolution of filaments over a wide range of parameters including the initial density, initial temperature, strength of the dissociating UV background, and metallicity. In low-redshift, $Z \approx 0.1 Z_\odot $ filaments, the evolution is determined completely by the initial cooling time. If this is sufficiently short, the center of the filament always collapses to form dense, cold core containing a substantial fraction of molecules. In high-redshift, $Z=10^{-3} Z_\odot$ filaments, the collapse proceeds much more slowly. This is due mostly to the lower initial temperatures, which leads to a much more modest increase in density before the atomic cooling limit is reached, making subsequent molecular cooling much less efficient. Finally, we study how the gravitational potential from a nearby dwarf galaxy affects the collapse of the filament and compare this to NGC 5253, a nearby starbusting dwarf galaxy thought to be fueled by the accretion of filament gas. In contrast to our fiducial case, a substantial density peak forms at the center of the potential. This peak evolves faster than the rest of the filament due to the increased rate at which chemical species form and cooling occur. We find that we achieve similar accretion rates as NGC 5253 but our two-dimensional simulations do not recover the formation of the giant molecular clouds that are seen in radio observations.

The interstellar medium is a multiphase gas in which turbulent support is as important as thermal pressure. Sustaining this configuration requires both continuous turbulent stirring and continuous radiative cooling to match the decay of turbulent energy. While this equilibrium can persist for small turbulent velocities, if the one-dimensional velocity dispersion is larger than approximately 35 km/s, the gas moves into an unstable regime that leads to rapid heating. I study the implications of this turbulent runaway, showing that it causes a hot gas outflow to form in all galaxies with a gas surface density above approximately 50 solar masses/pc^2 corresponding to a star formation rate per unit area of 0.1$ solar masses/yr/kpc^2. For galaxies with escape velocities above 200 km/s, the sonic point of this hot outflow should lie interior to the region containing cold gas and stars, while for galaxies with smaller escape velocities, the sonic point should lie outside this region. This leads to efficient cold cloud acceleration in higher mass galaxies, while in lower mass galaxies, clouds may be ejected by random turbulent motions rather than accelerated by the wind. Finally, I show that energy balance cannot be achieved at all for turbulent media above a surface density of approximately 10^5 solar masses/pc^2.

We report on deep near-infrared F125W (J) and F160W (H) Hubble Space Telescope Wide Field Camera 3 images of the z=6.42 quasar J1148+5251 to attempt to detect rest-frame near-ultraviolet emission from the host galaxy. These observations included contemporaneous observations of a nearby star of similar near-infrared colors to measure temporal variations in the telescope and instrument point spread function (PSF). We subtract the quasar point source using both this direct PSF and a model PSF. Using direct subtraction, we measure an upper limit for the quasar host galaxy of m_J>22.8, m_H>23.0 AB mag (2 sigma). After subtracting our best model PSF, we measure a limiting surface brightness from 0.3"-0.5" radius of mu_J > 23.5, mu_H > 23.7 AB magarc (2 sigma). We test the ability of the model subtraction method to recover the host galaxy flux by simulating host galaxies with varying integrated magnitude, effective radius, and S\'ersic index, and conducting the same analysis. These models indicate that the surface brightness limit (mu_J > 23.5 AB magarc) corresponds to an integrated upper limit of m_J > 22 - 23 AB mag, consistent with the direct subtraction method. Combined with existing far-infrared observations, this gives an infrared excess log(IRX) > 1.0 and corresponding ultraviolet spectral slope beta > -1.2\pm0.2. These values match those of most local luminous infrared galaxies, but are redder than those of almost all local star-forming galaxies and z~6 Lyman break galaxies.

Several lines of evidence, from isotopic analyses of meteorites to studies of the Sun's elemental and isotopic composition, indicate that the solar system was contaminated early in its evolution by ejecta from a nearby supernova (SN). Previous models have invoked SN material being injected into an extant protoplanetary disk, or isotropically expanding ejecta sweeping over a distant (>10 pc) cloud core, simultaneously enriching it and triggering its collapse. Here we consider a new astrophysical setting: the injection of clumpy SN ejecta, as observed in the Cas A SN remnant, into the molecular gas at the periphery of an HII region created by the SN's progenitor star. To track these interactions we have conducted a suite of high-resolution (1500^3 effective) 3D simulations that follow the evolution of individual clumps as they move into molecular gas. Even at these high resolutions, our simulations do not quite achieve numerical convergence, due to the challenge of properly resolving the small-scale mixing of ejecta and molecular gas, although they do allow some robust conclusions to be drawn. Isotropically exploding ejecta do not penetrate into the molecular cloud, but, if cooling is properly accounted for, clumpy ejecta penetrate to distances ~10^18 cm and mix effectively with star-forming molecular gas. The ~2 M_\odot high-metallicity ejecta from a core-collapse SN is likely to mix with ~2 \times 10^4 M_\odot of molecular gas. Thus all stars forming late (~5 Myr) in the evolution of an HII region may be contaminated by SN ejecta at a level ~10^-4. This level of contamination is consistent with the abundances of short-lived radionuclides and possibly some stable isotopic shifts in the early solar system, and is potentially consistent with the observed variability in stellar elemental abundances. SN contamination of forming planetary systems may be a common, universal process.

The first generation of stars had very different properties than later stellar generations, as they formed from a "pristine" gas that was free of heavy elements. Normal star formation took place only after the first stars polluted the surrounding turbulent interstellar gas, increasing its local heavy element concentration, Z, beyond a critical value, Z_c (10^-8 < Z_c <10^-5). Motivated by this astrophysical problem, we investigate the fundamental physics of the pollution of pristine fluid elements in isotropic compressible turbulence. Turbulence stretches the pollutants, produces concentration structures at small scales, and brings the pollutants and the unpolluted flow in closer contact. Our theoretical approach employs the probability distribution function (PDF) method for turbulent mixing. We adopt three PDF closure models and derive evolution equations for the pristine fraction from the models. To test and constrain the theoretical models, we conduct numerical simulations for decaying passive scalars in isothermal turbulent flows with Mach numbers of 0.9 and 6.2, and compute the mass fraction, P(Z_c, t), of the flow with Z < Z_c. In the Mach 0.9 flow, the evolution of P(Z_c, t)$ is well described by a continuous convolution model and dP(Z_c, t)/dt = P(Z_c, t) ln[P(Z_c, t)]/tau_con, if the mass fraction of the polluted flow is larger than ~ 0.1. If the initial pollutant fraction is smaller than ~ 0.1, an early phase exists during which the pristine fraction follows an equation from a nonlinear integral model: dP(Z_c, t)/dt = P(Z_c, t) [P(Z_c, t)-1]/tau_int. The timescales tau_con and tau_int are measured from our simulations. When normalized to the flow dynamical time, the decay of P(Z_ c, t) in the Mach 6.2 flow is slower than at Mach 0.9, and we show that P(Z_c, t) in the Mach 6.2 flow can be well fit using a formula from a generalized version of the self-convolution model.

We carry out a comprehensive smooth particle hydrodynamics simulation survey of double-degenerate white dwarf binary mergers of varying mass combinations in order to establish correspondence between initial conditions and remnant configurations. We find that all but one of our simulation remnants share general properties such as a cold, degenerate core surrounded by a hot disk, while our least massive pair of stars forms only a hot disk. We characterize our remnant configurations by the core mass, the rotational velocity of the core, and the half-mass radius of the disk. We also find that some of our simulations with very massive constituent stars exhibit helium detonations on the surface of the primary star before complete disruption of the secondary. However, these helium detonations are insufficiently energetic to ignite carbon, and so do not lead to prompt carbon detonations.

The interstellar medium in star-forming galaxies is a multiphase gas in which turbulent support is at least as important as thermal pressure. Sustaining this configuration requires continuous radiative cooling, such that the overall average cooling rate matches the decay rate of turbulent energy into the medium. Here we carry out a set of numerical simulations of a stratified, turbulently stirred, radiatively cooled medium, which uncover a fundamental transition at a critical one-dimensional turbulent velocity of ~ 35 km/s. At turbulent velocities below ~35 km/s, corresponding to temperatures below 300,000 K, the medium is stable, as the time for gas to cool is roughly constant as a function of temperature. On the other hand, at turbulent velocities above the critical value, the gas is shocked into an unstable regime in which the cooling time increases strongly with temperature, meaning that a substantial fraction of the interstellar medium is unable to cool on a turbulent dissipation timescale. This naturally leads to runaway heating and ejection of gas from any stratified medium with a one-dimensional turbulent velocity above ~35 km/s, a result that has implications for galaxy evolution at all redshifts.

Rapid mass assembly, likely from mergers or smooth accretion, has been predicted to play a vital role in star-formation in high-redshift Lyman-alpha (Lya) emitters. Here we predict the major merger, minor merger, and smooth accreting Lya emitter fraction from z~3 to z~7 using a large dark matter simulation, and a simple physical model that is successful in reproducing many observations over this large redshift range. The central tenet of this model, different from many of the earlier models, is that the star-formation in Lya emitters is proportional to the mass accretion rate rather than the total halo mass. We find that at z~3, nearly 35% of the Lya emitters accrete their mass through major (3:1) mergers, and this fraction increases to about 50% at z~7. This imply that the star-formation in a large fraction of high-redshift Lya emitters is driven by mergers. While there is discrepancy between the model predictions and observed merger fractions, some of this difference (~15%) can be attributed to the mass-ratio used to define a merger in the simulation. We predict that future, deeper observations which use a 3:1 definition of major mergers will find >30% major merger fraction of Lya emitters at redshifts >3.

The early universe hosted a large population of low-mass virialized "minihalos," that were not massive enough to form stars on their own. While most minihalos were photoevaporated by ionizing photons from star-forming galaxies, these galaxies also drove large outflows, which in some cases would have reached the minihalos in advance of ionization fronts. In the previous papers in this series, we carried out high-resolution, three-dimensional adaptive mesh refinement simulations of outflow-minihalo interactions that included non-equilibrium chemistry, radiative cooling, and turbulent mixing. We found that, for a fiducial set of parameters, minihalos were transformed into dense, chemically homogenous stellar clusters. Here we conduct a suite of simulations that follow these interactions over a wide range of parameters including minihalo mass, minihalo formation redshift, outflow energy, outflow redshift, distance, concentration, and spin. In almost all cases, the shocked minihalos form molecules through nonequillibrium reactions and then cool rapidly to become compact, chemically-homogenous stellar clusters. Furthermore, we show that the unique properties of these clusters make them a prime target for direct study with the next generation of telescopes, and that there are many reasons to suspect that their low-redshift counterparts are the observed population of halo globular clusters.

In the primordial universe, low mass structures with virial temperatures less than 10$^{4}$ K were unable to cool by atomic line transitions, leading to a strong suppression of star formation. On the other hand, these "minihalos" were highly prone to triggered star formation by interactions from nearby galaxy outflows. In Gray & Scannapieco (2010), we explored the impact of nonequilibrium chemistry on these interactions. Here we turn our attention to the role of metals, carrying out a series of high-resolution three-dimensional adaptive mesh refinement simulations that include both metal cooling and a subgrid turbulent mixing model. Despite the presence of an additional coolant, we again we find that outflow-minihalo interactions produce a distribution of dense, massive stellar clusters. We also find that these clusters are evenly enriched with metals to a final abundance of Z $\approx$ 10$^{-2}$ Z$_{\odot}$. As in our previous simulations, all of these properties suggest that these interactions may have given rise to present-day halo globular clusters.

We conduct a systematic numerical study of passive scalar structures in supersonic turbulent flows. We find that the degree of intermittency in the scalar structures increases only slightly as the flow changes from transonic to highly supersonic, while the velocity structures become significantly more intermittent. This difference is due to the absence of shock-like discontinuities in the scalar field. The structure functions of the scalar field are well described by the intermittency model of She and L\'{e}v\^{e}que [Phys. Rev. Lett. 72, 336 (1994)], and the most intense scalar structures are found to be sheet-like at all Mach numbers.

The temperature profile of hot gas in galaxies and galaxy clusters is largely determined by the depth of the total gravitational potential and thereby by the dark matter (DM) distribution. We use high-resolution hydrodynamical simulations of galaxy formation to derive a surprisingly simple relation between the gas temperature and DM properties. We show that this relation holds not just for galaxy clusters but also for equilibrated and relaxed galaxies at radii beyond the central stellar-dominated region of typically a few kpc. It is then clarified how a measurement of the temperature and density of the hot gas component can lead to an indirect measurement of the DM velocity anisotropy in galaxies. We also study the temperature relation for galaxy clusters in the presence of self-regulated, recurrent active galactic nuclei (AGN), and demonstrate that this temperature relation even holds outside the inner region of 30 kpc in clusters with an active AGN.

We present a comprehensive study of white dwarf collisions as an avenue for creating type Ia supernovae. Using a smooth particle hydrodynamics code with a 13-isotope, {\alpha}-chain nuclear network, we examine the resulting 56Ni yield as a function of total mass, mass ratio, and impact parameter. We show that several combinations of white dwarf masses and impact parameters are able to produce sufficient quantities of 56Ni to be observable at cosmological distances. We find the 56Ni production in double-degenerate white dwarf collisions ranges from sub-luminous to the super-luminous, depending on the parameters of the collision. For all mass pairs, collisions with small impact parameters have the highest likelihood of detonating, but 56Ni production is insensitive to this parameter in high-mass combinations, which significantly increases their likelihood of detection. We also find that the 56Ni dependence on total mass and mass ratio is not linear, with larger mass primaries producing disproportionately more 56Ni than their lower mass secondary counterparts, and symmetric pairs of masses producing more 56Ni than asymmetric pairs.

We present a simple physical model for populating dark matter halos with Lyman Alpha Emiiters(LAEs) and predict the physical properties of LAEs at z~3-7. The central tenet of this model is that the Ly-alpha luminosity is proportional to the star formation rate (SFR) which is directly related to the halo mass accretion rate. The only free parameter in our model is then the star-formation efficiency (SFE). An efficiency of 2.5% provides the best-fit to the Ly-alpha luminosity function (LF) at redshift z=3.1, and we use this SFE to construct Ly-alpha LFs at other redshifts. Our model reproduce the Ly-alpha LFs, stellar ages, SFR ~1-10; Msun/yr, stellar masses ~ 10^7-10^8 Msun and the clustering properties of LAEs at z~3-7. We find the spatial correlation lengths ro ~ 3-6 Mpc/h, in agreement with the observations. Finally, we estimate the field-to-field variation ~ 30% for current volume and flux limited surveys, again consistent with observations. Our results suggest that the star formation, and hence Ly-alpha emission in LAEs is powered by the accretion of new material, and that the physical properties of LAEs do not evolve significantly over a wide range of redshifts. Relating the accreted mass, rather than the total mass of halos, to the Ly-alpha luminosity of LAEs naturally gives rise to the duty cycle of LAEs.

In many astrophysical environments, mixing of heavy elements occurs in the presence of a supersonic turbulent velocity field. Here we carry out the first systematic numerical study of such passive scalar mixing in isothermal supersonic turbulence. Our simulations show that the ratio of the scalar mixing timescale, $\tau_{\rm c}$, to the flow dynamical time, $\tau_{\rm dyn}$ (defined as the flow driving scale divided by the rms velocity), increases with the Mach number, $M$, for $M \lsim3$, and becomes essentially constant for $M \gsim3.$ This trend suggests that compressible modes are less efficient in enhancing mixing than solenoidal modes. However, since the majority of kinetic energy is contained in solenoidal modes at all Mach numbers, the overall change in $\tau_{\rm c}/\tau_{\rm dyn}$ is less than 20\% over the range $1 \lsim M \lsim 6$. At all Mach numbers, if pollutants are injected at around the flow driving scale, $\tau_{\rm c}$ is close to $\tau_{\rm dyn}.$ This suggests that scalar mixing is driven by a cascade process similar to that of the velocity field. The dependence of $\tau_{\rm c}$ on the length scale at which pollutants are injected into flow is also consistent with this cascade picture. Similar behavior is found for the variance decay timescales for scalars without continuing sources. Extension of the scalar cascade picture to the supersonic regime predicts a relation between the scaling exponents of the velocity and the scalar structure functions, with the scalar structure function becoming flatter as the velocity scaling steepens with Mach number. Our measurements of the volume-weighted velocity and scalar structure functions confirm this relation for $M\lsim 2,$ but show discrepancies at $M \gsim 3$.

We present tomography of the circum-galactic metal distribution at redshift 1.7 to 4.5 derived from echellete spectroscopy of binary quasars. We find CIV systems at similar redshifts in paired sightlines more often than expected for sightline-independent redshifts. As the separation of the sightlines increases from 36 kpc to 907 kpc, the amplitude of this clustering decreases. At the largest separations, the CIV systems cluster similar to Lyman-break galaxies (Adelberger et al. 2005a). The CIV systems are significantly less correlated than these galaxies, however, at separations less than R_1 ~ 0.42 +/- 0.15 h-1 comoving Mpc. Measured in real space, i.e., transverse to the sightlines, this length scale is significantly smaller than the break scale estimated from the line-of-sight correlation function in redshift space (Scannapieco et al. 2006a). Using a simple model, we interpret the new real-space measurement as an indication of the typical physical size of enriched regions. We adopt this size for enriched regions and fit the redshift-space distortion in the line-of-sight correlation function. The fitted velocity kick is consistent with the peculiar velocity of galaxies as determined by the underlying mass distribution and places an upper limit on the outflow (or inflow) speed of metals. The implied time scale for dispersing metals is larger than the typical stellar ages of Lyman-break galaxies (Shapley et al. 2001), and we argue that enrichment by galaxies at z > 4.3 played a greater role in dispersing metals. To further constrain the growth of enriched regions, we discuss empirical constraints on the evolution of the CIV correlation function with cosmic time. This study demonstrates the potential of tomography for measuring the metal enrichment history of the circum-galactic medium.

We use high-resolution three-dimensional adaptive mesh refinement simulations to investigate the interaction of high-redshift galaxy outflows with low-mass virialized clouds of primordial composition. While atomic cooling allows star formation in objects with virial temperatures above $10^4$ K, "minihaloes" below this threshold are generally unable to form stars by themselves. However, these objects are highly susceptible to triggered star formation, induced by outflows from neighboring high-redshift starburst galaxies. Here we conduct a study of these interactions, focusing on cooling through non-equilibrium molecular hydrogen (H$_2$) and hydrogen deuteride (HD) formation. Tracking the non-equilibrium chemistry and cooling of 14 species and including the presence of a dissociating background, we show that shock interactions can transform minihaloes into extremely compact clusters of coeval stars. Furthermore, these clusters are all less than $\approx 10^6 M_\odot,$ and they are ejected from their parent dark matter halos: properties that are remarkably similar to those of the old population of globular clusters.

During the early evolution of an AM CVn system, helium is accreted onto the surface of a white dwarf under conditions suitable for unstable thermonuclear ignition. The turbulent motions induced by the convective burning phase in the He envelope become strong enough to influence the propagation of burning fronts and may result in the onset of a detonation. Such an outcome would yield radioactive isotopes and a faint rapidly rising thermonuclear ".Ia" supernova. In this paper, we present hydrodynamic explosion models and observable outcomes of these He shell detonations for a range of initial core and envelope masses. The peak UVOIR bolometric luminosities range by a factor of 10 (from 5e41 - 5e42 erg/s), and the R-band peak varies from M_R,peak = -15 to -18. The rise times in all bands are very rapid (<10 d), but the decline rate is slower in the red than the blue due to a secondary near-IR brightening. The nucleosynthesis primarily yields heavy alpha-chain elements (40Ca through 56Ni) and unburnt He. Thus, the spectra around peak light lack signs of intermediate mass elements and are dominated by CaII and TiII features, with the caveat that our radiative transfer code does not include the non-thermal effects necessary to produce He features.

We present three-dimensional, adaptive mesh simulations of dwarf galaxy out- flows driven by supersonic turbulence. Here we develop a subgrid model to track not only the thermal and bulk velocities of the gas, but also its turbulent velocities and length scales. This allows us to deposit energy from supernovae directly into supersonic turbulence, which acts on scales much larger than a particle mean free path, but much smaller than resolved large-scale flows. Unlike previous approaches, we are able to simulate a starbursting galaxy modeled after NGC 1569, with realistic radiative cooling throughout the simulation. Pockets of hot, diffuse gas around individual OB associations sweep up thick shells of material that persist for long times due to the cooling instability. The overlapping of high-pressure, rarefied regions leads to a collective central outflow that escapes the galaxy by eating away at the exterior gas through turbulent mixing, rather than gathering it into a thin, unstable shell. Supersonic, turbulent gas naturally avoids dense regions where turbulence decays quickly and cooling times are short, and this further enhances density contrasts throughout the galaxy- leading to a complex, chaotic distribution of bubbles, loops and filaments as observed in NGC 1569 and other outflowing starbursts.

We study the age and metallicity distribution function (MDF) of metal-poor stars in the Milky Way halo as a function of galactocentric radius by combining N-body simulations and semi-analytical methods. We find that the oldest stars populate the innermost region, while extremely metal-poor stars are more concentrated within r < 60 kpc. The MDF of [Fe/H] < -2 stars varies only very weakly within the central 50 kpc, while the relative contribution of [Fe/H] < -2 stars strongly increases with r, varying from 16% within 7 kpc < r < 20 kpc up to > 40% for r > 20 kpc. This is due to the faster descent of the spatial distribution (as seen from Earth) of the more enriched population. This implies that the outer halo < 40 kpc is the best region to search for very metal-poor stars. Beyond ~ 60 kpc the density of [Fe/H] < -2 stars is maximum within dwarf galaxies. All these features are imprinted by a combination of (i) the virialization epoch of the star-forming haloes, and (ii) the metal enrichment history of the Milky Way environment.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

The time delay between the formation of a population of stars and the onset of type Ia supernovae (SNe Ia) sets important limits on the masses and nature of SN Ia progenitors. Here we use a new observational technique to measure this time delay by comparing the spatial distributions of SNe Ia to their local environments. Previous work attempted such analyses encompassing the entire host of each SN Ia, yielding inconclusive results. Our approach confines the analysis only to the relevant portions of the hosts, allowing us to show that even so-called "prompt" SNe Ia that trace star-formation on cosmic timescales exhibit a significant delay time of 200-500 million years. This implies that either the majority of Ia companion stars have main-sequence masses less than 3 solar masses, or that most SNe Ia arise from double-white dwarf binaries. Our results are also consistent with a SNe Ia rate that traces the white dwarf formation rate, scaled by a fixed efficiency factor.

The physics of X-ray cavities in galaxy clusters is constrained by their observed morphological evolution, which depends on such poorly-understood properties as the turbulent density field and magnetic fields. Here we combine numerical simulations that include subgrid turbulence and software that produces synthetic X-ray observations to examine the evolution of X-ray cavities in the the absence of magnetic fields. Our results reveal an anisotropic size evolution that is very different from simplified, analytical predictions. These differences highlight some of the key issues that must be accurately quantified when studying AGN-driven cavities, and help to explain why the inferred pV energy in these regions appears to be correlated with their distance from the cluster center. Interpreting X-ray observations will require detailed modeling of effects including mass-entrainment, distortion by drag forces, and projection. Current limitations do not allow a discrimination between purely hydrodynamic and magnetically-dominated models for X-ray cavities.

Cool cores of galaxy clusters are thought to be heated by low-power active galactic nuclei (AGN), whose accretion is regulated by feedback. However, the interaction between the hot gas ejected by the AGN and the ambient intracluster medium is extremely difficult to simulate, as it involves a wide range of spatial scales and gas that is Rayleigh-Taylor (RT) unstable. Here we use a subgrid model for RT-driven turbulence to overcome these problems and present the first observationally-consistent hydrodynamical simulations of AGN self-regulation in galaxy clusters. For a wide range of parameter choices the cluster in our three-dimensional simulations regulates itself for at least several Gyrs years. Heating balances cooling through a string of outbreaks with a typical recurrence time of approximately 80 Myrs, a timescale that depends only on the global cluster properties.

We explore collisions between two white dwarfs as a pathway for making Type Ia Supernovae (SNIa). White dwarf number densities in globular clusters allow 10-100 redshift <1 collisions per year, and observations by (Chomiuk et al. 2008) of globular clusters in the nearby S0 galaxy NGC 7457 have detected what is likely to be a SNIa remnant. We carry out simulations of the collision between two 0.6 solar mass white dwarfs at various impact parameters and mass resolutions. For impact parameters less than half the radius of the white dwarf, we find such collisions produce approximately 0.4 solar masses of Ni56, making such events potential candidates for underluminous SNIa or a new class of transients between Novae and SNIa.

Astronomers have proposed a number of mechanisms to produce supernova explosions. Although many of these mechanisms are now not considered primary engines behind supernovae, they do produce transients that will be observed by upcoming ground-based surveys and NASA satellites. Here we present the first radiation-hydrodynamics calculations of the spectra and light curves from three of these "failed" supernovae: supernovae with considerable fallback, accretion induced collapse of white dwarfs, and energetic helium flashes (also known as type .Ia supernovae).

The first objects to arise in a cold dark matter universe present a daunting challenge for models of structure formation. In the ultra small-scale limit, CDM structures form nearly simultaneously across a wide range of scales. Hierarchical clustering no longer provides a guiding principle for theoretical analyses and the computation time required to carry out credible simulations becomes prohibitively high. To gain insight into this problem, we perform high-resolution (N=720^3 - 1584^3) simulations of an Einstein-de Sitter cosmology where the initial power spectrum is P(k) propto k^n, with -2.5 < n < -1. Self-similar scaling is established for n=-1 and n=-2 more convincingly than in previous, lower-resolution simulations and for the first time, self-similar scaling is established for an n=-2.25 simulation. However, finite box-size effects induce departures from self-similar scaling in our n=-2.5 simulation. We compare our results with the predictions for the power spectrum from (one-loop) perturbation theory and demonstrate that the renormalization group approach suggested by McDonald improves perturbation theory's ability to predict the power spectrum in the quasilinear regime. In the nonlinear regime, our power spectra differ significantly from the widely used fitting formulae of Peacock & Dodds and Smith et al. and a new fitting formula is presented. Implications of our results for the stable clustering hypothesis vs. halo model debate are discussed. Our power spectra are inconsistent with predictions of the stable clustering hypothesis in the high-k limit and lend credence to the halo model. Nevertheless, the fitting formula advocated in this paper is purely empirical and not derived from a specific formulation of the halo model.

We study the statistical properties of the cosmological 21 cm signal from both the intergalactic medium (IGM) and minihalos, using a reionization simulation that includes a self--consistent treatment of minihalo photoevaporation. We consider two models for minihalo formation and three typical thermal states of the IGM -- heating purely by ionization, heating from both ionizing and $\rm Ly\alpha$ photons, and a maximal "strong heating" model. We find that the signal from the IGM is almost always dominant over that from minihalos. In our calculation, the differential brightness temperature, ${\delta}T_b,$ of minihalos is never larger than 2 mK. Although there are indeed some differences in the signals from the minihalos and from the IGM, even with the planned generation of radio telescopes it will be unfeasible to detect them. However, minihalos significantly affect the ionization state of the IGM and the corresponding 21 cm flux.

Cool cores of galaxy clusters are thought to be heated by low-power active galactic nuclei (AGN), whose accretion is regulated by feedback. However, the interaction between the hot gas ejected by the AGN and the ambient intracluster medium is extremely difficult to simulate as it involves a wide range of spatial scales and gas that is Rayleigh-Taylor (RT) unstable. Here we present a series of three-dimensional hydrodynamical simulations of a self-regulating AGN in a galaxy cluster. Our adaptive-mesh simulations include prescriptions for radiative cooling, AGN heating and a subgrid model for RT-driven turbulence, which is crucial to simulate this evolution. AGN heating is taken to be proportional to the rest-mass energy that is accreted onto the central region of the cluster. For a wide range of feedback efficiencies, the cluster regulates itself for at least several $10^9$ years. Heating balances cooling through a string of outbursts with typical recurrence times of around 80 Myrs, a timescale that depends only on global cluster properties. Under certain conditions we find central dips in the metallicity of the intracluster medium. Provided the sub-grid model used here captures all its key properties, turbulence plays an essential role in the AGN self-regulation in cluster cores.

We explore the dependence of the subhalo mass function on the spectral index n of the linear matter power spectrum using scale-free Einstein-de Sitter simulations with n=-1 and n=-2.5. We carefully consider finite volume effects that may call into question previous simulations of n<-2 power spectra. Subhaloes are found using a 6D friends-of-friends algorithm in all haloes originating from high-sigma peaks. For n=-1, we find that the cumulative subhalo mass function is independent of the parameters used in the subhalo finding algorithm and is consistent with the subhalo mass function found in LCDM simulations. In particular, the subhalo mass function is well fit by a power-law with an index of alpha=-0.9, that is the mass function has roughly equal mass in subhaloes per logarithmic interval in subhalo mass. Conversely, for n=-2.5, the algorithm parameters affect the subhalo mass function since subhaloes are more triaxial with less well defined boundaries. We find that the index alpha is generally larger with alpha>=-0.75. We infer that although the subhalo mass function appears to be independent of n so long as n>=-2, it begins to flatten as n->-3. Thus, the common practice of using alpha=-1.0 may greatly overestimate the number of subhaloes at the smallest scales in the CDM hierarchy.

Fundamental mysteries remain regarding the physics of Type Ia supernovae (SNIa) and their stellar progenitors. We argue here that important clues to these questions may emerge by the identification of those SNIa that occur in extragalactic globular clusters--stellar systems with well defined ages and metallicities. We estimate an all-sky rate of approximately 0.1 eta (D/100 Mpc)^3 per year for SNIa in globular clusters within a distance D, where eta is the rate enhancement per unit mass as a result of dynamical production channels that are inaccessible in the galactic field. If eta is approximately 2-10, as suggested by observations and theory, the combined efforts of accurate supernova astrometry and deep follow-up imaging should identify the > 1% of nearby (D < 100 Mpc) SNIa that occur in globular clusters.

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

A wide range of recent observations have shown that AGN-driven cavities may provide the energy source that balances the cooling observed in the centres of cool-core galaxy clusters. One tool for better understanding the physics of these cavities is their observed morphological evolution, which is dependent on such poorly-understood properties as the turbulent density field and the impact of magnetic fields. Here we combine numerical simulations that include subgrid turbulence and software that produces synthetic X-ray observations to examine the evolution of X-ray cavities in the absence of magnetic fields. Our results reveal an anisotropic size evolution of that is dramatically different from simplified, analytical predictions. These differences highlight some of the key issues that must be accurately quantified when studying AGN-driven cavities, and help to explain why the inferred pV energy in these regions appears to be correlated with their distance from the cluster center. Interpreting X- ray observations will require detailed modeling of effects including mass-entrainment, distortion by drag forces, and pro jection. Current limitations do not allow a discrimination between purely hydrodynamic and magnetically-dominated models for X-ray cavities.

We show that current clustering observations of quasars and luminous AGN can be explained by a merger model augmented by feedback from outflows. Using numerical simulations large enough to study clustering out to 25 comoving h^{-1} Mpc, we calculate correlation functions, biases, and correlation lengths as a function of AGN redshift and optical and X-ray luminosity. At optical wavelengths, our results match a wide range of current observations and generate predictions for future data sets. We reproduce the weak luminosity dependence of clustering over the currently well-measured range, and predict a much stronger dependence at higher luminosities. The increase in the amplitude of binary quasar clustering observed in the Sloan Digital Sky Survey (SDSS) is also reproduced and is predicted to occur at higher redshift, an effect that is due to the one halo term in the correlation function. On the other hand, our results do not match the rapid evolution of the correlation length observed in the SDSS at z\simeq 3, a discrepancy that is at least partially due to differences in the scales probed by our simulation versus this survey. In fact, we show that changing the distances sampled from our simulations can produce changes as large as 40% in the fitted correlation lengths. Finally, in the X-ray, our simulations produce correlation lengths similar to that observed in the Chandra Deep Field (CDF) North, but not the significantly larger correlation length observed in the CDF South.

We carry out a comprehensive theoretical examination of the relationship between the spatial distribution of optical transients and the properties of their progenitor stars. By constructing analytic models of star-forming galaxies and the evolution of stellar populations within them, we are able to place constraints on candidate progenitors for core-collapse supernovae (SNe), long-duration gamma ray bursts, and supernovae Ia. In particular we first construct models of spiral galaxies that reproduce observations of core-collapse SNe, and we use these models to constrain the minimum mass for SNe Ic progenitors to approximately 25 solar masses. Secondly, we lay out the parameters of a dwarf irregular galaxy model, which we use to show that the progenitors of long-duration gamma-ray bursts are likely to have masses above approximately 43 solar masses. Finally, we introduce a new method for constraining the time scale associated with SNe Ia and apply it to our spiral galaxy models to show how observations can better be analyzed to discriminate between the leading progenitor models for these objects.

Hot, underdense bubbles powered by active galactic nuclei (AGN) are likely to play a key role in halting catastrophic cooling in the centers of cool-core galaxy clusters. We present three-dimensional simulations that capture the evolution of such bubbles, using an adaptive-mesh hydrodynamic code, FLASH3, to which we have added a subgrid model of turbulence and mixing. While pure-hydro simulations indicate that AGN bubbles are disrupted into resolution-dependent pockets of underdense gas, proper modeling of subgrid turbulence indicates that this a poor approximation to a turbulent cascade that continues far beyond the resolution limit. Instead, Rayleigh-Taylor instabilities act to effectively mix the heated region with its surroundings, while at the same time preserving it as a coherent structure, consistent with observations. Thus bubbles are transformed into hot clouds of mixed material as they move outwards in the hydrostatic intracluster medium (ICM), much as large airbursts lead to a distinctive ``mushroom cloud'' structure as they rise in the hydrostatic atmosphere of Earth. Properly capturing the evolution of such clouds has important implications for many ICM properties. In particular, it significantly changes the impact of AGN on the distribution of entropy and metals in cool-core clusters such as Perseus.

One of the most important and poorly-understood issues in structure formation is the role of outflows driven by active galactic nuclei (AGN). Using large-scale cosmological simulations, we compute the impact of such outflows on the small-scale distribution of the cosmic microwave background (CMB). Like gravitationally-heated structures, AGN outflows induce CMB distortions both through thermal motions and peculiar velocities, by processes known as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects, respectively. For AGN outflows the thermal SZ effect is dominant, doubling the angular power spectrum on arcminute scales. But the most distinct imprint of AGN feedback is a substantial increase in the thermal SZ distortions around elliptical galaxies, post-starburst ellipticals, and quasars, which is linearly proportional to the outflow energy. While point source subtraction is difficult for quasars, we show that by appropriately stacking microwave measurements around early-type galaxies, the new generation of small-scale microwave telescopes will be able to directly measure AGN feedback at the level important for current theoretical models.

The difference in density profiles of the contributions from different density peaks to dark matter halos results in certain expectations about the Milky Way's stellar halo. We cut our simulated halo stars into two populations: those forming before/during the last major merger, and those accreted after the last major merger. The former population are more centrally located at z=0, while stars forming in low mass late forming proto-galaxies are spread through the halo. A difference in observed binding energy distinguishes these two populations. We look at possible chemical abundance signatures of the two populations. We also show that galaxies forming in isolated low sigma peaks will form from primordial material. Thus, even though the oldest stars are centrally concentrated as they originated in the early collapsing, densest regions, primordial stars would be found distributed throughout the halo. Thus, the lack of observed metal free stars can be taken as directly constraining the Population III IMF, and the lowest metallicity observed stars can be interpreted as holding clues to the chemical yields of Pop III stars.

We use cosmological, chemo-dynamical, smoothed particle hydrodynamical simulations of Milky-Way-analogue galaxies to find the expected present-day distributions of both metal-free stars that formed from primordial gas and the oldest star populations. We find that metal-free stars continue to form until z~4 in halos that are chemically isolated and located far away from the biggest progenitor of the final system. As a result, if the Population III initial mass function allows stars with low enough mass to survive until z=0 (< 0.8 Msol), they would be distributed throughout the Galactic halo. On the other hand, the oldest stars form in halos that collapsed close to the highest density peak of the final system, and at z=0 they are located preferentially in the central region of the Galaxy, i.e., in the bulge. According to our models, these trends are not sensitive to the merger histories of the disk galaxies or the implementation of supernova feedback. Furthermore, these full hydrodynamics results are consistent with our N-body results in Paper I, and lend further weight to the conclusion that surveys of low-metallicity stars in the Galactic halo can be used to directly constrain the properties of primordial stars. In particular, they suggest that the current lack of detections of metal-free stars implies that their lifetimes were shorter than a Hubble time, placing constraints on the metal-free initial mass function.

We study the spatial distribution of Galactic metal-free stars by combining an extremely high-resolution (7.8 X 10^5 solar masses per particle) Cold Dark Matter N-body simulation of the Milky-Way with a semi-analytic model of metal enrichment. This approach allows us to resolve halos with virial temperatures down to the 10^4K atomic cooling limit, and it is sufficiently flexible to make a number of robust conclusions, despite the extremely uncertain properties of the first stars. Galactic metal-free stars are formed over a large redshift range, which peaks at z~10, but continues down to z~5, contributing stars at wide range of Galactocentric radii. Stars containing only metals generated by primordial stars are similarly widespread. Neither changing the efficiency of metal dispersal by two orders of magnitude, nor drastically changing the approximations in our semi-analytical model can affect these result. Thus, if they have sufficiently long lifetimes, a significant number of stars formed in initially primordial star clusters should be found in the nearby Galactic halo regardless of the specifics of metal-free star formation. Observations of metal abundances in Galactic halo stars should be taken as directly constraining the properties of primordial stars, and the lack of metal-free halo stars today should be taken as strongly suggesting a 0.8 solar mass lower limit on the primordial initial mass function.

I review the physical properties of pair-production supernovae (PPSNe) as well as the prospects for them to be constrained observationally. In very massive (140-260 solar mass) stars, much of the pressure support comes from the radiation field, meaning that they are loosely bound, with an adiabatic coefficient that is close to the minimum stable value. Near the end of C/O burning, the central temperature increases to the point that photons begin to be converted into electron-positron pairs, softening gamma below this critical value. The result is a runaway collapse, followed by explosive burning that completely obliterates the star. While these explosions can be up to 100 times more energetic that core collapse and Type Ia supernovae, their peak luminosities are only slightly greater. However, due both to copious Ni-56 production and hydrogen recombination, they are brighter much longer, and remain observable for ~ 1 year. Since metal enrichment is a local process, PPSNe should occur in pockets of metal-free gas over a broad range of redshifts, greatly enhancing their detectability, and distributing their nucleosyntehtic products about the Milky Way. This means that measurements of the abundances of metal-free stars should be thought of as directly constraining these objects. It also means that ongoing supernova searches, already provide weak constraints for PPSN models. A survey with the NIRCam instrument on JWST, on the other hand, would be able to extend these limits to z ~ 10. Observing a 0.3 deg^2 patch of sky for one week per year for three consecutive years, such a program would either detect or rule out the existence of these remarkable objects.

(Abridged) We explore the idea that the anti-hierarchical turn-off observed in the quasar population arises from self-regulating feedback, via an outflow mechanism. Using a detailed hydrodynamic simulation we calculate the luminosity function of quasars down to a redshift of z=1 in a large, cosmologically representative volume. Outflows are included explicitly by tracking halo mergers and driving shocks into the surrounding intergalactic medium. Our results are in excellent agreement with measurements of the spatial distribution of quasars, and we detect an intriguing excess of galaxy-quasar pairs at very short separations. We also reproduce the anti-hierarchical turnoff in the quasar luminosity function, however, the magnitude of the turn-off falls short of that observed as well as that predicted by analogous semi-analytic models. The difference can be traced to the treatment of gas heating within galaxies. The simulated galaxy cluster L_X-T relationship is close to that observed for z~1 clusters, but the simulated galaxy groups at z=1 are significantly perturbed by quasar outflows, suggesting that measurements of X-ray emission in high-redshift groups could well be a "smoking gun" for the AGN heating hypothesis.

We derive the equations for the propagation of relativistic ionization fronts in both static and moving gases. We focus on the supersonic R-type phase that occurs right after a source turns on, and we compare the nonrelativistic and relativistic solutions for several important cases. Relativistic corrections can be significant up until the light-crossing time of the equilibrium Stromgren sphere. For a static medium, we obtain exact analytical solutions and apply them to the illustrative problems of an O star in a molecular cloud and a starburst in a high-redshift cosmological halo. Relativistic corrections can be important at early times when the H II regions are small, as well as at later times, if a density gradient causes the I-front to accelerate. For the cosmologically-expanding IGM, we derive an analytical solution in the case of a steady source and a constant clumping factor. Here relativistic corrections are significant for short-lived, highly-luminous sources like QSOs at the end of reionization, but negligible for weaker or higher-redshift sources. Finally, we numerically calculate the evolution of relativistic I-fronts in the presence of small-scale structure and infall, for a large galaxy undergoing a starburst and a luminous, high-redshift QSO. For such strong and short-lived sources, the relativistic corrections are quite significant, and small-scale structure can decrease the size of the H II region by up to an additional ~25%. (abridged)

One of the most debated issues in the theoretical modeling of cosmic reionization is the impact of small-mass gravitationally-bound structures. We carry out the first numerical investigation of the role of such sterile `minihaloes', which serve as self-shielding screens of ionizing photons. Minihaloes are too small to be properly resolved in current large-scale cosmological simulations, and thus we estimate their effects using a sub-grid model, considering two cases that bracket their effect within this framework. In the `extreme suppression' case in which minihalo formation ceases once a region is partially ionized, their effect on cosmic reionization is modest, reducing the volume-averaged ionization fraction by an overall factor of less than 15%. In the other extreme, in which minihalo formation is never suppressed, they delay complete reionization as much as Delta z~2, in rough agreement with the results from a previous semi-analytical study by the authors. Thus, depending on the details of the minihalo formation process, their effect on the overall progress of reionization can range from modest to significant, but the minihalo photon consumption is by itself insufficient to force an extended reionization epoch.

We study the impact of outflows driven by active galactic nuclei (AGN) on galaxy formation. Outflows move into the surrounding intergalactic medium (IGM) and heat it sufficiently to prevent it from condensing onto galaxies. In the dense, high-redshift IGM, such feedback requires highly energetic outflows, driven by a large AGN. However, in the more tenuous low-redshift IGM, equivalently strong feedback can be achieved by less energetic winds (and thus smaller galaxies). Using a simple analytic model, we show that this leads to the anti-hierarchical quenching of star-formation in large galaxies, consistent with current observations. At redshifts prior to the formation of large AGN, galaxy formation is hierarchical and follows the growth of dark-matter halos. The transition between the two regimes lies at the z ~ 2 peak of AGN activity.

The majority of the heavy elements produced by stars 2 billion years after the Big Bang (redshift z~3) are presently undetected at those epochs. We propose a solution to this cosmic `missing metals' problem in which such elements are stored in gaseous halos produced by supernova explosions around star-forming galaxies. By using data from the ESO/VLT Large Program, we find that:(i) only 5%-9% of the produced metals reside in the cold phase, the rest being found in the hot (log T=5.8-6.4) phase; (ii) 1%-6% (3%-30%) of the observed CIV (OVI) is in the hot phase. We conclude that at z~3 more than 90% of the metals produced during the star forming history can be placed in a hot phase of the IGM, without violating any observational constraint. The observed galaxy mass-metallicity relation, and the intergalactic medium and intracluster medium metallicity evolution are also naturally explained by this hypothesis.

We study the clustering properties of metals in the intergalactic medium (IGM) as traced by 619 CIV, 81 SiIV, N >= 10^12 cm^-2 and 316 MgII, and 82 FeII N >= 10^11.5 cm^-2 absorption components in 19 high signal-to-noise (60-100 per pixel), high resolution (R = 45000) quasar spectra. Over the redshift range probed (1.5-3.0), CIV and SiIV trace each other closely and their line-of-sight correlation functions exhibit a steep decline at large separations and a flatter profile below ~ 150 km s^-1, with a large overall bias. These features do not depend on column depth. Carrying out a detailed SPH simulation (2 X 320^3, 57 Mpc^3 comoving), we show that this behavior can not be reproduced by models in which the IGM metallicity is constant or a local function of density. However, the CIV correlation function is consistent with a model in which metals are confined within bubbles with a typical radius Rs = 2 comoving Mpc about sources of mass >= Ms = 10^12 solar masses at z=3. Our lower redshift (0.5-2) measurements of the MgII and FeII correlation functions also uncover a steep decline at large separations and a flatter profile at small separations, but the clustering is even higher, and the turn-over is shifted to ~ 75 km s^-1. Again these features do not change with column depth. We describe an analytical bubble model for these species, which come from regions that are too compact to be simulated numerically, deriving best-fit values of R_s ~ 2.4 Mpc and M_s ~ 10^12 solar masses. Equally good fits to all four species are found in a similarly biased high-redshift enrichment model in which metals are placed within 2.4 comoving Mpc of 3 x 10^9 solar mass sources at z = 7.5.

We study the detectability of primordial metal-free stars. Cosmological enrichment is a local process that takes place over an extended redshift range. While the duration of this transition depends on several unknown factors, in all cases late-forming metal-free stars are found in ~ 10^8 solar mass objects, which are just large enough to cool, but small enough to not be clustered near areas of previous star formation. We discuss the observational properties of these objects, some of which may have already been detected in ongoing Lyman-alpha surveys. If metal-free stars have masses between 140 and 260 solar masses, they are expected to end their lives as pair-production supernovae (PPSNe). Using the KEPLER code, we calculate a set of light curves for these objects. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia supernovae, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches are able to place stringent limits on the fraction of very massive stars out to redshift 2. The planned Joint Dark Energy Mission satellite will be able to extend these constraints out to redshift 6

We explore the idea that the Type Ia supernovae (SNe Ia) rate consists of two components: a prompt piece that is proportional to the star formation rate (SFR) and an extended piece that is proportional to the total stellar mass. We fit the parameters of this model to the local observations of Mannucci and collaborators and then study its impact on three important problems. On cosmic scales, the model reproduces the observed SNe Ia rate density below z=1, and predicts that it will track the measured SFR density at higher redshift, reaching a value of 1-3.5 X 10^-4 per yr per Mpc^3 at z=2. In galaxy clusters, a large prompt contribution helps explain the iron content of the intracluster medium. Within the Galaxy, the model reproduces the observed stellar [O/Fe] abundance ratios if we allow a short (approximately 0.7 Gyr) delay in the prompt component. Ongoing medium-redshift SN surveys will yield more accurate parameters for our model

Nonrotating, zero metallicity stars with initial masses 140 < M < 260 solar masses are expected to end their lives as pair-production supernovae (PPSNe), in which an electron-positron pair-production instability triggers explosive nuclear burning. Interest in such stars has been rekindled by recent theoretical studies that suggest primordial molecular clouds preferentially form stars with these masses. Since metal enrichment is a local process, the resulting PPSNe could occur over a broad range of redshifts, in pockets of metal-free gas. Using the implicit hydrodynamics code KEPLER, we have calculated a set of PPSN light curves that addresses the theoretical uncertainties and allows us to assess observational strategies for finding these objects at intermediate redshifts. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches may soon be able to limit the contribution of these very massive stars to < 1% of the total star formation rate density out to z=2 which already provides useful constraints for theoretical models. The planned Joint Dark Energy Mission satellite will be able to extend these limits out to z=6.

The propagation of cosmological ionization fronts (I-fronts) during reionization is strongly influenced by small-scale structure. Here we summarize our recent attemps to understand the effect of this small-scale structure. We present high resolution cosmological N-body simulations at high-z (z>6) which resolve a wide range of halo mass, from minihalos to clusters of large, rare halos. We also study how minihalos affect I-fronts, through simulations of minihalo photoevaporation. Furthermore, we modify the I-front propagation equations to account for evolving small-scale structure, and incorporate these results into a semi-analytical reionization model. When intergalactic medium clumping and minihalo clustering around sources are included, small-scale structure affects reionization by slowing it down and extending it in time. This helps to explain observations by the Wilkinson Microwave Anisotropy Probe, which imply an early and extended reionization epoch. We also study how source clustering affects the evolution and size of H II regions, finding, in agreement with simulations, that H II regions usually expand, rarely shrinking. Hence, "relic H II regions'' are an exception, rather than the rule. When the suppression of small-mass sources in already-ionized regions by Jeans-mass filtering is accounted for, H II regions are smaller, delaying overlap. We also present a new numerical method for radiative transfer which is fast, efficient, and easily coupled to hydrodynamics and N-body codes, along with sample tests and applications.

I study the relationship between the spatial distribution of intergalactic metals and the masses and ejection energies of the sources that produced them. Over a wide range of models, metal enrichment is dominated by the smallest efficient sources, as the enriched volume scales roughly as E^{3/5} ~ M^{3/5} while the number density of sources goes as 1/M. In all cases, the earliest sources have the biggest impact, because fixed comoving distances correspond to smaller physical distances at higher redshifts. This means that most of the enriched volume is found around rare peaks, and intergalactic metals are naturally highly clustered. Furthermore, this clustering is so strong as to lead to a large overlap between individual bubbles. Thus the typical radius of enriched z ~ 3 regions should be interpreted as a constraint on groupings of sources rather than the ejection radius of a typical source. Similarly, the clustering of enriched regions should be taken as a measurement of source bias rather than mass.

The propagation of cosmological ionization fronts during the reionization of the universe is strongly influenced by small-scale gas inhomogeneities due to structure formation. These inhomogeneities include both collapsed minihalos, which are generally self-shielding, and lower-density structures, which are not. The minihalos are dense and sufficiently optically-thick to trap intergalactic ionization fronts, blocking their path and robbing them of ionizing photons until the minihalo gas is expelled as an evaporative wind. The lower-density structures do not trap these fronts, but they can slow them down by increasing the overall recombination rate in the intergalactic medium. In this paper we study the effects of both types of inhomogeneities, including nonlinear clustering effects, and we find that both IGM clumping and collapsed minihalos have significant yet qualitatively different impacts on reionization. While the number density of minihalos on average increases strongly with time, the density of minihalos inside H II regions around ionizing sources is largely constant. Thus the impact of minihalos is essentially to decrease the number of ionizing photons available to the IGM at all epochs, which is equivalent to a reduction in the luminosity of each source. On the other hand, the effect of IGM clumping increases strongly with time, slowing down reionization and extending it. Thus while the impact of minihalos is largely degenerate with the unknown source efficiency, IGM clumping can help significantly in reconciling the recent observations of cosmic microwave background polarization with quasar absorption spectra at z~6, which together point to an early but extended reionization epoch.

We present the results of an extensive spectroscopic survey of RX J0152.7-1357, one of the most massive distant clusters of galaxies known. Multi-object spectroscopy, carried out with FORS1 and FORS2 on the ESO Very Large Telescope (VLT), has allowed us to measure more than 200 redshifts in the cluster field and to confirm 102 galaxies as cluster members. The mean redshift of the cluster is $z=0.837 \pm 0.001$ and we estimate the velocity dispersion of the overall cluster galaxy distribution to be $\sim 1600 \mathrm{km \ s^{-1}}$. The distribution of cluster members is clearly irregular, with two main clumps that follow the X-ray cluster emission mapped by Chandra. A third clump of galaxies to the east of the central structure and at the cluster redshift has also been identified. The two main clumps have velocity dispersions of $\sim919$ and $\sim737 \mathrm{km s^{-1}}$ respectively, and the peculiar velocity of the two clumps suggests that they will merge into a single more massive cluster. A segregation in the star formation activity of the member galaxies is observed. All star forming galaxies are located outside the high-density peaks, which are populated only by passive galaxies. A population of red galaxies (belonging to the cluster red sequence) with clear post-starburst spectral features and [OII] ($\lambda$3727) emission lines is observed in the outskirts of the cluster. Two AGNs, which were previously confused with the diffuse X-ray emission from the intracluster medium in ROSAT and BeppoSAX observations, are found to be cluster members.

We carry out a numerical study of the effects of supernova-driven shocks on galaxy formation at z=9. These "cosmic explosions" can have a significant impact on galaxies forming nearby. We study such interactions in two key cases. In the first case in which the forming galaxy has already virialized, the impinging shock has only a small effect (< 1% of the gas is removed) and star formation continues relatively unimpeded. However, in the second case in which the nearby forming galaxy is at the more diffuse turn-around stage, a large fraction (~70%) of the gas is stripped away from the host dark-matter halo and ejected into the intergalactic medium. As the time spent near turn-around is much longer than the interval from virialization to galaxy formation due to strong radiative losses, we expect the second case to be more representative of the majority of outflow-galaxy interactions. Thus SN-driven pregalactic outflows may be an efficient mechanism for inhibiting the formation of neighbouring galaxies at high redshift. We briefly outline the possible cosmological consequences of this effect.

We carry out a detailed numerical investigation of the spatial correlation function of the initial positions of cosmological dark matter halos. In this Lagrangian coordinate system, which is especially useful for analytic studies of cosmological feedback, we are able to construct cross-correlation functions of objects with varying masses and formation redshifts and compare them with a variety of analytical approaches. For the case in which both formation redshifts are equal, we find good agreement between our numerical results and the bivariate model of Scannapieco & Barkana (2002; SB02) at all masses, redshifts, and separations, while the model of Porciani et al. (1998) does well for all parameters except for objects with different masses at small separations. We find that the standard mapping between Lagrangian and Eulerian bias performs well for rare objects at all separations, but fails if the objects are highly-nonlinear (low-sigma) peaks. In the Lagrangian case in which the formation redshifts differ, the SB02 model does well for all separations and combinations of masses, apart from a discrepancy at small separations in situations in which the smaller object is formed earlier and the difference between redshifts or masses is large. As this same limitation arises in the standard approach to the single-point progenitor distribution developed by Lacey & Cole (1993), we conclude that a more complete understanding of the progenitor distribution is the most important outstanding issue in the analytic modeling of Lagrangian bias.

We investigate the interactions of high-redshift galaxy outflows with low-mass virialized (Tvir < 10,000K) clouds of primordial composition. While atomic cooling allows star formation in larger primordial objects, such "minihalos" are generally unable to form stars by themselves. However, the large population of high-redshift starburst galaxies may have induced widespread star formation in these objects, via shocks that caused intense cooling both through nonequilibrium H2 formation and metal-line emission. Using a simple analytic model, we show that the resulting star clusters naturally reproduce three key features of the observed population of halo globular clusters (GCs). First, the 10,000 K maximum virial temperature corresponds to the ~ 10^6 solar mass upper limit on the stellar mass of GCs. Secondly, the momentum imparted in such interactions is sufficient to strip the gas from its associated dark matter halo, explaining why GCs do not reside in dark matter potential wells. Finally, the mixing of ejected metals into the primordial gas is able to explain the ~ 0.1 dex homogeneity of stellar metallicities within a given GC, while at the same time allowing for a large spread in metallicity between different clusters. To study this possibility in detail, we use a simple 1D numerical model of turbulence transport to simulate mixing in cloud-outflow interactions. We find that as the shock shears across the side of the cloud, Kelvin-Helmholtz instabilities arise, which cause mixing of enriched material into > 20% of the cloud. Such estimates ignore the likely presence of large-scale vortices, however, which would further enhance turbulence generation. Thus quantitative mixing predictions must await more detailed numerical studies.

We consider the impact of quasar outflows on structure formation. Such outflows are potentially more important than galactic winds, which appear insufficient to produce the level of preheating inferred from X-ray observations of galaxy clusters. Using a simple analytical model for the distribution of quasars with redshift, coupled with a one-dimensional Sedov-Taylor model for outflows, we are able to make robust statements about their impact of on structure formation. As large regions of the IGM are heated above a critical entropy of approximately 100 keV cm^2, cooling become impossible within them, regardless of changes in density. On quasar scales, this has the effect of inhibiting further formation, resulting in the observed fall-off in their number densities below z = 2. On galaxy scales, quasar feedback fixes the turn-over scale in the galaxy luminosity function (L_*) as the nonlinear scale at the redshift of strong feedback. The galaxy luminosity function then remains largely fixed after this epoch, consistent with recent observations and in contrast to the strong evolution predicted in more standard galaxy-formation models. Finally, strong quasar feedback explains why the intracluster medium is observed to have been pre-heated to entropy levels just above S_crit, the minimum excess that would not have been erased by cooling. The presence of such outflows is completely consistent with the observed properties of the Lyman-alpha forest at z ~ 2, but is expected to have a substantial and detectable impact on Compton distortions observed in the microwave background and the multiphase properties of the "warm-hot" (z=0) circumgalactic medium.

It is a well-established fact that massive cosmological objects exhibit a ``geometrical bias'' that boosts their spatial correlations with respect to the underlying mass distribution. Although this geometrical bias is a simple function of mass, this is only half of the story. We show using numerical simulations that objects that are in the midst of accreting material also exhibit a ``temporal bias,'' which further boosts their clustering far above geometrical bias levels. These results may help to resolve a discrepancy between spectroscopic and clustering mass estimates of Lyman Break Galaxies, a population of high-redshift galaxies that are caught in the act of forming large numbers of new stars.

We measure the spatial clustering of metals in the intergalactic medium from z = 1.7 to 3.0, as traced by 643 CIV and 104 SiIV N >= 10^12 cm^-2 absorption systems in 19 high signal-to-noise (40-80) and high resolution (R = 45000) quasar spectra. The number densities and two-point correlation functions of both these species are largely constant with redshift, suggesting the bulk of metal ejection occurred at z >= 3. However, at z <= 1.9 some additional signature appears in the CIV correlation function at 500 km/s, associated with four strong and peculiar systems. At all redshifts, the CIV and SiIV correlation functions exhibit a steep rise at large separations and a flatter profile at small separations, with an elbow occurring at ~150 km/s. We show that these properties are consistent with metals confined within bubbles with a typical radius R_s about sources of mass >= M_s, and use numerical simulations to derive best-fit values of R_s ~ 2 comoving Mpc and M_s ~ 5x10^11 solar masses at z = 3. This does not exclude that metals could have been produced at higher redshifts in smaller, but equally rare, objects. At the level of detection of this survey, IGM enrichment is likely to be incomplete and inhomogeneous, with a filling factor ~ 10%.

Motivated by the observed discrepancy between the strong spatial correlations of Lyman break galaxies (LBGs) and their velocity dispersions, we consider a theoretical model in which these starbursting galaxies are associated with dark matter halos that experience appreciable infall of material. We show using numerical simulation that selecting halos that substantially increase in mass within a fixed time interval introduces a ``temporal bias'' which boosts their clustering above that of the underlying population. If time intervals consistent with the observed LBGs star formation rates of 50 solar masses per year are chosen, then spatial correlations are enhanced by up to a factor of two. These values roughly correspond to the geometrical bias of objects three times as massive. Thus, it is clear that temporal biasing must be taken into account when interpreting the properties of Lyman break galaxies.

In the standard Cold Dark Matter (CDM) theory of structure formation, virialized minihalos (with T_vir = 10,000 K or less) form in abundance at high redshift (z > 6), during the cosmic "dark ages." The hydrogen in these minihalos, the first nonlinear baryonic structures to form in the universe, is mostly neutral and sufficiently hot and dense to emit strongly at the 21-cm line. We calculate the emission from individual minihalos and the radiation background contributed by their combined effect. Minihalos create a "21-cm forest" of emission lines. We predict that the angular fluctuations in this 21-cm background should be detectable with the planned LOFAR and SKA radio arrays, thus providing a direct probe of structure formation during the "dark ages." Such a detection will serve to confirm the basic CDM paradigm while constraining the background cosmology parameters, the shape of the power-spectrum of primordial density fluctuations, the onset and duration of the reionization epoch, and the conditions which led to the first stars and quasars. We present results here for the currently-favored, flat LCDM model, for different tilts of the primordial power spectrum. These minihalos will also cause a "21-cm forest" of absorption lines, as well, in the spectrum of radio continuum sources at high redshift, if the latter came into existence before the end of reionization.

We conduct a comprehensive investigation of the detectability of the first stars and their enrichment signatures in galaxy clusters. We show that the mean metallicity of outflows from objects containing these Population III (PopIII) stars is well above the critical transition metallicity (Z_cr \sim 10^-4) that marks the formation of normal stars. Thus the fraction of PopIII objects formed as a function of redshift is heavily dependent on the distribution of metals and fairly independent of the precise value of Z_cr. Using an analytical model of inhomogenous structure formation, we study the evolution of PopIII objects as a function of the star formation efficiency, IMF, and efficiency of outflow generation. For all models, PopIII objects tend to be in the 10^6.5-10^7.0 solar mass range, just large enough to cool within a Hubble time, but small enough that they are not clustered near areas of previous star formation. Although the mean metallicity exceeds Z_cr at a redshift of 15 in all models, the peak of PopIII star formation occurs at z \sim 10, and such stars continue to form well into the observable range. We discuss the observational properties of these objects, some of which may have already been detected in ongoing surveys of high-redshift Lyman-alpha emitters. Finally, we combine our PopIII distributions with the yield models of Heger and Woosley (2002) to study their impact on the intracluster medium (ICM) in galaxy clusters. We find that PopIII stars can contribute no more than 20% of the iron observed in the ICM, although their peculiar elemental yields help to reconcile theoretical models with the observed Fe and Si/Fe abundances. However, these stars tend to overproduce S/Fe and their associated SN heating falls far short of the observed level of 1 keV per ICM gas particle.

Hydrogen atoms inside virialized minihalos generate a radiation background from redshifted 21-cm line emission whose angular fluctuations reflect clustering during before and during reionization. We have shown elsewhere that this emission may be detectable with the planned Low Frequency Array (LOFAR) and Square Kilometer Array (SKA) in a flat LambdaCDM. This is a direct probe of structure during the ``Dark Ages'' at redshifts z>6 and down to smaller scales than have previously been constrained. In our original calculation, we used a standard approximation known as the ``linear bias'' [e.g. Mo & White (1996)]. Here we improve upon that treatment by considering the effect of nonlinear clustering. To accomplish this, we develop a new analytical method for calculating the nonlinear Eulerian bias of halos, which should be useful for other applications as well. Predictions of this method are compared with the results of LambdaCDM N-body simulations, showing significantly better agreement than the standard linear bias approximation. When applied to the 21-cm background from minihalos, our formalism predicts fluctuations that differ from our original predictions by up to 30% at low frequencies and small scales. However, within the range of frequencies and angular scales at which the signal could be observable by LOFAR and SKA as currently planned, the differences are small and our original predictions prove robust. Our results indicate that while a smaller frequency bandwidth of observation leads to a higher signal that is more sensitive to nonlinear effects, this effect is counteracted by the lowered sensitivity of the radio arrays. We calculate the best frequency bandwidth for these observations to be Delta nu~2 MHz.(abbreviated)

A wide range of observational and theoretical arguments suggest that the universe experienced a period of heating and metal enrichment, most likely from starbursting dwarf galaxies. Using a hydrodynamic simulation we have conducted a uniquely detailed theoretical investigation of this epoch at the end of the cosmological ``dark ages''. Outflows strip baryons from pre-viralized halos with total masses $\lesssim10{}^{10}$ M${}_\odot$, reducing their number density and the overall star formation rate, while pushing these quantities toward their observed values. We show that the metallicity of $\lesssim10{}^{10}$ M${}_\odot$ objects increases with size, but with a large scatter, reproducing the metallicity-luminosity relation of dwarf galaxies. Galaxies $\gtrsim10{}^{10}$ M${}_\odot$ form with a roughly constant initial metallicity of 10% solar, explaining the observed lack of metal-poor disk stars in these objects. Outflows enrich roughly 20% of the simulation volume, yielding a mean metallicity of 0.3% solar, in agreement with observations of CIV in QSO absorption-line systems.

We present a phenomenological model of feedback in early-type galaxies that tracks the evolution of the interstellar medium gas mass, metallicity, and temperature. Modeling the star formation rate as a Schmidt law with a temperature-dependent efficiency, we find that intermittent episodes of star formation are common in moderate-size ellipticals. Our model is applicable in the case in which the thermalization time from SN is sufficiently long that spatial variations are relatively unimportant, an appropriate assumption for the empirical parameters adopted here. The departure from a standard scenario of passive evolution implies significantly younger luminosity-weighted ages for the stellar populations of low-mass galaxies at moderate redshifts, even though the more physically meaningful mass-weighted ages are changed only slightly. Secondary bursts of star formation also lead to a natural explanation of the large scatter in the NUV-optical relation observed in clusters at moderate redshift and account for the population of E+A galaxies that display a spheroidal morphology. As the late-time formation of stars in our model is due to the gradual cooling of the interstellar medium, which is heated to temperatures ~1 keV by the initial burst of supernovae, our conclusions do not rely on any environmental effects or external mechanisms. Furthermore, a simple estimate of the X-ray emission from this supernova heated gas leads to an L_X vs L_B correlation that is in good agreement with observed values. Thus feedback processes may be essential to understanding the observed properties of early-type galaxies from the optical to the X-ray.

We develop an analytical formalism that is suitable for studying inhomogeneous structure formation, by studying the joint statistics of dark matter halos forming at two points. Extending the Bond et al. (1991) derivation of the mass function of virialized halos, based on excursion sets, we derive an approximate analytical expression for the ``bivariate'' mass function of halos forming at two redshifts and separated by a fixed comoving Lagrangian distance. Our approach also leads to a self-consistent expression for the nonlinear biasing and correlation function of halos, generalizing a number of previous results including those by Kaiser (1984) and Mo & White (1996). We compare our approximate solutions to exact numerical results within the excursion-set framework and find them to be consistent to within 2% over a wide range of parameters. Our formalism can be used to study various feedback effects during galaxy formation analytically, as well as to simply construct observable quantities dependent on the spatial distribution of objects. A code that implements our method is publicly available at http://www.arcetri.astro.it/~evan/Gemini

Supernova-driven outflows from early galaxies may have had a large impact on the kinetic and chemical structure of the intergalactic medium (IGM). We use three-dimensional Monte Carlo cosmological realizations of a simple linear peaks model to track the time evolution of such metal-enriched outflows and their feedback on galaxy formation. We find that at most 30% of the IGM by volume is enriched to values above 10^-3 solar in models that only include objects that cool by atomic transitions. The majority of enrichment occurs relatively early (5 < z < 12) and resulting in a mass-averaged cosmological metallicity between 10^-3 and 10^-1.5 solar. The inclusion of Population III objects that cool through H2 line emission has only a minor impact on these results: increasing the mean metallicity and filling factor by at most a factor of 1.4, and moving the dawn of the enrichment epoch to a redshift of approximately 14 at the earliest. Thus enrichment by outflowing galaxies is likely to have been incomplete and inhomogeneous, biased to the areas near the starbursting galaxies themselves. Models with a 10% star formation efficiency can satisfactorily reproduce the nearly constant (2 < z < 5, Z approximately 3.5 x 10^-4 solar) metallicity of the low column density Ly-alpha forest derived by Songaila (2001), an effect of the decreasing efficiency of metal loss from larger galaxies. Finally, we show that IGM enrichment is intimately tied to the ram-pressure stripping of baryons from neighboring perturbations. This results in the suppression of at least 20% of the dwarf galaxies in the mass range 10^8.5 to 10^9.5 solar, in all models with filling factors greater than 2%, and an overall suppression of approximately 50% of dwarf galaxies in the most observationally-favored model.

We study the effect of a prolonged epoch of reionization on the angular power spectrum of the Cosmic Microwave Background. Typically reionization studies assume a sudden phase transition, with the intergalactic gas moving from a fully neutral to a fully ionized state at a fixed redshift. Such models are at odds, however, with detailed investigations of reionization, which favor a more extended transition. We have modified the code CMBFAST to allow the treatment of more realistic reionization histories and applied it to data obtained from numerical simulations of reionization. We show that the prompt reionization assumed by CMBFAST in its original form heavily contaminates any constraint derived on the reionization redshift. We find, however, that prompt reionization models give a reasonable estimate of the epoch at which the mean cosmic ionization fraction was ~50%, and provide a very good measure of the overall Thomson optical depth. The overall differences in the temperature (polarization) angular power spectra between prompt and extended models with equal optical depths are less than 1% (10%).

We examine the effects of galaxy outflows on the formation of dwarf galaxies in numerical simulations of the high-redshift Universe. Using a Smoothed Particle Hydrodynamic code, we conduct two detailed simulations of a (5.2 Mpc/h)^3 comoving volume of the Universe. In both simulations we implement simple, well-motivated models of galaxy identification and star formation, while our second simulation also includes a simple ``blow-out'' model of galaxy outflows in which supernova driven winds from newly formed disk galaxies punch-out and shock the intergalactic medium while leaving the host galaxies intact. A direct comparison between these simulations suggests that there are two major mechanisms by which outflows affect dwarf formation. Firstly, the formation of an outflow slows down the further accretion of gas onto a galaxy, causing an overall decrease of approximately 50% in the total gas mass accreted by the objects in our simulations. Additionally, our simulations uncover a significant population of approximately 10^9 solar mass objects whose formation is suppressed by the mechanism of ``baryonic stripping,'' in which outflows from early galaxies strip the gas out of nearby overdense regions that would have otherwise later formed into dwarf galaxies. This mechanism may be important in explaining the observed discrepancy between the number of dwarf galaxies predicted and observed in the local group and provide a natural explanation for the formation of empty halos which may be required by the existence of the extremely gas-poor extra-galactic High-Velocity Clouds.

We compare the metallicity distribution of Galactic Halo stars with 3D realizations of hierarchical galaxy formation. Outflows from dwarf galaxies enrich the intergalactic medium inhomogeneously, at a rate depending on the local galaxy density. Consequently, the first stars created in small early-forming galaxies are less metal-rich that the first stars formed in more massive galaxies which typically form later. As most halo stars are likely to originate in accreted dwarfs, while disk stars formed out of outflow-enriched gas, this scenario naturally generates a ``metallicity floor'' for old disk stars, which we find to be roughly coincident with the higher end of our predicted metallicity distribution of halo stars, in agreement with observations. The broad and centrally peaked distribution of halo star metallicities is well reproduced in our models, with a natural dispersion depending on the exact accretion history. Our modeling includes the important ``baryonic stripping'' effect of early outflows, which brush away the tenuously held gas in neighboring pre-virialized density perturbations. This stripping process does not significantly modify the predicted shape of the halo star metal distribution but inhibits star-formation and hence the number of accreted stars, helping to reproduce the observed total Galactic halo luminosity and also the lack of low-luminosity local dwarf galaxies relative to N-body predictions.

We show that the winds identified with high-redshift low-mass galaxies may strongly affect the formation of stars in more massive galaxies that form later. With 3D realizations of a simple linear growth model we track gas shocking, metal enrichment, and cooling, together with dark halo formation. We show that outflows typically strip baryonic material out of collapsing intermediate mass halos, suppressing star formation. More massive halos can trap the heated gas but collapse late, leading to a broad bimodal redshift distribution, with a larger characteristic mass associated with the lower redshift peak. This scenario accounts for the observed bell-shaped luminosity function of early type galaxies, explains the small number of Milky Way satellite galaxies relative to Cold Dark Matter models predictions, and provides a possible explanation for the lack of metal poor G-dwarfs in the solar neighborhood and the more general lack of low-metallicity stars in massive galaxies relative to ``closed-box'' models of chemical enrichment. Intergalactic medium heating from outflows should produce spectral distortions in the cosmic microwave background that will be measurable with the next generation of experiments.

We show that the gas in growing density perturbations is vulnerable to the influence of winds outflowing from nearby collapsed galaxies that have already formed stars. This suggests that the formation of nearby galaxies with masses less than 10^9 solar masses is likely to be suppressed, irrespective of the details of galaxy formation. An impinging wind may shock heat the gas of a nearby perturbation to above the virial temperature, thereby mechanically evaporating the gas, or the baryons may be stripped from the perturbation entirely if they are accelerated to above the escape velocity. We show that baryonic stripping is the most effective of these two processes, because shock-heated clouds that are too large to be stripped are able to radiatively cool within a sound-crossing time, limiting evaporation. The intergalactic medium temperatures and star-formation rates required for outflows to have a significant influence on the formation of low-mass galaxies are consistent with current observations, but may soon be examined directly via associated distortions in the cosmic microwave background, and with near-infrared observations from the Next Generation Space Telescope, which may detect the supernovae from early-forming stars.

We examine the gravitational redshift of radiation emitted from within the potential of a cluster. Spectral lines from the intracluster medium (ICM) are redshifted in proportion to the emission-weighted mean potential along the line of sight, amounting to approximately 50 km/s at a radius of 100 kpc/h, for a cluster dispersion of 1200 km/s. We show that the relative redshifts of different ionization states of metals in the ICM provide a unique probe of the three-dimensional matter distribution. An examination of the reported peculiar velocities of cD galaxies in well studied Abell clusters reveals they are typically redshifted by an average of $\sim +200$ km/s. This can be achieved by gravity with the addition of a steep central potential associated with the cD galaxy. Note that in general gravitational redshifts cause a small overestimate of the recessional velocities of clusters by an average of $\sim$ 20 km/s.

The dominant linear contribution to cosmic microwave background (CMB) fluctuations at small angular scales (less than one arcsec) is a second-order contribution known as the Vishniac or Ostriker-Vishniac effect. This effect is caused by the scattering of CMB photons off free electrons after the universe has been reionized, and is dominated by linear perturbations near the R_V =2 Mpc/(h Gamma/0.2) scale in the Cold Dark Matter cosmogony. As the reionization of the universe requires that nonlinear objects exist on some scale, however, one can compare the scale responsible for reionization to R_V and ask if a linear treatment is even feasible in different scenarios of reionization. For an Omega_0 = 1 cosmology normalized to cluster abundances, only about 65% of the linear integral is valid if reionization is due to quasars in halos of mass 10^9 solar, while 75% of the integral is valid if reionization was caused by stars in 10^6 solar mass halos. In lambda or open cosmologies, both the redshift of reionization and z_V are pushed further back, but still only 75% to 85% of the linear integral is valid, independent of the ionization scenario. We point out that all odd higher-order moments from Vishniac fluctuations are zero while even moments are non-zero, regardless of the gaussianity of the density perturbations. This provides a defining characteristic of the Vishniac effect that differentiates it from other secondary perturbations and may be helpful in separating them.

We describe a measurement of the angular power spectrum of anisotropies in the Cosmic Microwave Background (CMB) from 0.3 degrees to ~10 degrees from the North American test flight of the BOOMERANG experiment. BOOMERANG is a balloon-borne telescope with a bolometric receiver designed to map CMB anisotropies on a Long Duration Balloon flight. During a 6-hour test flight of a prototype system in 1997, we mapped > 200 square degrees at high galactic latitudes in two bands centered at 90 and 150 GHz with a resolution of 26 and 16.6 arcmin FWHM respectively. Analysis of the maps gives a power spectrum with a peak at angular scales of ~1 degree with an amplitude ~70 uK.

We consider the effect of lensing by galaxy clusters on the angular distribution of submillimeter wavelength objects. While lensing does not change the total flux and number counts of submillimeter sources, it can affect the number counts and fluxes of flux-limited samples. Therefore imposing a flux cut on point sources not only reduces the overall Poisson noise, but imprints the correlations between lensing clusters on the unresolved flux distribution. Using a simple model, we quantify the lensing anisotropy induced in flux-limited samples and compare this to Poisson noise. We find that while the level of induced anisotropies on the scale of the cluster angular correlation length is comparable to Poisson noise for a slowly evolving cluster model, it is negligible for more realistic models of cluster evolution. Thus the removal of point sources is not expected to induce measurable structure in the microwave or far-infrared backgrounds.

The Wavelength-Oriented Microwave Background Analysis Team (WOMBAT) is constructing microwave maps which will be more realistic than previous simulations. Our foreground models represent a considerable improvement: where spatial templates are available for a given foreground, we predict the flux and spectral index of that component at each place on the sky and estimate uncertainties. We will produce maps containing simulated CMB anisotropy combined with expected foregrounds. The simulated maps will be provided to the community as the WOMBAT Challenge, so such maps can be analyzed to extract cosmological parameters by scientists who are unaware of their input values. This will test the efficacy of foreground subtraction, power spectrum analysis, and parameter estimation techniques and help identify the areas most in need of progress. These maps are also part of the FORECAST project, which allows web-based access to the known foreground maps for the planning of CMB missions.

We study the effect of a finite topology on the temperature correlations of the cosmic microwave background in a flat universe. Analytic expressions for the angular power spectrum are given for all possible finite flat models. We examine the angular correlation function itself, pointing out visible and discrete features that arise from topology. While observations of the power spectrum on large angular scales can be used to place bounds on the minimum topology length, cosmic variance generally restricts us from differentiating one flat topology from another. Schemes that acknowledge anisotropic structures, such as searches for ghosts, circles or geometric patterns, will be needed to further probe topology.

The Wavelength-Oriented Microwave Background Analysis Team (WOMBAT) is constructing microwave skymaps which will be more realistic than previous simulations. Our foreground models represent a considerable improvement: where spatial templates are available for a given foreground, we predict the flux and spectral index of that component at each place on the sky and estimate the uncertainties in these quantities. We will produce maps containing simulated Cosmic Microwave Background anisotropies combined with all major expected foreground components. The simulated maps will be provided to the cosmology community as the WOMBAT Challenge, a "hounds and hares" exercise where such maps can be analyzed to extract cosmological parameters by scientists who are unaware of their input values. This exercise will test the efficacy of current foreground subtraction, power spectrum analysis, and parameter estimation techniques and will help identify the areas most in need of progress.

The universe displays a three-dimensional pattern of hot and cold spots in the radiation remnant from the big bang. The global geometry of the universe can be revealed in the spatial distribution of these spots. In a topologically compact universe, distinctive patterns are especially prominent in spatial correlations of the radiation temperature. Whereas these patterns are usually washed out in statistical averages, we propose a scheme which uses the universe's spots to observe global geometry in a manner analogous to the use of multiple images of a gravitationally lensed quasar to study the geometry of the lens. To demonstrate how the geometry of space forms patterns in observations of the microwave sky, we develop a simple real-space approximation to estimate temperature correlations for any set of cosmological parameters and any global geometry. We present correlated spheres which clearly show geometric pattern formation for compact flat universes as well as for the compact negatively curved space introduced by Weeks and another discovered by Best. These examples illustrate how future satellite-based observations of the microwave background can determine the full geometry of the universe.

We analyze the effects of the detector response time on bolometric measurements of the anisotropy of the cosmic microwave background (CMB). We quantify the effect in terms of a single dimensionless parameter $L$ defined as the ratio between the time the beam sweeps its own size and the bolometer response time. As $L$ decreases below $\sim 2.5$ the point source response of the experiment becomes elongated. We introduce a window function matrix based on the timestream data to assess the effects of the elongated beam. We find that the values of the window function matrix elements decrease slowly as a function $L$. Our analysis and results apply to other cases of beam asymmetry. For the High Frequency Instrument on board the Planck Surveyor satellite we show that for a broad range of $L$ the ability of the experiment to extract the cosmological parameters is not degraded. Our analysis enhances the flexibility in tuning the design parameters of CMB anisotropy experiments.

Clues as to the geometry of the universe are encoded in the cosmic background radiation. Hot and cold spots in the primordial radiation may be randomly distributed in an infinite universe while in a universe with compact topology distinctive patterns can be generated. With improved vision, we could actually see if the universe is wrapped into a hexagonal prism or a hyperbolic horn. We discuss the search for such geometric patterns in predictive maps of the microwave sky.

The global geometry of the universe is in principle as observable an attribute as local curvature. Previous studies have established that if the universe is wrapped into a flat hypertorus, the simplest compact space, then the fundamental domain must be at least 0.4 times the diameter of the observable universe. Despite a standard lore that the other five compact, orientable flat spaces are more weakly constrained, we find the same bound holds for all. Our analysis provides the first limits on compact cosmologies built from the identifications of hexagonal prisms.

We propose a new method for constraining a primordial homogeneous magnetic field with the cosmic microwave background. Such a field will induce an observable parity odd cross correlation between the polarization anisotropy and the temperature anisotropy by Faraday rotation. We analyze the necessary experimental features to match, and improve, current constraints of such a field by measuring this correlation.