Little is known about the mass-metallicity relation (MZR) in galaxies at cosmic dawn. Studying the first appearance of the MZR is one of the keys to understanding the formation and evolution of the first galaxies. In order to lay the groundwork for upcoming observational campaigns, we analyze 290 galaxies in halos spanning Mh = 10^9-10^11 Msun selected from the FirstLight (FL) cosmological zoom simulations to predict the MZR at z=5-8. Over this interval, the metallicity of FL galaxies with stellar mass Ms= 10^8 Msun declines by < 0.2 dex. This contrasts with the observed tendency for metallicities to increase at lower redshifts, and reflects weakly-evolving or even increasing gas fractions. We assess the use of the R3 strong-line diagnostic as a metallicity indicator, finding that it is informative for 12+log(O/H) < 8 but saturates to R3=3 at higher metallicities owing to a cancellation between enrichment and spectral softening. Nonetheless, campaigns with JWST should be able to detect a clear trend between R3 and stellar mass for Ms > 10^7.5 Msun. We caution that, at fixed metallicity, galaxies with higher specific star formation show higher R3 owing to their more intense radiation fields, indicating a potential for selection biases.

We explore the expected galaxy environments of CIV absorbers at z>5 using the Technicolor Dawn simulations. These simulations reproduce the observed history of reionization, the z~6 galaxy stellar mass function, the Ly$\alpha$ forest transmission at z>5, and the SiIV column density distribution (CDD) at z~5.5. Nonetheless, the CIV CDD remains underproduced. Comparison with observed CII/SiII equivalent width ratios and the CII line incidence suggests that a low carbon yield accounts for some, but not all, of the CIV discrepancy. Alternatively, a density-bounded escape scenario could harden the metagalactic ionizing background more dramatically even than binary stellar evolution, boosting the CIV CDD into near-agreement with observations. In this case galaxies ionize more efficiently and fewer are required to host a given high-ionization absorber. Absorbers' environments therefore constrain ionizing escape. Regardless of the escape scenario, galaxies correlate with CIV absorbers out to 300 proper kpc (pkpc). The correlation strengthens independently with galaxy luminosity and CIV column density. Around strong systems (log(N$_{\rm CIV}$/cm$^{-2}$)>14)), the overdensity of galaxies with M$_{\rm UV}$<-18 or log($L_{{\rm Ly}\alpha}$/erg s$^{-1}$) > 41.9 declines from 200-300 within 100 pkpc to 40-60 within 250 pkpc. The previously-suggested association between strong CIV absorbers and Ly$\alpha$ emitters at z>5 is not expected. It may arise if both populations inhabit large-scale voids, but for different reasons. Although most neighboring galaxies are too faint for HST, JWST will, with a single pointing, identify ~10 neighboring galaxies per strong CIV absorber at z > 5. Ground-based tests of these predictions are possible via deep surveys for Ly$\alpha$ emission using integral field units.

We study six luminous Lyman-alpha emitters (LAEs) with very blue rest-frame UV continua at $5.7\le z \le 6.6$. These LAEs have previous HST and Spitzer IRAC observations. Combining our newly acquired HST images, we find that their UV-continuum slopes $\beta$ are in a range of $-3.4\le \beta \le -2.6$. Unlike previous, tentative detections of $\beta \simeq -3$ in photometrically selected, low-luminosity galaxies, our LAEs are spectroscopically confirmed and luminous ($M_{\rm UV}<-20$ mag). We model their broadband spectral energy distributions (SEDs), and find that two $\beta\simeq-2.6\pm0.2$ galaxies can be well fitted with young and dust-free stellar populations. However, it becomes increasingly difficult to fit bluer galaxies. We explore further interpretations by including non-zero LyC escape fraction $f_{\rm esc}$, very low metallicities, and/or AGN contributions. Assuming $f_{\rm esc}\simeq0.2$, we achieve the bluest slopes $\beta\simeq-2.7$ when nebular emission is considered. This can nearly explain the SEDs of two galaxies with $\beta\simeq-2.8$ and --2.9 ($\sigma_{\beta}=0.15$). Larger $f_{\rm esc}$ values and very low metallicities are not favored by the strong nebular line emission (evidenced by the IRAC flux) or the observed (IRAC 1 - IRAC 2) color. Finally, we find that the $\beta\simeq-2.9$ galaxy can potentially be well explained by the combination of a very young population with a high $f_{\rm esc}$ ($\ge0.5$) and an old, dusty population. We are not able to produce two $\beta \simeq -3.4 \pm0.4$ galaxies. Future deep spectroscopic observations are needed to fully understand these galaxies.

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 use synthetic sightlines drawn through snapshots of the Technicolor Dawn simulations to explore how the statistics of neutral oxygen OI absorbers respond to hydrogen reionization. The ionization state of the circumgalactic medium (CGM) initially roughly tracks that of the intergalactic medium, but beginning at $z=8$ the CGM grows systematically more neutral owing to self-shielding. Weak absorbers trace diffuse gas that lies farther from halos, hence they are ionized first, whereas stronger systems are less sensitive to reionization. The overall OI covering fraction decreases slowly with time owing to competition between ongoing enrichment and gradual encroachment of ionization fronts into increasingly overdense gas. While the declining covering fraction is partially offset by continued formation of new halos, the ionization of the diffuse gas causes the predicted line-of-sight incidence rate of OI absorbers to decline abruptly at the overlap epoch, in qualitative agreement with observations. In comparison to the recently-observed equivalent width (EW) distribution at $z\approx6$, the simulations underproduce systems with $EW \geq 0.1 \unicode{x212B}$, although they reproduce weaker systems with $EW \geq 0.05 \unicode{x212B}$. By $z\approx5$, the incidence of $EW < 0.1 \unicode{x212B}$ systems are overproduced, consistent with previous indications that the simulated ionizing background is too weak at $z<6$. The summed column densities of SiII and SiIV trace the total oxygen column, and hence the ratio of the OI and SiII + SiIV comoving mass densities traces the progress of reionization. This probe may prove particularly useful in the regime where $x_{HI} > 10\%$

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

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.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

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.

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.

We test the galactic outflow model by probing associated galaxies of four strong intergalactic CIV absorbers at $z=5$--6 using the Hubble Space Telescope (HST) ACS ramp narrowband filters. The four strong CIV absorbers reside at $z=5.74$, $5.52$, $4.95$, and $4.87$, with column densities ranging from $N_{\rm{CIV}}=10^{13.8}$ cm$^{-2}$ to $10^{14.8}$ cm$^{-2}$. At $z=5.74$, we detect an i-dropout Ly$\alpha$ emitter (LAE) candidate with a projected impact parameter of 42 physical kpc from the CIV absorber. This LAE candidate has a Ly$\alpha$-based star formation rate (SFR$_{\rm{Ly\alpha}}$) of 2 $M_\odot$ yr$^{-1}$ and a UV-based SFR of 4 $M_\odot$ yr$^{-1}$. Although we cannot completely rule out that this $i$-dropout emitter may be an [OII] interloper, its measured properties are consistent with the CIV powering galaxy at $z=5.74$. For CIV absorbers at $z=4.95$ and $z=4.87$, although we detect two LAE candidates with impact parameters of 160 kpc and 200 kpc, such distances are larger than that predicted from the simulations. Therefore we treat them as non-detections. For the system at $z=5.52$, we do not detect LAE candidates, placing a 3-$\sigma$ upper limit of SFR$_{\rm{Ly\alpha}}\approx 1.5\ M_\odot$ yr$^{-1}$. In summary, in these four cases, we only detect one plausible CIV source at $z=5.74$. Combining the modest SFR of the one detection and the three non-detections, our HST observations strongly support that smaller galaxies (SFR$_{\rm{Ly\alpha}} \lesssim 2\ M_\odot$ yr$^{-1}$) are main sources of intergalactic CIV absorbers, and such small galaxies play a major role in the metal enrichment of the intergalactic medium at $z\gtrsim5$.

The recent low value of Planck (2016) integrated optical depth to Thomson scattering suggests that the reionization occurred fairly suddenly, disfavoring extended reionization scenarios. This will have a significant impact on the 21cm power spectrum. Using a semi-numerical framework, we improve our model from Hassan et al. (2016) to include time-integrated ionisation and recombination effects, and find that this leads to more sudden reionisation. It also yields larger HII bubbles which leads to an order of magnitude more 21cm power on large scales, while suppressing the small scale ionization power. Local fluctuations in the neutral hydrogen density play the dominant role in boosting the 21cm power spectrum on large scales, while recombinations are subdominant. We use a Monte Carlo Markov Chain approach to constrain our model to observations of the star formation rate functions at z = 6,7,8 from Bouwens et al. (2015), the Planck (2016) optical depth measurements, and the Becker & Bolton (2013) ionising emissivity data at z~5. We then use this constrained model to perform 21cm forecasting for LOFAR, HERA, and SKA in order to determine how well such data can characterise the sources driving reionisation. We find that the 21cm power spectrum alone can somewhat constrain the halo mass dependence of ionising sources, the photon escape fraction and ionising amplitude, but combining the 21cm data with other current observations enables us to separately constrain all these parameters. Our framework illustrates how 21cm data can play a key role in understanding the sources and topology of reionisation as observations improve.

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {\AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

This chapter reviews how galactic inflows influence galaxy metallicity. The goal is to discuss predictions from theoretical models, but particular emphasis is placed on the insights that result from using models to interpret observations. Even as the classical G-dwarf problem endures in the latest round of observational confirmation, a rich and tantalizing new phenomenology of relationships between $M_*$, $Z$, SFR, and gas fraction is emerging both in observations and in theoretical models. A consensus interpretation is emerging in which star-forming galaxies do most of their growing in a quiescent way that balances gas inflows and gas processing, and metal dilution with enrichment. Models that explicitly invoke this idea via equilibrium conditions can be used to infer inflow rates from observations, while models that do not assume equilibrium growth tend to recover it self-consistently. Mergers are an overall subdominant mechanism for delivering fresh gas to galaxies, but they trigger radial flows of previously-accreted gas that flatten radial gas-phase metallicity gradients and temporarily suppress central metallicities. Radial gradients are generically expected to be steep at early times and then flattened by mergers and enriched inflows of recycled gas at late times. However, further theoretical work is required in order to understand how to interpret observations. Likewise, more observational work is needed in order to understand how metallicity gradients evolve to high redshifts.

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.

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.

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.

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.

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$\alpha$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$\alpha$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$\alpha$ emission are indistinguishable in terms of age, stellar mass, and SFR.

From a large sample of $\approx 170,000$ local SDSS galaxies, we find that the Fundamental Metallicity Relation (FMR) has an overabundance of outliers, compared to what would be expected from a Gaussian distribution of residuals, with significantly lower metallicities than predicted from their stellar mass and star formation rate (SFR). This low-metallicity population has lower stellar masses, bimodial specific SFRs with enhanced star formation within the aperture and smaller half-light radii than the general sample, and is hence a physically distinct population. We show that they are consistent with being galaxies that are merging or have recently merged with a satellite galaxy. In this scenario, low-metallicity gas flows in from large radii, diluting the metallicity of star-forming regions and enhancing the specific SFR until the inflowing gas is processed and the metallicity has recovered. We introduce a simple model in which mergers with a mass ratio larger than a minimum dilute the central galaxy's metallicity by an amount that is proportional to the stellar mass ratio for a constant time, and show that it provides an excellent fit to the distribution of FMR residuals. We find the dilution time-scale to be $\tau=1.568_{-0.027}^{+0.029}$ Gyr, the average metallicity depression caused by a 1:1 merger to be $\alpha=0.2480_{-0.0020}^{+0.0017}$ dex and the minimum mass ratio merger that can be discerned from the intrinsic Gaussian scatter in the FMR to be $\xi_\text{min}=0.2030_{-0.0095}^{+0.0127}$ (these are statistical errors only). From this we derive that the average metallicity depression caused by a merger with mass ratio between 1:5 and 1:1 is 0.114 dex.

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.

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

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.

We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

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 results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 \AA. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 \AA) that could be the dominant contribution to the total ionizing flux at z >= 6.

We present abundance measurements from 26 emission-line selected galaxies at z~0.6-0.7. By reaching stellar masses as low as 10^8 M_{\sun}, these observations provide the first measurement of the intermediate redshift mass-metallicity (MZ) relation below 10^9 M_{\sun} For the portion of our sample above M > 10^9 M_{\sun} (8/26 galaxies), we find good agreement with previous measurements of the intermediate redshift MZ relation. Compared to the local relation, we measure an evolution that corresponds to a 0.12 dex decrease in oxygen abundances at intermediate redshifts. This result confirms the trend that metallicity evolution becomes more significant towards lower stellar masses, in keeping with a downsizing scenario where low mass galaxies evolve onto the local MZ relation at later cosmic times. We show that these galaxies follow the local fundamental metallicity relation, where objects with higher specific (mass-normalized) star formation rates (SFRs) have lower metallicities. Furthermore, we show that the galaxies in our sample lie on an extrapolation of the SFR-M_{*} relation (the star-forming main sequence). Leveraging the MZ relation and star-forming main sequence (and combining our data with higher mass measurements from the literature), we test models that assume an equilibrium between mass inflow, outflow and star formation. We find that outflows are required to describe the data. By comparing different outflow prescriptions, we show that momentum driven winds can describe the MZ relation; however, this model under-predicts the amount of star formation in low mass galaxies. This disagreement may indicate that preventive feedback from gas-heating has been overestimated, or it may signify a more fundamental deviation from the equilibrium assumption.

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.

Observations of the Lyman-alpha forest and of high-redshift galaxies at z~5-10 imply that there were just enough photons to maintain the universe in an ionized state at z~5-6, indicating a "photon-starved" end to reionization. The ionizing emissivity must have been larger at earlier times in order to yield the extended reionization history implied by the electron scattering optical depth constraint from WMAP. Here we address the possibility that a faint population of galaxies with host halo masses of ~1e8-1e9 Msun dominated the ionizing photon budget at redshifts of about z>9, due to their much higher escape fractions. Such faint, early galaxies, would not have formed in ionized regions due to suppression by heating from the UV background (UVB), and would therefore not contribute to the ionizing background at z<6, after reionization is complete. Our model matches: (1) the low escape fractions observed for high-redshift galaxies, (2) the WMAP constraint of tau_es~0.09, (3) the low values for the UVB at z<6, and (4) the observed star formation rate density inferred from Lyman-break galaxies. A top heavy IMF from Pop III stars is not required in this scenario. We compare our model to recent ones in the literature that were forced to introduce an escape fraction that increases strongly towards high redshift, and show that a similar evolution occurs naturally if low mass galaxies possess high escape fractions.

Recent studies have found a dramatic difference between the observed number density evolution of low mass galaxies and that predicted by semi-analytic models. While models accurately reproduce the z=0 number density, they require that the evolution occurs rapidly at early times, which is incompatible with the strong late evolution found observationally. We report here the same discrepancy in two state-of-the-art cosmological hydrodynamical simulations, which is evidence that the problem is fundamental. We search for the underlying cause of this problem using two complementary methods. Firstly, we look for evidence of a different history of today's low mass galaxies in models and observations and we find that the models yield too few young, strongly star-forming galaxies. Secondly, we construct a toy model to link the observed evolution of specific star formation rates (sSFR) with the evolution of the galaxy stellar mass function. We infer from this model that a key problem in both semi-analytic and hydrodynamical models is the presence of a positive instead of a negative correlation between sSFR and stellar mass. A similar positive correlation is found between the specific dark matter halo accretion rate and the halo mass, indicating that model galaxies are growing in a way that follows the growth of their host haloes too closely. It therefore appears necessary to find a mechanism that decouples the growth of low mass galaxies, which occurs primarily at late times, from the growth of their host haloes, which occurs primarily at early times. We argue that the current form of star-formation driven feedback implemented in most galaxy formation models is unlikely to achieve this goal, owing to its fundamental dependence on host halo mass and time. [Abridged]

We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synthetic stellar population models, and measure the value of the UV spectral slope (beta) from the best-fit model spectrum. The median value of beta evolves significantly from -1.82 (+0.00,-0.04) at z = 4, to -2.37 (+0.26,-0.06) at z = 7. Additionally, we find that faint galaxies at z = 7 have beta = -2.68 (+0.39,-0.24) (~ -2.4 after correcting for observational bias); this is redder than previous claims in the literature, and does not require "exotic" stellar populations to explain their colors. This evolution can be explained by an increase in dust extinction, with the timescale consistent with low-mass AGB stars forming the bulk of the dust. We find no significant (< 2-sigma) correlation between beta and M_UV when measuring M_UV at a consistent rest-frame wavelength of 1500 A. This is particularly true at bright magnitudes, though our results do show evidence for a weak correlation at faint magnitudes when galaxies in the HUDF are considered separately, hinting that dynamic range in sample luminosities may play a role. We do find a strong correlation between beta and the stellar mass at all redshifts, in that more massive galaxies exhibit redder colors. The most massive galaxies in our sample have red colors at each redshift, implying that dust can build up quickly in massive galaxies, and that feedback is likely removing dust from low-mass galaxies at z > 7. Thus the stellar-mass - metallicity relation, previously observed up to z ~ 3, may extend out to z = 7 - 8.

We present measurements of the specific ultraviolet luminosity density from a sample of 483 galaxies at 6<z<8. These galaxies were selected from new deep near-infrared HST imaging from the CANDELS, HUDF09 and ERS programs. In contrast to the majority of previous analyses, which assume that the distribution of galaxy ultraviolet (UV) luminosities follows a Schechter distribution, and that the distribution continues to luminosities far below our observable limit, we investigate the contribution to reionization from galaxies which we can observe, free from these assumptions. We find that the observable population of galaxies can sustain a fully reionized IGM at z=6, if the average ionizing photon escape fraction (f_esc) is ~30%. A number of previous studies have measured UV luminosity densities at these redshifts that vary by 5X, with many concluding that galaxies could not complete reionization by z=6 unless a large population of galaxies fainter than the detection limit were invoked, or extremely high values of f_esc were present. The observed UV luminosity density from our observed galaxy samples at z=7-8 is not sufficient to maintain a fully reionized IGM unless f_esc>50%. Combining our observations with constraints on the emission rate of ionizing photons from Ly-alpha forest observations at z=6, we can constrain f_esc<34% (2-sigma) if the observed galaxies are the only contributors to reionization, or <13% (2-sigma) if the luminosity function extends to M_UV = -13. These escape fractions are sufficient to complete reionization by z=6. These constraints imply that the volume ionized fraction of the IGM becomes less than unity at z>7, consistent with a number of complementary reionization probes. If faint galaxies dominate reionization, future JWST observations will probe deep enough to see them, providing an indirect constraint on the ionizing photon escape fraction [abridged].

I discuss recent advances in the study of hydrogen reionization, focusing on progress that was achieved during the years 2010-2011. First, I discuss recent measurements of the progress of reionization. Next, I discuss recent observational constraints on the nature and abundance of the dominant ionizing sources. Finally, I discuss recent progress in modeling reionization. This review is written for an audience of astronomers who do not specialize in the high-redshift Universe.

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.

The stars in the dwarf spheroidal satellite galaxies (dSphs) of the Milky Way are significantly more metal-poor than would be expected from a closed box model of chemical evolution. Gas outflows likely carried away most of the metals produced by the dSphs. Based on previous Keck/DEIMOS observations and models, we calculate the mass in Mg, Si, Ca, and Fe expelled from each of eight dSphs. Essentially, these masses are the differences between the observed amount of metals present in the dSphs' stars today and the inferred amount of metals produced by supernovae. We conclude that the dSphs lost 96% to >99% of the metals their stars manufactured. We apply the observed mass function of Milky Way dSphs to the ejected mass function to determine that a single large dSph, like Fornax, lost more metals over 10 Gyr than all smaller dSphs combined. Therefore, small galaxies like dSphs are not significant contributors to the metal content of the intergalactic medium. Finally, we compare our ejected mass function to previous X-ray measurements of the metal content of the winds from the post-starburst dwarf irregular galaxy NGC 1569. Remarkably, the most recent starburst in that galaxy falls exactly on the ejected mass-stellar mass relation defined by the Milky Way dSphs.

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

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

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

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.

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)

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 study the topology of reionization using accurate three-dimensional radiative transfer calculations post-processed on outputs from cosmological hydrodynamic simulations. In our simulations, reionization begins in overdense regions and then "leaks" directly into voids, with filaments reionizing last owing to their combination of high recombination rate and low emissivity. This result depends on the uniquely-biased emissivity field predicted by our prescriptions for star formation and feedback, which have previously been shown to account for a wide array of measurements of the post-reionization Universe. It is qualitatively robust to our choice of simulation volume, ionizing escape fraction, and spatial resolution (in fact it grows stronger at higher spatial resolution) even though the exact overlap redshift is sensitive to each of these. However, it weakens slightly as the escape fraction is increased owing to the reduced density contrast at higher redshift. We also explore whether our results are sensitive to commonly-employed approximations such as using optically-thin Eddington tensors or substantially altering the speed of light. Such approximations do not qualitatively change the topology of reionization. However, they can systematically shift the overlap redshift by up to $\Delta z\sim 0.5$, indicating that accurate radiative transfer is essential for computing reionization. Our model cannot simultaneously reproduce the observed optical depth to Thomson scattering and ionization rate per hydrogen atom at $z=6$, which could owe to numerical effects and/or missing early sources of ionization.

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.

We introduce a new code for computing time-dependent continuum radiative transfer and non-equilibrium ionization states in static density fields with periodic boundaries. Our code solves the moments of the radiative transfer equation, closed by an Eddingtion tensor computed using a long characteristics method. We show that pure (i.e., not source-centered) short characteristics and the optically-thin approximation are inappropriate for computing Eddington factors for the problem of cosmological reionization. We evolve the non-equilibrium ionization field via an efficient and accurate (errors <1%) technique that switches between fully implicit or explicit finite-differencing depending on whether the local timescales are long or short compared to the timestep. We tailor our code for the problem of cosmological reionization. In tests, the code conserves photons, accurately treats cosmological effects, and reproduces analytic Stromgren sphere solutions. Its chief weakness is that the computation time for the long characteristics calculation scales relatively poorly compared to other techniques (t_{LC} \propto N_{cells}^1.5); however, we mitigate this by only recomputing the Eddington tensor when the radiation field changes substantially. Our technique makes almost no physical approximations, so it provides a way to benchmark faster but more approximate techniques. It can readily be extended to evolve multiple frequencies, though we do not do so here. Finally, we note that our method is generally applicable to any problem involving the transfer of continuum radiation through a periodic volume.

(Abridged) Using cosmological hydrodynamic simulations in combination with analytic modeling, we show that the galaxy stellar mass-metallicity relation (MZR) provides strong constraints on galactic outflows across cosmic time. We compare three outflow models: No outflows, a "constant wind" (cw) model that emulates the popular Dekel & Silk (1986) scenario, and a ``momentum-driven wind" (vzw) model. We find that only the vzw scaling reproduces the observed z~2 MZR's slope, amplitude, and scatter. Comparing our fully three-dimensional simulations with a simple one-zone chemical evolution model, we find that the MZR can be understood in terms of three parameters: (1) The equilibrium metallicity Z_eq=y*SFR/ACC (where y=true yield), reflecting the enrichment balance between star formation rate SFR and gas accretion rate ACC; (2) the dilution time t_d=M_g/ACC, representing the timescale for a galaxy to return to Z_eq after a metallicity-perturbing interaction; and (3) the blowout mass M_blowout, which is the galaxy stellar mass above which winds can escape its halo. When outflows with mass loading factor MLF are present, galaxies below M_blowout obey Z_eq = y/(1+MLF), while above M_blowout, Z_eq->y. Our cw model has M_blowout ~ 10^10 M_sun, which yields a sharp upturn in the MZR above this scale and a flat MZR with large scatter below it, in strong disagreement with observations. Our vzw model reproduces the observed Z_g \propto M_*^0.3 because Z_eq \propto MLF^-1 \propto M_*^1/3 when MLF >> 1 (i.e. at low masses). The flattening of the MZR at M_* > 10^10.5 M_sun observed by Tremonti et al. (2004) reflects the mass scale where MLF~1, rather than a characteristic wind speed. The tight observed MZR scatter is ensured when t_d<1 dynamical time, which is only satisified at all masses and epochs in our momentum-driven wind model.

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

We use 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 use cosmological hydrodynamic simulations with enriched galactic outflows to compare predictions for the galaxy mass-metallicity (M*-Z) relation versus observations at z~2 from Erb et al. (2006). With no outflows, galaxies are over-enriched, indicating that outflows are required not only to suppress star formation and enrich the IGM but also to lower galaxy metal content. The observed M*-Z slope is matched both in our model without winds as well as our favored outflow model where the outflow velocity scales as the escape velocity, but is too steep in a model with constant outflow speeds. If outflows are too widespread at early times, the IGM out of which smaller galaxies form can become pre-polluted, resulting in a low-mass flattening of the M*-Z relation that is inconsistent with data. Remarkably, the same momentum-driven wind model that provides the best agreement with IGM enrichment data also yields the best agreement with the z~2 M*-Z relation, showing the proper outflow scaling and strength to match the observed slope and amplitude. In this model, the M*-Z relation evolves slowly from z=6-2; an (admittedly uncertain) extrapolation to z=0 broadly matches local M*-Z observations. Overall, the M*-Z relation provides critical constraints on galactic outflow processes during the heydey of star formation in the Universe.

(Abridged) We present the evolution of the volume averaged properties of the rest-frame optically luminous galaxy population to z~3, determined from four disjoint deep fields with optical to near-infrared wavelength coverage. We select galaxies above a rest-frame V-band luminosity of 3x10^10 Lsol and characterize their rest-frame UV through optical properties via the mean spectral energy distribution (SED). To measure evolution we apply the same selection criteria to a sample of galaxies from the Sloan Digital Sky Survey and COMBO-17. The mean rest-frame 2200Ang through V-band SED becomes steadily bluer with increasing redshift but at z<3 the mean SED falls within the range defined by ``normal'' galaxies in the nearby Universe. We measure stellar mass-to-light ratios (Mstar/L) by fitting models to the rest-frame UV-optical SEDs and derive the stellar mass density. The stellar mass density in luminous galaxies has increased by a factor of 3.5-7.9 from z=3 to z=0.1, including field-to-field variance uncertainties. After correcting to total, the measured mass densities at z<2 lie below the integral of the star formation rate (SFR) density as a function of redshift as derived from UV selected samples. This may indicate a systematic error in the mass densities or SFR(z) estimates. We find large discrepancies between recent model predictions for the evolution of the mass density and our results, even when our observational selection is applied to the models. Finally we determine that Distant Red Galaxies (selected to have J_s - K_s>2.3) in our LV selected samples contribute 30% and 64% of the stellar mass budget at z~2 and z~ 2.8 respectively. These galaxies are largely absent from UV surveys and this result highlights the need for mass selection of high redshift galaxies.

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)

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 study the physical and photometric properties of galaxies at z=4 in cosmological hydrodynamic simulations of a lambda-CDM universe. We focus on galaxies satisfying the GOODS "B-dropout" criteria. Our goals are: (1) to study the nature of high-redshift galaxies; (2) to test the simulations against published measurements of high-redshift galaxies; (3) to find relations between photometric measurements by HST/ACS (0.4 -- 1 micron) and Spitzer/IRAC (3.6 -- 8 micron) and the intrinsic physical properties of GOODS "B-dropouts" such as stellar mass, stellar age, dust reddening, and star-formation rate; and (4) to assess how representative the GOODS survey is at this epoch. Our simulations predict that high-redshift galaxies show strong correlations in star formation rate versus stellar mass, and weaker correlations versus environment and age, such that GOODS galaxies are predicted to be the most massive, most rapidly star-forming galaxies at that epoch, living preferentially in dense regions. The simulated rest-frame UV luminosity function (LF) and integrated luminosity density are in broad agreement with observations at z~4. The predicted rest-frame optical (observed 3.6 micron) LF is similar to the rest-frame UV LF, shifted roughly one magnitude brighter. We predict that GOODS detects less than 50% of the total stellar mass density formed in galaxies more massive than 10^8.7 M_sun by z=4, mainly because of brightness limits in the HST/ACS bands. The most rapidly star forming galaxies in our simulations have rates exceeding 1000 M_sun yr^-1, similar to observed sub-mm galaxies. The star formation rates of these galaxies show at most a mild excess (2--3x) over the rates that would be expected for their stellar mass. Whether these bright galaxies would be observable as LBGs depends on the uncertain effects of dust reddening.

We positionally match sources observed by the Sloan Digital Sky Survey (SDSS), the Two Micron All Sky Survey (2MASS), and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. Practically all 2MASS sources are matched to an SDSS source within 2 arcsec; ~11% of them are optically resolved galaxies and the rest are dominated by stars. About 1/3 of FIRST sources are matched to an SDSS source within 2 arcsec; ~80% of these are galaxies and the rest are dominated by quasars. Based on these results, we project that by the completion of these surveys the matched samples will include about 10^7 stars and 10^6 galaxies observed by both SDSS and 2MASS, and about 250,000 galaxies and 50,000 quasars observed by both SDSS and FIRST. Here we present a preliminary analysis of the optical, infrared and radio properties for the extragalactic sources from the matched samples. In particular, we find that the fraction of quasars with stellar colors missed by the SDSS spectroscopic survey is probably not larger than ~10%, and that the optical colors of radio-loud quasars are ~0.05 mag. redder (with 4-sigma significance) than the colors of radio-quiet quasars.

We discuss optical and infrared photometric properties of stars matched in the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) commissioning data for 50 deg2 of sky. About 98% (63,000) of objects listed in the 2MASS Point Source Catalog in the analyzed area are matched within 2 arcsec to an SDSS source. The matched sources represent 8% of the about 800,000 SDSS sources in this area, and 15% of them are resolved in SDSS imaging data, although they are detected as point sources in 2MASS data. For about 14,000 stars with the smallest photometric errors (< 10%) in both surveys, we present optical and infrared color-magnitude and color-color diagrams. We use optical (SDSS) colors to identify the stellar spectral sequence and show that stars of different spectral types can have similar infrared colors, thus making the classification of stars based on only 2MASS data very difficult. However, a broad separation into ``early'' and ``late'' spectral types (relative to type K0) is possible with a reliability of about 95% even with 2MASS colors alone. The distributions of matched sources in color-magnitude and color-color diagrams are compared to the predictions of a stellar population synthesis code. We find that the models are in fair overall agreement with the data. The most significant discrepancies are found for the number ratio of ``early'' to ``late'' type stars (by about a factor of 2) and in the colors of M stars (up to 0.2 mag).

We present a sample of 148 candidate RR Lyrae stars selected from Sloan Digital Sky Survey (SDSS) commissioning data for about 100 deg2 of sky surveyed twice with Delta(t) = 1.9946 days. Although the faint magnitude limit of the SDSS allows us to detect RR Lyrae stars to large galactocentric distances (about 100 kpc, or r = 21), we find no candidates fainter than r = 20, i.e. further than about 65 kpc from the Galactic center. On the assumption that all 148 candidates are indeed RR Lyrae stars (contamination by other species of variable star is probably less than 10%), we find that their volume density has roughly a power-law dependence on galactocentric radius, 1/R^(2.7+-0.2), between 10 and 50 kpc, and drops abruptly at R = 50-60 kpc, possibly indicating a sharp edge to the stellar halo as traced by RR Lyrae stars. The Galactic distribution of stars in this sample is very inhomogeneous and shows a clump of over 70 stars at about 45 kpc from the Galactic center. This clump is also detected in the distribution of nonvariable objects with RR Lyrae star colors. When sources in the clump are excluded, the best power-law fit becomes consistent with the 1/R^3 distribution found from surveys of bright RR Lyrae stars. These results imply that the halo contains clumpy overdensities inhomogeneously distributed within a smooth 1/R^3 background, with a possible cutoff at about 50 kpc.