We examine the thermodynamic state and cooling of the low-$z$ Circum-Galactic Medium (CGM) in five FIRE-2 galaxy formation simulations of Milky Way-mass galaxies. We find that the CGM in these simulations is generally multiphase and dynamic, with a wide spectrum of largely nonlinear density perturbations sourced by the accretion of gas from the IGM and winds from both the central and satellite galaxies. We investigate the origin of the multiphase structure of the CGM with a particle tracking analysis and find that most of the low entropy gas has cooled from the hot halo as a result of thermal instability triggered by these perturbations. The ratio of cooling to free-fall timescales $t_{\rm cool}/t_{\rm ff}$ in the hot component of the CGM spans a wide range $\sim 1-100$ at a given radius, but exhibits approximately constant median values $\sim 5-20$ at all radii $0.1 R_{\rm vir} < r < R_{\rm vir}$. These are similar to the $\approx 10-20$ value typically adopted as the thermal instability threshold in "precipitation" models of the ICM. Consequently, a one-dimensional model based on the assumption of a constant $t_{\rm cool}/t_{\rm ff}$ and hydrostatic equilibrium approximately reproduces the simulation number density and entropy profiles, but only if it assumes the metallicity profiles taken directly from the simulations. We explicitly show that the $t_{\rm cool}/t_{\rm ff}$ value of a gas parcel in the hot component of the CGM does not predict its probability of cooling and subsequently accreting onto the central galaxy. This suggests that the value of $t_{\rm cool}/t_{\rm ff}$ is a poor predictor of thermal stability in gaseous halos in which large-amplitude density perturbations are prevalent.

Galaxy evolution is thought to be driven in large part by the flow of gas between galaxies and the circumgalactic medium (CGM), a halo of metal-enriched gas extending out to $\gtrsim100$ kpc from each galaxy. Studying the spatial structure of the CGM holds promise for understanding these gas flow mechanisms; however, the common method using background quasar sightlines provides minimal spatial information. Recent works have shown the utility of extended background sources such as giant gravitationally lensed arcs. Using background lensed arcs from the CSWA 38 lens system, we continuously probed, at a resolution element of about 15 kpc$^2$, the spatial and kinematic distribution of MgII absorption in a star-forming galaxy at $z=0.77$ (stellar mass $\approx 10^{9.7}$ M$_\odot$, star formation rate $\approx 10$ M$_\odot$ yr$^{-1}$) at impact parameters $D=5-40$ kpc. Our results present an anisotropic, optically thick medium whose absorption strength decreases with increasing impact parameter, in agreement with the statistics towards quasars and other gravitational arcs. Furthermore, we detect low line-of-sight velocities ($v\approx-25-80$ km s$^{-1}$) and relatively high velocity dispersion ($\sigma\approx50\pm30$ km s$^{-1}$) in the MgII gas. These measures provide evidence of a mainly pressure-supported, metal-enriched diffuse gas recycling through the CGM rather than an active galactic outflow.

We present optical light curves from the Transiting Exoplanet Survey Satellite (TESS) for the archetypical dwarf AGN in the nearby galaxy NGC 4395 hosting a $\sim 10^5\,M_\odot$ supermassive black hole. Significant variability is detected on timescales from weeks to hours before reaching the background noise level. The $\sim$month-long, 30-minute-cadence, high-precision TESS light curve can be well fit by a simple damped random walk model, with the damping timescale $\tau_{\rm DRW}$ constrained to be $2.3_{-0.7}^{+1.8}$~days ($1\sigma$). \obj\ lies almost exactly on the extrapolation of the $\tau_{\rm DRW}-M_{\rm BH}$ relation measured for AGN with BH masses that are more than three orders of magnitude larger. The optical variability periodogram can be well fit by a broken power-law with the high-frequency slope ($-1.88\pm0.15$) and the characteristic timescale ($\tau_{\rm br}\equiv 1/(2\pi f_{\rm br})=1.4_{-0.5}^{+1.9}\,$days) consistent with the DRW model within 1$\sigma$. This work demonstrates the power of TESS light curves in identifying low-mass accreting SMBHs with optical variability, and a potential global $\tau_{\rm DRW}-M_{\rm BH}$ relation that can be used to estimate SMBH masses with optical variability measurements.

Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $\gamma$-ray constraints. Here, we compare the effects of these models relative to simpler 'diffusion+streaming' models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie 'in between' models which ignore CRs entirely and models which treat CRs with constant diffusivity.

While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and distribution of metals in the CGM of a single galaxy using ions that probe a wide range of gas phases (Si II, Si III, Si IV, C II, C IV, and O VI, the latter being from the Far Ultraviolet Spectroscopic Explorer). We find that Si III and O VI have near unity covering factor maintained all the way out to 1.2Rvir and 1.9Rvir, respectively. We show that Si III is the dominant ion over Si II and Si IV at any R. While we do not find that the properties of the CGM of M31 depend strongly on the azimuth, we show that they change remarkably around 0.3-0.5Rvir, conveying that the inner regions of the CGM of M31 are more dynamic and have more complicated multi-phase gas-structures than at R>0.5Rvir. We estimate the metal mass of the CGM within Rvir as probed by Si II, Si III, and Si IV is 2x10^7 Msun and by O VI is >8x10^7 Msun, while the baryon mass of the 10^4-10^5.5 K gas is ~4x10^10 (Z/0.3 Zsun)^(-1) Msun within Rvir. We show that different zoom-in cosmological simulations of L* galaxies better reproduce the column density profile of O VI with R than Si III or the other studied ions. We find that observations of the M31 CGM and zoom-in simulations of L* galaxies have both lower ions showing higher column density dispersion and dependence on R than higher ions, indicating that the higher ionization structures are larger and/or more broadly distributed.

The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{\ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $\gamma$-rays ($L_{\gamma}$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_{\gamma}$), as transport is rate-limited by the ionized WIM and 'inner CGM' gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_{\gamma}$ can be dominated by trapping in small 'patches.' Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $\kappa$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $\kappa \gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $\kappa$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $\gamma$-ray emission. But around $\kappa\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.

We present an updated model of the cosmic ionizing background from the UV to the X-rays. Relative to our previous model (Faucher-Giguere et al. 2009), the new model provides a better match to a large number of up-to-date empirical constraints, including: 1) new galaxy and AGN luminosity functions; 2) stellar spectra including binary stars; 3) obscured and unobscured AGN; 4) a measurement of the non-ionizing UV background; 5) measurements of the intergalactic HI and HeII photoionization rates at z~0-6; 6) the local X-ray background; and 7) improved measurements of the intergalactic opacity. In this model, AGN dominate the HI ionizing background at z<~3 and star-forming galaxies dominate it at higher redshifts. Combined with the steeply declining AGN luminosity function beyond z~2, the slow evolution of the HI ionization rate inferred from the high-redshift HI Lya forest requires an escape fraction from star-forming galaxies that increases with redshift (a population-averaged escape fraction of ~1% suffices to ionize the intergalactic medium at z=3 when including the contribution from AGN). We provide effective photoionization and photoheating rates calibrated to match the Planck 2018 reionization optical depth and recent constraints from the HeII Lya forest in hydrodynamic simulations.

We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. Commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds.

Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: when did the main progenitor of a MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We investigate the formation times and progenitor galaxies of MW/M31-mass galaxies using the FIRE-2 cosmological simulations, including 6 isolated MW/M31-mass galaxies and 6 galaxies in Local Group (LG)-like pairs at z = 0. We examine main progenitor "formation" based on two metrics: (1) transition from primarily ex-situ to in-situ stellar mass growth and (2) mass dominance compared to other progenitors. We find that the main progenitor of a MW/M31-mass galaxy emerged typically at z ~ 3-4 (11.6-12.2 Gyr ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z < 5 (< 12.6 Gyr ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (\Delta z ~ 2, \Delta t ~ 1.6 Gyr), with ~ 4x higher stellar mass at all z > 4 (> 12.2 Gyr ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. Overall, about 100 galaxies with M_star > 10^5 M_sun formed a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ~ 5x) of the progenitor population.

We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as its location at z=0 (its "ultimate" fate). Of the CGM at z=2, about half is found in the ISM or stars of the central galaxy by z=0 in Mh(z=2) ~ 5e11 Msun halos, but most of the CGM in lower-mass halos ends up in the IGM. This is so even though most of the CGM in M_h(z=2) ~ 5e10 Msun halos first accretes onto the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z=0.25 remains in the CGM by z=0 at all halo masses analyzed. Of the CGM gas that subsequently accretes onto the central galaxy in the progenitors of Mh(z=0) ~10^12 Msun halos, most of it is cool (T~10^4 K) at z=2 but hot (~Tvir) at z=0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z=0.25 and z=2 >~80% of the cool gas in Mh >~ 10^11 Msun halos will accrete onto a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM.

We present the implementation and the first results of cosmic ray (CR) feedback in the Feedback In Realistic Environments (FIRE) simulations. We investigate CR feedback in non-cosmological simulations of dwarf, sub-$L\star$ starburst, and $L\star$ galaxies with different propagation models, including advection, isotropic and anisotropic diffusion, and streaming along field lines with different transport coefficients. We simulate CR diffusion and streaming simultaneously in galaxies with high resolution, using a two moment method. We forward-model and compare to observations of $\gamma$-ray emission from nearby and starburst galaxies. We reproduce the $\gamma$-ray observations of dwarf and $L\star$ galaxies with constant isotropic diffusion coefficient $\kappa \sim 3\times 10^{29}\,{\rm cm^{2}\,s^{-1}}$. Advection-only and streaming-only models produce order-of-magnitude too large $\gamma$-ray luminosities in dwarf and $L\star$ galaxies. We show that in models that match the $\gamma$-ray observations, most CRs escape low-gas-density galaxies (e.g.\ dwarfs) before significant collisional losses, while starburst galaxies are CR proton calorimeters. While adiabatic losses can be significant, they occur only after CRs escape galaxies, so they are only of secondary importance for $\gamma$-ray emissivities. Models where CRs are ``trapped'' in the star-forming disk have lower star formation efficiency, but these models are ruled out by $\gamma$-ray observations. For models with constant $\kappa$ that match the $\gamma$-ray observations, CRs form extended halos with scale heights of several kpc to several tens of kpc.

We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range log(Mh/Msun) ~ 10-12 , and we focus on z=0.25 and z=2. Owing to strong stellar feedback, only ~L* halos retain a baryon mass >~50% of their cosmic budget. Metals are more efficiently retained by halos, with a retention fraction >~50%. Across all masses and redshifts analyzed >~60% of the CGM mass originates as IGM accretion (some of which is associated with infalling halos). Overall, the second most important contribution is wind from the central galaxy, though gas ejected or stripped from satellites can contribute a comparable mass in ~L* halos. Gas can persist in the CGM for billions of years, resulting in well-mixed halo gas. Sight lines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ~L* halos, cool gas (T<10^4.7 K) is distributed on average preferentially along the galaxy plane, however with strong halo-to-halo variability. The metallicity of IGM accretion is systematically lower than the metallicity of winds (typically by >~1 dex), although CGM and IGM metallicities depend significantly on the treatment of subgrid metal diffusion. Our results highlight the multiple physical mechanisms that contribute to the CGM and will inform observational efforts to develop a cohesive picture.

We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of "satellite vs. central" and "LG vs. individual MWvs. isolated dwarf central." Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the "near field" have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated ("true field") dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_\ast = 10^6 - 10^7\,M_\odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG.

We study a suite of extremely high-resolution cosmological FIRE simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}$$M_{\odot}$), run to $z=0$ with $30 M_{\odot}$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6} M_{\odot}$ is populated by a resolved {\em stellar} galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast}<10^{5}\,M_{\odot}$) have their star formation truncated early ($z\gtrsim2$), likely by reionization, while classical dwarfs ($M_{\ast}>10^{5} M_{\odot}$) continue forming stars to $z<0.5$. The systems have bursty star formation (SF) histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with $M_{\ast}/M_{\rm halo} > 10^{-4}$ to form a dark matter core $>200$pc, while lower-mass UFDs exhibit cusps down to $\lesssim100$pc, as expected from energetic arguments. Our dwarfs with $M_{\ast}>10^{4}\,M_{\odot}$ have half-mass radii ($R_{\rm 1/2}$) in agreement with Local Group (LG) dwarfs; dynamical mass vs. $R_{1/2}$ and the degree of rotational support also resemble observations. The lowest-mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).

We study how jets driven by active galactic nuclei influence the cooling flow in Perseus-like galaxy cluster cores with idealised, non-relativistic, hydrodynamical simulations performed with the Eulerian code ATHENA using high-resolution Godunov methods with low numerical diffusion. We use novel analysis methods to measure the cooling rate, the heating rate associated to multiple mechanisms, and the power associated with adiabatic compression/expansion. A significant reduction of the cooling rate and cooling flow within 20 kpc from the centre can be achieved with kinetic jets. However, at larger scales and away from the jet axis, the system relaxes to a cooling flow configuration. Jet feedback is anisotropic and is mostly distributed along the jet axis, where the cooling rate is reduced and a significant fraction of the jet power is converted into kinetic power of heated outflowing gas. Away from the jet axis weak shock heating represents the dominant heating source. Turbulent heating is significant only near the cluster centre, but it becomes inefficient at 50 kpc scales where it only represents a few percent of the total heating rate. Several details of the simulations depend on the choice made for the hydro solver, a consequence of the difficulty of achieving proper numerical convergence for this problem: current physics implementations and resolutions do not properly capture multi-phase gas that develops as a consequence of thermal instability. These processes happen at the grid scale and leave numerical solutions sensitive to the properties of the chosen hydro solver.

We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This algorithm ensures manifest conservation of mass, energy, and momentum, and avoids imprinting 'preferred directions' on the ejecta. We show that it is critical to incorporate both momentum and thermal energy of mechanical ejecta in a self-consistent manner, accounting for SNe cooling radii when they are not resolved. Using idealized simulations of single SN explosions, we show that the FIRE-2 algorithm, independent of resolution, reproduces converged solutions in both energy and momentum. In contrast, common 'fully-thermal' (energy-dump) or 'fully-kinetic' (particle-kicking) schemes in the literature depend strongly on resolution: when applied at mass resolution >100 solar masses, they diverge by orders-of-magnitude from the converged solution. In galaxy-formation simulations, this divergence leads to orders-of-magnitude differences in galaxy properties, unless those models are adjusted in a resolution-dependent way. We show that all models that individually time-resolve SNe converge to the FIRE-2 solution at sufficiently high resolution. However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments (FIRE) physics to both Local Group (LG)-like and isolated Milky Way (MW)-like volumes (ten host systems in total with baryonic particle mass $\simeq 3,500-7,000\,M_\odot$). We study the stellar mass functions, circular velocity or mass profiles, and velocity dispersions of the dwarf galaxy populations. The simulations reproduce the stellar mass function and central densities of MW satellite dwarfs for $M_\ast \geq10^{5.5}\,M_\odot$ and predict the existence of $\sim3$ unidentified galaxies with $M_\ast\sim10^5\,M_\odot$ within $300$ kpc of the MW. Overall, we find no evidence for the classical missing satellites or too-big-to-fail (TBTF) problems for satellite galaxies in our sample. Among the satellites, TBTF is resolved primarily by subhalo disruption and overall mass loss; central density profiles of subhalos are of secondary importance. For non-satellite galaxies, our LG-like simulations predict as many as $\sim10$ as-of-yet unseen galaxies at distances $0.3-1$ Mpc from both hosts, with $M_\ast\simeq10^{5-6}\,M_\odot$ (in halos with $V_\mathrm{max}\sim 20~\mathrm{km\,s}^{-1}$), albeit with large halo-to-halo variance. None of our simulations produces a compact, baryon-dominated, high-density dwarf elliptical-type galaxy (with $V_\mathrm{circ} \gtrsim 35~\mathrm{km\,s}^{-1}$ at $r<1\,$kpc), of which six may appear in the LG (but none in the MW). It may therefore remain a challenge to reproduce the full diversity of the dwarf population, including both the highest and lowest density systems.

We recently used hydro-chemical simulations to demonstrate that molecular outflows observed in luminous quasars can be explained by molecule formation within the AGN wind. However, these simulations cover a limited parameter space, due to their computational cost. We have therefore developed an analytic model to follow cooling in the shocked ISM layer of an AGN wind. We explore different ambient densities ($1-10^{4} \, \rm{cm}^{-3}$), density profile slopes ($0-1.5$), AGN luminosities ($10^{44}-10^{47} \, \rm{erg} \, \rm{s}^{-1}$), and metallicities ($0.1-3 \rm{Z}_{\odot}$). The swept up gas mostly cools within ~1 Myr. Based on our previous simulations, we predict that this gas would produce observable molecular outflows. The instantaneous momentum boost initially increases as the outflow decelerates. However, it reaches a maximum of $\approx$20, due to work done against the gravitational potential. The predicted time-averaged observational estimate of the molecular outflow momentum boost reaches a maximum of $\approx1-2$, partly due to our assumed molecular fraction, 0.2, but also because the instantaneous and observational, time-averaged definitions are not equivalent. Thus recent observational estimates of order unity momentum boosts do not necessarily rule out energy-driven outflows. Finally, we find that dust grains are likely to re-form by accretion of metals after the shocked ISM layer has cooled, assuming that a small fraction of dust grains swept up after this layer has cooled are able to mix into the cool phase, and assuming that grain growth remains efficient in the presence of the strong AGN radiation field. This would enable rapid molecule formation, as assumed in our models.

Galaxy formation simulations are an essential part of the modern toolkit of astrophysicists and cosmologists alike. Astrophysicists use the simulations to study the emergence of galaxy populations from the Big Bang, as well as problems including the formation of stars and supermassive black holes. For cosmologists, galaxy formation simulations are needed to understand how baryonic processes affect measurements of dark matter and dark energy. Owing to the extreme dynamic range of galaxy formation, advances are driven by novel approaches using simulations with different tradeoffs between volume and resolution. Large-volume but low-resolution simulations provide the best statistics, while higher resolution simulations of smaller cosmic volumes can be evolved with more self-consistent physics and reveal important emergent phenomena. I summarize recent progress in galaxy formation simulations, including major developments in the past five years, and highlight some key areas likely to drive further advances over the next decade.

We study the impact of a warm dark matter (WDM) cosmology on dwarf galaxy formation through a suite of cosmological hydrodynamical zoom-in simulations of $M_{\rm halo} \approx10^{10}\,M_{\odot}$ dark matter halos as part of the Feedback in Realistic Environments (FIRE) project. A main focus of this paper is to evaluate the combined effects of dark matter physics and stellar feedback on the well-known small-scale issues found in cold dark matter (CDM) models. We find that the $z=0$ stellar mass of a galaxy is strongly correlated with the central density of its host dark matter halo at the time of formation, $z_{\rm f}$, in both CDM and WDM models. WDM halos follow the same $M_{\star}(z=0)-V_{\rm max}(z_{\rm f})$ relation as in CDM, but they form later, are less centrally dense, and therefore contain galaxies that are less massive than their CDM counterparts. As a result, the impact of baryonic effects on the central gravitational potential is typically diminished relative to CDM. However, the combination of delayed formation in WDM and energy input from stellar feedback results in dark matter profiles with lower overall densities. The WDM galaxies studied here have a wider diversity of star formation histories (SFHs) than the same systems simulated in CDM, and the two lowest $M_{\star}$ WDM galaxies form all of their stars at late times. The discovery of young ultra-faint dwarf galaxies with no ancient star formation -- which do not exist in our CDM simulations -- would therefore provide evidence in support of WDM.

We explore the origin of fast molecular outflows that have been observed in Active Galactic Nuclei (AGN). Previous numerical studies have shown that it is difficult to create such an outflow by accelerating existing molecular clouds in the host galaxy, as the clouds will be destroyed before they can reach the high velocities that are observed. In this work, we consider an alternative scenario where molecules form in-situ within the AGN outflow. We present a series of hydro-chemical simulations of an isotropic AGN wind interacting with a uniform medium. We follow the time-dependent chemistry of 157 species, including 20 molecules, to determine whether molecules can form rapidly enough to produce the observed molecular outflows. We find H$_2$ outflow rates up to 140 M$_\odot$ yr$^{-1}$, which is sensitive to density, AGN luminosity, and metallicity. We compute emission and absorption lines of CO, OH and warm (a few hundred K) H$_2$ from the simulations in post-processing. The CO-derived outflow rates and OH absorption strengths at solar metallicity agree with observations, although the maximum line of sight velocities from the model CO spectra are a factor $\approx$2 lower than is observed. We derive a CO (1-0) to H$_2$ conversion factor of $\alpha_{\rm{CO} (1-0)}$ = 0.13 M$_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, 6 times lower than is commonly assumed in observations of such systems. We find strong emission from the mid-infrared lines of H$_2$. The mass of H$_2$ traced by this infrared emission is within a few per cent of the total H$_2$ mass. This H$_2$ emission may be observable by JWST.

We propose a simple analytic model to understand when star formation is time-steady versus bursty in galaxies. Recent models explain the observed Kennicutt-Schmidt relation between star formation rate and gas surface densities in galaxies as resulting from a balance between stellar feedback and gravity. We argue that bursty star formation occurs when such an equilibrium cannot be stably sustained, and identify two regimes in which galaxy-scale star formation should be bursty: i) at high redshift (z>~1) for galaxies of all masses, and ii) at low masses (depending on gas fraction) for galaxies at any redshift. At high redshift, characteristic galactic dynamical timescales become too short for supernova feedback to effectively respond to gravitational collapse in galactic discs (an effect recently identified for galactic nuclei), whereas in dwarf galaxies star formation occurs in too few bright star-forming regions to effectively average out. Burstiness is also enhanced at high redshift owing to elevated gas fractions in the early Universe. Our model can thus explain the bursty star formation rates predicted in these regimes by recent high-resolution galaxy formation simulations, as well as the bursty star formation histories observationally-inferred in both local dwarf and high-redshift galaxies. In our model, bursty star formation is associated with particularly strong spatio-temporal clustering of supernovae. Such clustering can promote the formation of galactic winds and our model may thus also explain the much higher wind mass loading factors inferred in high-redshift massive galaxies relative to their z~0 counterparts.

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between halo assembly, galaxy formation, and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution. The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way ($M_\star/{\rm M}_\odot\approx 10^5-10^7$). This non-negligible spread in stellar mass at $z=0$ in halos within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity $V_{\rm max}$. Much of this dependence of $M_\star$ on a second parameter (beyond $M_{\rm halo}$) is a direct consequence of the $M_{\rm halo}\sim10^{10}\,{\rm M}_\odot$ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming halos remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Halos that form galaxies with $M_\star\gtrsim2\times10^{6}\,{\rm M}_\odot$ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius $r_{1/2}$. This apparent stellar mass threshold of $M_\star \approx 2\times 10^{6} \approx 2\times 10^{-4} \,M_{\rm halo}$ is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in LCDM.

We present a deep search for HI 21-cm emission from the gaseous halo of Messier 31 as part of Project AMIGA, a large program Hubble Space Telescope program to study the circumgalactic medium of the Andromeda galaxy. Our observations with the Robert C. Byrd Green Bank Telesope target sight lines to 48 background AGNs, more than half of which have been observed in the ultraviolet with the Cosmic Origins Spectrograph, with impact parameters $25 \lesssim \rho \lesssim 330$ kpc ($0.1 \lesssim \rho / R_{\rm vir} \lesssim 1.1$). We do not detect any 21-cm emission toward these AGNs to limits of $N({\rm HI}) \approx 4 \times10^{17}$ cm$^{-2}$ ($5\sigma$, per 2 kpc diameter beam). This column density corresponds to an optical depth of $\sim2.5$ at the Lyman limit, thus our observations overlap with absorption line studies of Lyman limit systems at higher redshift. Our non-detections place a limit on the covering factor of such optically-thick gas around M31 to $f_c < 0.051$ (at 90\% confidence) for $\rho \leq R_{\rm vir}$. While individual clouds have previously been found in the region between M31 and M33, the covering factor of strongly optically-thick gas is quite small. Our upper limits on the covering factor are consistent with expectations from recent cosmological "zoom" simulations. Recent COS-Halos ultraviolet measurements of \HI\ absorption about an ensemble of galaxies at $z \approx 0.2$ show significantly higher covering factors within $\rho \lesssim 0.5 R_{\rm vir}$ at the same $N({\rm H I})$, although the metal ion-to-H I ratios appear to be consistent with those seen in M31.

Among the most important goals in cosmology is detecting and quantifying small ($M_{\rm halo}\simeq10^{6-9}~\mathrm{M}_\odot$) dark matter (DM) subhalos. Current probes around the Milky Way (MW) are most sensitive to such substructure within $\sim20$ kpc of the halo center, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in $\Lambda$CDM using cosmological zoom-in baryonic simulations of MW-mass halos from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two halos run using (1) DM-only (DMO), (2) full baryonic physics, and (3) DM with an embedded disk potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain $\sim2\times$ as many subhalos within 100 kpc of the halo center; this excess is $\gtrsim5\times$ within 25 kpc. At $z=0$, the baryonic simulations are completely devoid of subhalos down to $3\times10^6~\mathrm{M}_\odot$ within $15$ kpc of the MW-mass galaxy, and fewer than 20 surviving subhalos have orbital pericenters <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disk potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhalos on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving subhalo population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a disk potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.

Using high-resolution simulations with explicit treatment of stellar feedback physics based on the FIRE (Feedback in Realistic Environments) project, we study how galaxy formation and the interstellar medium (ISM) are affected by magnetic fields, anisotropic Spitzer-Braginskii conduction and viscosity, and sub-grid metal diffusion from unresolved turbulence. We consider controlled simulations of isolated (non-cosmological) galaxies but also a limited set of cosmological "zoom-in" simulations. Although simulations have shown significant effects from these physics with weak or absent stellar feedback, the effects are much weaker than those of stellar feedback when the latter is modeled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs) . In contrast, removing stellar feedback leads to SFRs being over-predicted by factors of $\sim 10 -100$. Without feedback, neither galactic winds nor volume filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scale-heights with unphysically dense clumps congregating at the galactic center. With stellar feedback, a multi-phase, turbulent medium with galactic fountains and winds is established. At currently achievable resolutions and for the investigated halo mass range $10^{10}-10^{13} M_{\odot}$, the additional physics investigated here (MHD, conduction, viscosity, metal diffusion) have only weak ($\sim10\%$-level) effects on regulating SFR and altering the balance of phases, outflows, or the energy in ISM turbulence, consistent with simple equipartition arguments. We conclude that galactic star formation and the ISM are primarily governed by a combination of turbulence, gravitational instabilities, and feedback. We add the caveat that AGN feedback is not included in the present work.

We use cosmological hydrodynamic simulations with stellar feedback from the FIRE project to study the physical nature of Lyman limit systems (LLSs) at z<1. At these low redshifts, LLSs are closely associated with dense gas structures surrounding galaxies, such as galactic winds, dwarf satellites, and cool inflows from the intergalactic medium. Our analysis is based on 14 zoom-in simulations covering the halo mass range M_h~10^9-10^13 Msun at z=0, which we convolve with the dark matter halo mass function to produce cosmological statistics. We find that the majority of cosmologically-selected LLSs are associated with halos in the mass range 10^10 < M_h < 10^12 Msun. The incidence and HI column density distribution of simulated absorbers with columns 10^16.2 < N_HI < 2x10^20 cm^-2 are consistent with observations. High-velocity outflows (with radial velocity exceeding the halo circular velocity by a factor >~2) tend to have higher metallicities ([X/H] ~ -0.5) while very low metallicity ([X/H] < -2) LLSs are typically associated with gas infalling from the intergalactic medium. However, most LLSs occupy an intermediate region in metallicity-radial velocity space, for which there is no clear trend between metallicity and radial kinematics. Metal-enriched inflows arise in the FIRE simulations as a result of galactic winds that fall back onto galaxies at low redshift. The overall simulated LLS metallicity distribution has a mean (standard deviation) [X/H] = -0.9 (0.4) and does not show significant evidence for bimodality, in contrast to recent observational studies but consistent with LLSs arising from halos with a broad range of masses and metallicities.

We present an analysis of the flow of metals through the circumgalactic medium (CGM) in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to $L*$ galaxies. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that reach $0.25 R_{\rm vir}$. When measured at $0.25 R_{\rm vir}$ the metallicity of outflows is slightly higher than the interstellar medium (ISM) metallicity. Many metals thus reside in the CGM. Cooling and recycling from this reservoir determine the metal budget in the ISM. The outflowing metal flux decreases by a factor of $\sim2-5$ between $0.25 R_{\rm vir}$ and $R_{\rm vir}$. Furthermore, outflow metallicity is typically lower at $R_{\rm vir}$ owing to dilution of the remaining outflow by metal-poor material swept up from the CGM. The inflow metallicity at $R_{\rm vir}$ is generally low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, although they preferentially retain metal ejecta. In all but the least massive galaxy considered, a majority of the metals are within the halo at $z=0$. We measure the fraction of metals in CGM, ISM and stars, and quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated $L*$ galaxies is consistent with estimates from the COS halos survey, while for the other two it appears to be lower.

We use idealized three-dimensional hydrodynamic simulations of global galactic discs to study the launching of galactic winds by supernovae (SNe). The simulations resolve the cooling radii of the majority of supernova remnants (SNRs) and thus self-consistently capture how SNe drive galactic winds. We find that SNe launch highly supersonic winds with properties that agree reasonably well with expectations from analytic models. The energy loading ($\eta_E = \dot{E}_{\rm wind} / \dot{E}_{\rm SN}$) of the winds in our simulations are well converged with spatial resolution while the wind mass loading ($\eta_M = \dot{M}_{\rm wind} / \dot{M}_\star$) decreases with resolution at the resolutions we achieve. We present a simple analytic model based on the concept that SNRs with cooling radii greater than the local scale height breakout of the disc and power the wind. This model successfully explains the dependence (or lack thereof) of $\eta_E$ (and by extension $\eta_M$) on the gas surface density, star formation efficiency, disc radius, and the clustering of SNe. The winds in the majority of our simulations are weaker than expected in reality, likely due to the fact that we seed SNe preferentially at density peaks. Clustering SNe in time and space substantially increases the wind power.

We study the structure, age and metallicity gradients, and dynamical evolution using a cosmological zoom-in simulation of a Milky Way-mass galaxy from the Feedback in Realistic Environments project. In the simulation, stars older than 6 Gyr were formed in a chaotic, bursty mode and have the largest vertical scale heights (1.5-2.5 kpc) by z=0, while stars younger than 6 Gyr were formed in a relatively calm, stable disk. The vertical scale height increases with stellar age at all radii, because (1) stars that formed earlier were thicker "at birth", and (2) stars were kinematically heated to an even thicker distribution after formation. Stars of the same age are thicker in the outer disk than in the inner disk (flaring). These lead to positive vertical age gradients and negative radial age gradients. The radial metallicity gradient is neg- ative at the mid-plane, flattens at larger disk height |Z|, and turns positive above |Z|~1.5kpc. The vertical metallicity gradient is negative at all radii, but is steeper at smaller radii. These trends broadly agree with observations in the Milky Way and can be naturally understood from the age gradients. The vertical stellar density profile can be well-described by two components, with scale heights 200-500 pc and 1-1.5 kpc, respectively. The thick component is a mix of stars older than 4 Gyr which formed through a combination of several mechanisms. Our results also demonstrate that it is possible to form a thin disk in cosmological simulations even with strong stellar feedback.

Galactic accretion interacts in complex ways with gaseous halos, including galactic winds. As a result, observational diagnostics typically probe a range of intertwined physical phenomena. Because of this complexity, cosmological hydrodynamic simulations have played a key role in developing observational diagnostics of galactic accretion. In this chapter, we review the status of different observational diagnostics of circumgalactic gas flows, in both absorption (galaxy pair and down-the-barrel observations in neutral hydrogen and metals; kinematic and azimuthal angle diagnostics; the cosmological column density distribution; and metallicity) and emission (Lya; UV metal lines; and diffuse X-rays). We conclude that there is no simple and robust way to identify galactic accretion in individual measurements. Rather, progress in testing galactic accretion models is likely to come from systematic, statistical comparisons of simulation predictions with observations. We discuss specific areas where progress is likely to be particularly fruitful over the next few years.

Recent spatially resolved observations of galaxies at z=0.6-3 reveal that high-redshift galaxies show complex kinematics and a broad distribution of gas-phase metallicity gradients. To understand these results, we use a suite of high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, which include physically motivated models of the multi-phase ISM, star formation, and stellar feedback. Our simulations reproduce the observed diversity of kinematic properties and metallicity gradients, broadly consistent with observations at z=0-3. Strong negative metallicity gradients only appear in galaxies with a rotating disk, but not all rotationally supported galaxies have significant gradients. Strongly perturbed galaxies with little rotation always have flat gradients. The kinematic properties and metallicity gradient of a high-redshift galaxy can vary significantly on short time-scales, associated with starburst episodes. Feedback from a starburst can destroy the gas disk, drive strong outflows, and flatten a pre-existing negative metallicity gradient. The time variability of a single galaxy is statistically similar to the entire simulated sample, indicating that the observed metallicity gradients in high-redshift galaxies reflect the instantaneous state of the galaxy rather than the accretion and growth history on cosmological time-scales. We find weak dependence of metallicity gradient on stellar mass and specific star formation rate (sSFR). Low-mass galaxies and galaxies with high sSFR tend to have flat gradients, likely due to the fact that feedback is more efficient in these galaxies. We argue that it is important to resolve feedback on small scales in order to produce the diverse metallicity gradients observed.

The morphology of massive star-forming galaxies at high redshift is often dominated by giant clumps of mass ~10^8-10^9 Msun and size ~100-1000 pc. Previous studies have proposed that giant clumps might have an important role in the evolution of their host galaxy, particularly in building the central bulge. However, this depends on whether clumps live long enough to migrate from their original location in the disc or whether they get disrupted by their own stellar feedback before reaching the centre of the galaxy. We use cosmological hydrodynamical simulations from the FIRE (Feedback in Realistic Environments) project that implement explicit treatments of stellar feedback and ISM physics to study the properties of these clumps. We follow the evolution of giant clumps in a massive (stellar mass ~10^10.8 Msun at z=1), discy, gas-rich galaxy from redshift z>2 to z=1. Even though the clumpy phase of this galaxy lasts over a gigayear, individual gas clumps are short-lived, with mean lifetime of massive clumps of ~20 Myr. During that time, they turn between 0.1% and 20% of their gas into stars before being disrupted, similar to local GMCs. Clumps with M>10^7 Msun account for ~20% of the total star formation in the galaxy during the clumpy phase, producing ~10^10 Msun of stars. We do not find evidence for net inward migration of clumps within the galaxy. The number of giant clumps and their mass decrease at lower redshifts, following the decrease in the overall gas fraction and star-formation rate.

We use cosmological simulations from the Feedback In Realistic Environments (FIRE) project, which implement a comprehensive set of stellar feedback processes, to study ultra-violet (UV) metal line emission from the circum-galactic medium of high-redshift (z=2-4) galaxies. Our simulations cover the halo mass range Mh ~ 2x10^11 - 8.5x10^12 Msun at z=2, representative of Lyman break galaxies. Of the transitions we analyze, the low-ionization C III (977 A) and Si III (1207 A) emission lines are the most luminous, with C IV (1548 A) and Si IV (1394 A) also showing interesting spatially-extended structures. The more massive halos are on average more UV-luminous. The UV metal line emission from galactic halos in our simulations arises primarily from collisionally ionized gas and is strongly time variable, with peak-to-trough variations of up to ~2 dex. The peaks of UV metal line luminosity correspond closely to massive and energetic mass outflow events, which follow bursts of star formation and inject sufficient energy into galactic halos to power the metal line emission. The strong time variability implies that even some relatively low-mass halos may be detectable. Conversely, flux-limited samples will be biased toward halos whose central galaxy has recently experienced a strong burst of star formation. Spatially-extended UV metal line emission around high-redshift galaxies should be detectable by current and upcoming integral field spectrographs such as the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope and Keck Cosmic Web Imager (KCWI).

Low-mass "dwarf" galaxies represent the most significant challenges to the cold dark matter (CDM) model of cosmological structure formation. Because these faint galaxies are (best) observed within the Local Group (LG) of the Milky Way (MW) and Andromeda (M31), understanding their formation in such an environment is critical. We present first results from the Latte Project: the Milky Way on FIRE (Feedback in Realistic Environments). This simulation models the formation of a MW-mass galaxy to z = 0 within LCDM cosmology, including dark matter, gas, and stars at unprecedented resolution: baryon particle mass of 7070 Msun with gas kernel/softening that adapts down to 1 pc (with a median of 25 - 60 pc at z = 0). Latte was simulated using the GIZMO code with a mesh-free method for accurate hydrodynamics and the FIRE-2 model for star formation and explicit feedback within a multi-phase interstellar medium. For the first time, Latte self-consistently resolves the spatial scales corresponding to half-light radii of dwarf galaxies that form around a MW-mass host down to Mstar > 10^5 Msun. Latte's population of dwarf galaxies agrees with the LG across a broad range of properties: (1) distributions of stellar masses and stellar velocity dispersions (dynamical masses), including their joint relation; (2) the mass-metallicity relation; and (3) a diverse range of star-formation histories, including their mass dependence. Thus, Latte produces a realistic population of dwarf galaxies at Mstar > 10^5 Msun that does not suffer from the "missing satellites" or "too big to fail" problems of small-scale structure formation. We conclude that baryonic physics can reconcile observed dwarf galaxies with standard LCDM cosmology.

Observations reveal that quasar host halos at z~2 have large covering fractions of cool dense gas (>~60% for Lyman limit systems within a projected virial radius). Most simulations have so far have failed to explain these large observed covering fractions. We analyze a new set of 15 simulated massive halos with explicit stellar feedback from the FIRE project, covering the halo mass range M_h~2x10^12-10^13 Msun at z=2. This extends our previous analysis of the circum-galactic medium of high-redshift galaxies to more massive halos. AGN feedback is not included in these simulations. We find Lyman limit system covering fractions consistent with those observed around quasars. The large HI covering fractions arise from star formation-driven galactic winds, including winds from low-mass satellite galaxies that interact with cosmological filaments. We show that it is necessary to resolve these satellite galaxies and their winds to reproduce the large Lyman limit system covering fractions observed in quasar-mass halos. Our simulations predict that galaxies occupying dark matter halos of mass similar to quasars but without a luminous AGN should have Lyman limit system covering fractions comparable to quasars.

We use local Cartesian simulations with a vertical gravitational potential to study how supernova (SN) feedback in stratified galactic discs drives turbulence and launches galactic winds. Our analysis includes three disc models with gas surface densities ranging from Milky Way-like galaxies to gas-rich ultra-luminous infrared galaxies (ULIRGs), and two different SN driving schemes (random and correlated with local gas density). In order to isolate the physics of SN feedback, we do not include additional feedback processes. We find that, in these local box calculations, SN feedback excites relatively low mass-weighted gas turbulent velocity dispersions ~3-7 km/s and low wind mass loading factors < 1 in all the cases we study. The low turbulent velocities and wind mass loading factors predicted by our local box calculations are significantly below those suggested by observations of gas-rich and rapidly star-forming galaxies; they are also in tension with global simulations of disc galaxies regulated by stellar feedback. Using a combination of numerical tests and analytic arguments, we argue that local Cartesian boxes cannot predict the properties of galactic winds because they do not capture the correct global geometry and gravitational potential of galaxies. The wind mass loading factors are in fact not well-defined in local simulations because they decline significantly with increasing box height. More physically realistic calculations (e.g., including a global galactic potential and disc rotation) will likely be needed to fully understand disc turbulence and galactic outflows, even for the idealized case of feedback by SNe alone.

Empirical constraints on reionization require galactic ionizing photon escape fractions fesc>20%, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ~1-5%. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the BPASS model). Binary mass transfer and mergers enhance the population of massive stars at late times (>3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the ISM, and so efficiently leak out of galaxies. As a result, the time-averaged "effective" escape fraction (ratio of escaped ionizing photons to observed 1500 A photons) increases by factors 4-10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

We study the interaction of feedback from active galactic nuclei (AGN) and a multi-phase interstellar medium (ISM), in simulations including explicit stellar feedback, multi-phase cooling, accretion-disk winds, and Compton heating. We examine radii ~0.1-100 pc around a black hole (BH), where the accretion rate onto the BH is determined and where AGN-powered winds and radiation couple to the ISM. We conclude: (1) The BH accretion rate is determined by exchange of angular momentum between gas and stars in gravitational instabilities. This produces accretion rates ~0.03-1 Msun/yr, sufficient to power luminous AGN. (2) The gas disk in the galactic nucleus undergoes an initial burst of star formation followed by several Myrs where stellar feedback suppresses the star formation rate (SFR). (3) AGN winds injected at small radii with momentum fluxes ~L/c couple efficiently to the ISM and have dramatic effects on ISM properties within ~100 pc. AGN winds suppress the nuclear SFR by factors ~10-30 and BH accretion rate by factors ~3-30. They increase the outflow rate from the nucleus by factors ~10, consistent with observational evidence for galaxy-scale AGN-driven outflows. (4) With AGN feedback, the predicted column density distribution to the BH is consistent with observations. Absent AGN feedback, the BH is isotropically obscured and there are not enough optically-thin sightlines to explain Type-I AGN. A 'torus-like' geometry arises self-consistently as AGN feedback evacuates gas in polar regions.

Stars in star clusters are thought to form in a single burst from a common progenitor cloud of molecular gas. However, massive, old globular clusters -- with ages greater than 10 billion years and masses of several hundred thousand solar masses -- often harbour multiple stellar populations, indicating that more than one star-forming event occurred during their lifetimes. Colliding stellar winds from late-stage, asymptotic-giant-branch stars are often invoked as second-generation star-formation trigger. The initial cluster masses should be at least 10 times more massive than they are today for this to work. However, large populations of clusters with masses greater than a few million solar masses are not found in the local Universe. Here we report on three 1-2 billion-year-old, massive star clusters in the Magellanic Clouds, which show clear evidence of burst-like star formation that occurred a few hundred million years after their initial formation era. We show that such clusters could accrete sufficient gas reservoirs to form new stars if the clusters orbited in their host galaxies' gaseous discs throughout the period between their initial formation and the more recent bursts of star formation. This may eventually give rise to the ubiquitous multiple stellar populations in globular clusters.

Quasar feedback models often predict an expanding hot gas bubble which drives a galaxy-scale outflow. In many circumstances this hot gas radiates inefficiently and is therefore difficult to observe directly. We present an indirect method to detect the presence of a hot bubble using hydrostatic photoionization calculations of the cold (~10^4 K) line-emitting gas. We compare our calculations with observations of the broad line region, the inner face of the torus, the narrow line region (NLR), and the extended NLR, and thus constrain the hot gas pressure at distances 0.1 pc -- 10 kpc from the center. We find that emission line ratios observed in the average quasar spectrum are consistent with radiation-pressure-dominated models on all scales. On scales <40 pc a dynamically significant hot gas pressure is ruled out, while on larger scales the hot gas pressure cannot exceed six times the local radiation pressure. In individual quasars, ~25% of quasars exhibit NLR ratios that are inconsistent with radiation-pressure-dominated models, though in these objects the hot gas pressure is also unlikely to exceed the radiation pressure by an order of magnitude or more. The derived upper limits on the hot gas pressure imply that the instantaneous gas pressure force acting on galaxy-scale outflows falls short of the time-averaged force needed to explain the large momentum fluxes \dot{p} >> L_AGN/c inferred for galaxy-scale outflows. This apparent discrepancy can be reconciled if optical quasars previously experienced a buried, fully-obscured phase during which the hot gas bubble was more effectively confined and during which galactic wind acceleration occurred.

The cosmic noon (z~1.5-3) marked a period of vigorous star formation for most galaxies. However, about a third of the more massive galaxies at those times were quiescent in the sense that their observed stellar populations are inconsistent with rapid star formation. The reduced star formation activity is often attributed to gaseous outflows driven by feedback from supermassive black holes, but the impact of black hole feedback on galaxies in the young Universe is not yet definitively established. We analyze the origin of quiescent galaxies with the help of ultra-high resolution, cosmological simulations that include feedback from stars but do not model the uncertain consequences of black hole feedback. We show that dark matter halos with specific accretion rates below ~0.25-0.4 per Gyr preferentially host galaxies with reduced star formation rates and red broad-band colors. The fraction of such halos in large dark matter only simulations matches the observed fraction of massive quiescent galaxies (~10^10-10^11 Msun). This strongly suggests that halo accretion rate is the key parameter determining which massive galaxies at z~1.5-3 become quiescent. Empirical models that connect galaxy and halo evolution, such as halo occupation distribution or abundance matching models, assume a tight link between galaxy properties and the masses of their parent halos. These models will benefit from adding the specific accretion rate of halos as a second model parameter.

Supernova (SN) feedback is one of the key processes shaping the interstellar medium (ISM) of galaxies. SNe contribute to (and in some cases may dominate) driving turbulence in the ISM and accelerating galactic winds. Modern cosmological simulations have sufficient resolution to capture the main structures in the ISM of galaxies, but are typically still not capable of explicitly resolving all of the small-scale stellar feedback processes, including the expansion of supernova remnants (SNRs). We perform a series of controlled three-dimensional hydrodynamic (adaptive mesh refinement) simulations of single SNRs expanding in an inhomogeneous density field with statistics motivated by those of the turbulent ISM. We use these to quantify the momentum and thermal energy injection from SNe as a function of spatial scale and the density, metallicity, and structure of the ambient medium. We develop a series of analytic formulae that we fit to the simulations. These formulae can be used as a basis for a more predictive sub-resolution model for SN feedback for galaxy formation simulations. We then use simulations of multiple, stochastically driven SNe that resolve the key phases of SNRs to test the sub-resolution model, and show that it accurately captures the turbulent kinetic energy and thermal energy in the ISM. By contrast, proposed SN feedback models in the literature based on `delayed cooling' significantly overpredict the late-time thermal energy and momentum in SNRs.

We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simulations resolves formation of star-forming giant molecular clouds to $z=0$, and features an explicit stellar feedback model on small scales. Our simulations reveal that high redshift galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows that eject much of the ISM and temporarily suppress star formation. At low redshift, however, sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch into a continuous and quiescent mode of star formation that does not drive outflows far into the halo. Mass-loading factors for winds in L*-progenitors are $\eta \approx 10$ at high redshift, but decrease to $\eta \ll 1$ at low redshift. Although lower values of $\eta$ are expected as halos grow in mass over time, we show that the strong suppression of outflows with decreasing redshift cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable and broadly distributed, but typically range between one and three times the circular velocity of the halo. Much of the ejected material builds a reservoir of enriched gas within the circumgalactic medium, some of which could be later recycled to fuel further star formation. However, a fraction of the gas that leaves the virial radius through galactic winds is never regained, causing most halos with mass $M_h \le 10^{12} M_{\odot}$ to be deficient in baryons compared to the cosmic mean by $z=0$.

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo=10^9-10^13 Msun and stellar mass Mstar=10^4-10^11 Msun at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass-metallicity relations at z=0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z=0-6 as log(Zgas/Zsun)=12+log(O/H)-9.0=0.35[log(Mstar/Msun)-10]+0.93 exp(-0.43 z)-1.05 and log(Zstar/Zsun)=[Fe/H]-0.2=0.40[log(Mstar/Msun)-10]+0.67 exp(-0.50 z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above Mstar=10^6 Msun are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between "sub-grid" wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.

Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).

We present a series of high-resolution (20-2000 Msun, 0.1-4 pc) cosmological zoom-in simulations at z~6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 10^9-10^11 Msun and rest-frame ultraviolet magnitude Muv = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01%-20%), while the time-averaged fesc over long time-scales generally remains ~5%, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3-10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g., models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that "sub-grid" star formation models, which do not explicitly resolve star formation in dense clouds with n >> 1 cm^-3, will dramatically over-predict fesc.

Massive early-type galaxies commonly have gas discs which are kinematically misaligned with the stellar component. These discs feel a torque from the stars and the angular momentum vectors are expected to align quickly. We present results on the evolution of a misaligned gas disc in a cosmological simulation of a massive early-type galaxy from the Feedback In Realistic Environments project. This galaxy experiences a merger which, together with a strong galactic wind, removes most of the original gas disc. The galaxy subsequently reforms a gas disc through accretion of cold gas, but it is initially 120 degrees misaligned with the stellar rotation axis. This misalignment persists for about 2 Gyr before the gas-star misalignment angle drops below 20 degrees. The time it takes for the gaseous and stellar components to align is much longer than previously thought, because the gas disc is accreting a significant amount of mass for about 1.5 Gyr after the merger, during which the angular momentum change induced by accreted gas dominates over that induced by stellar torques. Once the gas accretion rate has decreased sufficiently, the gas disc decouples from the surrounding halo gas and realigns with the stellar component in about 6 dynamical times. During the late evolution of the misaligned gas disc, the centre aligns faster than the outskirts, resulting in a warped disc. We discuss the observational consequences of the long survival of our misaligned gas disc and how our results can be used to calibrate merger rate estimates from observed gas misalignments.

We use high-resolution cosmological zoom-in simulations from the FIRE project to make predictions for the covering fractions of neutral hydrogen around galaxies at z=2-4. These simulations resolve the interstellar medium of galaxies and explicitly implement a comprehensive set of stellar feedback mechanisms. Our simulation sample consists of 16 main halos covering the mass range M_h~10^9-6x10^12 Msun at z=2, including 12 halos in the mass range M_h~10^11-10^12 Msun corresponding to Lyman break galaxies (LBGs). We process our simulations with a ray tracing method to compute the ionization state of the gas. Galactic winds increase the HI covering fractions in galaxy halos by direct ejection of cool gas from galaxies and through interactions with gas inflowing from the intergalactic medium. Our simulations predict HI covering fractions for Lyman limit systems (LLSs) consistent with measurements around z~2-2.5 LBGs; these covering fractions are a factor ~2 higher than our previous calculations without galactic winds. The fractions of HI absorbers arising in inflows and in outflows are on average ~50% but exhibit significant time variability, ranging from ~10% to ~90%. For our most massive halos, we find a factor ~3 deficit in the LLS covering fraction relative to what is measured around quasars at z~2, suggesting that the presence of a quasar may affect the properties of halo gas on ~100 kpc scales. The predicted covering fractions, which decrease with time, peak at M_h~10^11-10^12 Msun, near the peak of the star formation efficiency in dark matter halos. In our simulations, star formation and galactic outflows are highly time dependent; HI covering fractions are also time variable but less so because they represent averages over large areas.

We quantify the stellar abundances of neutron-rich r-process nuclei in cosmological zoom-in simulations of a Milky Way-mass galaxy from the Feedback In Realistic Environments project. The galaxy is enriched with r-process elements by binary neutron star (NS) mergers and with iron and other metals by supernovae. These calculations include key hydrodynamic mixing processes not present in standard semi-analytic chemical evolution models, such as galactic winds and hydrodynamic flows associated with structure formation. We explore a range of models for the rate and delay time of NS mergers, intended to roughly bracket the wide range of models consistent with current observational constraints. We show that NS mergers can produce [r-process/Fe] abundance ratios and scatter that appear reasonably consistent with observational constraints. At low metallicity, [Fe/H]<-2, we predict there is a wide range of stellar r-process abundance ratios, with both supersolar and subsolar abundances. Low-metallicity stars or stars that are outliers in their r-process abundance ratios are, on average, formed at high redshift and located at large galactocentric radius. Because NS mergers are rare, our results are not fully converged with respect to resolution, particularly at low metallicity. However, the uncertain rate and delay time distribution of NS mergers introduces an uncertainty in the r-process abundances comparable to that due to finite numerical resolution. Overall, our results are consistent with NS mergers being the source of most of the r-process nuclei in the Universe.

We present a series of high-resolution cosmological simulations of galaxy formation to z=0, spanning halo masses ~10^8-10^13 M_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of stellar feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_halo~10^12 M_sun (including dwarfs, satellites, MW-mass disks, and small groups). By extension, this leads to reasonable agreement with the stellar mass function for M_star<10^11 M_sun. We predict weak redshift evolution in the M_star-M_halo relation, consistent with current constraints to z>6. We find that the M_star-M_halo relation is insensitive to numerical details, but is sensitive to the feedback physics. Simulations with only supernova feedback fail to reproduce the observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to disrupt GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars in the simulation, and is instead determined by self-regulation via stellar feedback. Feedback reduces star formation rates considerably and produces a reservoir of gas that leads to rising late-time star formation histories significantly different from the halo accretion history. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. Many of these properties are not captured by common 'sub-grid' galactic wind models.

We predict the observational signatures of galaxy scale outflows powered by active galactic nuclei (AGN). Most of the emission is produced by the forward shock driven into the ambient interstellar medium (ISM) rather than by the reverse shock. AGN powered galactic winds with energetics suggested by phenomenological feedback arguments should produce spatially extended 1-10 keV X-ray emission of 10^(41-44) erg/s, significantly in excess of the spatially extended X-ray emission associated with normal star forming galaxies. The presence of such emission is a direct test of whether AGN outflows significantly interact with the ISM of their host galaxy. We further show that even radio quiet quasars should have a radio luminosity comparable to or in excess of the far infrared-radio correlation of normal star forming galaxies. This radio emission directly constrains the total kinetic energy flux in AGN-powered galactic winds. Radio emission from AGN wind shocks can also explain the recently highlighted correlations between radio luminosity and the kinematics of AGN narrow-line regions in radio quiet quasars.

We present a two-zone theory for feedback-regulated star formation in galactic discs, consistently connecting the galaxy-averaged star formation law with star formation proceeding in giant molecular clouds (GMCs). Our focus is on galaxies with gas surface density Sigma_g>~100 Msun pc^-2. In our theory, the galactic disc consists of Toomre-mass GMCs embedded in a volume-filling ISM. Radiation pressure on dust disperses GMCs and most supernovae explode in the volume-filling medium. A galaxy-averaged star formation law is derived by balancing the momentum input from supernova feedback with the gravitational weight of the disc gas. This star formation law is in good agreement with observations for a CO conversion factor depending continuously on Sigma_g. We argue that the galaxy-averaged star formation efficiency per free fall time, epsilon_ff^gal, is only a weak function of the efficiency with which GMCs convert their gas into stars. This is possible because the rate limiting step for star formation is the rate at which GMCs form: for large efficiency of star formation in GMCs, the Toomre Q parameter obtains a value slightly above unity so that the GMC formation rate is consistent with the galaxy-averaged star formation law. We contrast our results with other theories of turbulence-regulated star formation and discuss predictions of our model. Using a compilation of data from the literature, we show that the galaxy-averaged star formation efficiency per free fall time is non-universal and increases with increasing gas fraction, as predicted by our model. We also predict that the fraction of the disc gas mass in bound GMCs decreases for increasing values of the GMC star formation efficiency. This is qualitatively consistent with the smooth molecular gas distribution inferred in local ultra-luminous infrared galaxies and the small mass fraction in giant clumps in high-redshift galaxies.

Active galactic nuclei (AGN) drive fast winds in the interstellar medium of their host galaxies. It is commonly assumed that the high ambient densities and intense radiation fields in galactic nuclei imply short cooling times, thus making the outflows momentum-conserving. We show that cooling of high-velocity, shocked winds in AGN is in fact inefficient in a wide range of circumstances, including conditions relevant to ultra-luminous infrared galaxies (ULIRGs), resulting in energy-conserving outflows. We further show that fast energy-conserving outflows can tolerate a large amount of mixing with cooler gas before radiative losses become important. For winds with initial velocity v_in>~10,000 km s^-1, as observed in ultra-violet and X-ray absorption, the shocked wind develops a two-temperature structure. While most of the thermal pressure support is provided by the protons, the cooling processes operate directly only on the electrons. This significantly slows down inverse Compton cooling, while free free cooling is negligible. Slower winds with v_in~1,000 km s^-1, such as may be driven by radiation pressure on dust, can also experience energy-conserving phases but under more restrictive conditions. During the energy-conserving phase, the momentum flux of an outflow is boosted by a factor ~v_in/2v_s by work done by the hot post-shock gas, where v_s is the velocity of the swept-up material. Energy-conserving outflows driven by fast AGN winds (v_in~0.1c) may therefore explain the momentum fluxes Pdot>>L_AGN/c of galaxy-scale outflows recently measured in luminous quasars and ULIRGs. Shocked wind bubbles expanding normal to galactic disks may also explain the large-scale bipolar structures observed in some systems, including around the Galactic Center, and can produce significant radio, X-ray, and gamma-ray emission. [Abridged]

Recent observations have constrained the galaxy UV luminosity function up to z~10. However, these observations alone allow for a wide range of reionization scenarios due to uncertainties in the abundance of faint galaxies and the escape fraction of ionizing photons. We show that requiring continuity with post-reionization (z<6) measurements, where the Lya forest provides a complete probe of the cosmological emissivity of ionizing photons, significantly reduces the permitted parameter space. Models that are simultaneously consistent with the measured UV luminosity function, the Thomson optical depth to the CMB, and the Lya forest data require either: 1) extrapolation of the galaxy luminosity function down to very faint UV magnitudes M_lim ~ -10, corresponding roughly to the UV background suppression scale; 2) an increase of f_esc by a factor > ~10 from z=4 (where the best fit is 4%) to z=9; or 3) more likely, a hybrid solution in which undetected galaxies contribute significantly and f_esc increases more modestly. Models in which star formation is strongly suppressed in low-mass, reionization-epoch haloes of mass up to ~10^10 M_sun (e.g., owing to a metallicity dependence) are only allowed for extreme assumptions for the evolution of f_esc. However, variants of such models in which the suppression mass is reduced (e.g., assuming an earlier or higher metallicity floor) are in better agreement with the data. Concordance scenarios satisfying the available data predict a consistent redshift of 50% ionized fraction z_reion(50%) ~ 10. On the other hand, the duration of reionization is sensitive to the relative contribution of bright versus faint galaxies, with scenarios dominated by faint galaxies predicting a more extended reionization event. Scenarios relying too heavily on high-z dwarfs are disfavored by kinetic Sunyaev-Zeldovich measurements, which prefer a short reionization history.

Photoionization modeling of certain low-ionization broad absorption lines in quasars implies very compact (Delta R~0.01 pc), galaxy-scale (R kpc) absorbers blueshifted by several 1000 km s^-1. While these are likely signatures of quasar outflows, the lifetimes of such compact absorbers are too short for them to be direct ejecta from a nuclear wind. Instead, I argue that the absorbing clouds must be transient and created in situ. Following arguments detailed by Faucher-Giguere, Quataert, & Murray (2011), I show that a model in which the cool absorbers form in radiative shocks arising when a quasar blast wave impacts an interstellar cloud along the line of sight successfully explains the key observed properties. Using this radiative shock model, the outflow kinetic luminosities for three luminous quasars are estimated to be Edot,k~2-5% L_AGN (with corresponding momentum fluxes Pdot~2-15 L_AGN/c), consistent with feedback models of the M-sigma relation. These energetics are similar to those recently inferred of molecular outflows in local ultra-luminous infrared galaxies and in post-starburt winds, suggesting that active galactic nuclei (AGN) are capable of driving such outflows. Radiative shocks probably affect the multiphase structure of outflows in a range of other systems, including narrower and higher-ionization quasar absorption lines, and compact intergalactic absorbers ejected by star formation and AGN activity.

Photoionization modeling of the low-ionization broad absorption lines of certain quasars, known as FeLoBALs, has recently revealed the number density of the wind absorbers and their distance from the central supermassive black hole. From these, the feedback efficiency of the quasars can in principle be derived. The implied properties of the FeLoBALs are, however, surprising, with the thickness of the absorbers relative to their distance from the black hole, Delta R/R, as small as ~10^-5. Such absorbers are unlikely to survive the journey from the supermassive black hole to their inferred location. We show that the observed FeLoBAL properties are readily explained if they are formed in situ in radiative shocks produced when a quasar blast wave impacts a moderately dense interstellar clump along the line of sight. This physical picture differs significantly from the thin shell approximation often assumed, and implies outflow rates, kinetic luminosities and momentum fluxes that differ correspondingly, in some cases at the order of magnitude level. Using the radiative shock model, we estimate the ratio of the outflow kinetic luminosity to bolometric luminosity for three bright FeLoBAL quasars in the literature. We find Edot/Lbol~2-5% (and corresponding momentum fluxes Pdot~2-15 Lbol/c), similar to what is adopted in models reproducing the M-sigma relation. These outflow properties are also comparable to those recently inferred for molecular outflows in local ultra-luminous infrared galaxies, suggesting that active galactic nuclei are capable of driving such outflows.

We study the properties of self-shielding intergalactic absorption systems and their implications for the ionizing background. We find that cosmological simulations post-processed with detailed radiative transfer calculations generally are able to reproduce the observed abundance of Lyman-limit systems, and we highlight possible discrepancies between the observations and simulations. This comparison tests cosmological simulations at overdensities of ~100. Furthermore, we show that the properties of Lyman-limit systems in these simulations, in simple semi-analytic arguments, and as suggested by recent observations indicate that a small change in the ionizing emissivity of the sources would have resulted in a much larger change in the amplitude of the intergalactic HI-ionizing background (with this scaling strengthening with increasing redshift). This strong scaling could explain the rapid evolution in the Lyman-alpha forest transmission observed at z = 6. Our calculations agree with the suggestion of simpler models that the comoving ionizing emissivity was constant or even increasing from z = 3 to 6. Our calculations also provide a more rigorous estimate than in previous studies for the clumping factor of intergalactic gas after reionization, which we estimate was 2 - 3 at z = 6.

We use a suite of cosmological hydrodynamic simulations to quantify the accretion rates of baryons into dark matter halos and the resulting baryon mass fractions, as a function of halo mass, redshift, and baryon type (including cold and hot gas). We find that the net baryonic accretion rates through the virial radius are sensitive to galactic outflows and explore a range of outflow parameters to illustrate the effects. We show that the cold gas accretion rate is in general not a simple universal factor of the dark matter accretion rate, and that galactic winds can cause star formation rates to deviate significantly from the external gas accretion rates, both via gas ejection and re-accretion. Furthermore, galactic winds can inject enough energy and momentum in the surrounding medium to slow down accretion altogether, especially in low-mass halos and at low redshift. By resolving the accretion rates versus radius from the halo centers, we show how cold streams penetrate the hot atmospheres of massive halos at z>2, but gradually disappear at lower redshift. The total baryon mass fraction is also strongly suppressed by outflows in low-mass halos, but is nearly universal in the absence of feedback in halos above the UV background suppression scale. The transition halo mass, at which the gas mass in halos is equal for the cold and hot components, is roughly constant at ~10^11.5 Msun and does not depend sensitively on the wind prescription. We provide simple fitting formulae for the cold gas accretion rate into halos in the no-wind case. Finally, we show that cold accretion is broadly consistent with driving the bulk of the highly star-forming galaxies observed at z~2, but that the more intense star formers likely sample the high end of the accretion rate distribution, and may be additionally fueled by a combination of gas recycling, gas re-accretion, hot mode cooling, and mergers.

Binaries consisting of a pulsar and a black hole (BH) are a holy grail of astrophysics, both for their significance for stellar evolution and for their potential application as probes of strong gravity. In spite of extensive surveys of our Galaxy and its system of globular clusters, no pulsar-black hole (PSR-BH) binary has been found to date. Clues as to where such systems might exist are therefore highly desirable. We show that if the central parsec around Sgr A* harbors a cluster of ~25,000 stellar BHs (as predicted by mass segregation arguments) and if it is also rich in recycled pulsar binaries (by analogy with globular clusters), then 3-body exchange interactions should produce PSR-BHs in the Galactic center. Simple estimates of the formation rate and survival time of these binaries suggest that a few PSR-BHs should be present in the central parsec today. The proposed formation mechanism makes unique predictions for the PSR-BH properties: 1) the binary would reside within ~1 pc of Sgr A*; 2) the pulsar would be recycled, with a period of ~1 to a few tens of milliseconds, and a low magnetic field B<~10^10 G; 3) the binary would have high eccentricity, e~0.8, but with a large scatter; and 4) the binary would be relatively wide, with semi-major axis a_b~0.1 - >~3 AU. The potential discovery of a PSR-BH binary therefore provides a strong motivation for deep, high-frequency radio searches for recycled pulsars toward the Galactic center.

Theoretical models of galaxy formation predict that galaxies acquire most of their baryons via cold mode accretion. Observations of high-redshift galaxies, while showing ubiquitous outflows, have so far not revealed convincing traces of the predicted cold streams, which has been interpreted as a challenge for the current models. Using high-resolution, zoom-in smooth particle hydrodynamics simulations of Lyman break galaxy (LBG) halos combined with ionizing radiative transfer, we quantify the covering factor of the cold streams at z=2-4. We focus specifically on Lyman limit systems (LLSs) and damped Ly-alpha absorbers (DLAs), which can be probed by absorption spectroscopy using a background galaxy or quasar sightline, and which are closely related to low-ionization metal absorbers. We show that the covering factor of these systems is relatively small and decreases with time. At z=2, the covering factor of DLAs within the virial radius of the simulated galaxies is ~3% (~1% within twice this projected distance), and arises principally from the galaxy itself. The corresponding values for LLSs are ~10% and ~4%. Because of their small covering factor compared to the order unity covering fraction expected for galactic winds, the cold streams are naturally dominated by outflows in stacked spectra. We conclude that the existing observations are consistent with the predictions of cold mode accretion, and outline promising kinematic and chemical diagnostics to separate out the signatures of galactic accretion and feedback.

Recent studies have shown that galaxies accrete most of their baryons via the cold mode, from streams with temperatures T~10^4-10^5 K. At these temperatures, the streams should radiate primarily in the Lya line and have therefore been proposed as a model to power the extended, high-redshift objects known as Lya blobs and other high-redshift Lya sources. We introduce a new Lya radiative transfer code, aRT, and apply it to cosmological hydrodynamical simulations. We address physical and numerical issues that are critical to making accurate predictions for the cooling luminosity, but that have been mostly neglected or treated simplistically so far. We highlight the importance of self-shielding and of properly treating sub-resolution models in simulations. Most existing simulations do not self-consistently incorporate these effects, which can lead to order-of-magnitude errors in the predicted cooling luminosity. Using a combination of post-processing ionizing radiative transfer and re-simulation techniques, we develop an approximation to the consistent evolution of the self-shielded gas. We quantify the dependence of the Lya cooling luminosity on halo mass at z=3 for the simplified problem of pure gas accretion. While cooling in massive halos (without additional energy input from star formation and AGN) is in principle sufficient to produce L_alpha~10^43-10^44 erg s^-1 blobs, this requires including energy released in gas of density sufficient to form stars, but which is kept 100% gaseous in our optimistic estimates. Excluding emission from such dense gas yields lower luminosities by up to one to two orders of magnitude at high masses, making it difficult to explain the observed Lya blobs with pure cooling. Resonant scattering produces diffuse Lya halos, even for centrally concentrated emission, and broad double peaked line profiles. [Abridged]

When a spinning system experiences a transient gravitational encounter with an external perturber, a quasi-resonance occurs if the spin frequency of the victim matches the peak orbital frequency of the perturber. Such encounters are responsible for the formation of long tails and bridges of stars during galaxy collisions. For high-speed encounters, the resulting velocity perturbations can be described within the impulse approximation. The traditional impulse approximation, however, does not distinguish between prograde and retrograde encounters, and therefore completely misses the resonant response. Here, using perturbation theory, we compute the effects of quasi-resonant phenomena on stars orbiting within a disk. Explicit expressions are derived for the velocity and energy change to the stars induced by tidal forces from an external gravitational perturber passing either on a straight line or parabolic orbit. Comparisons with numerical restricted three-body calculations illustrate the applicability of our analysis.

We estimate the contribution of Galactic pulsars, both ordinary and millisecond pulsars (MSPs), to the high-energy (>100 MeV) gamma-ray background. We pay particular attention to the high-latitude part of the background that could be confused with an extragalactic component in existing analyses that subtract a Galactic cosmic-ray model. Our pulsar population models are calibrated to the results of large-scale radio surveys and we employ a simple empirical gamma-ray luminosity calibration to the spin-down rate that provides a good fit to existing data. We find that while ordinary pulsars are expected to contribute only a fraction ~10^-3 of the high-latitude gamma-ray intensity (I_X~1x10^-5 ph s^-1 cm^-2 sr^-1), MSPs could provide a much larger contribution and even potentially overproduce it, depending on the model parameters. We explore these dependences using a range of MSP models as a guide to how gamma-ray measurements can usefully constrain the MSP population. Existing gamma-ray background measurements and source counts already rule out several models. Finally, we show how fluctuations in the gamma-ray sky can be used to distinguish between different sources of the background.

The amplitude of fluctuations in the Ly-a forest on small spatial scales is sensitive to the temperature of the IGM and its spatial fluctuations. The temperature of the IGM and its spatial variations contain important information about hydrogen and helium reionization. We present a new measurement of the small-scale structure in the Ly-a forest from 40 high resolution, high signal-to-noise, VLT spectra at z=2.2-4.2. We convolve each Ly-a forest spectrum with a suitably chosen wavelet filter, which allows us to extract the amount of small-scale structure in the forest as a function of position across each spectrum. We compare these measurements with high resolution hydrodynamic simulations of the Ly-a forest which track more than 2 billion particles. This comparison suggests that the IGM temperature close to the cosmic mean density (T_0) peaks near z=3.4, at which point it is greater than 20,000 K at 2-sigma confidence. The temperature at lower redshift is consistent with the fall-off expected from adiabatic cooling ($T_0 \propto (1+z)^2$), after the peak temperature is reached near z=3.4. At z=4.2 our results favor a temperature of T_0 = 15-20,000 K. However, owing mostly to uncertainties in the mean transmitted flux at this redshift, a cooler IGM model with T_0 = 10,000 K is only disfavored at the 2-sigma level here, although such cool IGM models are strongly discrepant with the z ~ 3-3.4 measurement. We do not detect large spatial fluctuations in the IGM temperature at any redshift covered by our data set. The simplest interpretation of our measurements is that HeII reionization completes sometime near z ~ 3.4, although statistical uncertainties are still large [Abridged].

The ionizing background determines the ionization balance and the thermodynamics of the cosmic gas. It is therefore a fundamental ingredient to theoretical and empirical studies of both the IGM and galaxy formation. We present here a new calculation of its spectrum that satisfies the empirical constraints we recently obtained by combining state-of-the-art luminosity functions and intergalactic opacity measurements. In our preferred model, star-forming galaxies and quasars each contribute substantially to the HI ionizing field at z<3, with galaxies rapidly overtaking quasars at higher redshifts as quasars become rarer. In addition to our fiducial model, we explore the physical dependences of the calculated background and clarify how recombination emission contributes to the ionization rates. We find that recombinations do not simply boost the ionization rates by the number of reemitted ionizing photons as many of these rapidly redshift below the ionization edges and have a distribution of energies. A simple analytic model that captures the main effects seen in our numerical radiative transfer calculations is given. Finally, we discuss the effects of HeII reionization by quasars on both the spectrum of the ionizing background and on the thermal history of the IGM. In regions that have yet to be reionized, the spectrum is expected to be almost completely suppressed immediately above 54.4 eV while a background of higher-energy (>~0.5 keV) photons permeates the entire universe owing to the frequency-dependence of the photoionization cross section. We provide an analytic model of the heat input during HeII reionization and its effects on the temperature-density relation.

Observations of the intergalactic medium (IGM) suggest that quasars reionize HeII in the IGM at z ~ 3. We have run a set of 190 and 430 comoving Mpc simulations of HeII being reionized by quasars to develop an understanding of the nature of HeII reionization and its potential impact on observables. We find that HeII reionization heats regions in the IGM by as much as 25,000 K above the temperature that is expected otherwise, with the volume-averaged temperature increasing by ~ 12,000 K and with large temperature fluctuations on ~ 50 Mpc scales. Much of the heating occurs far from QSOs by hard photons. We find a temperature-density equation of state of gamma -1 ~ 0.3 during HeII reionization, but with a wide dispersion in this relation having sigma ~ 10^4 K. HeII reionization by the observed population of quasars cannot produce an inverted relation (gamma - 1 < 0). Our simulations are consistent with the observed evolution in the mean transmission of the HeII Ly-alpha forest. We argue that the heat input due to HeII reionization is unable to cause the observed depression at z = 3.2 in the HI Ly-alpha forest opacity as has been suggested. We investigate how uncertainties in the properties of QSOs and of HeII Lyman-limit systems influence our predictions.

We investigate the implications of the intergalactic opacity for the evolution of the cosmic UV luminosity density and its sources. Our main constraint is our measurement of the Lya forest opacity at redshifts 2<z<4.2 from 86 high-resolution quasar spectra. In addition, we impose the requirements that HI must be reionized by z=6 and HeII by z~3, and consider estimates of the hardness of the ionizing background from HI to HeII column density ratios. The derived hydrogen photoionization rate is remarkably flat over the Lya forest redshift range covered. Because the quasar luminosity function is strongly peaked near z~2, the lack of redshift evolution indicates that star-forming galaxies likely dominate the photoionization rate at z>~3. Combined with direct measurements of the galaxy UV luminosity function, this requires only a small fraction f_esc~0.5% of galactic hydrogen ionizing photons to escape their source for galaxies to solely account for the entire ionizing background. Under the assumption that the galactic UV emissivity traces the star formation rate, current state-of-the-art observational estimates of the star formation rate density appear to underestimate the total photoionization rate at z~4 by a factor ~4, are in tension with recent determinations of the UV luminosity function, and fail to reionize the Universe by z~6 if extrapolated to arbitrarily high redshift. A theoretical star formation history peaking earlier fits the Lya forest photoionization rate well, reionizes the Universe in time, and is in better agreement with the rate of z~4 gamma-ray bursts observed by Swift. Quasars suffice to doubly ionize helium by z~3 and likely contribute a non-negligible and perhaps dominant fraction of the hydrogen ionizing background at their z~2 peak. [Abridged]

We directly measure the evolution of the intergalactic Lya effective optical depth, tau_eff, over the redshift range 2<z<4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar spectra obtained with the ESI and HIRES spectrographs on Keck, and with the MIKE spectrograph on Magellan. This represents an improvement over previous analyses of the Lya forest from high-resolution spectra in this redshift interval of a factor of two in the size of the data set alone. We pay particular attention to robust error estimation and extensively test for systematic effects. We find that our estimates of the quasar continuum levels in the Lya forest obtained by spline fitting are systematically biased low, with the magnitude of the bias increasing with redshift, but that this bias can be accounted for using mock spectra. The mean fractional error is <1% at z=2, 4% at z=3, and 12% at z=4. Previous measurements of tau_eff at z>~3 based on directly fitting the quasar continua in the Lya forest, which have generally neglected this effect, are therefore likely biased low. We provide estimates of the level of absorption arising from metals in the Lya forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement, attaining 3% in redshift bins of width Delta z=0.2 around z=3, indicates significant departures from the best-fit power-law redshift evolution (tau_eff=0.0018(1+z)^3.92, when metals are left in), particularly near z=3.2. The observed downward departure is statistically consistent with a similar feature detected in a precision statistical measurement using Sloan Digital Sky Survey spectra by Bernardi and coworkers using an independent approach.

We investigate the implications of our measurement of the Lyman-alpha forest opacity at redshifts 2<z<4.2 from a sample of 86 high-resolution quasar spectra for the evolution of the cosmic ultraviolet luminosity density and its sources. The derived hydrogen photoionization rate is remarkably flat over this redshift range, implying an increasing comoving ionizing emissivity with redshift. Because the quasar luminosity function is strongly peaked near z~2, star-forming galaxies likely dominate the ionizing emissivity at z>~3. Our measurement argues against a star formation rate density declining beyond z~3, in contrast with existing state-of-the-art determinations of the cosmic star formation history from direct galaxy counts. Stellar emission from galaxies therefore likely reionized the Universe.

Binary pulsar systems are superb probes of stellar and binary evolution and the physics of extreme environments. In a survey with the Arecibo telescope, we have found PSR J1903+0327, a radio pulsar with a rotational period of 2.15 ms in a highly eccentric (e = 0.44) 95-day orbit around a solar mass companion. Infrared observations identify a possible main-sequence companion star. Conventional binary stellar evolution models predict neither large orbital eccentricities nor main-sequence companions around millisecond pulsars. Alternative formation scenarios involve recycling a neutron star in a globular cluster then ejecting it into the Galactic disk or membership in a hierarchical triple system. A relativistic analysis of timing observations of the pulsar finds its mass to be 1.74+/-0.04 Msun, an unusually high value.

The universe is permeated by a network of filaments, sheets, and knots collectively forming a "cosmic web.'' The discovery of the cosmic web, especially through its signature of absorption of light from distant sources by neutral hydrogen in the intergalactic medium, exemplifies the interplay between theory and experiment that drives science, and is one of the great examples in which numerical simulations have played a key and decisive role. We recount the milestones in our understanding of cosmic structure, summarize its impact on astronomy, cosmology, and physics, and look ahead by outlining the challenges faced as we prepare to probe the cosmic web at new wavelengths.

We investigate the birth and evolution of isolated radio pulsars using a population synthesis method, modeling the birth properties of the pulsars, their time evolution, and their detection in the Parkes and Swinburne Multibeam (MB) surveys. Together, the Parkes and Swinburne MB surveys have detected nearly 2/3 of the known pulsars and provide a remarkably homogeneous sample to compare with simulations. New proper motion measurements and an improved model of the distribution of free electrons in the interstellar medium, NE2001, also make revisiting these issues particularly worthwhile. We present a simple population model that reproduces the actual observations well, and consider others that fail. We conclude that: pulsars are born in the spiral arms, with the birthrate of 2.8+/-0.5 pulsars/century peaking at a distance ~3 kpc from the Galactic centre, and with mean initial speed of 380^{+40}_{-60} km/s; the birth spin period distribution extends to several hundred milliseconds, with no evidence of multimodality, implying that characteristic ages overestimate the true ages of the pulsars by a median factor >2 for true ages <30,000 yr; models in which the radio luminosities of the pulsars are random generically fail to reproduce the observed P-Pdot diagram, suggesting a relation between intrinsic radio luminosity and (P, Pdot); radio luminosities L Edot^0.5 provide a good match to the observed P-Pdot diagram; for this favored radio luminosity model, we find no evidence for significant magnetic field decay over the lifetime of the pulsars as radio sources ~100 Myr.

We directly measure the evolution of the intergalactic Lyman-alpha effective optical depth, tau_eff, over the redshift range 2<z<4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar spectra obtained with Keck/ESI, Keck/HIRES, and Magellan/MIKE. We find that our estimates of the quasar continuum levels in the Ly-alpha forest obtained by spline fitting are systematically biased low, but that this bias can be accounted for using mock spectra. The mean fractional error <Delta C/C_true> is <1% at z=2, 4% at z=3, and 12% at z=4. We provide estimates of the level of absorption arising from metals in the Ly-alpha forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement indicates significant departures from the best-fit power-law redshift evolution, particularly near z=3.2.

We investigate the birth and evolution of Galactic isolated radio pulsars. We begin by estimating their birth space velocity distribution from proper motion measurements of Brisken et al. (2002, 2003). We find no evidence for multimodality of the distribution and favor one in which the absolute one-dimensional velocity components are exponentially distributed and with a three-dimensional mean velocity of 380^{+40}_{-60} km s^-1. We then proceed with a Monte Carlo-based population synthesis, modelling the birth properties of the pulsars, their time evolution, and their detection in the Parkes and Swinburne Multibeam surveys. We present a population model that appears generally consistent with the observations. Our results suggest that pulsars are born in the spiral arms, with a Galactocentric radial distribution that is well described by the functional form proposed by Yusifov & Kucuk (2004), in which the pulsar surface density peaks at radius ~3 kpc. The birth spin period distribution extends to several hundred milliseconds, with no evidence of multimodality. Models which assume the radio luminosities of pulsars to be independent of the spin periods and period derivatives are inadequate, as they lead to the detection of too many old simulated pulsars in our simulations. Dithered radio luminosities proportional to the square root of the spin-down luminosity accommodate the observations well and provide a natural mechanism for the pulsars to dim uniformly as they approach the death line, avoiding an observed pile-up on the latter. There is no evidence for significant torque decay (due to magnetic field decay or otherwise) over the lifetime of the pulsars as radio sources (~100 Myr). Finally, we estimate the pulsar birthrate and total number of pulsars in the Galaxy.

We show that the Lyman-alpha optical depth statistics in the proximity regions of quasar spectra depend on the mass of the dark matter halos hosting the quasars. This is owing to both the overdensity around the quasars and the associated infall of gas toward them. For a fiducial quasar host halo mass of (3.0+/-1.6) h^-1 x 10^12 Msun, as inferred by Croom et al. from clustering in the 2dF QSO Redshift Survey, we show that estimates of the ionizing background (Gamma^bkg) from proximity effect measurements could be biased high by a factor of ~2.5 at z=3 owing to neglecting these effects alone. The clustering of galaxies and other active galactic nuclei around the proximity effect quasars enhances the local background, but is not expected to skew measurements by more than a few percent. Assuming the measurements of Gamma^bkg based on the mean flux decrement in the Ly-alpha forest to be free of bias, we demonstrate how the proximity effect analysis can be inverted to measure the mass of the dark matter halos hosting quasars. In ideal conditions, such a measurement could be made with a precision comparable to the best clustering constraints to date from a modest sample of only about 100 spectra. We discuss observational difficulties, including continuum flux estimation, quasar systematic redshift determination, and quasar variability, which make accurate proximity effect measurements challenging in practice. These are also likely to contribute to the discrepancies between existing proximity effect and flux decrement measurements of Gamma^bkg.

We report results from the initial stage of a long-term pulsar survey of the Galactic plane using the Arecibo L-band Feed Array (ALFA), a seven-beam receiver operating at 1.4 GHz with 0.3 GHz bandwidth. The search targets Galactic latitudes |b| < 5 deg in the longitude ranges 32 deg < l < 77 deg and 168 deg < l < 77 deg. Data discussed here were collected over a 100 MHz passband centered on 1.42 GHz using a spectrometer that recorded 256 channels every 64 microsec. In a preliminary, standard period-DM analysis, we have detected 29 previously known pulsars and discovered 11 new ones. One of these, with a period of 69 ms and a low characteristic age of 82 kyr, is a plausible candidate for association with the unidentified EGRET source 3EG J1928+1733. Another is a non-recycled pulsar in a relativistic binary with orbital period of 3.98 hr. We also search the data for isolated dispersed pulses, a technique that yielded discovery of an extremely sporadic radio emitter with a spin period of 1.2 s. Simulations we have carried out indicate that about 1000 new pulsars will be found in the ALFA survey. In addition to providing a large sample for use in population analyses and for probing the magnetoionic interstellar medium, the survey maximizes the chances of finding rapidly spinning millisecond pulsars and pulsars in compact binary systems. Our search algorithms will exploit the multiple data streams from ALFA to discriminate between radio frequency interference and celestial signals, including pulsars and possibly new classes of transient radio sources.