Source author record

Kristian Finlator

Kristian Finlator appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

ALMA reveals extended cool gas and hot ionized outflows in a typical star-forming galaxy at $z=7.13$

We present spatially-resolved morphological properties of [CII] 158 $μ$m, [OIII] 88 $μ$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is extended up to a radius of $r \sim 12$ kpc. Using multi-band rest-frame far-infrared (FIR) continuum data ranging from 52-400 $μ$m, we find an average dust temperature and emissivity index of $T_{\rm dust} = 41^{+17}_{-14}$ K and $β= 1.7^{+1.1}_{-0.7}$, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star-formation rate (SFR) via IR and UV luminosities and determine a total SFR of $37\pm 1~M_\odot~{\rm yr}^{-1}$ with an obscured fraction of $87\%$. While the [OIII] line is a good tracer of the SFR, the [CII] line shows deviation from the local $L_{\rm [CII]}$-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [CII] and [OIII], with significant residuals ($\sim 5σ$) in the [OIII] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [CII] structure from the cooling of hot ionized outflows. The extended [CII] and high-velocity [OIII] emission may both contribute in part to the high $L_{\rm [OIII]}$/$L_{\rm [CII]}$ ratios recently reported in $z>6$ galaxies.

preprint2022arXiv

Assuming Ionization Equilibrium and the Impact on the Lyman-$\mathrmα$ Forest Power Spectrum during the End of Reionization at $8 \geq z \geq 5$

We explore how the assumption of ionization equilibrium modulates the modeled intergalactic medium (IGM) at the end of the hydrogen Epoch of Reionization using the cosmological radiation hydrodynamic \textsc{Technicolor Dawn} simulation. In neutral and partially-ionized regions where the metagalactic ultraviolet background (UVB) is weak, the ionization timescale $t_\mathrm{ion}\equiv Γ^{-1}$ exceeds the Hubble time. Assuming photoionization equilibrium in such regions artificially boosts the ionization rate, accelerating reionization. By contrast, the recombination time $t_\mathrm{rec} < t_\mathrm{ion}$ in photoionized regions, with the result that assuming photoionization equilibrium artificially increases the neutral hydrogen fraction. Using snapshots between $8 \geq z \geq 5$, we compare the predicted Lyman-$α$ forest flux power spectrum with and without the assumption of ionization equilibrium. Small scales ($k > 0.1$ rad s km$^{-1}$) exhibit reduced power from $7 \leq z \leq 5.5$ in the ionization equilibrium case while larger scales are unaffected. This occurs for the same reasons: ionization equilibrium artificially suppresses the neutral fraction in self-shielded gas and boosts ionizations in voids, suppressing small-scale fluctuations in the ionization field. When the volume-averaged neutral fraction drops below $10^{-4}$, the signature of non-equilibrium ionizations on the Lyman-$α$ forest (LAF) disappears. Comparing with recent observations indicates that these non-equilibrium effects are not yet observable in the LAF flux power spectrum.

preprint2022arXiv

SED Analysis of 13 Spectroscopically Confirmed Galaxies at z$\simeq$6 to Constrain UV-Slope, Model Dust Attenuation and Escape Fractions

The reionization of the hydrogen in the Universe is thought to have completed by redshift $z\simeq5.5-6$. To probe this era, galaxy observations in the Subaru Deep Field (SDF) have identified more than 100 galaxies at $z\simeq6$, many spectroscopically confirmed through follow-up observations. We model the spectral energy distributions (SEDs) of 13 SDF galaxies with the CIGALE and Dense Basis codes using available optical/IR data. Modeling deep IR photometry has the potential to constrain the galaxy's Lyman continuum (LyC) escape fraction (\fesc). We use the modeled nebular emission lines and find that the implied escape fractions ranges from \textbf{0 to 0.8 with a median of $\sim$0.35 for Dense Basis and $\sim$0.55 for CIGALE.} Significant uncertainties in the data exist, so that fitting results in a large range of \fesc\ for individual objects. The implied median \fesc-values may be high enough for galaxies to finish reionization by $z\sim6$. Furthermore, we find no strong trends between the UV-slope $β$ or \EBminV with model \fesc. If true, the lack of trends suggest that other factors besides nebular emission or dust extinction could have led to LyC escaping, such as the presence of holes in the ISM with sufficiently wide opening angles from outflows of supernovae and/or weak AGN, resulting in a range of implied \fesc-values depending on the viewing angle of each galaxy. The current \textit{HST, Spitzer} and ground-based photometric and model errors for the galaxies remain large, so IR spectroscopic observations with the \textit{James Webb Space Telescope} are needed to constrain this possibility.

preprint2021arXiv

CLEAR: Boosted Ly$α$ Transmission of the Intergalactic Medium in UV bright Galaxies

Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalactic medium (IGM) around overdense regions of galaxies. Recent Lyman-alpha (Ly$α$) studies during the epoch of reionization show growing evidence that ionized bubbles formed earlier around brighter galaxies, suggesting higher IGM transmission of Ly$α$ from these galaxies. We investigate this problem using IR slitless spectroscopy from the Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) G102 grism observations of 148 galaxies selected via photometric redshifts at $6.0<z<8.2$. These galaxies have spectra extracted from the CANDELS Ly$α$ Emission at Reionization (CLEAR) survey. We combine the CLEAR data for 275 galaxies with the Keck/DEIMOS+MOSFIRE dataset from the Texas Spectroscopic Search for Ly$α$ Emission at the End of Reionization Survey. We then constrain the Ly$α$ equivalent-width (EW) distribution at $6.0<z<8.2$, which is described by an exponential form, $dN/d\text{EW}\propto\text{exp(-EW)}/W_0$, with the characteristic $e$-folding scale width ($W_0$). We confirm a significant drop of the Ly$α$ strength (or $W_0$) at $z>6$. Furthermore, we compare the redshift evolution of $W_0$ between galaxies at different UV luminosities. The UV-bright ($M_{\text{UV}}<-21$, or $L_{\text{UV}}>L^{*}$) galaxies show weaker evolution with a decrease of 0.4 ($\pm$0.2) dex in $W_0$ at $z>6$ while UV-faint ($M_{\text{UV}}>-21$, or $L_{\text{UV}}<L^{*}$) galaxies exhibit a significant drop by a factor of 0.7-0.8 ($\pm0.2$) dex in $W_0$ from $z<6$ to $z>6$. Our results add to the accumulating evidence that UV-bright galaxies exhibit boosted Ly$α$ transmission in the IGM, suggesting that reionization completes sooner in regions proximate to galaxies of higher UV luminosity.

preprint2020arXiv

ARTIST: Fast radiative transfer for large-scale simulations of the epoch of reionisation

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.

preprint2020arXiv

Luminous Lyman-alpha Emitters with Very Blue UV-continuum Slopes at Redshift 5.7 <= z <= 6.6

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 $β$ are in a range of $-3.4\le β\le -2.6$. Unlike previous, tentative detections of $β\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 $β\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 $β\simeq-2.7$ when nebular emission is considered. This can nearly explain the SEDs of two galaxies with $β\simeq-2.8$ and --2.9 ($σ_β=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 $β\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 $β\simeq -3.4 \pm0.4$ galaxies. Future deep spectroscopic observations are needed to fully understand these galaxies.

preprint2020arXiv

Testing Galaxy Formation Simulations with Damped Lyman-$α$ Abundance and Metallicity Evolution

We examine the properties of damped Lyman-$α$ 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-$α$ 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, $Ω_{\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.

preprint2020arXiv

The Faint Host Galaxies of C IV Absorbers at z > 5

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$α$ 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}α}$/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$α$ 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$α$ emission using integral field units.

preprint2020arXiv

Weak Evolution of the Mass-metallicity Relation at Cosmic Dawn in the FirstLight Simulations

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.

preprint2016arXiv

Equilibrium model prediction for the scatter in the star-forming main sequence

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.

preprint2016arXiv

Simulating the 21-cm signal from reionisation including non-linear ionisations and inhomogeneous recombinations

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.

preprint2015arXiv

A Critical Assessment of Stellar Mass Measurement Methods

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, $Δ\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.

preprint2015arXiv

Equilibrium Model Constraints on Baryon Cycling Across Cosmic Time

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.

preprint2015arXiv

Merging galaxies produce outliers from the Fundamental Metallicity Relation

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 $τ=1.568_{-0.027}^{+0.029}$ Gyr, the average metallicity depression caused by a 1:1 merger to be $α=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 $ξ_\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.

preprint2015arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at $z\ge6$. III. Stellar Populations from SED Modeling with Secure Ly$α$ Emission and Redshifts

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$α$ 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$α$ 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$α$ emission are indistinguishable in terms of age, stellar mass, and SFR.

preprint2015arXiv

The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS

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, $σ(\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.

preprint2013arXiv

A Comparison between Semi-Analytic Model Predictions for the CANDELS Survey

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

preprint2013arXiv

A Critical Assessment of Photometric Redshift Methods: A CANDELS Investigation

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.

preprint2013arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. I. Basic Characteristics of the Rest-Frame UV Continuum and Lyman-alpha emission

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 Å. 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 Å) that could be the dominant contribution to the total ionizing flux at z >= 6.

preprint2013arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. II. Morphology of the Rest-Frame UV Continuum and Lyman-alpha Emission

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.

preprint2013arXiv

The metallicity evolution of low mass galaxies: New constraints at intermediate redshift

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.

preprint2012arXiv

A fundamental problem in our understanding of low mass galaxy evolution

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]

preprint2012arXiv

CANDELS: The Contribution of the Observed Galaxy Population to Cosmic Reionization

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

preprint2012arXiv

CANDELS: The Evolution of Galaxy Rest-Frame Ultraviolet Colors from z = 8 to 4

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.

preprint2012arXiv

Constraints on the Ionizing Efficiency of the First Galaxies

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.

preprint2012arXiv

Recent Advances in Cosmological Hydrogen Reionization

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.

preprint2011arXiv

An Analytic Model for the Evolution of the Stellar, Gas, and Metal Content of Galaxies

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.

preprint2011arXiv

Galaxy Evolution in Cosmological Simulations With Outflows I: Stellar Masses and Star Formation Rates

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.

preprint2011arXiv

Galaxy Evolution in Cosmological Simulations with Outflows II: Metallicities and Gas Fractions

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.

preprint2011arXiv

Metals Removed by Outflows from Milky Way Dwarf Spheroidal Galaxies

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.

preprint2010arXiv

Smoothly-Rising Star Formation Histories During the Reionization Epoch

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)

preprint2010arXiv

The Nature of Sub-millimetre Galaxies in Cosmological Hydrodynamic Simulations

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.

preprint2009arXiv

Tracing the Reionization-Epoch Intergalactic Medium with Metal Absorption Lines

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.