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

Matching the number counts of high-$z$ sub-millimetre-selected galaxies (SMGs) has been a long standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the sub-mm emission from galaxies in the SIMBA cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850$\mathrm{\mu m}$-selected sources. We find unprecedented agreement with the integrated 850$\mathrm{\mu m}$ luminosity function, along with good agreement in the redshift distribution of bright SMGs. The excellent agreement is driven primarily by SIMBA's good match to infrared measures of the star formation rate (SFR) function between $z = 2-4$ at high SFRs. Also important is the self-consistent on-the-fly dust model in SIMBA, which predicts, on average, higher dust masses (by up to a factor of 7) compared to using a fixed dust-to-metals ratio of 0.3. We construct a lightcone to investigate the effect of far-field blending, and find minimal contribution to the shape and normalisation of the luminosity function. We provide new fits to the 850$\mathrm{\mu m}$ luminosity as a function of SFR and dust mass. Our results demonstrate that exotic solutions to the discrepancy between sub-mm counts in simulations and observations, such as a top-heavy IMF, are unnecessary, and that sub-millimetre-bright phases are a natural consequence of massive galaxy evolution.

We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ mass functions (HIMF, H2MF), CO(1-0) luminosity functions (COLF), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density $\mu_*$, with some quantitative differences. To compare to the COLF, we apply an H$_2$-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in $\alpha_{\rm CO}$ and significant differences between models. Using this, predicted $z=0$ COLFs agree better with data than predicted H2MFs. Out to $z\sim 2$, EAGLE's and Simba's HIMF and COLF strongly increase, while TNG's HIMF declines and COLF evolves slowly. EAGLE and Simba reproduce high $L_{\rm CO1-0}$ galaxies at $z\sim 1-2$ as observed, owing partly to a median $\alpha_{\rm CO}(z=2)\sim 1$ versus $\alpha_{\rm CO}(z=0)\sim 3$. Examining \HI, H$_2$, and CO scaling relations, their trends with $M_*$ are broadly reproduced in all models, but EAGLE yields too little HI in green valley galaxies, TNG and Simba overproduce cold gas in massive galaxies, and Simba overproduces molecular gas in small systems. Using Simba variants that exclude individual AGN feedback modules, we find that Simba's AGN jet feedback is primarily responsible by lowering cold gas contents from $z\sim 1\to0$ by suppressing cold gas in $M_*> 10^{10}{\rm M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-$\mu_*$ systems.

Galaxies forming in low-mass halos are thought to be primarily responsible for reionizing the Universe during the first billion years after the Big Bang. Yet, these halos are extremely inefficient at forming stars in the nearby Universe. In this work, we address this apparent tension, and ask whether a physically motivated model of galaxy formation that reproduces the observed abundance of faint galaxies in the nearby Universe is also consistent with available observational constraints on the reionization history. By interfacing the Santa Cruz semi-analytic model for galaxy formation with an analytic reionization model, we constructed a computationally efficient pipeline that connects `ground-level' galaxy formation physics to `top-level' cosmological-scale observables. Based on photometric properties of the galaxy populations predicted up to $z=15$, we compute the reionization history of intergalactic hydrogen. We quantify the three degenerate quantities that influence the total ionizing photon budget, including the abundance of galaxies, the intrinsic production rate of ionizing photons, and the LyC escape fraction. We explore covariances between these quantities using a Markov chain Monte Carlo method. We find that our locally calibrated model is consistent with all currently available constraints on the reionization history, under reasonable assumptions about the LyC escape fraction. We quantify the fraction of ionizing photons contributed by galaxies of different luminosities and find that the galaxies expected to be detected in JWST NIRCam wide and deep surveys are responsible for producing $\sim 40$-$80\%$ of ionizing photons throughout the EoR. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.

The primary method for inferring the stellar mass ($M_*$) of a galaxy is through spectral energy distribution (SED) modeling. However, the technique rests on assumptions such as the galaxy star formation history and dust attenuation law that can severely impact the accuracy of derived physical properties from SED modeling. Here, we examine the effect that the assumed star formation history (SFH) has on the stellar properties inferred from SED fitting by ground truthing them against mock observations of high-resolution cosmological hydrodynamic galaxy formation simulations. Classically, SFHs are modeled with simplified parameterized functional forms, but these forms are unlikely to capture the true diversity of galaxy SFHs and may impose systematic biases with under-reported uncertainties on results. We demonstrate that flexible nonparametric star formation histories outperform traditional parametric forms in capturing variations in galaxy star formation histories, and as a result, lead to significantly improved stellar masses in SED fitting. We find a decrease in the average bias of 0.4 dex with a delayed-$\tau$ model to a bias of just under 0.05 dex for the nonparametric model. Similarly, using nonparametric star formation histories in SED fitting result in increased accuracy in recovered galaxy star formation rates (SFRs) and stellar ages.

We study specific star formation rate (sSFR) and gas profiles of star forming and green valley galaxies in the Simba cosmological hydrodynamic simulation. Star-forming galaxy half-light radii ($R_{\rm half}$) at $z=0$ and their evolution ($\propto(1+z)^{-0.78}$) agree with observations. Passive galaxy $R_{\rm half}$ agree with observations at high redshift, but by $z=0$ are too large, owing to numerical heating. We compare Simba $z=0$ sSFR radial profiles for star forming and green valley galaxies to observations. Simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in green valley galaxies and massive star-forming systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, Simba yields green valley galaxies with strongly dropping sSFR profiles beyond $\gtrsim R_{\rm half}$, regardless of AGN feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. Simba's satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At $z=2$ Simba does not show central depressions in massive star-forming galaxies, suggesting Simba's X-ray feedback should be more active at high-$z$. High resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in $z=0$ green valley galaxies represents a unique success of Simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms.

We present an analysis of the origin and properties of the circum-galactic medium (CGM) in a suite of 11 cosmological zoom simulations resembling present day spiral galaxies. On average the galaxies retain about 50\% of the cosmic fraction in baryons, almost equally divided into disc (interstellar medium) gas, cool CGM gas and warm-hot CGM gas. At radii smaller than 50 kpc the CGM is dominated by recycled warm-hot gas injected from the central galaxy, while at larger radii it is dominated by cool gas accreted onto the halo. The recycled gas typically accounts for one-third of the CGM mass. We introduce the novel publicly available analysis tool \textsc{pygad} to compute ion abundances and mock absorption spectra. For Lyman-${\alpha}$ absorption we find good agreement of the simulated equivalent width (EW) distribution and observations out to large radii. Disc galaxies with quiescent assembly histories show significantly more absorption along the disc major axis. By comparing the EW and HI column densities we find that CGM Lyman-${\alpha}$ absorbers are best represented by an effective line-width $b\approx 50 - 70$ km s$^{-1}$ that increases mildly with halo mass, larger than typically assumed.

We use the Simba cosmological hydrodynamic simulation suite to explore the impact of feedback on the circumgalactic medium (CGM) and intergalactic medium (IGM) around $2 \leq z \leq 3$ quasars. We identify quasars in Simba as the most rapidly-accreting black holes, and show that they are well-matched in bolometric luminosity and correlation strength to real quasars. We extract Lyman-alpha (Ly-a) absorption in spectra passing at different transverse distances (10 kpc $\lesssim b \lesssim$ 10 Mpc) around those quasars, and compare to observations of the mean Ly-a absorption profile. The observations are well reproduced, except within 100 kpc from the foreground quasar, where Simba overproduces absorption; this could potentially be mitigated by including ionisation from the quasar itself. By comparing runs with different feedback modules activated, we find that (mechanical) AGN feedback has little impact on the surrounding CGM even around these most highly luminous black holes, while stellar feedback has a significant impact. By further investigating thermodynamic and kinematic properties of CGM gas, we find that stellar feedback, and not AGN feedback, is the primary physical driver in determining the average properties of the CGM around $z\sim 2-3$ quasars. We also compare our results with previous works, and find that Simba predicts much more absorption within 100 kpc than the Nyx and Illustris simulations, showing that the Ly-a absorption profile can be a powerful constraint on simulations. Instruments such as VLT-MUSE and upcoming surveys (e.g., WEAVE and DESI) promise to further improve such constraints.

Only a handful of galaxies in the local Universe appear to be very young. We estimate the fraction of very young galaxies (VYGs), defined as those with more than half their stellar masses formed within the last Gyr. We fit non-parametric star formation histories (SFHs) to ~280 000 galaxy spectra from a flux- and volume-limited subsample of the Main Galaxy Sample (MGS) of the SDSS, which is also complete in mass-to-light ratio, thus properly accounting for passive galaxies of a given mass. The VYG fractions decrease with increasing galaxy stellar mass, from ~50% at $m = 10^8\,\rm M_\odot$ to ~0.1% at $m = 10^{11.5}\,\rm M_\odot$, with differences of up to 1 dex between the different spectral models used to estimate the SFH and on how we treat aperture effects. But old stellar populations may hide in our VYGs despite our conservative VYG sample built with galaxies that are globally bluer than within the region viewed by the SDSS fibre. The VYG fractions versus mass decrease more gradually compared to the Tweed et al. predictions using analytical and semi-analytical models of galaxy formation, but agree better with the SIMBA hydrodynamical simulation. These discrepancies highlight the usefulness of VYGs in constraining the strong uncertainties in both galaxy formation models and spectral modelling of galaxy SFHs. Given the lognormal cosmic SFH, these mean VYG fractions suggest that galaxies with $m > 10^8\,\rm M_\odot$ undergo at most 4 major starbursts on average.

We study the photometric properties and sizes of the reionization-epoch galaxies in high-resolution Simba cosmological hydrodynamical simulations with box sizes of $[25,50]~h^{-1}$Mpc. Assuming various attenuation laws, we compute photometry by extincting each star particle's spectrum using the line-of-sight gas metal column density. The predicted ultraviolet luminosity function (UVLF) generally agrees with observations at $z=6$, owing to a partial cancellation between the high metallicities of the simulated galaxies and lower dust-to-metal ratios. The simulated $z=8$ UVLF is low compared to observations, likely owing to excessive dust extinction. Simba predicts UV continuum slopes ($\beta$) in agreement with the $z=6$ observations, with the best agreement obtained using a Calzetti extinction law. Interestingly, the gas-phase mass-metallicity relation in Simba is higher at $z\sim 6$ than at $z\sim 2$, suggesting that rapid early enrichment (and dust growth) might be necessary to match the observed $\beta$. We find that $\beta$ is more sensitive to the dust extinction law than the UVLF. By generating mock James Webb Space Telescope (JWST) images and analysing in a manner similar to observations, we show that Simba's galaxy size-luminosity relation well reproduces the current $z=6$ Hubble observations. Unlike observations at lower redshifts, Simba predicts similar rest-UV and rest-optical sizes of $z=6$ galaxies, owing to weak age gradients and dust extinction in star-forming regions counteract each other to weaken the color gradients within galaxies. These predictions will be testable with JWST.

We examine X-ray scaling relations for massive halos ($M_{500}>10^{12.3}M_\odot$) in the Simba galaxy formation simulation. The X-ray weighted luminosity, $L_X$ vs. $M_{500}$ has power-law slopes $\approx\frac{5}{3}$ and $\approx\frac{8}{3}$ above and below $10^{13.5} M_{\odot}$, deviating from the self-similarity increasingly to low masses. $T_X-M_{500}$ is self-similar above this mass, and slightly shallower below it. Comparing Simba to observed $T_X$ scalings, we find that $L_X$, $L_X$-weighted [Fe/H], and entropies at $0.1 R_{200}$ ($S_{0.1}$) and $R_{500}$ ($S_{500}$) all match reasonably well. $S_{500}-T_X$ is consistent with self-similar expectations, but $S_{0.1}-T_X$ is shallower at lower $T_X$, suggesting the dominant form of heating moves from gravitational shocks in the outskirts to non-gravitational feedback in the cores of smaller groups. Simba matches observations of $L_X$ versus central galaxy stellar mass $M_*$, predicting the additional trend that star-forming galaxies have higher $L_X(M_*)$. The electron density profiles for $M_{500}>10^{14}M_\odot$ halos show a $\sim 0.1R_{200}$ core, but at lower masses the profiles are flat at all radii. $T_X$ profiles are flat for $M_{500}<10^{13}M_{\odot}$, and increasingly centrally peaked with mass. The resulting entropy profiles are centrally depressed for the highest and lowest mass halos but flatter for intermediate-mass halos, with $S_{\rm core}\approx200-400$ keV cm$^2$. Simba's [Fe/H] profile matches observations in the core but over-predicts at larger radii, suggesting overly widespread metal distribution. We show that Simba's bipolar jet AGN feedback is most responsible for increasingly evacuating lower-mass halos, thereby suppressing core densities and raising core entropies to improve agreement with X-ray scaling relations.

We investigate the Baryonic Tully-Fisher Relation (BTFR) in the $(100\,h^{-1}{\rm Mpc})^3$ Simba hydrodynamical galaxy formation simulation together with a higher-resolution $(25\,h^{-1}{\rm Mpc})^3$ Simba run, for over $10,000$ disk-dominated, HI-rich galaxies. We generate simulated galaxy rotation curves from the mass distribution, which we show yields similar results to using the gas rotational velocities. From this we measure the galaxy rotation velocity $V_{\rm circ}$ using four metrics: $V_{\rm max}, V_{\rm flat}, V_{2R_e},$ and $V_{\rm polyex}$. We compare the predicted BTFR to the SPARC observational sample and find broad agreement. In detail, however, Simba is biased towards higher $V_{\rm circ}$ by up to 0.1 dex. We find evidence for the flattening of the BTFR in $V_{\rm circ}>300$ km s$^{-1}$ galaxies, in agreement with recent observational findings. Simba's rotation curves are more peaked for lower mass galaxies, in contrast with observations, suggesting overly bugle-dominated dwarf galaxies in our sample. We present the residuals around the BTFR versus HI mass, stellar mass, gas fraction, and specific star formation rate, which provide testable predictions for upcoming BTFR surveys. Simba's BTFR shows sub-optimal resolution convergence, with the higher-resolution run lowering $V$ in better agreement with data.

We present the first results from our CGM at Cosmic Noon with KCWI program to study gas flows in the circumgalactic medium (CGM) at $z=2-3$. Combining the power of a high-resolution VLT/UVES quasar spectrum, an HST/ACS image, and integral field spectroscopy with Keck/KCWI, we detected Lya emission from a $3.7L_{\ast}$ galaxy at $z_{\rm gal}=2.0711$ associated with a Lyman limit system with weak MgII ($W_r(2796)=0.24$ Ang) in quasar field J143040$+$014939. The galaxy is best modeled as a star-forming (SFR$_{\rm FUV}=37.8$ M$_{\odot}$ yr$^{-1}$) clumpy edge-on disk ($i=85^{\circ}$). The background quasar probes the galaxy at an impact parameter of $D=66$ kpc along the projected galaxy minor axis ($\Phi=89^{\circ}$). From photoionization modeling of the absorption system, which includes Lya, MgII, SiII, SiIII, SiIV, and CIV, we infer a total line-of-sight CGM metallicity of ${\rm [Si/H]}=-1.5^{+0.4}_{-0.3}$. The absorption system is roughly kinematically symmetric about the galaxy redshift, with a full MgII velocity spread of $\sim210$ km s$^{-1}$. Given the galaxy--quasar orientation, CGM metallicity, and gas kinematics, we interpret this gas as outflowing from the host galaxy. By modeling the absorption as a polar outflow cone, we find the gas is decelerating with radial velocity $V_{\rm out}=100-350$ km s$^{-1}$ for half opening angles of $\theta_0=20^{\circ}-80^{\circ}$. Assuming a constant $V_{\rm out}$, it would take $\sim190-650$ Myr for the gas to reach 66 kpc. The outflow is energetic, with a mass outflow rate of $\dot{M}_{\rm out}=45-51$ M$_{\odot}$ yr$^{-1}$ and mass loading factor of $\eta\sim1.3$. We aim to build a sample of $\sim50$ MgII absorber--galaxy pairs at this epoch to better understand gas flows when they are most actively building galaxies.

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

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

We examine the properties of damped Lyman-$\alpha$ absorbers (DLAs) emerging from a single set of cosmological initial conditions in two state-of-the-art cosmological hydrodynamic simulations: {\sc Simba} and {\sc Technicolor Dawn}. The former includes star formation and black hole feedback treatments that yield a good match with low-redshift galaxy properties, while the latter uses multi-frequency radiative transfer to model an inhomogeneous ultraviolet background (UVB) self-consistently and is calibrated to match the Thomson scattering optical depth, UVB amplitude, and Ly-$\alpha$ forest mean transmission at $z>5$. Both simulations are in reasonable agreement with the measured stellar mass and star formation rate functions at $z\geq 3$, and both reproduce the observed neutral hydrogen cosmological mass density, $\Omega_{\rm HI}(z)$. However, the DLA abundance and metallicity distribution are sensitive to the galactic outflows' feedback and the UVB amplitude. Adopting a strong UVB and/or slow outflows under-produces the observed DLA abundance, but yields broad agreement with the observed DLA metallicity distribution. By contrast, faster outflows eject metals to larger distances, yielding more metal-rich DLAs whose observational selection may be more sensitive to dust bias. The DLA metallicity distribution in models adopting an ${\rm H}_2$-regulated star formation recipe includes a tail extending to $[M/H] \ll -3$, lower than any DLA observed to date, owing to curtailed star formation in low-metallicity galaxies. Our results show that DLA observations play an imporant role in constraining key physical ingredients in galaxy formation models, complementing traditional ensemble statistics such as the stellar mass and star formation rate functions.

We introduce the "Asymmetric Radiative Transfer In Shells Technique" (ARTIST), a new method for photon propagation on large scales that explicitly conserves photons, propagates photons at the speed of light, approximately accounts for photon directionality, and closely reproduces results of more detailed radiative transfer (RT) codes. Crucially, it is computationally fast enough to evolve the large cosmological volumes required to predict the 21cm power spectrum on scales that will be probed by future experiments targeting the Epoch of Reionisation (EoR). Most semi-numerical models aimed at predicting the EoR 21cm signal make use of an excursion set formalism (ESF) approach, which achieves computational viability by compromising on photon conservation, constraining ionised regions to be spherical by construction, and not accounting for light-travel time. By implementing our RT method within the semi-numerical code SimFast21, we show that ARTIST predicts a significantly different evolution for the EoR ionisation field compared to the code's native ESF. In particular, ARTIST predicts a more gradual evolution of the volume-averaged ionisation fraction, and up to an order-of-magnitude difference in the ionisation power, depending on the physical parameters assumed. Its application to large-scale EoR simulations will therefore allow more physically-motivated constraints to be obtained for key EoR parameters, such as the escape fraction.

We examine the impact of black hole jet feedback on the properties of the low-redshift intergalactic medium (IGM) in the Simba simulation, with a focus on the Ly$\alpha$ forest mean flux decrement $D_A$ and the inferred HI photo-ionisation rate $\Gamma_{HI}$. Without jet feedback, we confirm the Photon Underproduction Crisis (PUC; Kollmeier et al. 2014) in which $\Gamma_{HI}$ at $z=0$ must be increased by $\times 6$ over the Haardt & Madau (2012) value in order to match the observed $D_A$. Turning on jet feedback lowers this discrepancy to $\sim\times 2.5$, and additionally using the recent Faucher-Gigu\`ere 2019 background results in even better agreement, nearly solving the PUC. The PUC becomes apparent at late epochs ($z \sim < 1$) where the jet and no-jet simulations diverge; at higher redshifts Simba reproduces the observed $D_A$ with no adjustment, with or without jets. The main impact of jet feedback is to lower the cosmic baryon fraction in the diffuse IGM from 39% to 16% at $z=0$, while increasing the warm-hot intergalactic medium (WHIM) baryon fraction from 30% to 70%; the lowering of the diffuse IGM content directly translates into a lowering of $D_A$ by a similar factor. Comparing to the older Mufasa simulation that employs different quenching feedback but is otherwise similar to Simba, Mufasa matches $D_A$ less well than Simba, suggesting that low-redshift measurements of $D_A$ and $\Gamma_{HI}$ could provide constraints on feedback mechanisms. Our results suggest that widespread IGM heating at late times is a plausible solution to the PUC, and that Simba's jet AGN feedback model, constrained to reproduce quenched massive galaxies, approximately yields this required heating.

We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We firstly constrain the stellar mass vs stellar age relationship, finding a steep trend towards earlier average formation with increasing stellar mass of $1.48^{+0.34}_{-0.39}$ Gyr per decade in mass, although this shows signs of flattening at $M_* > 10^{11} M_\odot$. We show that this is consistent with other spectroscopic studies from $0 < z < 2$. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by Simba and IllustrisTNG agree well with observations at $z=0.1$, they are too shallow at $z=1$, predicting an evolution of $<0.5$ Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies ($M_* \sim 10^{10.5} M_\odot$) are consistent with formation in recent ($z<2$), intense starburst events, with timescales of $\lesssim500$ Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by $z=1$, they formed at $z\sim5$ and quenched by $z=3$. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity.

We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall).

Reionization leads to large spatial fluctuations in the intergalactic temperature that can persist well after its completion. We study the imprints of such fluctuations on the $z\sim5$ Ly$\alpha$ forest flux power spectrum using a set of radiation-hydrodynamic simulations that model different reionization scenarios. We find that large-scale coherent temperature fluctuations bring $\sim20-60\%$ extra power at $k\sim0.002$ s/km, with the largest enhancements in the models where reionization is extended or ends the latest. On smaller scales ($k\gtrsim0.1$ s/km), we find that temperature fluctuations suppress power by $\lesssim10\%$. We find that the shape of the power spectrum is mostly sensitive to the reionization midpoint rather than temperature fluctuations from reionization's patchiness. However, for all of our models with reionization midpoints of $z\le 8$ ($z\le 12$) the shape differences are $\lesssim20\%$ ($\lesssim40\%$) because of a surprisingly well-matched cancellation between thermal broadening and pressure smoothing that occurs for realistic thermal histories. We also consider fluctuations in the ultraviolet background, finding their impact on the power spectrum to be much smaller than temperature fluctuations at $k\gtrsim0.01$ s/km. Furthermore, we compare our models to power spectrum measurements, finding that none of our models with reionization midpoints of $z<8$ is strongly preferred over another and that all of our models with midpoints of $z\geq8$ are excluded at $2.5\sigma$. Future measurements may be able to distinguish between viable reionization models if they can be performed at lower $k$ or, alternatively, if the error bars on the high-$k$ power can be reduced by a factor of $1.5$.

The long-anticipated James Webb Space Telescope (JWST) will be able to directly detect large samples of galaxies at very high redshift. Using the well-established, computationally efficient Santa Cruz semi-analytic model, with recently implemented multiphase gas partitioning and H2-based star formation recipes, we make predictions for a wide variety of galaxy properties for galaxy populations at $z$ = 4-10. In this work, we provide forecasts for the physical properties of high-redshift galaxies and links to their photometric properties. With physical parameters calibrated only to $z\sim0$ observations, our model predictions are in good agreement with current observational constraints on stellar mass and star formation rate distribution functions up to $z \sim 8$. We also provide predictions representing wide, deep, and lensed JWST survey configurations. We study the redshift evolution of key galaxy properties and the scaling relations among them. Taking advantage of our models' high computational efficiency, we study the impact of systematically varying the model parameters. All distribution functions and scaling relations presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We present predictions for the evolution of the galaxy dust-to-gas (DGR) and dust-to-metal (DTM) ratios from z=0 to 6, using a model for the production, growth, and destruction of dust grains implemented into the \simba\ cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z=0, but modestly under-predicts the mass function by ~x3 at z ~ 1-2. The z=0 DGR vs metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR-metallicity relationship between z=0-6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion timescale, stellar mass and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM vs. the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.

We use the Simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\leq$ 4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass $M_*$ well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate sSFR$>t_{H}^{-1}$ to sSFR$<0.2t_{H}^{-1}$, where $t_{H}$ is the Hubble time. At $z\approx 0-3$, mergers show $\sim\times 2-3$ higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\leq 1\%$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the $H_2$ content, but for $M_*\geq 10^{10.5} M_{\odot}$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift $\propto (1+z)^{3.5}$, but the quenching rate evolves much more slowly, $\propto (1+z)^{0.9}$; there are insufficient mergers to explain the quenching rate at $z\leq 1.5$. Simba first quenches galaxies at $z\geq 3$, with a number density in good agreement with observations. The quenching timescales $\tau_q$ are strongly bimodal, with `slow' quenchings ($\tau_q \sim 0.1t_{H}$) dominating overall, but `fast' quenchings ($\tau_q\sim 0.01 t_H$) dominating in $M_*\sim 10^{10}-10^{10.5}M_{\odot}$ galaxies, likely induced by Simba's jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, Simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

We examine the co-evolution of galaxies and supermassive black holes in the Simba cosmological hydrodynamic simulation. Simba grows black holes via gravitational torque-limited accretion from cold gas and Bondi accretion from hot gas, while feedback from black holes is modeled in radiative and jet modes depending on the Eddington ratio ($f_{Edd}$). Simba shows generally good agreement with local studies of black hole properties, such as the black hole mass--stellar velocity dispersion ($M_{BH}-\sigma$) relation, 2 the black hole accretion rate vs. star formation rate (BHAR--SFR), and the black hole mass function. $M_{BH}-\sigma$ evolves such that galaxies at a given $M_{BH}$ have higher $\sigma$ at higher redshift, consistent with no evolution in $M_{BH}-M_*$. For $M_{BH}<\sim 10^8 M_\odot$, $f_{Edd}$ is anti-correlated with $M_{BH}$ since the BHAR is approximately independent of $M_{BH}$, while at higher masses $f_{Edd}-M_{BH}$ flattens and has a larger scatter. BHAR vs. SFR is invariant with redshift, but $f_{Edd}$ drops steadily with time at a given $M_{BH}$, such that all but the most massive black holes are accreting in a radiatively efficient mode at $z>\sim 2$. The black hole mass function amplitude decreases with redshift and is locally dominated by quiescent galaxies for $M_{BH}>10^{8}M_{\odot}$, but for $z>\sim 1$ star forming galaxies dominate at all $M_{BH}$. The $z=0$ $f_{Edd}$ distribution is roughly lognormal with a peak at $f_{Edd}<\sim 0.01$ as observed, shifting to higher $f_{Edd}$ at higher redshifts. Finally, we study the dependence of black hole properties with \HI\ content and find that the correlation between gas content and star formation rate is modulated by black hole properties, such that higher SFR galaxies at a given gas content have smaller black holes with higher $f_{Edd}$

We use machine learning to classify galaxies according to their HI content, based on both their optical photometry and environmental properties. The data used for our analyses are the outputs in the range $z = 0-1$ from MUFASA cosmological hydrodynamic simulation. In our previous paper, where we predicted the galaxy HI content using the same input features, HI rich galaxies were only selected for the training. In order for the predictions on real observation data to be more accurate, the classifiers built in this study will first establish if a galaxy is HI rich ($\rm{log(M_{HI}/M_{*})} > -2 $) before estimating its neutral hydrogen content using the regressors developed in the first paper. We resort to various machine learning algorithms and assess their performance with various metrics such as accuracy for instance. The performance of the classifiers gets better with increasing redshift and reaches their peak performance around $z = 1$. Random Forest method, the most robust among the classifiers when considering only the mock data for both training and test in this study, reaches an accuracy above $98.6 \%$ at $z = 0$ and above $99.0 \%$ at $z = 1$. We test our algorithms, trained with simulation data, on classification of the galaxies in RESOLVE, ALFALFA and GASS surveys. Interestingly, SVM algorithm, the best classifier for the tests, achieves a precision, the relevant metric for the tests, above $87.60\%$ and a specificity above $71.4\%$ with all the tests, indicating that the classifier is capable of learning from the simulated data to classify HI rich/HI poor galaxies from the real observation data. With the advent of large HI 21 cm surveys such as the SKA, this set of classifiers, together with the regressors developed in the first paper, will be part of a pipeline, a very useful tool, which is aimed at predicting HI content of galaxies.

We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($\eta=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $\eta \propto M_{\ast}^{\beta}$ with $\beta \simeq -0.4$.

We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Angl\'es-Alc\'azar et al. (2017) from cold gas and Bondi accretion from hot gas, along with black hole feedback via kinetic bipolar outflows and X-ray energy. Ejection velocities are taken to be ~10^3 km/s at high Eddington ratios, increasing to ~8000 km/s at Eddington ratios below 2%, with a constant momentum input of 20L/c. Simba further includes an on-the-fly dust production, growth, and destruction model. Our Simba run with (100 Mpc/h)^3 and 1024^3 gas elements reproduces numerous observables, including galaxy stellar mass functions at z=0-6, the stellar mass--star formation rate main sequence, HI and H2 fractions, the mass-metallicity relation at z=0 and z=2, star-forming galaxy sizes, hot gas fractions in massive halos, and z=0 galaxy dust properties. However, Simba also yields an insufficiently sharp truncation of the z=0 mass function, and too-large sizes for low-mass quenched galaxies. We show that Simba's jet feedback is primarily responsible for quenching massive galaxies.

One of the key questions in understanding the formation and evolution of galaxies is how starbursts affect the assembly of stellar populations in galaxies over time. We define a burst indicator ($\eta$), which compares a galaxy's star formation rates on short ($\sim10$ Myr) and long ($\sim100$ Myr) timescales. To estimate $\eta$, we apply the detailed time-luminosity relationship for H$\alpha$ and near-ultraviolet emission to simulated star formation histories (SFHs) from semi-analytic models and the Mufasa hydrodynamical cosmological simulations. The average of $\eta$ is not a good indicator of star formation stochasticity (burstiness); indeed, we show that this average should be close to zero unless the population has an average SFH which is rising or falling rapidly. Instead, the width of the $\eta$ distribution characterizes the burstiness of a galaxy population's recent star formation. We find this width to be robust to variations in stellar initial mass function and metallicity. We apply realistic noise and selection effects to the models to generate mock HST and JWST galaxy catalogs and compare these catalogs with 3D-HST observations of 956 galaxies at $0.65<z<1.5$ detected in H$\alpha$. Measurements of $\eta$ are unaffected by dust measurement errors under the assumption that $E(B-V)_\mathrm{stars}=0.44\,E(B-V)_\mathrm{gas}$ (i.e., $Q_\mathrm{sg}=0.44$). However, setting $Q_\mathrm{sg}=0.8^{+0.1}_{-0.2}$ removes an unexpected dependence of the average value of $\eta$ upon dust attenuation and stellar mass in the 3D-HST sample while also resolving disagreements in the distribution of star formation rates. However, even varying the dust law cannot resolve all discrepancies between the simulated and the observed galaxies.

Modelling the turbulent diffusion of thermal energy, momentum, and metals is required in all galaxy evolution simulations due to the ubiquity of turbulence in galactic environments. The most commonly employed diffusion model, the Smagorinsky model, is known to be over-diffusive due to its strong dependence on the fluid velocity shear. We present a method for dynamically calculating a more accurate, locally appropriate, turbulent diffusivity: the dynamic localised Smagorinsky model. We investigate a set of standard astrophysically-relevant hydrodynamical tests, and demonstrate that the dynamic model curbs over-diffusion in non-turbulent shear flows and improves the density contrast in our driven turbulence experiments. In galactic discs, we find that the dynamic model maintains the stability of the disc by preventing excessive angular momentum transport, and increases the metal-mixing timescale in the interstellar medium. In both our isolated Milky Way-like galaxies and cosmological simulations, we find that the interstellar and circumgalactic media are particularly sensitive to the treatment of turbulent diffusion. We also examined the global gas enrichment fractions in our cosmological simulations, to gauge the potential effect on the formation sites and population statistics of Population III stars and supermassive black holes, since they are theorised to be sensitive to the metallicity of the gas out of which they form. The dynamic model is, however, not for galaxy evolution studies only. It can be applied to all astrophysical hydrodynamics simulations, including those modelling stellar interiors, planetary formation, and star formation.

We investigate the connection between the HI~content, SFR and environment of galaxies using a hydrodynamic simulation that incorporates scaling relations for galactic wind and a heuristic halo mass-based quenching prescription. We run two zoom-in simulations of galaxy groups with $M_{halo}>10^{13}M_\odot$ at z=0, selected to have quiet merger histories. We track galaxies as they become satellites, and compute the delay time $\tau_{d}$ during which the satellites are similar to central galaxies at a given stellar mass, and a fading time $\tau_{f}$ during which satellites go from gas-rich and star-forming to gas-poor and quiescent. We find $0.7< \tau_{d}< 3$ Gyr at $z=0$, and depends inversely on the satellite halo mass at infall. At z ~ 1 we find ~$0.3< \tau_{d} < 2$ Gyr, broadly consistent with a positive correlation with the Hubble time. For a given halo mass, lower stellar mass galaxies at infall time have higher $\tau_{d}$. We generally find $\tau_{f}\ll \tau_{d}$, ranging between ~ 150 Myr at z~0 and ~ 80 Myr at z ~ 1 based on linear interpolation, with some uncertainty because they are smaller than our simulation output frequency ($200-300$ Myr). $\tau_{f}$ has no obvious dependency on infall halo mass. Both timescales show little difference between HI depletion and SF quenching, indicating that using up the gas reservoir by star formation without refilling is the main mechanism to transform satellite galaxies at these halo masses. At a given physical distance from the center of the main halo of interest, higher redshift galaxies have on average higher cold gas content, but the ratio of gas (HI or H$_2$) to star formation rate is similar, indicating that star formation is consistently fed through reservoirs of HI then H$_2$. For a given amount of HI, galaxies have shorter consumption times in more massive halo structures.

In anticipation of the upcoming deployment of the James Webb Space Telescope (JWST), we present high-redshift predictions by the well-established Santa Cruz semi-analytic model. We update the models by re-calibrating them after adopting cosmological parameters consistent with recent constraints from Planck. We provide predictions for rest-frame UV luminosity functions for galaxy populations over a wide range of $M_\text{UV}$ from $\sim-6$ to $\sim-24$ between $z = 4 - 10$. In addition, we present the corresponding predictions for observed-frame galaxy number counts in different redshift bins in the full set of NIRCam filters. We provide predictions of the quantitative effect on these observables of varying the physical recipes implemented in the models, such as the molecular gas depletion time (star formation efficiency) scalings or the scalings of outflow rates driven by stars and supernovae with galaxy circular velocity. Based on these results, we discuss what may be learned about the physical processes that shape galaxy formation from JWST observations of galaxy number densities at different intrinsic luminosities. All data tables for the results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulations result in lower mass, older, less metal rich and slower rotating systems with less disky isophotes - in general agreement with observations. Two-dimensional kinematic maps of in-situ and accreted stars show that these differences result from reduced in-situ star formation due to AGN feedback. A full analysis of stellar orbits indicates that galaxies simulated with AGN are typically more triaxial and have higher fractions of x-tubes and box orbits and lower fractions of z-tubes. This trend can also be explained by reduced late in-situ star formation. We introduce a global parameter, $\xi_3$ , to characterise the anti-correlation between the third-order kinematic moment $h_3$ and the line-of-sight velocity ($v_{los}/{\sigma}$), and compare to ATLAS$^{3D}$ observations. The kinematic asymmetry parameter $\xi_3$ might be a useful diagnostic for large integral field surveys as it is a kinematic indicator for intrinsic shape and orbital content.

We use a suite of cosmological zoom galaxy formation simulations and dust radiative transfer calculations to explore the use of the monochromatic $850~\mu m$ luminosity (L$_{\rm \nu,850}$) as a molecular gas mass (M$_{\rm mol}$) estimator in galaxies between $0 < z < 9.5$ for a broad range of masses. For our fiducial simulations, where we assume the dust mass is linearly related to the metal mass, we find that empirical L$_{\rm \nu,850}$-M$_{\rm mol}$ calibrations accurately recover the molecular gas mass of our model galaxies, and that the L$_{\rm \nu,850}$-dependent calibration is preferred. We argue the major driver of scatter in the L$_{\rm \nu,850}$-M$_{\rm mol}$ relation arises from variations in the molecular gas to dust mass ratio, rather than variations in the dust temperature, in agreement with the previous study of Liang et al. Emulating a realistic measurement strategy with ALMA observing bands that are dependent on the source redshift, we find that estimating S$_{\rm \nu,850}$ from continuum emission at a different frequency contributes $10-20\%$ scatter to the L$_{\rm \nu,850}$-M$_{\rm mol}$ relation. This additional scatter arises from a combination of mismatches in assumed T$_{dust}$ and $\beta$ values, as well as the fact that the SEDs are not single-temperature blackbodies.Finally we explore the impact of a dust prescription in which the dust-to-metals ratio varies with metallicity. Though the resulting mean dust temperatures are $\sim50\%$ higher, the dust mass is significantly decreased for low-metallicity halos. As a result, the observationally calibrated L$_{\rm \nu,850}$-M$_{\rm mol}$ relation holds for massive galaxies, independent of the dust model, but below L$_{\rm \nu,850}\lesssim10^{28}$ erg s$^{-1}$ (metallicities $\log_{10}({\rm Z}/{\rm Z}_{\odot})\lesssim -0.8$) we expect galaxies may deviate from literature observational calibrations by $\gtrsim0.5$ dex.

To achieve a fuller understanding of galaxy evolution, SED fitting can be used to recover quantities beyond stellar masses (M$_*$) and star formation rates (SFRs). We use Star Formation Histories (SFHs) reconstructed via the Dense Basis method of Iyer \& Gawiser (2017) for a sample of $17,873$ galaxies at $0.5<z<6$ in the CANDELS GOODS-S field to study the nature and evolution of the SFR-M$_*$ correlation. The reconstructed SFHs represent trajectories in SFR-M$_*$ space, enabling us to study galaxies at epochs earlier than observed by propagating them backwards in time along these trajectories. We study the SFR-M$_*$ correlation at $z=1,2,3,4,5,6$ using both direct fits to galaxies observed at those epochs and SFR-M$_*$ trajectories of galaxies observed at lower redshifts. The SFR-M$_*$ correlations obtained using the two approaches are found to be consistent with each other through a KS test. Validation tests using SFHs from semi-analytic models and cosmological hydrodynamical simulations confirm the sensitivity of the method to changes in the slope, normalization and shape of the SFR-M$_*$ correlation. This technique allows us to further probe the low-mass regime of the correlation at high-z by $\sim 1$ dex and over an effective volume of $\sim 10\times$ larger than possible with just direct fits. We find that the SFR-M$_*$ correlation is consistent with being linear down to M$_*\sim 10^7 M_\odot$ at $z>4$. The evolution of the correlation is well described by $\log SFR= (0.80\pm 0.029 - 0.017\pm 0.010\times t_{univ})\log M_*$ $- (6.487\pm 0.282-0.039\pm 0.008\times t_{univ})$, where $t_{univ}$ is the age of the universe in Gyr.

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.

We present Bayesian Analysis of Galaxies for Physical Inference and Parameter EStimation, or BAGPIPES, a new Python tool which can be used to rapidly generate complex model galaxy spectra and to fit these to arbitrary combinations of spectroscopic and photometric data using the MultiNest nested sampling algorithm. We extensively test our ability to recover realistic star-formation histories (SFHs) by fitting mock observations of quiescent galaxies from the MUFASA simulation. We then perform a detailed analysis of the SFHs of a sample of 9289 quiescent galaxies from UltraVISTA with stellar masses, $M_* > 10^{10}\ \mathrm{M_\odot}$ and redshifts $0.25 < z < 3.75$. The majority of our sample exhibit SFHs which rise gradually then quench relatively rapidly, over $1{-}2$ Gyr. This behaviour is consistent with recent cosmological hydrodynamic simulations, where AGN-driven feedback in the low-accretion (jet) mode is the dominant quenching mechanism. At $z > 1$ we also find a class of objects with SFHs which rise and fall very rapidly, with quenching timescales of $< 1$ Gyr, consistent with quasar-mode AGN feedback. Finally, at $z < 1$ we find a population with SFHs which quench more slowly than they rise, over $>3$ Gyr, which we speculate to be the result of diminishing overall cosmic gas supply. We confirm the mass-accelerated evolution (downsizing) trend, and a trend towards more rapid quenching at higher stellar masses. However, our results suggest that the latter is a natural consequence of mass-accelerated evolution, rather than a change in quenching physics with stellar mass. We find $61\pm8$ per cent of $z > 1.5$ massive quenched galaxies undergo significant further evolution by $z = 0.5$. BAGPIPES is available at https://bagpipes.readthedocs.io

We present the Technicolor Dawn simulations, a suite of cosmological radiation-hydrodynamic simulations of the first 1.2 billion years. By modeling a spatially-inhomogeneous UVB on-the-fly with 24 frequencies and resolving dark matter halos down to $10^8 M_\odot$ within 12 $h^{-1}$ Mpc volumes, our simulations unify observations of the intergalactic and circumgalactic media, galaxies, and reionization into a common framework. The only empirically-tuned parameter, the fraction $f_{\mathrm{esc,gal}}(z)$ of ionizing photons that escape the interstellar medium, is adjusted to match observations of the Lyman-$\alpha$ forest and the cosmic microwave background. With this single calibration, our simulations reproduce the history of reionization; the stellar mass-star formation rate relation of galaxies; the number density and metallicity of damped Lyman-$\alpha$ absorbers (DLAs) at $z\sim5$; the abundance of weak metal absorbers; the ultraviolet background (UVB) amplitude; and the Lyman-$\alpha$ flux power spectrum at $z=5.4$. The galaxy stellar mass and UV luminosity functions are underproduced by $\leq2\times$, suggesting an overly vigorous feedback model. The mean transmission in the Lyman-$\alpha$ forest is underproduced at $z<6$, indicating tension between measurements of the UVB amplitude and Lyman-$\alpha$ transmission. The observed SiIV column density distribution is reasonably well-reproduced ($\sim 1\sigma$ low). By contrast, CIV remains significantly underproduced despite being boosted by an intense $>4$ Ryd UVB. Solving this problem by increasing metal yields would overproduce both weak absorbers and DLA metallicities. Instead, the observed strength of high-ionization emission from high-redshift galaxies and absorption from their environments suggest that the ionizing flux from conventional stellar population models is too soft.

We investigate the nature of the relation among stellar mass, star-formation rate, and gas-phase metallicity (the M$_*$-SFR-Z relation) at high redshifts using a sample of 260 star-forming galaxies at $z\sim2.3$ from the MOSDEF survey. We present an analysis of the high-redshift M$_*$-SFR-Z relation based on several emission-line ratios for the first time. We show that a M$_*$-SFR-Z relation clearly exists at $z\sim2.3$. The strength of this relation is similar to predictions from cosmological hydrodynamical simulations. By performing a direct comparison of stacks of $z\sim0$ and $z\sim2.3$ galaxies, we find that $z\sim2.3$ galaxies have $\sim0.1$ dex lower metallicity at fixed M$_*$ and SFR. In the context of chemical evolution models, this evolution of the M$_*$-SFR-Z relation suggests an increase with redshift of the mass-loading factor at fixed M$_*$, as well as a decrease in the metallicity of infalling gas that is likely due to a lower importance of gas recycling relative to accretion from the intergalactic medium at high redshifts. Performing this analysis simultaneously with multiple metallicity-sensitive line ratios allows us to rule out the evolution in physical conditions (e.g., N/O ratio, ionization parameter, and hardness of the ionizing spectrum) at fixed metallicity as the source of the observed trends with redshift and with SFR at fixed M$_*$ at $z\sim2.3$. While this study highlights the promise of performing high-order tests of chemical evolution models at high redshifts, detailed quantitative comparisons ultimately await a full understanding of the evolution of metallicity calibrations with redshift.

We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as v_circ^(-2.2), with an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically about 1 Gyr, virtually independent of halo mass. Recycled material is re-accreted farther out in the disk and with typically about 2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most of cosmic star formation occurs.

We develop a machine learning (ML) framework to populate large dark matter-only simulations with baryonic galaxies. Our ML framework takes input halo properties including halo mass, environment, spin, and recent growth history, and outputs central galaxy and halo baryonic properties including stellar mass ($M_*$), star formation rate (SFR), metallicity ($Z$), neutral ($\rm HI$) and molecular ($\rm H_2$) hydrogen mass. We apply this to the MUFASA cosmological hydrodynamic simulation, and show that it recovers the mean trends of output quantities with halo mass highly accurately, including following the sharp drop in SFR and gas in quenched massive galaxies. However, the scatter around the mean relations is under-predicted. Examining galaxies individually, at $z=0$ the stellar mass and metallicity are accurately recovered ($\sigma\lesssim 0.2$~dex), but SFR and $\rm HI$ show larger scatter ($\sigma\gtrsim 0.3$~dex); these values improve somewhat at $z=1,2$. Remarkably, ML quantitatively recovers second parameter trends in galaxy properties, e.g. that galaxies with higher gas content and lower metallicity have higher SFR at a given $M_*$. Testing various ML algorithms, we find that none perform significantly better than the others, nor does ensembling improve performance, likely because none of the algorithms reproduce the large observed scatter around the mean properties. For the random forest algorithm, we find that halo mass and nearby ($\sim 200$~kpc) environment are the most important predictive variables followed by growth history, while halo spin and $\sim$Mpc scale environment are not important. Finally we study the impact of additionally inputting key baryonic properties $M_*$, SFR and $Z$, as would be available e.g. from an equilibrium model, and show that particularly providing the SFR enables $\rm HI$ to be recovered substantially more accurately.

We develop a machine learning-based framework to predict the HI content of galaxies using more straightforwardly observable quantities such as optical photometry and environmental parameters. We train the algorithm on z=0-2 outputs from the Mufasa cosmological hydrodynamic simulation, which includes star formation, feedback, and a heuristic model to quench massive galaxies that yields a reasonable match to a range of survey data including HI. We employ a variety of machine learning methods (regressors), and quantify their performance using the root mean square error ({\sc rmse}) and the Pearson correlation coefficient (r). Considering SDSS photometry, 3$^{rd}$ nearest neighbor environment and line of sight peculiar velocities as features, we obtain r $> 0.8$ accuracy of the HI-richness prediction, corresponding to {\sc rmse}$<0.3$. Adding near-IR photometry to the features yields some improvement to the prediction. Compared to all the regressors, random forest shows the best performance, with r $>0.9$ at $z=0$, followed by a Deep Neural Network with r $>0.85$. All regressors exhibit a declining performance with increasing redshift, which limits the utility of this approach to $z\la 1$, and they tend to somewhat over-predict the HI content of low-HI galaxies which might be due to Eddington bias in the training sample. We test our approach on the RESOLVE survey data. Training on a subset of RESOLVE data, we find that our machine learning method can reasonably well predict the HI-richness of the remaining RESOLVE data, with {\sc rmse}$\sim0.28$. When we train on mock data from Mufasa and test on RESOLVE, this increases to {\sc rmse}$\sim0.45$. Our method will be useful for making galaxy-by-galaxy survey predictions and incompleteness corrections for upcoming HI 21cm surveys such as the LADUMA and MIGHTEE surveys on MeerKAT, over regions where photometry is already available.

Non-parametric morphology measures are a powerful tool for identifying galaxy mergers at low redshifts. We employ cosmological zoom simulations using Gizmo with the Mufasa feedback scheme, post-processed using 3D dust radiative transfer into mock observations, to study whether common morphological measures Gini G, M20, concentration C, and asymmetry A are effective at identifying major galaxy mergers at z ~ 2 - 4, i.e. "Cosmic Noon". Our zoom suite covers galaxies with 10^8.6 < M_* < 10^11 M_sun at z ~ 2, and broadly reproduces key global galaxy observations. Our primary result is that these morphological measures are unable to robustly pick out galaxies currently undergoing mergers during Cosmic Noon, typically performing no better than a random guess. This improves only marginally if we consider whether galaxies have undergone a merger within the last Gyr. When also considering minor mergers, galaxies display no trend of moving towards the merger regime with increasing merger ratio. From z = 4 -> 2, galaxies move from the non-merger towards the merger regime in all statistics, but this is primarily an effect of mass: Above a given noise level, higher mass galaxies display a more complex outer morphology induced by their clustered environment. We conclude that during Cosmic Noon, these morphological statistics are of limited value in identifying galaxy mergers.

We present the extended GALEX Arecibo SDSS Survey (xGASS), a gas fraction-limited census of the atomic (HI) gas content of 1179 galaxies selected only by stellar mass ($M_\star =10^{9}-10^{11.5} M_\odot$) and redshift ($0.01<z<0.05$). This includes new Arecibo observations of 208 galaxies, for which we release catalogs and HI spectra. In addition to extending the GASS HI scaling relations by one decade in stellar mass, we quantify total (atomic+molecular) cold gas fractions and molecular-to-atomic gas mass ratios, $R_{mol}$, for the subset of 477 galaxies observed with the IRAM 30 m telescope. We find that atomic gas fractions keep increasing with decreasing stellar mass, with no sign of a plateau down to $\log M_\star/M_\odot = 9$. Total gas reservoirs remain HI-dominated across our full stellar mass range, hence total gas fraction scaling relations closely resemble atomic ones, but with a scatter that strongly correlates with $R_{mol}$, especially at fixed specific star formation rate. On average, $R_{mol}$ weakly increases with stellar mass and stellar surface density $\mu_\star$, but individual values vary by almost two orders of magnitude at fixed $M_\star$ or $\mu_\star$. We show that, for galaxies on the star-forming sequence, variations of $R_{mol}$ are mostly driven by changes of the HI reservoirs, with a clear dependence on $\mu_\star$. Establishing if galaxy mass or structure plays the most important role in regulating the cold gas content of galaxies requires an accurate separation of bulge and disk components for the study of gas scaling relations.

We use observations of spatially-aligned C II, C IV, Si II, Si IV, and O I absorbers to probe the slope and intensity of the ultraviolet background (UVB) at $z \sim 6$. We accomplish this by comparing observations with predictions from a cosmological hydrodynamic simulation using three trial UVBs applied in post-processing: a spectrally soft, fluctuating UVB calculated using multi-frequency radiative transfer; a soft, spatially-uniform UVB; and a hard, spatially-uniform "quasars-only" model. When considering our paired high-ionization absorbers (C IV/Si IV), the observed statistics strongly prefer the hard, spatially-uniform UVB. This echoes recent findings that cosmological simulations generically underproduce strong C IV absorbers at $z>5$. A single low/high ionization pair (Si II/Si IV), by contrast, shows a preference for the HM12 UVB, while two more (C II/C IV and O I/C IV) show no preference for any of the three UVBs. Despite this, future observations of specific absorbers, particularly Si IV/C IV, with next-generation telescopes probing to lower column densities should yield tighter constraints on the UVB.

We investigate galaxy conformity using the Mufasa cosmological hydrodynamical simulation. We show a bimodal distribution in galaxy colour with radius, albeit with too many low-mass quenched satellite galaxies compared to observations. Mufasa produces conformity in observed properties such as colour, sSFR, and HI content; i.e neighbouring galaxies have similar properties. We see analogous trends in other properties such as in environment, stellar age, H$_2$ content, and metallicity. We introduce quantifying conformity using $S(R)$, measuring the relative difference in upper and lower quartile properties of the neighbours. We show that low-mass and non-quenched haloes have weak conformity ($S(R)\leq 0.5$) extending to large projected radii $R$ in all properties, while high-mass and quenched haloes have strong conformity ($S(R)\sim 1$) that diminishes rapidly with $R$ and disappears at $R\geq 1$ Mpc. $S(R)$ is strongest for environment in low-mass haloes, and sSFR (or colour) in high-mass haloes, and is dominated by one-halo conformity with the exception of HI in small haloes. Metallicity shows a curious anti-conformity in massive haloes. Tracking the evolution of conformity for $z=0$ galaxies back in time shows that conformity broadly emerges as a late-time ($z\leq 1$) phenomenon. However, for fixed halo mass bins, conformity is fairly constant with redshift out to $z\geq 2$. These trends are consistent with the idea that strong conformity only emerges once haloes grow above Mufasa's quenching mass scale of $\sim 10^{12}M_\odot$. A quantitative measure of conformity in various properties, along with its evolution, thus represents a new and stringent test of the impact of quenching on environment within current galaxy formation models.

We introduce xCOLD GASS, a legacy survey providing a census of molecular gas in the local Universe. Building upon the original COLD GASS survey, we present here the full sample of 532 galaxies with CO(1-0) measurements from the IRAM-30m telescope. The sample is mass-selected in the redshift interval $0.01<z<0.05$ from SDSS, and therefore representative of the local galaxy population with M$_{\ast}>10^9$M$_{\odot}$. The CO(1-0) flux measurements are complemented by observations of the CO(2-1) line with both the IRAM-30m and APEX telescopes, HI observations from Arecibo, and photometry from SDSS, WISE and GALEX. Combining the IRAM and APEX data, we find that the CO(2-1) to CO(1-0) luminosity ratio for integrated measurements is $r_{21}=0.79\pm0.03$, with no systematic variations across the sample. The CO(1-0) luminosity function is constructed and best fit with a Schechter function with parameters {$L_{\mathrm{CO}}^* = (7.77\pm2.11) \times 10^9\,\mathrm{K\,km\,s^{-1}\, pc^{2}}$, $\phi^{*} = (9.84\pm5.41) \times 10^{-4} \, \mathrm{Mpc^{-3}}$ and $\alpha = -1.19\pm0.05$}. With the sample now complete down to stellar masses of $10^9$M$_{\odot}$, we are able to extend our study of gas scaling relations and confirm that both molecular gas fraction and depletion timescale vary with specific star formation rate (or offset from the star-formation main sequence) much more strongly than they depend on stellar mass. Comparing the xCOLD GASS results with outputs from hydrodynamic and semi-analytic models, we highlight the constraining power of cold gas scaling relations on models of galaxy formation.

Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionise the Universe. We assess the viability of this scenario using a semi-numerical framework for modeling reionisation, to which we add a quasar contribution by constructing a Quasar Halo Occupation Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an AGN-only model cannot simultaneously match observations of the optical depth $\tau_e$, neutral fraction, and ionising emissivity. Such a model predicts $\tau_e$ too low by $\sim 2\sigma$ relative to Planck constraints, and reionises the Universe at $z\lesssim 5$. Arbitrarily increasing the AGN emissivity to match these results yields a strong mismatch with the observed ionising emissivity at $z\sim 5$. If we instead assume a redshift-independent AGN luminosity function yielding an emissivity evolution like that assumed in Madau & Haardt model, then we can match $\tau_e$ albeit with late reionisation, however such evolution is inconsistent with observations at $z\sim 4-6$ and poorly motivated physically. These results arise because AGN are more biased towards massive halos than typical reionising galaxies, resulting in stronger clustering and later formation times. AGN-dominated models produce larger ionising bubbles that are reflected in $\sim\times 2$ more 21cm power on all scales. A model with equal parts galaxies and AGN contribution is still (barely) consistent with observations, but could be distinguished using next-generation 21cm experiments HERA and SKA-low. We conclude that, even with recent claims of more faint AGN than previously thought, AGN are highly unlikely to dominate the ionising photon budget for reionisation.

Of the almost 40 star forming galaxies at z>~5 (not counting QSOs) observed in [CII] to date, nearly half are either very faint in [CII], or not detected at all, and fall well below expectations based on locally derived relations between star formation rate (SFR) and [CII] luminosity. Combining cosmological zoom simulations of galaxies with SIGAME (SImulator of GAlaxy Millimeter/submillimeter Emission) we have modeled the multi-phased interstellar medium (ISM) and its emission in [CII], [OI] and [OIII], from 30 main sequence galaxies at z~6 with star formation rates ~3-23Msun/yr, stellar masses ~(0.7-8)x10^9Msun, and metallicities ~(0.1-0.4)xZsun. The simulations are able to reproduce the aforementioned [CII]-faintness at z>5, match two of the three existing z>~5 detections of [OIII], and are furthermore roughly consistent with the [OI] and [OIII] luminosity relations with SFR observed for local starburst galaxies. We find that the [CII] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ~66% and ~27%, respectively. The molecular gas, which constitutes only ~10% of the total gas mass is thus a more efficient emitter of [CII] than the ionized gas making up ~85% of the total gas mass. A principal component analysis shows that the [CII] luminosity correlates with the star formation activity as well as average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [CII]-faintness, and we suggest these factors may also be responsible for the [CII]-faint normal galaxies observed at these early epochs.

We examine the growth and evolution of quenched galaxies in the Mufasa cosmological hydrodynamic simulations that include an evolving halo mass-based quenching prescription, with galaxy colours computed accounting for line-of-sight extinction to individual star particles. Mufasa reproduces the observed present-day red sequence reasonably well, including its slope, amplitude, and scatter. In Mufasa, the red sequence slope is driven entirely by the steep stellar mass-stellar metallicity relation, which independently agrees with observations. High-mass star-forming galaxies blend smoothly onto the red sequence, indicating the lack of a well-defined green valley at M*>10^10.5 Mo. The most massive galaxies quench the earliest and then grow very little in mass via dry merging; they attain their high masses at earlier epochs when cold inflows more effectively penetrate hot halos. To higher redshifts, the red sequence becomes increasingly contaminated with massive dusty star-forming galaxies; UVJ selection subtly but effectively separates these populations. We then examine the evolution of the mass functions of central and satellite galaxies split into passive and star-forming via UVJ. Massive quenched systems show good agreement with observations out to z~2, despite not including a rapid early quenching mode associated with mergers. However, low-mass quenched galaxies are far too numerous at z<1 in Mufasa, indicating that Mufasa strongly over-quenches satellites. A challenge for hydrodynamic simulations is to devise a quenching model that produces enough early massive quenched galaxies and keeps them quenched to z=0, while not being so strong as to over-quench satellites; Mufasa's current scheme fails at the latter.

The recent discovery of gravitational radiation from merging black holes poses a challenge of how to organize the electromagnetic follow-up of gravitational-wave events as well as observed bursts of neutrinos. We propose a technique to select the galaxies that are most likely to host the event given some assumptions of whether the particular event is associated with recent star formation, low metallicity stars or simply proportional to the total stellar mass in the galaxy. We combine data from the 2-MASS Photometric Redshift Galaxy Catalogue with results from galaxy formation simulations to develop observing strategies that potentially reduce the area of sky to search by up to a factor of two relative to an unweighted search of galaxies, and a factor twenty to a search over the entire LIGO localization region.

We compare galaxy scaling relations as a function of environment at $z\sim2$ with our ZFIRE survey where we have measured H$\alpha$ fluxes for 90 star-forming galaxies selected from a mass-limited [$\log(M_{\star}/M_{\odot})>9$] sample based on ZFOURGE. The cluster galaxies (37) are part of a confirmed system at z=2.095 and the field galaxies (53) are at $1.9<z<2.4$; all are in the COSMOS legacy field. There is no statistical difference between H$\alpha$-emitting cluster and field populations when comparing their star formation rate (SFR), stellar mass ($M_{\star}$), galaxy size ($r_{eff}$), SFR surface density [$\Sigma$(H$\alpha_{star}$)], and stellar age distributions. The only difference is that at fixed stellar mass, the H$\alpha$-emitting cluster galaxies are $\log(r_{eff})\sim0.1$ larger than in the field. Approximately 19% of the H$\alpha$-emitters in the cluster and 26% in the field are IR-luminous ($L_{IR}>2\times10^{11} L_{\odot}$). Because the LIRGs in our combined sample are $\sim5$ times more massive than the low-IR galaxies, their radii are $\sim70$% larger. To track stellar growth, we separate galaxies into those that lie above, on, and below the H$\alpha$ star-forming main sequence (SFMS) using $\Delta$SFR$(M_{\star})=\pm0.2$ dex. Galaxies above the SFMS (starbursts) tend to have higher H$\alpha$ SFR surface densities and younger light-weighted stellar ages compared to galaxies below the SFMS. Our results indicate that starbursts (+SFMS) in the cluster and field at $z\sim2$ are growing their stellar cores. Lastly, we compare to the (SFR-$M_{\star}$) relation from RHAPSODY cluster simulations and find the predicted slope is nominally consistent with the observations. However, the predicted cluster SFRs tend to be too low by a factor of $\sim2$ which seems to be a common problem for simulations across environment.

We examine galaxy star formation rates (SFRs), metallicities, and gas contents predicted by the MUFASA cosmological hydrodynamic simulations, which employ meshless hydrodynamics and novel feedback prescriptions that yield a good match to observed galaxy stellar mass assembly. We combine 50, 25, and 12.5 Mpc/h boxes with a quarter billion particles each to show that MUFASA broadly reproduces a wide range of relevant observations, including SFR and specific SFR functions, the mass-metallicity relation, HI and H2 fractions, HI (21 cm) and CO luminosity functions, and cosmic gas density evolution. There are mild but significant discrepancies, such as too many high-SFR galaxies, overly metal-rich and HI-poor galaxies at M*>10^{10} Mo, and sSFRs that are too low at z~1-2. The HI mass function increases by x2 out to z~1 then steepens to higher redshifts, while the CO luminosity function computed using the Narayanan et al. conversion factor shows a rapid increase of CO-bright galaxies out to z~2 in accord with data. Omega_HI and Omega_H2 both scale roughly as (1+z)^0.7 out to z~3, comparable to the rise in HI and H2 fractions. MUFASA galaxies with high SFR at a given M* have lower metallicities and higher HI and H2 fractions, following observed trends; we make quantitative predictions for how fluctuations in the baryon cycle drive correlated scatter around galaxy scaling relations. Most of these trends are well converged with numerical resolution. These successes highlight MUFASA as a viable platform to study many facets of cosmological galaxy evolution.

The analytic "equilibrium model" for galaxy evolution using a mass balance equation is able to reproduce mean observed galaxy scaling relations between stellar mass, halo mass, star formation rate (SFR) and metallicity across the majority of cosmic time with a small number of parameters related to feedback. Here we aim to test this data-constrained model to quantify deviations from the mean relation between stellar mass and SFR, i.e. the star-forming galaxy main sequence (MS). We implement fluctuation in halo accretion rates parameterised from merger-based simulations, and quantify the intrinsic scatter introduced into the MS under the assumption that fluctuations in star formation follow baryonic inflow fluctuations. We predict the 1-sigma MS scatter to be ~ 0.2 - 0.25 dex over the stellar mass range 10^8 Mo to 10^11 Mo and a redshift range 0.5 < z < 3 for SFRs averaged over 100 Myr. The scatter increases modestly at z > 3, as well as by averaging over shorter timescales. The contribution from merger-induced star formation is generally small, around 5% today and 10 - 15% during the peak epoch of cosmic star formation. These results are generally consistent with available observations, suggesting that deviations from the MS primarily reflect stochasticity in the inflow rate owing to halo mergers.

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

Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection.

We present the MUFASA suite of cosmological hydrodynamic simulations, which employs the GIZMO meshless finite mass (MFM) code including H2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching. Our fiducial (50 Mpc/h)^3 volume is evolved to z=0 with a quarter billion particles, The predicted galaxy stellar mass functions (GSMF) reproduce observations from z=4-0 to <1.2sigma in cosmic variance, providing an unprecedented match to this key diagnostic. The cosmic star formation history and stellar mass growth show general agreement with data, with a strong archaeological downsizing trend such that dwarf galaxies form the majority of their stars after z~1. We run 25 Mpc/h and 12.5 Mpc/h volumes to z=2 with identical feedback prescriptions, the latter resolving all hydrogen-cooling halos, and the three runs display fair resolution convergence. The specific star formation rates broadly agree with data at z=0, but are underpredicted at z~2 by a factor of three, re-emphasizing a longstanding puzzle in galaxy evolution models. We compare runs using MFM and two flavours of Smoothed Particle Hydrodynamics, and show that the GSMF is sensitive to hydrodynamics methodology at the ~x2 level, which is sub-dominant to choices for parameterising feedback.

We perform a kinematic analysis of galaxies at $z\sim2$ in the COSMOS legacy field using near-infrared (NIR) spectroscopy from Keck/MOSFIRE as part of the ZFIRE survey. Our sample consists of 75 Ks-band selected star-forming galaxies from the ZFOURGE survey with stellar masses ranging from log(M$_{\star}$/M$_{\odot}$)$=9.0-11.0$, 28 of which are members of a known overdensity at $z=2.095$. We measure H$\alpha$ emission-line integrated velocity dispersions ($\sigma_{\rm int}$) from 50$-$230 km s$^{-1}$, consistent with other emission-line studies of $z\sim2$ field galaxies. From these data we estimate virial, stellar, and gas masses and derive correlations between these properties for cluster and field galaxies at $z\sim2$. We find evidence that baryons dominate within the central effective radius. However, we find no statistically significant differences between the cluster and the field, and conclude that the kinematics of star-forming galaxies at $z\sim2$ are not significantly different between the cluster and field environments.

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

We investigate black hole-host galaxy scaling relations in cosmological simulations with a self-consistent black hole growth and feedback model. The sub-grid accretion model captures the key scalings governing angular momentum transport from galactic scales down to parsec scales, while our kinetic feedback implementation enables the injection of outflows with properties chosen to match observed nuclear outflows. We show that "quasar mode" feedback can have a large impact on the thermal properties of the intergalactic medium and the growth of galaxies and massive black holes for kinetic feedback efficiencies as low as 0.1% relative to the bolometric luminosity. Nonetheless, our simulations suggest that the black hole-host scaling relations are only weakly dependent on the effects of black hole feedback on galactic scales, owing to feedback suppressing the growth of galaxies and massive black holes by a similar amount. In contrast, the rate at which gravitational torques feed the central black hole relative to the host galaxy star formation rate governs the slope and normalization of the black hole-host correlations. Our results suggest that a common gas supply regulated by gravitational torques is the primary driver of the observed co-evolution of black holes and galaxies.

The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): Bright quasars (QSOs) generate a hard UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller-scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensitive to different energies. Likewise, they probe spatial fluctuations because they originate in regions where a galaxy-driven UVB is harder and more intense. We take a first step towards studying the reionization-epoch UVB's slope and inhomogeneity by comparing observations of 12 metal absorbers at $z\sim6$ versus predictions from a cosmological hydrodynamic simulation using three different UVBs: a soft, spatially-inhomogeneous "galaxies+QSOs" UVB; a homogeneous "galaxies+QSOs" UVB (Haardt & Madau 2012); and a QSOs-only model. All UVBs reproduce the observed column density distributions of CII, SiIV, and CIV reasonably well although high-column, high-ionization absorbers are underproduced, reflecting numerical limitations. With upper limits treated as detections, only a soft, fluctuating UVB reproduces both the observed SiIV/CIV and CII/CIV distributions. The QSOs-only UVB overpredicts both CIV/CII and CIV/SiIV, indicating that it is too hard. The Haardt & Madau (2012) UVB underpredicts CIV/SiIV, suggesting that it lacks amplifications near galaxies. Hence current observations prefer a soft, fluctuating UVB as expected from a predominantly Population II background although they cannot rule out a harder one. Future observations probing a factor of two deeper in metal column density will distinguish between the soft, fluctuating and QSOs-only UVBs.

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

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

We explore the impact of incorporating physically motivated ionisation and recombination rates on the history and topology of cosmic reionisation, by incorporating inputs from small-volume hydrodynamic simulations into a semi-numerical code, SimFast21, that evolves reionisation on large scales. We employ radiative hydrodynamic simulations to parameterize the ionisation rate Rion and recombination rate Rrec as functions of halo mass, overdensity and redshift. We find that Rion is super-linearly dependent on halo mass (Rion ~ Mh^1.41), in contrast to previous assumptions. We implement these scalings into SimFast21 to identify the ionized regions. We tune our models to be consistent with recent observations of the optical depth, ionizing emissivity, and neutral fraction by the end of reionisation. We require an average photon escape fraction fesc=0.04 within ~ 0.5 cMpc cells, independent of halo mass or redshift, to simultaneously match these data. We present predictions for the 21cm power spectrum, and show that it is converged with respect to simulation volume. We find that introducing superlinearly mass-dependent ionisations increases the duration of reionisation and boosts the small-scale 21cm power by ~ 2-3 at intermediate phases of reionisation. Introducing inhomogeneous recombinations reduces ionised bubble sizes and suppresses large-scale 21cm power by ~ 2-3. Moreover, gas clumping on sub-cell scales has a minimal effect on the 21cm power, indicating that robust predictions do not depend on the behaviour of kpc-scale structures. The superlinear ionisations significantly increase the median halo mass scale for ionising photon output to >10^10 Mo, giving greater hope for detecting most of ionising sources with next-generation facilities. These results highlight the importance of more accurately treating ionising sources and recombinations for modeling reionisation and its 21cm signal.

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

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

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

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

The Bullet Cluster has provided some of the best evidence for the $\Lambda$ cold dark matter ($\Lambda\rm{CDM}$) model via direct empirical proof of the existence of collisionless dark matter, while posing a serious challenge owing to the unusually high inferred pairwise velocities of its progenitor clusters. Here we investigate the probability of finding such a high-velocity pair in large-volume N-body simulations, particularly focusing on differences between halo finding algorithms. We find that algorithms that do not account for the kinematics of infalling groups yield vastly different statistics and probabilities. When employing the ROCKSTAR halo finder that considers particle velocities, we find numerous Bullet-like pair candidates that closely match not only the high pairwise velocity, but also the mass, mass ratio, separation distance, and collision angle of the initial conditions that have been shown to produce the Bullet Cluster in non-cosmological hydrodynamic simulations. The probability of finding a high pairwise velocity pair among haloes with $M_{\rm halo}\geq10^{14} M_{\odot}$ is $4.6\times 10^{-4}$ using ROCKSTAR, while it is $\approx 34\times$ lower using a friends-of-friends (FOF) based approach as in previous studies. This is because the typical spatial extent of Bullet progenitors is such that FOF tends to group them into a single halo despite clearly distinct kinematics. Further requiring an appropriately high average mass among the two progenitors, we find the comoving number density of potential Bullet-like candidates to be on the order of $\approx10^{-10} \,{\rm Mpc}^{-3}$. Our findings suggest that $\Lambda\rm{CDM}$ straightforwardly produces massive, high relative velocity halo pairs analogous to Bullet Cluster progenitors, and hence the Bullet Cluster does not present a challenge to the $\Lambda\rm{CDM}$ model.

Galaxies strongly self-regulate their growth via energetic feedback from stars, supernovae, and black holes, but these processes are among the least understood aspects of galaxy formation theory. We present an analytic galaxy evolution model that directly constrains such feedback processes from observed galaxy scaling relations. The equilibrium model, which is broadly valid for star-forming central galaxies that dominate cosmic star formation, is based on the ansatz that galaxies live in a slowly-evolving equilibrium between inflows, outflows, and star formation. Using a Bayesian Monte Carlo Markov chain approach, we constrain our model to match observed galaxy scaling relations between stellar mass and halo mass, star formation rate, and metallicity from 0<z<2. A good fit (chi^2~1.6) is achieved with eight free parameters. We further show that constraining our model to any two of the three data sets also produces a fit to the third that is within reasonable systematic uncertainties. The resulting best-fit parameters that describe baryon cycling suggest galactic outflow scalings intermediate between energy and momentum-driven winds, a weak dependence of wind recycling time on mass, and a quenching mass scale that evolves modestly upwards with redshift. This model further predicts a stellar mass-star formation rate relation that is in good agreement with observations to z~6. Our results suggest that this simple analytic framework captures the basic physical processes required to model the mean evolution of stars and metals in galaxies, despite not incorporating many canonical ingredients of galaxy formation models such as merging or disk formation.

Recent studies have shown that the local mass-metallicity (M-Z) relation depends on the specific star formation rate (SSFR). Whether such a dependence exists at higher redshifts, and whether the resulting M-Z-SFR relation is redshift invariant, is debated. We re-examine these issues by applying the non-parametric techniques of Salim et al. (2014) to ~130 $z\sim2.3$ galaxies with N2 and O3 measurements from KBSS (Steidel et al. 2014). We find that the KBSS M-Z relation depends on SSFR at intermediate masses, where such dependence exists locally. KBSS and SDSS galaxies of the same mass and SSFR ("local analogs") are similarly offset in the BPT diagram relative to the bulk of local star-forming galaxies, and thus we posit that metallicities can be compared self-consistently at different redshifts as long as the masses and SSFRs of the galaxies are similar. We find that the M-Z-SFR relation of $z\sim2$ galaxies is consistent with the local one at $\log M_*<10$, but is offset up to -0.25 dex at higher masses, so it is altogether not redshift invariant. This high-mass offset could arise from a bias that high-redshift spectroscopic surveys have against high-metallicity galaxies, but additional evidence disfavors this possibility. We identify three causes for the reported discrepancy between N2 and O3N2 metallicities at $z\sim2$: (1) a smaller offset that is also present for SDSS galaxies, which we remove with new N2 calibration, (2) a genuine offset due to differing ISM condition, which is also present in local analogs, (3) an additional offset due to unrecognized AGN contamination.

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

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

We investigate the evolution of stellar population gradients from $z=2$ to $z=0$ in massive galaxies at large radii ($r > 2R_{\mathrm{eff}}$) using ten cosmological zoom simulations of halos with $6 \times 10^{12} M_{\odot} < M_{\mathrm{halo}} < 2 \times 10^{13}M_{\odot}$. The simulations follow metal cooling and enrichment from SNII, SNIa and AGB winds. We explore the differential impact of an empirical model for galactic winds that reproduces the mass-metallicity relation and its evolution with redshift. At larger radii the galaxies, for both models, become more dominated by stars accreted from satellite galaxies in major and minor mergers. In the wind model, fewer stars are accreted, but they are significantly more metal poor resulting in steep global metallicity ($\langle \nabla Z_{\mathrm{stars}} \rangle= -0.35$ dex/dex) and color (e.g. $\langle \nabla g-r \rangle = -0.13$ dex/dex) gradients in agreement with observations. In contrast, colour and metallicity gradients of the models without winds are inconsistent with observations. Age gradients are in general mildly positive at $z=0$ ($\langle \nabla Age_{\mathrm{stars}} \rangle= 0.04$ dex/dex) with significant differences between the models at higher redshift. We demonstrate that for the wind model, stellar accretion is steepening existing in-situ metallicity gradients by about 0.2 dex by the present day and helps to match observed gradients of massive early-type galaxies at large radii. Colour and metallicity gradients are significantly steeper for systems which have accreted stars in minor mergers, while galaxies with major mergers have relatively flat gradients, confirming previous results. This study highlights the importance of stellar accretion for stellar population properties of massive galaxies at large radii, which can provide important constraints for formation models.

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $\sigma(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within 2-4 times that of Haardt & Madau (2012), which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman-alpha forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially-uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when HII regions overlap globally even at column densites of 10^{12} cm^{-2}. We find, consistent with observations, that the CII mass fraction drops to low redshift while CIV rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z=6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially-neutral Universe at z=6.

It has been proposed that the mass-metallicity relation of galaxies exhibits a secondary dependence on star formation rate (SFR), and that the resulting M-Z-SFR relation may be redshift-invariant, i.e., "fundamental." However, conflicting results on the character of the SFR dependence, and whether it exists, have been reported. To gain insight into the origins of the conflicting results, we (a) devise a non-parametric, astrophysically motivated analysis framework based on the offset from the star-forming ("main") sequence at a given stellar mass (relative specific SFR), (b) apply this methodology and perform a comprehensive re-analysis of the local M-Z-SFR relation, based on SDSS, GALEX, and WISE data, and (c) study the impact of sample selection, and of using different metallicity and SFR indicators. We show that metallicity is anti-correlated with specific SFR regardless of the indicators used. We do not find that the relation is spurious due to correlations arising from biased metallicity measurements, or fiber aperture effects. We emphasize that the dependence is weak/absent for massive galaxies ($\log M_*>10.5$), and that the overall scatter in the M-Z-SFR relation does not greatly decrease from the M-Z relation. We find that the dependence is stronger for the highest SSFR galaxies above the star-forming sequence. This two-mode behavior can be described with a broken linear fit in 12+log(O/H) vs. log (SFR$/M_*$), at a given $M_*$. Previous parameterizations used for comparative analysis with higher redshift samples that do not account for the more detailed behavior of the local M-Z-SFR relation may incorrectly lead to the conclusion that those samples follow a different relationship.

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

Observations indicate that galaxies with high stellar masses or in dense environments have low specific star formation rates, i.e. they are quenched. Based on cosmological hydrodynamic simulations that include a prescription where quenching occurs in regions dominated by hot (>10^5 K) gas, we argue that this hot gas quenching in halos >10^12 Msun drives both mass quenching (i.e. central quenching) and environment quenching (i.e. satellite quenching). These simulations reproduce a broad range of locally observed trends among quenching, halo mass, stellar mass, environment, and distance to halo center. We show that mass quenching is independent of environment because 10^12-10^13 Msun "quenching halos" -- those where most mass quenching occurs -- inhabit a large range of environments. On the other hand, environment quenching is independent of stellar mass because galaxies of all stellar masses may live in dense environments as satellites of groups and clusters. Furthermore, satellite galaxies show signs of mass quenching independent of halo mass because massive satellites at z=0 have typically been mass quenched as centrals in their own hot halos at higher z -- a kind of pre-processing. As in observations, the fraction of quenched satellites increases with halo mass and decreases with distance to the center of the group or cluster. We investigate quenched centrals in low-mass halos (<10^12 Msun), and show that most of these are ejected former satellites of groups or clusters, while about 20 per cent were never satellites but are enveloped in hot gas that extends up to 3 Rvir from the centers of clusters. The agreement of our model with key observational trends suggests that hot gas in massive halos plays a leading role in quenching low-redshift galaxies.

Modeling galaxy formation in a cosmological context presents one of the greatest challenges in astrophysics today, due to the vast range of scales and numerous physical processes involved. Here we review the current status of models that employ two leading techniques to simulate the physics of galaxy formation: semi-analytic models and numerical hydrodynamic simulations. We focus on a set of observational targets that describe the evolution of the global and structural properties of galaxies from roughly Cosmic High Noon ($z\sim 2-3$) to the present. Although minor discrepancies remain, overall, models show remarkable convergence between different methods and make predictions that are in qualitative agreement with observations. Modelers seem to have converged on a core set of physical processes that are critical for shaping galaxy properties. This core set includes cosmological accretion, strong stellar-driven winds that are more efficient at low masses, black hole feedback that preferentially suppresses star formation at high masses, and structural and morphological evolution through merging and environmental processes. However, all cosmological models currently adopt phenomenological implementations of many of these core processes, which must be tuned to observations. Many details of how these diverse processes interact within a hierarchical structure formation setting remain poorly understood. Emerging multi-scale simulations are helping to bridge the gap between stellar and cosmological scales, placing models on a firmer, more physically grounded footing. Concurrently, upcoming telescope facilities will provide new challenges and constraints for models, particularly by directly constraining inflows and outflows through observations of gas in and around galaxies.

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

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($\alpha \approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} \rho_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($\beta$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

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

[Abridged] With VLT/X-shooter, we obtain optical and NIR spectra of six Ly-alpha blobs at z~2.3. Using three measures --- the velocity offset between the Lya line and the non-resonant [OIII] or H-alpha line (Dv_Lya), the offset of stacked interstellar metal absorption lines, and the spectrally-resolved [OIII] line profile --- we study the kinematics of gas along the line of sight to galaxies within each blob center. These three indicators generally agree in velocity and direction, and are consistent with a simple picture in which the gas is stationary or slowly outflowing at a few hundred km/s from the embedded galaxies. The absence of stronger outflows is not a projection effect: the covering fraction for our sample is limited to <1/8 (13%). The outflow velocities exclude models in which star formation or AGN produce "super" or "hyper" winds of up to ~1000km/s. The Dv_Lya offsets here are smaller than typical of LBGs, but similar to those of compact LAEs. The latter suggests that outflow speed cannot be a dominant factor in driving extended Lya emission. For one Lya blob (CDFS-LAB14), whose Lya profile and metal absorption line offsets suggest no significant bulk motion, we use a simple radiative transfer model to make the first column density measurement of gas in an embedded galaxy, finding it consistent with a DLA system. Overall, the absence of clear inflow signatures suggests that the channeling of gravitational cooling radiation into Lya is not significant over the radii probed here. However, one peculiar system (CDFS-LAB10) has a blueshifted Lya component that is not obviously associated with any galaxy, suggesting either displaced gas arising from tidal interactions among blob galaxies or gas flowing into the blob center. The former is expected in these overdense regions, and the latter might signify the predicted but elusive cold gas accretion along filaments.

We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$\sigma$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.

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

We use high-resolution cosmological zoom simulations with ~200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ~ 10^11--2x10^12 Msun at z = 2. We present a detailed comparison to spatially and spectrally resolved H{\alpha} and other observations of z ~ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics - among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in H{\alpha} but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth.

Recent models of black hole growth in a cosmological context have forwarded a paradigm in which the growth is self-regulated by feedback from the black hole itself. Here we use cosmological zoom simulations of galaxy formation down to z = 2 to show that such strong self-regulation is required in the popular spherical Bondi accretion model, but that a plausible alternative model in which black hole growth is limited by galaxy-scale torques does not require self-regulation. Instead, this torque-limited accretion model yields black holes and galaxies evolving on average along the observed scaling relations by relying only on a fixed, 5% mass retention rate onto the black hole from the radius at which the accretion flow is fed. Feedback from the black hole may (and likely does) occur, but does not need to couple to galaxy-scale gas in order to regulate black hole growth. We show that this result is insensitive to variations in the initial black hole mass, stellar feedback, or other implementation details. The torque-limited model allows for high accretion rates at very early epochs (unlike the Bondi case), which if viable can help explain the rapid early growth of black holes, while by z = 2 it yields Eddington factors of 1%-10%. This model also yields a less direct correspondence between major merger events and rapid phases of black hole growth. Instead, growth is more closely tied to cosmological disk feeding, which may help explain observational studies showing that, at least at z > 1, active galaxies do not preferentially show merger signatures.

We examine the cosmic evolution of a stellar initial mass function (IMF) in galaxies that varies with the Jeans mass in the interstellar medium, paying particular attention to the K-band stellar mass to light ratio (M/L_K) of present-epoch massive galaxies. We calculate the typical Jeans mass using high-resolution hydrodynamic simulations coupled with a fully radiative model for the ISM, which yields a parameterisation of the IMF characteristic mass as a function of galaxy star formation rate (SFR). We then calculate the star formation histories of galaxies utilising an equilibrium galaxy growth model coupled with constraints on the star formation histories set by abundance matching models. We find that at early times, energetic coupling between dust and gas drive warm conditions in the ISM, yielding bottom-light/top- heavy IMFs associated with large ISM Jeans masses for massive star-forming galaxies. Owing to the remnants of massive stars that formed during the top-heavy phases at early times, the resultant M/L_K(sigma) in massive galaxies at the present epoch is increased relative to the non- varying IMF case. At late times, lower cosmic ray fluxes allow for cooler ISM temperatures in massive galaxies, and hence newly formed clusters will exhibit bottom-heavy IMFs, further increasing M/L_K(sigma). Our central result is hence that a given massive galaxy may go through both top-heavy and bottom-heavy IMF phases during its lifetime, though the bulk of the stars form during a top-heavy phase. Qualitatively, the variations in M/L_K(sigma) with galaxy mass are in agreement with observations, however, our model may not be able to account for bottom-heavy mass functions as indicated by stellar absorption features.

Studies of QSO absorber-galaxy connections are often hindered by inadequate information on whether faint/dwarf galaxies are located near the QSO sight lines. To investigate the contribution of faint galaxies to QSO absorber populations, we are conducting a deep galaxy redshift survey near low-z C IV absorbers. Here we report a blindly-detected C IV absorption system (z(abs) = 0.00348) in the spectrum of PG1148+549 that appears to be associated either with an edge-on dwarf galaxy with an obvious disk (UGC 6894, z(gal) = 0.00283) at an impact parameter of rho = 190 kpc or with a very faint dwarf irregular galaxy at rho = 23 kpc, which is closer to the sightline but has a larger redshift difference (z(gal) = 0.00107, i.e., dv = 724 km/s). We consider various gas/galaxy associations, including infall and outflows. Based on current theoretical models, we conclude that the absorber is most likely tracing (1) the remnants of an outflow from a previous epoch, a so-called 'ancient outflow', or (2) intergalactic gas accreting onto UGC 6894, 'cold mode' accretion. The latter scenario is supported by H I synthesis imaging data that shows the rotation curve of the disk being codirectional with the velocity offset between UGC 6894 and the absorber, which is located almost directly along the major axis of the edge-on disk.

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

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

We use a radiation hydrodynamic simulation of the hydrogen reionization epoch to study OI absorbers at z~6. The intergalactic medium (IGM) is reionized before it is enriched, hence OI absorption originates within dark matter halos. The predicted abundance of OI absorbers is in reasonable agreement with observations. At z=10, roughly 70% of sightlines through atomically-cooled halos encounter a visible (N_OI > 10^14 cm^-2) column. Reionization ionizes and removes gas from halos less massive than 10^8.4 M_0, but 20% of sightlines through more massive halos encounter visible columns even at z=5. The mass scale of absorber host halos is 10-100 times smaller than the halos of Lyman break galaxies and Lyman-alpha emitters, hence absorption probes the dominant ionizing sources more directly. OI absorbers have neutral hydrogen columns of 10^19-10^21 cm^-2, suggesting a close resemblance between objects selected in OI and HI absorption. Finally, the absorption in the foreground of the z=7.085 quasar ULASJ1120+0641 cannot originate in a dark matter halo because halo gas at the observed HI column density is enriched enough to violate the upper limits on the OI column. By contrast, gas at less than one third the cosmic mean density satisfies the constraints. Hence the foreground absorption likely originates in the IGM.

We present evidence for a skewed distribution of UV FeII emission in quasars within candidate overdense regions spanning spatial scales of ~ 50 Mpc at 1.11 < z < 1.67, compared to quasars in field environments at comparable redshifts. The overdense regions have an excess of high equivalent width sources (W2400 > 42 \AA), and a dearth of low equivalent width sources. There are various possible explanations for this effect, including dust, Ly\alpha fluorescence, microturbulence, and iron abundance. We find that the most plausible of these is enhanced iron abundance in the overdense regions, consistent with an enhanced star formation rate in the overdense regions compared to the field.

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

We have updated our radially-resolved SAMs of galaxy formation, which track both the atomic and molecular gas phases of the ISM. The models are adapted from those of Guo et al. using similar methodology as in Fu et al. and are run on halo merger trees from the MS and MS II with the following main changes: (1) We adopt a simple star formation law where \Sigma_SFR \propto \Sigma_H2 (2) We inject the heavy elements produced by supernovae directly into the halo hot gas, instead of first mixing them with the cold gas in the disk. (3) We include radial gas inflows in disks using a model of the form v_inflow = \alpha r. The models are used to study the radial profiles of star formation rate and gas-phase metallicity in present-day galaxies. The \Sigma_H2 profiles in L* galaxies place strong constraints on inflow velocities, favouring models where v_inflow~7km/s at a galactocentric radius of 10kpc. Radial gas inflow has little influence on gas-phase and stellar metallicity gradients, which are affected much more strongly by the fraction of metals that are directly injected into the halo gas, rather than mixed with the cold gas. Metals ejected out of the galaxy in early epochs result in late infall of pre-enriched gas and flatter present-day gas-phase metallicity gradients. A prescription in which 80% of the metals are injected into the halo gas results in good fits to the flat observed metallicity gradients in galaxies with stellar masses greater than 10^10 M_sun, as well as the relations between gas-phase metallicity and sSFR in the outer parts of galactic disks. We examine the correlation between gas-phase metallicity gradient and global galaxy properties, finding that it is most strongly correlated with the bulge-to-total ratio of the galaxy. This is because gas is consumed when the bulge forms during galaxy mergers, and the gas-phase metallicity gradient is then set by newly-accreted gas.

We use long-slit spectroscopy from Moran et al. to study the radial dependence of the recent star formation histories of nearby galaxies with stellar masses greater than 10^10M_sun. We fit stellar population models to the combination of SSFR, D4000 and Hdelta_A and show that many galaxies have Balmer absorption line equivalent widths that require recent short-lived episodes or bursts of star formation. The fraction of galaxies that have experienced episodic rather than continuous star formation is highest for late-type galaxies with low stellar masses. In these systems, bursts occur both in the inner and outer regions of the galaxy. The fraction of stars formed in a single burst episode is typically around 15% of the total stellar mass in the inner regions of the galaxy and around 5% of the mass in the outer regions. When we average over the population, we find that such bursts contribute around a half of the total mass in stars formed in the last 2 Gyr. In massive galaxies, bursts occur predominantly in the outer disk. Around a third of all massive, bulge-dominated galaxies have experienced recent star formation episodes that are fully confined to their outer (R > 0.7R_90) regions. The fraction of stars formed in a single episode is only 2 - 3 % of the underlying stellar mass, but such bursts contribute nearly all the stellar mass formed in the last 2 Gyr. Recent star formation in outer disks is strongly correlated with the global atomic gas fraction of the galaxy, but not its global molecular gas fraction. We suggest that outer episodic star formation is triggered by gas accretion events.

We analyze the galactic H I content and nebular log(O/H) for 60 spiral galaxies in the Moustakas et al. (2006) spectral catalog. After correcting for the mass-metallicity relationship, we show that the spirals in cluster environments show a positive correlation for log(O/H) on DEF, the galactic H I deficiency parameter, extending the results of previous analyses of the Virgo and Pegasus I clusters. Additionally, we show for the first time that galaxies in the field obey a similar dependence. The observed relationship between H I deficiency and galactic metallicity resembles similar trends shown by cosmological simulations of galaxy formation including inflows and outflows. These results indicate the previously observed metallicity-DEF correlation has a more universal interpretation than simply a cluster's effects on its member galaxies. Rather, we observe in all environments the stochastic effects of metal-poor infall as minor mergers and accretion help to build giant spirals.

We have exploited the HST CANDELS WFC3/IR imaging to study the properties of (sub-)mm galaxies in GOODS-South. After using the deep radio and Spitzer imaging to identify galaxy counterparts for the (sub-)mm sources, we have used the new CANDELS data in two ways. First, we have derived improved photometric redshifts and stellar masses, confirming that the (sub-)mm galaxies are massive (<M*>=2.2x10^11 M_solar) galaxies at z=1-3. Second, we have exploited the depth and resolution of the WFC3/IR imaging to determine the sizes and morphologies of the galaxies at rest-frame optical wavelengths, fitting two-dimensional axi-symmetric Sersic models. Crucially, the WFC3/IR H-band imaging enables modelling of the mass-dominant galaxy, rather than the blue high-surface brightness features which often dominate optical (rest-frame UV) images of (sub-)mm galaxies, and can confuse visual morphological classification. As a result of this analysis we find that >95% of the rest-frame optical light in almost all of the (sub-)mm galaxies is well-described by either a single exponential disk, or a multiple-component system in which the dominant constituent is disk-like. We demonstrate that this conclusion is consistent with the results of high-quality ground-based K-band imaging, and explain why. The massive disk galaxies which host luminous (sub-)mm emission are reasonably extended (r_e=4 kpc), consistent with the sizes of other massive star-forming disks at z~2. In many cases we find evidence of blue clumps within the sources, with the mass-dominant disk becoming more significant at longer wavelengths. Finally, only a minority of the sources show evidence for a major galaxy-galaxy interaction. Taken together, these results support the view that most (sub-)mm galaxies at z~2 are simply the most extreme examples of normal star-forming galaxies at that era.

Measuring rest-frame ultraviolet rotational transitions from the Lyman and Werner bands in absorption against a bright background continuum is one of the few ways to directly measure molecular hydrogen (H2). Here we report the detection of Lyman-Werner absorption from H2 at z=0.56 in a sub-damped Ly-alpha system with neutral hydrogen column density N(HI) = 10^(19.5 +/- 0.2) cm^-2. This is the first H2 system analysed at a redshift < 1.5 beyond the Milky Way halo. It has a surprisingly high molecular fraction: log f(H2) > -1.93 +/- 0.36 based on modelling the line profiles, with a robust model-independent lower limit of f(H2) > 10^-3. This is higher than f(H2) values seen along sightlines with similar N(HI) through the Milky Way disk and the Magellanic clouds. The metallicity of the absorber is 0.19 +0.21 -0.10 solar, with a dust-to-gas ratio < 0.36 times the value in the solar neighbourhood. Absorption from associated low-ionisation metal transitions such as OI and FeII is observed in addition to OVI. Using Cloudy models we show that there are three phases present; a ~100 K phase giving rise to H2, a ~10^4 K phase where most of the low-ionisation metal absorption is produced; and a hotter phase associated with OVI. Based on similarities to high velocity clouds in the Milky Way halo showing H2 and the presence of two nearby galaxy candidates with impact parameters of ~10 kpc, we suggest that the absorber may be produced by a tidally-stripped structure similar to the Magellanic Stream.

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

We use a suite of cosmological hydrodynamic simulations including a self-consistent treatment for inhomogeneous reionisation to study the impact of galactic outflows and photoionisation heating on the volume-averaged recombination rate of the intergalactic medium (IGM). By incorporating an evolving ionising escape fraction and a treatment for self-shielding within Lyman limit systems, we have run the first simulations of "photon-starved" reionisation scenarios that simultaneously reproduce observations of the abundance of galaxies, the optical depth to electron scattering of cosmic microwave background photons \tau, and the effective optical depth to Lyman\alpha absorption at z=5. We confirm that an ionising background reduces the clumping factor C by more than 50% by smoothing moderately-overdense (\Delta=1--100) regions. Meanwhile, outflows increase clumping only modestly. The clumping factor of ionised gas is much lower than the overall baryonic clumping factor because the most overdense gas is self-shielded. Photoionisation heating further suppresses recombinations if reionisation heats gas above the canonical 10,000 K. Accounting for both effects within our most realistic simulation, C rises from <1 at z>10 to 3.3 at z=6. We show that incorporating temperature- and ionisation-corrected clumping factors into an analytical reionisation model reproduces the numerical simulation's \tau to within 10%. Finally, we explore how many ionising photons are absorbed during the process of heating filaments by considering the overall photon cost of reionisation in analytical models that assume that the IGM is heated at different redshifts. For reionisation redshifts of 9--10, cold filaments boost the reionisation photon budget by ~1 photon per hydrogen atom.

We examine the cosmic growth of the red sequence in a cosmological hydrodynamic simulation that includes a heuristic prescription for quenching star formation that yields a realistic passive galaxy population today. In this prescription, halos dominated by hot gas are continually heated to prevent their coronae from fueling new star formation. Hot coronae primarily form in halos above \sim10^12 M\odot, so that galaxies with stellar masses \sim10^10.5 M\odot are the first to be quenched and move onto the red sequence at z > 2. The red sequence is concurrently populated at low masses by satellite galaxies in large halos that are starved of new fuel, resulting in a dip in passive galaxy number densities around \sim10^10 M\odot. Stellar mass growth continues for galaxies even after joining the red sequence, primarily through minor mergers with a typical mass ratio \sim1:5. For the most massive systems, the size growth implied by the distribution of merger mass ratios is typically \sim2\times the corresponding mass growth, consistent with observations. This model reproduces mass-density and colour-density trends in the local universe, with essentially no evolution to z = 1, with the hint that such relations may be washed out by z \sim 2. Simulated galaxies are increasingly likely to be red at high masses or high local overdensities. In our model, the presence of surrounding hot gas drives the trends with both mass and environment.

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

Observations of star-forming galaxies at high-z have suggested discrepancies in the inferred star formation rates (SFRs) either between data and models, or between complementary measures of the SFR. These putative discrepancies could all be alleviated if the stellar IMF is systematically weighted toward more high-mass star formation in rapidly star-forming galaxies. Here, we explore how the IMF might vary under the central assumption that the turnover mass in the IMF, Mc, scales with the Jeans mass in giant molecular clouds (GMCs), M_J. We employ hydrodynamic and radiative transfer simulations of galaxies to predict how the typical GMC Jeans mass, and hence the IMF, varies with galaxy property. We then study the impact of such an IMF on the star formation law, the SFR-M* relation, submillimetre galaxies (SMGs), and the cosmic SFR density. Our main results are: The H2 mass-weighted Jeans mass in a galaxy scales with the SFR when the SFR is greater a few M_sun/yr. SPS modeling shows that this results in a nonlinear relation between SFR and Lbol, such that SFR Lbol^0.88. Using this model relation, the inferred SFR of local ULIRGs decreases by ~2, and that of high-z SMGs decreases by ~3-5. At z 2, this results in a lowered normalisation of the SFR-M* relation in better agreement with models, a reduced discrepancy between the observed cosmic SFR density and stellar mass density evolution, and SMG SFRs that are easier to accommodate in current hierarchical structure formation models. It further results in a Schmidt relation with slope of ~1.6 when utilising a physically motivated form for the CO-H2 conversion factor. While each of the discrepancies considered here could be alleviated without appealing to a varying IMF, the modest variation implied by assuming Mc M_J is a plausible solution that simultaneously addresses numerous thorny issues regarding the SFRs of high-z galaxies.

In this investigation we quantify the metallicities of low mass galaxies by constructing the most comprehensive census to date. We use galaxies from the SDSS and DEEP2 survey and estimate metallicities from their optical emission lines. We also use two smaller samples from the literature which have metallicities determined by the direct method using the temperature sensitive [OIII]4363 line. We examine the scatter in the local mass-metallicity (MZ) relation determined from ~20,000 star-forming galaxies in the SDSS and show that it is larger at lower stellar masses, consistent with the theoretical scatter in the MZ relation determined from hydrodynamical simulations. We determine a lower limit for the scatter in metallicities of galaxies down to stellar masses of ~10^7 M_solar that is only slightly smaller than the expected scatter inferred from the SDSS MZ relation and significantly larger than what is previously established in the literature. The average metallicity of star-forming galaxies increases with stellar mass. By examining the scatter in the SDSS MZ relation, we show that this is mostly due to the lowest metallicity galaxies. The population of low mass, metal-rich galaxies have properties which are consistent with previously identified galaxies that may be transitional objects between gas-rich dwarf irregulars and gas-poor dwarf spheroidals and ellipticals.

We have obtained long-slit spectra of 174 star-forming galaxies with stellar masses greater than 10^10 M_\odot from the GALEX Arecibo SDSS (GASS) survey. These galaxies have both HI and H_2 mass measurements. The average metallicity profile is strikingly flat out to R_90, the radius enclosing 90% of the r-band light. Metallicity profiles which decline steadily with radius are found primarily for galaxies in our sample with low stellar mass (Log(M_*)<10.2), concentration, and/or mean stellar mass density. Beyond ~R_90, however, around 10 percent of the galaxies in our sample exhibit a sharp downturn in metallicity. Remarkably, we find that the magnitude of the outer metallicity drop is well correlated with the total HI content of the galaxy (measured as f_HI=M_HI/M_*). We examine the radial profiles of stellar population ages and star formation rate densities, and conclude that the galaxies with largest outer metallicity drops are actively growing their stellar disks, with mass doubling times across the whole disk only one third as long as a typical GASS galaxy. We also describe a correlation between local stellar mass density and metallicity, which is valid across all galaxies in our sample. We argue that much of the recent stellar mass growth at the edges of these galaxies can be linked to the accretion or radial transport of relatively pristine gas from beyond the galaxies' stellar disks.

We present an analytic formalism that describes the evolution of the stellar, gas, and metal content of galaxies. It is based on the idea, inspired by hydrodynamic simulations, that galaxies live in a slowly-evolving equilibrium between inflow, outflow, and star formation. We argue that this formalism broadly captures the behavior of galaxy properties evolving in simulations. The resulting equilibrium equations for the star formation rate, gas fraction, and metallicity depend on three key free parameters that represent ejective feedback, preventive feedback, and re-accretion of ejected material. We schematically describe how these parameters are constrained by models and observations. Galaxies perturbed off the equilibrium relations owing to inflow stochasticity tend to be driven back towards equilibrium, such that deviations in star formation rate at a given mass are correlated with gas fraction and anti-correlated with metallicity. After an early gas accumulation epoch, quiescently star-forming galaxies are expected to be in equilibrium over most of cosmic time. The equilibrium model provides a simple intuitive framework for understanding the cosmic evolution of galaxy properties, and centrally features the cycle of baryons between galaxies and surrounding gas as the driver of galaxy growth.

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

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We carry out a new suite of cosmological radiation hydrodynamic simulations and explore the relative impacts on reionization-epoch star formation of galactic outflows and photoionization heating. By itself, an extragalactic ultraviolet background (EUVB) suppresses the luminosity function by less than 50% at z=6, overproducing the observed galaxy abundance by a factor of 3-5. Galactic outflows restore agreement with observations without preventing Population II star formation from reionizing the Universe by z=6. The resulting EUVB suppresses star formation in halos with virial temperatures below 10^5K but has a weaker impact in more massive halos. Nonetheless, the low-mass halos contribute up to 50% of all ionizing photons owing to the EUVB's inhomogeneity. Overall, star formation rate scales as halo mass M_h to the 1.3-1.4 in halos with $M_h=10^{8.2--10.2}\msun$. This is a steeper dependence than is often assumed in reionization models, boosting the expected power spectrum of 21 centimeter fluctuations on large scales. The luminosity function rises steeply to at least M_1600=-13, indicating that reionization was driven by faint galaxies (M_1600 >= -15) that have not yet been observed. Our models cannot simultaneously explain observations of galaxies, the cosmic microwave background, and the intergalactic medium. Increased dynamic range will alleviate the existing discrepancies, but observations may still require additional physics such as a variable ionizing escape fraction (abridged).

Many observed massive star-forming z\approx2 galaxies are large disks that exhibit irregular morphologies, with \sim1kpc, \sim10^(8-10)Msun clumps. We present the largest sample to date of high-resolution cosmological SPH simulations that zoom-in on the formation of individual M*\sim10^(10.5)Msun galaxies in \sim10^(12)Msun halos at z\approx2. Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model reproduces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, smooth and monotonic velocity gradients, high gas fractions (f_g\sim50%) and high specific star-formation rates (\gtrsim1Gyr^(-1)). In accord with recent models, giant clumps (Mclump\sim(5x10^8-10^9)Msun) form in-situ via gravitational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations reproduce well new observational constraints on clump kinematics and clump ages. In particular, the observation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of inside-out formation of the disks rather than inward clump migration.

We use cosmological hydrodynamic simulations to investigate how inflows, star formation, and outflows govern the the gaseous and metal content of galaxies. In our simulations, galaxy metallicities are established by a balance between inflows and outflows as governed by the mass outflow rate, implying that the mass-metallicity relation reflects how the outflow rate varies with stellar mass (M*). Gas content is set by a competition between inflow into and gas consumption within the ISM, the latter being governed by the SF law, while the former is impacted by both wind recycling and preventive feedback. Stochastic variations in the inflow rate move galaxies off the equilibrium M*-Z and Z*-fgas relations in a manner correlated with star formation rate, and the scatter is set by the timescale to re-equilibrate. The evolution of both relations from z=3-0 is slow, as individual galaxies tend to evolve mostly along the relations. Gas fractions at a given M* slowly decrease with time because the cosmic inflow rate diminishes faster than the consumption rate, while metallicities slowly increase as infalling gas becomes more enriched. Observations from z~3-0 are better matched by simulations employing momentum-driven wind scalings rather than constant wind speeds, but all models predict too low gas fractions at low masses and too high metallicities at high M*. All our models reproduce observed second-parameter trends of the mass-metallicity relation with star formation rate and environment, indicating that these are a consequence of equilibrium and not feedback. Overall, the analytical framework of our equilibrium scenario broadly captures the relevant physics establishing the galaxy gas and metal content in simulations, which suggests that the cycle of baryonic inflows and outflows centrally governs the cosmic evolution of these properties in typical star-forming galaxies.

Exploring the origin of Ly-alpha nebulae (blobs) requires measurements of their gas kinematics that are impossible with only the resonant, optically-thick LyA line. To define gas motions relative to the systemic velocity of the blob, the LyA line must be compared with an optically-thin line like Halpha, which is not much altered by radiative transfer effects. We obtain optical and NIR spectra of the two brightest LyA blobs from Yang et al. sample using the Magellan/MagE and VLT/SINFONI. Both the LyA and Halpha lines confirm that these blobs lie at the survey redshift, z~2.3. Within each blob, we detect several Halpha sources, which roughly correspond to galaxies seen in HST images. The Halpha detections show that these galaxies have large internal velocity dispersions (130 - 190km/s) and that, in the one system (LAB01), their velocity difference is ~440 km/s. The presence of multiple galaxies within the blobs, and those galaxies' large velocity dispersions and large relative motion, is consistent with our previous finding that LyA blobs inhabit massive dark matter halos that will evolve into those typical of rich clusters today. To determine whether the gas near the embedded galaxies is predominantly infalling or outflowing, we compare the LyA and Halpha line centers, finding that LyA is not offset (Delta LyA = +0km/s) in LAB01 and redshifted by only +230 km/s in LAB02. These offsets are small compared to those of Lyman break galaxies, which average +450 km/s and extend to about +700 km/s. We test and rule out the simplest infall models and those outflow models with super/hyper-winds, which require large outflow velocities. Because of the unknown geometry of the gas distribution and the possibility of multiple sources of LyA emission embedded in the blobs, a larger sample and more sophisticated models are required to test more complex or a wider range of infall and outflow scenarios.

We examine the growth of the stellar content of galaxies from z=3-0 in cosmological hydrodynamic simulations incorporating parameterised galactic outflows. Without outflows, galaxies overproduce stellar masses (M*) and star formation rates (SFRs) compared to observations. Winds introduce a three-tier form for the galaxy stellar mass and star formation rate functions, where the middle tier depends on differential (i.e. mass-dependent) recycling of ejected wind material back into galaxies. A tight M*-SFR relation is a generic outcome of all these simulations, and its evolution is well-described as being powered by cold accretion, although current observations at z>2 suggest that star formation in small early galaxies must be highly suppressed. Roughly one-third of z=0 galaxies at masses below M^* are satellites, and star formation in satellites is not much burstier than in centrals. All models fail to suppress star formation and stellar mass growth in massive galaxies at z<2, indicating the need for an external quenching mechanism such as black hole feedback. All models also fail to produce dwarfs as young and rapidly star-forming as observed. An outflow model following scalings expected for momentum-driven winds broadly matches observed galaxy evolution around M^* from z=0-3, which is a significant success since these galaxies dominate cosmic star formation, but the failures at higher and lower masses highlight the challenges still faced by this class of models. We argue that central star-forming galaxies are well-described as living in a slowly-evolving equilibrium between inflows from gravity and recycled winds, star formation, and strong and ubiquitous outflows that regulate how much inflow forms into stars. Star-forming galaxy evolution is thus primarily governed by the continual cycling of baryons between galaxies and intergalactic gas.

Herschel has opened new windows into studying the evolution of rapidly star-forming galaxies out to high redshifts. Today's massive starbursts are characterized by star formation rates (SFRs) of 100+ Mo/yr and display a chaotic morphology and nucleated star formation indicative of a major merger. At z~2, galaxies of similar mass and SFR are characterized by ordered rotation and distributed star formation. The emerging cold accretion paradigm provides an intuitive understanding for such differences. In it, halo accretion rates govern the supply of gas into star-forming regions, modulated by strong outflows. The high accretion rates at high-z drive more rapid star formation, while also making disks thicker and clumpier; the clumps are expected to be short-lived in the presence of strong galactic outflows as observed. Hence equivalently rapid star-formers at high redshift are not analogous to local merger-driven starbursts, but rather to local disks with highly enhanced accretion rates.

We report the first results from the Hubble Infrared Pure Parallel Imaging Extragalactic Survey, which utilizes the pure parallel orbits of the Hubble Space Telescope to do deep imaging along a large number of random sightlines. To date, our analysis includes 26 widely separated fields observed by the Wide Field Camera 3, which amounts to 122.8 sq.arcmin in total area. We have found three bright Y098-dropouts, which are candidate galaxies at z >~ 7.4. One of these objects shows an indication of peculiar variability and its nature is uncertain. The other two objects are among the brightest candidate galaxies at these redshifts known to date L>2L*. Such very luminous objects could be the progenitors of the high-mass Lyman break galaxis (LBGs) observed at lower redshifts (up to z~5). While our sample is still limited in size, it is much less subject to the uncertainty caused by "cosmic variance" than other samples because it is derived using fields along many random sightlines. We find that the existence of the brightest candidate at z~7.4 is not well explained by the current luminosity function (LF) estimates at z~8. However, its inferred surface density could be explained by the prediction from the LFs at z~7 if it belongs to the high-redshift tail of the galaxy population at z~7.

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

It is difficult to reconcile the observed evolution of the star formation rate versus stellar mass (SFR-M*) relation with expectations from current hierarchical galaxy formation models. The observed SFR-M* relation shows a rapid rise in SFR(M*) from z=0-2, and then a surprisingly lack of amplitude evolution out to z~6+. Hierarchical models of galaxy formation match this trend qualitatively but not quantitatively, with a maximum discrepancy of ~x3 in SFR at z~2. One explanation, albeit radical, is that the IMF becomes modestly weighted towards massive stars out to z~2, and then evolves back towards its present-day form by z~4 or so. We observe that this redshift trend mimics that of the cosmic fraction of obscured star formation, perhaps hinting at a physical connection. Such IMF evolution would concurrently go towards explaining persistent discrepancies between integrated measures of star formation and present-day stellar mass or cosmic colors.

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

Cosmological hydrodynamic simulations robustly predict that high-redshift galaxy star formation histories (SFHs) are smoothly-rising and vary with mass only by a scale factor. We use our latest simulations to test whether this scenario can account for recent observations at z>=6 from WFC3/IR, NICMOS, and IRAC. Our simulations broadly reproduce the observed ultraviolet (UV) luminosity functions and stellar mass densities and their evolution at z=6-8, all of which are nontrivial tests of the mean SFH. In agreement with observations, simulated galaxies possess blue UV continua owing to young ages (50-150 Myr), low metallicities (0.1-0.5 Zsun), and low dust columns (E(B-V) <= 0.05). Observations imply a near-unity slope in the stellar mass--star formation rate relation at all z=6-8, confirming the prediction that SFH shapes are invariant. Current surveys detect the majority of galaxies with stellar masses exceeding 10^9 Msun and few galaxies less massive than 10^{8.5} Msun, implying that they probe no more than the brightest 30% of the complete star formation and stellar mass densities at z>=6. Finally, we demonstrate that there is no conflict between smoothly-rising SFHs and recent clustering observations. This is because momentum-driven outflows suppress star formation in low-mass halos, leading to overall occupancies of 0.2-0.4 even though the star formation duty cycle is one. This leads to many interesting predictions at z>=4, among them that (1) optically-selected and UV-selected samples largely overlap; (2) few galaxies exhibit significantly suppressed specific star formation rates; and (3) occupancy is constant or increasing with decreasing luminosity. These predictions are in tentative agreement with current observations, but further analysis of existing and upcoming data sets is required in order to test them more thoroughly. (abridged)

Understanding the nature of distant Ly-alpha nebulae ("blobs") and connecting them to their present-day descendants requires constraining their number density, clustering, and large-scale environment. To measure these basic quantities, we conduct a deep narrowband imaging survey in four different fields, Chandra Deep Field South (CDFS), Chandra Deep Field North, and two COSMOS subfields, for a total survey area of 1.2deg^2. We discover 25 blobs at z=2.3 with Ly-alpha luminosities of 0.7-8x10^43 erg/s and isophotal areas of Aiso = 10-60 arcsec^2. The transition from compact Ly-alpha emitters (Aiso ~ a few arcsec^2) to extended blobs (Aiso > 10 arcsec^2) is continuous, suggesting a single family perhaps governed by similar emission mechanisms. Surprisingly, most blobs (16/25) are in one survey field, the CDFS. The six brightest, largest blobs with L > 1.5x10^43 erg/s and Aiso > 16 arcsec^2 lie only in the CDFS. These large, bright blobs have a field-to-field variance of sigma_v >~ 1.5 (150%) about their number density n ~ 1.0x10^-5 Mpc^-3. This variance is large, significantly higher than that of unresolved LAEs (sigma_v ~ 0.3 or 30%), and can adversely affect comparisons of blob number densities and luminosity functions among different surveys. We compare the statistics of our blobs with dark matter halos in a 1 Gpc/h cosmological N-body simulation. At z=2.3, the number density (n) implies that each bright, large blob could occupy a halo of M_halo > 10^13 Msun if most halos have detectable blobs. The predicted variance in n is consistent with that observed and corresponds to a bias of ~7. Blob halos lie at the high end of the halo mass distribution at z=2.3 and are likely to evolve into the ~10^14 Msun halos typical of galaxy clusters today. On larger scales of ~10 co-moving Mpc, blobs cluster where compact LAEs do, indicating that blobs lie in coherent, highly overdense structures.

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

The bimodality in observed present-day galaxy colours has long been a challenge for hierarchical galaxy formation models, as it requires some physical process to quench (and keep quenched) star formation in massive galaxies. Here we examine phenomenological models of quenching by post-processing the star formation histories of galaxies from cosmological hydrodynamic simulations that reproduce observations of star-forming galaxies reasonably well. We consider recipes for quenching based on major mergers, halo mass thresholds, gas temperature thresholds, and variants thereof. We compare the resulting simulated star formation histories to observed g-r colour-magnitude diagrams and red and blue luminosity functions from SDSS. The merger and halo mass quenching scenarios each yield a distinct red sequence and blue cloud of galaxies that are in broad agreement with data, albeit only under rather extreme assumptions. In detail, however, the simulated red sequence slope and amplitude in each scenario is somewhat discrepant, perhaps traceable to low metallicities in simulated galaxies. Merger quenching produces more massive blue galaxies, earlier quenching, and more frosting of young stars; comparing to relevant data tends to favor merger over halo mass quenching. Although physically-motivated quenching models can produce a red sequence, interesting generic discrepancies remain that indicate that additional physics is required to reproduce the star formation and enrichment histories of red and dead galaxies.

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

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

We present an analysis of the Ly-a forest toward 3C 273 from the Space Telescope Imaging Spectrograph at ~7 km/s resolution, along with re-processed data from the Far Ultraviolet Spectroscopic Explorer. The high UV flux of 3C 273 allows us to probe the weak, low z absorbers. The main sample consists of 21 HI absorbers that we could discriminate to a sensitivity of log NHI~ 12.5. The redshift density for absorbers with 13.1<log NHI<14.0 is ~1.5 sigma below the mean for other lines of sight; for log NHI >= 12.5, it is consistent with numerical model predictions. The Doppler parameter distribution is consistent with other low z samples. We find no evidence for a break in the column density power-law distribution to log NHI=12.3. A broad Ly-a absorber (BLA) is within Delta v =< 50 km/s and 1.3 local frame Mpc of two ~0.5L* galaxies, with an OVI absorber ~700 km/s away, similarly close to three galaxies and indicating overdense environments. We detect clustering on the Delta v<1000 km/s scale at 3.4 sigma significance for log NHI >= 12.6, consistent with the level predicted from hydrodynamical simulations, and indication for a Ly-a forest void at 0.09<z<0.12. We find at least two components for the z=0.0053 Virgo absorber, but the total NHI column is not significantly changed.

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

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

We post-process galaxy star formation histories in cosmological hydrodynamics simulations to test quenching mechanisms associated with AGN. By comparing simulation results to color-magnitude diagrams and luminosity functions of SDSS galaxies, we examine whether "quasar mode" or "radio mode" AGN feedback can yield a realistic red sequence. Both cases yield red sequences distinct from the blue cloud, decent matches to the luminosity function, and galaxies that are too blue by about 0.1 magnitudes in g-r. Our merger-based prescription for quasar mode feedback, however, yields a red sequence build-up inconsistent with observations: the luminosity function lacks a characteristic knee, and the brightest galaxies include a small number of young stars.

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

Astronomical software is now a fact of daily life for all hands-on members of our community. Purpose-built software for data reduction and modeling tasks becomes ever more critical as we handle larger amounts of data and simulations. However, the writing of astronomical software is unglamorous, the rewards are not always clear, and there are structural disincentives to releasing software publicly and to embedding it in the scientific literature, which can lead to significant duplication of effort and an incomplete scientific record. We identify some of these structural disincentives and suggest a variety of approaches to address them, with the goals of raising the quality of astronomical software, improving the lot of scientist-authors, and providing benefits to the entire community, analogous to the benefits provided by open access to large survey and simulation datasets. Our aim is to open a conversation on how to move forward. We advocate that: (1) the astronomical community consider software as an integral and fundable part of facility construction and science programs; (2) that software release be considered as integral to the open and reproducible scientific process as are publication and data release; (3) that we adopt technologies and repositories for releasing and collaboration on software that have worked for open-source software; (4) that we seek structural incentives to make the release of software and related publications easier for scientist-authors; (5) that we consider new ways of funding the development of grass-roots software; (6) and that we rethink our values to acknowledge that astronomical software development is not just a technical endeavor, but a fundamental part of our scientific practice.

IGM metal absorption lines observed in z>6 spectra offer the opportunity to probe early feedback processes, the nature of enriching sources, and the topology of reionization. We run high-resolution cosmological simulations including galactic outflows to study the observability and physical properties of 5 ions (C II, C IV, O I, Si II, Si IV) in absorption between z=8->5. We apply three cases for ionization conditions: Fully neutral, fully reionized, and a patchy model based on the flux from the nearest galaxy. We find that our simulations broadly fit available z~5-6 IGM metal-line data, although all observations cannot be accommodated with a single ionization condition. Variations in O I absorbers among sight lines seen by Becker et al. (2006) suggest significant neutral IGM patches down to z~6. Strong C IV absorbers at z~6 may be the result of ionization by their parent galaxy. Our outflows have typical speeds of ~200 km/s and mass loading factors of ~6. Such high mass loading is critical for enriching the IGM to the observed levels while curtailing star formation to match the observed z~6 rest-frame UV luminosity function. The volume filling factor of metals increases during this epoch, but only reaches ~1% for Z>10^(-3) Zsolar by z=5. C IV is an ideal tracer of IGM metals at z~5-6, with dropping global ionization fractions to either higher or lower redshifts. This results in a strongly increasing global Omega(C IV) from z=8->5, in contrast to its relative constancy from z=5->2. Our simulations do not support widespread early IGM enrichment from e.g. Pop III stars. High-z absorbers arise from metals on their first outward journey from galaxies, at distances less than 50 kpc. The galaxies responsible for early IGM enrichment have typical M*=10^(7.0-8.5) Msolar.

The O VI ion observed in quasar absorption line spectra is the most accessible tracer of the cosmic metal distribution in the low redshift (z<0.5) intergalactic medium (IGM). We explore the nature and origin of O VI absorbers using cosmological hydrodynamic simulations including galactic outflows. We consider the effects of ionization background variations, non-equilibrium ionization and cooling, uniform metallicity, and small-scale (sub-resolution) turbulence. Our main results are 1) IGM O VI is predominantly photo-ionized with T= 10^(4.2+/-0.2) K. A key reason for this is that O VI absorbers preferentially trace over-enriched regions of the IGM at a given density, which enhances metal-line cooling such that absorbers can cool within a Hubble time. As such, O VI is not a good tracer of the WHIM. 2) The predicted O VI properties fit observables only if sub-resolution turbulence is added. The required turbulence increases with O VI absorber strength such that stronger absorbers arise from more recent outflows with turbulence dissipating on the order of a Hubble time. The amount of turbulence is consistent with other examples of turbulence observed in the IGM and galactic halos. 3) Metals traced by O VI and H I do not trace exactly the same baryons, but reside in the same large-scale structure. Observed alignment statistics are reproduced in our simulations. 4) Photo-ionized O VI traces gas in a variety of environments, and is not directly associated with the nearest galaxy, though is typically nearest to ~0.1L* galaxies. Weaker O VI components trace some of the oldest cosmic metals. 5) Very strong absorbers are more likely to be collisionally ionized, tracing more recent enrichment (<2 Gyr) within or near galactic halos.

We present predictions for galactic halo baryon fractions from cosmological hydrodynamic simulations with a well-constrained model for galactic outflows. Without outflows, halos contain roughly the cosmic fraction of baryons, slightly lowered at high masses owing to pressure support from hot gas. The star formation efficiency is large and increases monotonically to low masses, in disagreement with data. With outflows, the baryon fraction is increasingly suppressed in halos to lower masses. A Milky Way-sized halo at z=0 has about 60% of the cosmic fraction of baryons, so "missing" halo baryons have largely been evacuated, rather than existing in some hidden form. Large halos (>10^13 Mo) contain 85% of their cosmic share of baryons, which explains the mild missing baryon problem seen in clusters. By comparing results at z=3 and z=0, we show that most of the baryon removal occurs at early epochs in larger halos, while smaller halos lose baryons more recently. Star formation efficiency is maximized in halos of ~10^13 Mo, dropping significantly to lower masses, which helps reconcile the sub-L* slope of the observed stellar and halo mass functions. These trends are predominantly driven by differential wind recycling, namely, that wind material takes longer to return to low-mass galaxies than high-mass galaxies. The hot gas content of halos is mostly unaffected by outflows, showing that outflows tend to blow holes and escape rather than deposit their energy into halo gas.

To obtain an unbiased sample of bright LyA blobs [L(LyA) > 10^43 ergs/s], we have undertaken a blind, wide-field, narrow-band imaging survey in the NOAO Deep Wide Field Survey Bootes field with the Steward Bok-2.3m telescope. After searching over 4.82 sq. degrees at z=2.3, we discover four LyA blobs with L(LyA) = 1.6-5.3 x 10^43 ergs/s, isophotal areas of 28-57 sq. arcsec, and broad LyA line profiles (FWHM = 900-1250 km/s). In contrast with the extended Lyman alpha halos associated with high-z radio galaxies, none of our four blobs are radio-loud. The X-ray luminosities and optical spectra of these blobs are diverse. Two blobs (3 and 4) are X-ray-detected with L_X(2-7 keV) = 2-4 x 10^44 ergs/s and have broad optical emission lines (C IV) characteristic of AGN, implying that 50% of our sample blobs are associated with strong AGN. The other 50% of blobs (1 and 2) are not X-ray or optically-detected as AGN down to similar limits. The number density of the four blobs is ~3 x 10^{-6} Mpc^{-3}, comparable to that of galaxy clusters at similar redshifts and 3x lower than that found in the SSA22 proto-cluster at z=3.1, even after accounting for the over-density of that region. The two X-ray undetected blobs are separated by only 70" (550 kpc) and have almost identical redshifts (< 360 kpc along the line-of-sight), suggesting that they are part of the same system. Given the rarity of the blobs and our discovery of a close pair, we speculate that blobs occupy the highest density regions and thus may be precursors of today's rich cluster galaxies.

We examine metal and entropy content in galaxy groups having T_X~0.5-2 keV in cosmological hydrodynamic simulations. Our simulations include a well-constrained prescription for galactic outflows following momentum-driven wind scalings, and a sophisticated chemical evolution model. Our simulation with no outflows reproduces observed iron abundances in X-ray emitting gas, but the oxygen abundance is too low; including outflows yields iron and oxygen abundances in good agreement with data. X-ray measures of [O/Fe] primarily reflect metal distribution mechanisms into hot gas, not the ratio of Type Ia to Type II supernovae within the group. Iron abundance increases by x2 from z=1-0 independent of group size, consistent with that seen in clusters, while [O/Fe] drops by ~30%. Core entropy versus temperature is elevated over self-similar predictions regardless of outflows due to radiative cooling removing low-entropy gas, but outflows provide an additional entropy contribution below 1 keV. This results in a noticeable break in the L_X-T_X relation below 1 keV, as observed. Importantly, outflows serve to reduce the stellar content of groups to observed levels. Radial profiles from simulations are in broad agreement with observations, but there remain non-trivial discrepancies that may reflect an excess of late-time star formation in central group galaxies in our simulations. Our model with outflows suggests a connection between physical processes of galaxy formation and both pre-heating and enrichment in intragroup gas, though more definitive conclusions must await a model that simultaneously suppresses cooling flows as observed.

Using Gadget-2 cosmological hydrodynamic simulations including an observationally-constrained model for galactic outflows, we investigate how feedback from star formation distributes mass, metals, and energy on cosmic scales from z=6->0. We include instantaneous enrichment from Type II SNe, delayed enrichment from Type Ia SNe and stellar (AGB) mass loss, and we individually track C, O, Si, and Fe. Following on the successes of the momentum-driven wind scalings, we improve our implementation with an on-the-fly galaxy finder to derive wind properties based on host galaxy masses. By tracking wind particles in a suite of simulations, we find: (1) Wind material reaccretes onto a galaxy on a recycling timescale that varies inversely with galaxy mass. Hence metals driven into the IGM by galactic superwinds cannot be assumed to leave their galaxy forever. Wind material is typically recycled several times; the median number of ejections for a given wind particle is 3, so by z=0 the total mass ejected in winds exceeds 0.5\Omega_b. (2) The physical distance winds travel is fairly independent of redshift and galaxy mass (60-100 physical kpc). For sizable galaxies at later epochs, winds typically do not escape the galaxy halo, and rain back down in a halo fountain. High-z galaxies enrich a significantly larger comoving volume of the IGM, with metals migrating back into galaxies to lower z. (3) The energy imparted into winds scales with M_{gal}^{1/3}, and energy from another source besides SNe (such as photons from young stars) may be required to distributed cosmic metals as observed. (4) The production of all 4 metals tracked is globally dominated by Type II SNe at all epochs. However, intracluster gas iron content triples as a result of non-Type II sources, and the low-z IGM carbon content is boosted by AGB feedback.

The evolution of the galaxy stellar mass--star formation rate relationship (M*-SFR) provides key constraints on the stellar mass assembly histories of galaxies. For star-forming galaxies, M*-SFR is observed to be fairly tight with a slope close to unity from z~0-2. Simulations of galaxy formation reproduce these trends owing to the generic dominance of smooth and steady cold accretion in these systems. In contrast, the amplitude of the M*-SFR relation evolves markedly differently than in models. Stated in terms of a star formation activity parameter alpha=(M*/SFR)/(t_H-1 Gyr), models predict a constant alpha~1 out to redshifts z=4+, while the observed M*-SFR relation indicates that alpha increases by X3 from z~2 until today. The low alpha at high-z not only conflicts with models, but is also difficult to reconcile with other observations of high-z galaxies. Systematic biases could significantly affect measurements of M* and SFR, but detailed considerations suggest that none are obvious candidates to reconcile the discrepancy. A speculative solution is considered in which the stellar initial mass function (IMF) evolves towards more high-mass star formation at earlier epochs. Following Larson, a model is investigated in which the characteristic mass Mhat where the IMF turns over increases with redshift. The observed and predicted M*-SFR evolution may be brought into agreement if Mhat=0.5(1+z)^2 Mo out to z~2. Such evolution broadly matches recent observations of cosmic stellar mass growth, and the resulting z=0 cumulative IMF is similar to the paunchy IMF favored by Fardal et al to reconcile the observed cosmic star formation history with present-day fossil light measures. [abridged]

We use cosmological hydrodynamic simulations including galactic feedback based on observations of local starbursts to find a self-consistent evolutionary model capable of fitting the observations of the intergalactic metallicity history as traced by C IV between z=6.0->1.5. Our main finding is that despite the relative invariance in the measurement of Omega(C IV) as well as the column density and linewidth distributions over this range, continual feedback from star formation-driven winds are able to reproduce the observations, while an early enrichment scenario where a majority of the metals are injected into the IGM at z>6 is disfavored. The constancy of the C IV observations results from a rising IGM metallicity content balanced by a declining C IV ionization fraction due to a 1) decreasing physical densities, 2) increasing ionization background strength, and 3) metals becoming more shock-heated at lower redshift. Our models predict that ~20x more metals are injected into the IGM between z=6->2 than at z>6. We show that the median C IV absorber at z=2 traces metals injected 1 Gyr earlier indicating that the typical metals traced by C IV are neither from very early times nor from very recent feedback.

We present predictions for the cosmic metal budget in various phases of baryons from redshift z=6-0, taken from a cosmological hydrodynamic simulation that includes a well-constrained model for enriched galactic outflows. We find that substantial amounts of metals are found in every baryonic phase at all epochs, with diffuse intergalactic gas dominating the metal budget at early epochs and stars and halo gas dominating at recent epochs. We provide a full accounting of metals in the context of the missing metals problem at z~2.5, showing that ~40% of the metals are in galaxies, and the remainder is divided between diffuse IGM gas and shocked gas in halos and filamentary structures. Comparisons with available observations of metallicity and metal mass fraction evolution show broad agreement. We predict stars have a mean metallicity of one-tenth solar already at z=6, which increases slowly to one-half solar today, while stars just forming today have typically solar metallicity. Our HI column density-weighted mean metallicity (comparable to Damped Ly-alpha system metallicities) slowly increases from one-tenth to one-third solar from z=6-1, then falls to one-quarter solar at z=0. The global mean metallicity of the universe tracks ~50% higher than that of the diffuse phase down to z~1, and by z=0 it has a value around one-tenth solar. Metals move towards higher densities and temperatures with time, peaking around the mean cosmic density at z=2 and an overdensity of 100 at z=0. We study how carbon and oxygen ions trace the path of metals in phase space, and show that OIII-OVII lines provide the most practical option for constraining intergalactic medium metals at z<2.

We investigate models of self-consistent chemical enrichment of the intergalactic medium (IGM) from z=6 to 1.5, based on hydrodynamic simulations of structure formation that incorporate galactic outflows. Our main result is that outflow parameterizations based on momentum-driven winds as seen in local starburst galaxies provide the best agreement with observations of CIV absorption at z~2-5. Such models sufficiently enrich the high-z IGM to produce a global mass density of CIV absorbers that is relatively invariant from z=5.5 to 1.5, in agreement with observations. This occurs despite an increase in the volume-averaged metallicity by x5-10 over this redshift range, because energy input from outflows causes a drop in the global ionization fraction of CIV. Comparisons to observed CIV column density and linewidth distributions and CIV-based pixel optical depth ratios provide significant constraints on wind models. Our best-fitting models show mean IGM temperatures only slightly above our no-outflow case, metal filling factors of just a few % with volume-weighted metallicities around 0.001 at z~3, significant amounts of collisionally-ionized CIV absorption, and a metallicity-density relationship that rises rapidly at low overdensities and flattens at higher ones. In general, we find that outflow speeds must be high enough to enrich the low-density IGM at early times but low enough not to overheat it, and concurrently must significantly suppress early star formation while still producing enough early metals. It is therefore non-trivial that locally-calibrated momentum-driven wind scenarios naturally yield the desired strength and evolution of outflows, and suggest that such models represent a significant step towards understanding the impact of galactic outflows on galaxies and the IGM across cosmic time.

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

We present predictions drawn from cosmological hydrodynamic simulations for the physical, photometric and emission line properties of galaxies present during the latter stages of reionization from z=9-6. We find significant numbers of galaxies that have stellar masses exceeding 10^8 Mo during this epoch, with metallicities exceeding one-thirtieth solar. Far from primeval "first-star" objects,these objects are likely to have reionized their infall regions prior to z=9, are dominated by atomic rather than molecular cooling, and are not expected to be forming Population III stars. Galaxies exhibit a slowly evolving comoving autocorrelation length from z=9-6, continuing a trend seen at lower redshifts in which the rapidly dropping bias counteracts the rapidly increasing matter clustering. These sources can be marginally detected using current instruments, but modest increases in sensitivity or survey area would yield significantly increased samples. We compare to current observations of the z~6 rest-UV and Ly-alpha line luminosity functions, and find good agreement. We also compare with the z~7 object observed by Egami et al., and find that such systems are ubiquitous in our simulations. The intrinsic Ly-alpha luminosity function evolves slowly from z=9-6, meaning that it should also be possible to detect these objects with upcoming narrow band surveys such as DAzLE. We make predictions for near-IR surveys with JWST, and show that while a high density of sources will be found, Population III objects may remain elusive. We present and compare simulations with several recipes for superwind feedback, and show that while our broad conclusions are insensitive to this choice, a feedback model based on momentum-driven winds is favored in comparisons with available data. (abridged)

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

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

We present the predictions for the photometric and emission line properties of galaxies present during the latter stages of reionization from z=8 to 6. These preliminary predictions are made from cosmological hydrodynamic simulations that include star formation and feedback, but not the effects of radiative transfer. We find significant numbers of galaxies that have stellar masses exceeding 10^8 Mo by z=8, with metallicities in the range of one-tenth solar. These galaxies are just beyond the reach of current near-infrared surveys, but should be found in large numbers by next-generation programs. The Lyman alpha luminosity function does not evolve much from z=6 to z=8, meaning that it should also be possible to detect these objects in significant numbers with upcoming narrow band surveys, unless the escape fraction of Ly-alpha evolves significantly between those epochs.

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

We present a new implicit numerical algorithm for the calculation of the time dependent non-Local Thermodynamic Equilibrium of a gas in an external radiation field that is accurate, fast and unconditionally stable for all spatial and temporal increments. The method is presented as a backward difference scheme in 1-D but can be readily generalised to 3-D. We apply the method for calculating the evolution of ionisation domains in a hydrogen plasma with plane-parallel Gaussian density enhancements illuminated by sources of UV radiation. We calculate the speed of propagation of ionising fronts through different ambient densities and the interaction of such ionising fronts with density enhancements. We show that for a typical UV source that may be present in the early universe, the introduction of a density enhancement of a factor ~10 above an ambient density 10^{-4} atoms/cm^3 could delay the outward propagation of an ionisation front by millions of years. Our calculations show that within the lifetime of a single source (~a few million years), and for ambient intergalactic densities appropriate to redshifts z~6-20, neutral fractions of of ~10^{-3}-10^{-5} can be achieved within its zone of influence. We also present calculations which demonstrate that once started, ionisation will proceed very efficiently as multiple sources are subsequently introduced, even if the time between the appearence of such sources may be much longer than their lifetimes.

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

I present a review by epoch of baryons in the intergalactic medium (IGM), from the first star until today. Recent observations indicate a protracted period of reionization, suggesting multiple populations of reionizers; detection of these z>6 sources is a key goal that is now coming within reach. The optical Lyman alpha forest (2<z<4) is well-described by the Fluctuating Gunn-Peterson Approximation, but puzzling observations of associated metal lines and nearby galaxies may provide insights into galactic feedback processes. IGM studies at z<1.5 are progressing rapidly thanks to ultraviolet absorption line studies of both Ly-alpha absorbers and the warm-hot intergalactic medium. Cosmological hydrodynamic simulations have played an integral part in these advances, and have helped to shape our understanding of the IGM at each epoch. Working jointly, observations and theory continue to expand our knowledge of the IGM as the earliest stages of galaxy formation and the dominant reservoir of baryons at all redshifts.

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

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

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

We present a preliminary study of the metallicity of low-redshift Lyman alpha absorbers from the HST/STIS spectrum of 3C273. Using a pixel-based shift-and-coadd technique, we compare observations to carefully-constructed mock quasar spectra from a cosmological hydrodynamic simulation. We place an upper limit of [C/H]<-1.5 at >90% confidence using two Ly-alpha absorbers with N_HI>10^14 cm^-2, from the fact that we see little absorption near the expected CIV positions, whereas the artificial spectra predict significant absorption. We assume a Haardt & Madau (quasar-based) ionizing background spectrum, but a softer spectrum results in a stronger constraint on the metallicity. With only two absorbers and various modeling uncertainties, it is too early to make a definitive determination of the diffuse IGM metallicity, but upcoming observations should improve constraints.

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

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

We examine the Ly-alpha absorber population at z<0.3 detected in spectra of the QSOs PG0953+415 and H1821+643 taken with the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope. We compare their statistical properties to those in carefully-constructed mock quasar spectra drawn from a cosmological hydrodynamic simulation of a Lambda-CDM universe. We find very good agreement in the column density and b-parameter distributions, down to the smallest observable absorbers with N_HI\approx 10^{12.3} cm^-2. The observed absorber population is complete for N_HI\ga 10^{13} cm^-2, with a column density distribution slope of \beta=2.04\pm 0.23 and a median b-parameter of 21 km/s above this limit. The intergalactic gas giving rise to these weak absorbers is analogous to that at high redshift, located in diffuse large-scale structures that are highly photoionized by the metagalactic UV flux, though a greater number arise within shock-heated warm gas. The density, temperature, and column density of these absorbers follow similar relationships to those at high redshift, though with substantially larger scatter due to the shock-heated gas. The b-parameters typically have a significant contribution from thermal broadening, which facilitates a measurement of the low-z IGM temperature as traced by Ly-alpha absorbers. From our simulation we estimate T_IGM\sim 5000 K, with an upper limit of 10^4 K, at the mean density. The agreement in predicted and observed amplitude of the column density distributions allows us to measure the HI photoionization rate at z=0.17 to be 10^{-13.3\pm 0.7} s^{-1} (estimated modeling uncertainty), close to predictions based on quasar properties.

I review results obtained from studies of the high-redshift Lyman alpha forest, and present new results from HST/STIS spectra of low-redshift quasars in comparison with cosmological hydrodynamic simulations. The evolution of the Lyman alpha forest from redshift 3 to 0 is well-described by current structure formation models, in which Ly-a forest absorbing gas at all redshifts traces moderate-overdensity large-scale structures. I describe some of the insights provided by hydrodynamic simulations into the observed statistical trends of Ly-a absorbers and the physical state of the absorbing gas.

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

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

We present a comparison of halo properties in cosmological simulations of collisionless cold dark matter (CDM) and self-interacting dark matter (SIDM) for a range of dark matter cross sections. We find, in agreement with various authors, that CDM yields cuspy halos that are too centrally concentrated as compared to observations. Conversely, SIDM simulations using a Monte Carlo N-body technique produce halos with significantly reduced central densities and flatter cores with increasing cross section. We introduce a concentration parameter based on enclosed mass that we expect will be straightforward to determine observationally, unlike that of Navarro, Frenk & White, and provide predictions for SIDM and CDM. SIDM also produces more spherical halos than CDM, providing possibly the strongest observational test of SIDM. We discuss our findings in relation to various relevant observations as well as SIDM simulations of other groups. Taking proper account of simulation limitations, we find that a dark matter cross section per unit mass of sigma_DM ~= 10^{-23}-10^{-24} cm^2/GeV is consistent with all current observational constraints.

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

High resolution N-body simulations of cold dark matter (CDM) models predict that galaxies and clusters have cuspy halos with excessive substructure. Observations reveal smooth halos with central density cores. One possible resolution of this conflict is that the dark matter is warm (WDM); this will suppress the power spectrum on small scales. The Lyman-alpha forest is a powerful probe of the linear power spectrum on these scales. We use collisionless N-body simulations to follow the evolution of structure in WDM models, and analyze artificial Lyman-alpha forest spectra extracted from them. By requiring that there is enough small-scale power in the linear power spectrum to reproduce the observed properties of the Lyman-alpha forest in quasar spectra, we derive a lower limit to the mass of the WDM particle of 750 eV. This limit is robust to reasonable uncertainties in our assumption about the temperature of the mean density gas (T0) at z=3. We argue that any model that suppresses the CDM linear theory power spectrum more severely than a 750 eV WDM particle cannot produce the Lyman-alpha forest.

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

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

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

A new class of geometric statistics for analyzing galaxy catalogs is presented. Filament statistics quantify filamentarity and planarity in large scale structure in a manner consistent with catalog visualizations. These statistics are based on sequences of spatial links which follow local high-density structures. From these link sequences we compute the discrete curvature, planarity, and torsion. Filament statistics are applied to CDM and CHDM ($\Omega_\nu = 0.3$) simulations of Klypin \etal (1996), the CfA1-like mock redshift catalogs of Nolthenius, Klypin and Primack (1994, 1996), and the CfA1 catalog. We also apply the moment-based shape statistics developed by Babul \& Starkman (1992), Luo \& Vishniac (1995), and Robinson \& Albrecht (1996) to these same catalogs, and compare their robustness and discriminatory power versus filament statistics. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find discrimination of $\sim 4\sigma$ (where $\sigma$ represents resampling errors) between CHDM and CDM for selected filament statistics and shape statistics, including variations in the galaxy identification scheme. Comparing the CfA1 data versus the models does not yield a conclusively favored model; no model is excluded at more than a $\sim 2\sigma$ level for any statistic, not including cosmic variance which could further degrade the discriminatory power. We find that CfA1 discriminates between models poorly mainly due to its sparseness and small number of galaxies, not due to redshift distortion, magnitude limiting, or geometrical effects. We anticipate that the proliferation of large redshift surveys and simulations will enable the statistics presented here to provide robust discrimination between large-scale structure in various cosmological models.

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

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

A new class of geometric statistics for analyzing galaxy catalogs is presented. {\it Filament statistics} quantify filamentarity and planarity in large scale structure in a manner consistent with catalog visualizations. These statistics are based on sequences of spatial links which follow local high-density structures. From these link sequences we compute the discrete curvature, planarity, and torsion. Filament statistics are applied to CDM and CHDM ($\Omega_\nu = 0.3$) simulations of Klypin \etal (1995), the CfA1-like sky catalogs of Nolthenius, Klypin and Primack (1994, 1995), and the CfA1 catalog. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find robust discrimination of over 4$\sigma$ (where $\sigma$ represents resampling errors) between CHDM and CDM. The {\it reduced filament statistics} show that CfA1 data is intermediate between CHDM and CDM, but more consistent with the CHDM models. Filament statistics provide robust and discriminatory shape statistics with which to test cosmological simulations of various models against present and future redshift survey data.