The Reionization of Carbon

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.Comment: 16 pages, 14 figures, accepted to MNRA

The impact of nebular emission on the broadband fluxes of high-redshift galaxies

A substantial fraction of the light emitted from young or star-forming galaxies at ultraviolet to near-infrared wavelengths comes from the ionized interstellar medium in the form of emission lines and a nebular continuum. At high redshifts, star formation rates are on average higher and stellar populations younger than in the local Universe. Both of these effects act to boost the impact of nebular emission on the overall spectrum of galaxies. Even so, the broadband fluxes and colours of high-redshift galaxies are routinely analyzed under the assumption that the light observed originates directly from stars. Here, we assess the impact of nebular emission on broadband fluxes in Johnson/Cousins BVRIJHK, Sloan Digital Sky Survey griz and Spitzer IRAC/MIPS filters as a function of observed redshift (up to z=15) for galaxies with different star formation histories. We find that nebular emission may account for a non-negligible fraction of the light received from high-redshift galaxies. The ages and masses inferred for such objects through the use of spectral evolutionary models that omit the nebular contribution are therefore likely to contain systematic errors. We argue that a careful treatment of the nebular component will be essential for the interpretation of the rest-frame ultraviolet-to-infrared properties of the first galaxies formed, like the ones expected to be detected with the James Webb Space Telescope.Comment: 5 pages, 3 figures, accepted for publication in ApJ

A Link to the Past: Using Markov Chain Monte Carlo Fitting to Constrain Fundamental Parameters of High-Redshift Galaxies

We have a developed a new method for fitting spectral energy distributions (SEDs) to identify and constrain the physical properties of high-redshift (4 < z < 8) galaxies. Our approach uses an implementation of Bayesian based Markov Chain Monte Carlo (PiMC^2) that allows us to compare observations to arbitrarily complex models and to compute 95% credible intervals that provide robust constraints for the model parameters. The work is presented in 2 sections. In the first, we test PiMC^2 using simulated SEDs to not only confirm the recovery of the known inputs but to assess the limitations of the method and identify potential hazards of SED fitting when applied specifically to high redshift (z>4) galaxies. Our tests reveal five critical results: 1) the ability to confidently constrain metallicity, population ages, and Av all require photometric accuracy better than what is currently achievable (i.e. less than a few percent); 2) the ability to confidently constrain stellar masses (within a factor of two) can be achieved without the need for high-precision photometry; 3) the addition of IRAC photometry does not guarantee that tighter constraints of the stellar masses and ages can be defined; 4) different assumptions about the star formation history can lead to significant biases in mass and age estimates; and 5) we are able to constrain stellar age and Av of objects that are both young and relatively dust free. In the second part of the paper we apply PiMC^2 to 17 4<z<8 objects, including the GRAPES Ly alpha sample (4<z<6), supplemented by HST/WFC3 near-IR observations, and several broad band selected z>6 galaxies. Using PiMC^2, we are able to constrain the stellar mass of these objects and in some cases their stellar age and find no evidence that any of these sources formed at a redshift much larger than z_f=8, a time when the Universe was ~ 0.6 Gyr old.Comment: Submitted to ApJ (Full abstract, 47 pages, 17 figures, 7 tables

A Link to the Past: Using Markov Chain Monte Carlo Fitting to Constrain Fundamental Parameters of High-Redshift Galaxies

We have a developed a new method for fitting spectral energy distributions (SEDs) to identify and constrain the physical properties of high-redshift (4 < z < 8) galaxies. Our approach uses an implementation of Bayesian based Markov Chain Monte Carlo (PiMC^2) that allows us to compare observations to arbitrarily complex models and to compute 95% credible intervals that provide robust constraints for the model parameters. The work is presented in 2 sections. In the first, we test PiMC^2 using simulated SEDs to not only confirm the recovery of the known inputs but to assess the limitations of the method and identify potential hazards of SED fitting when applied specifically to high redshift (z>4) galaxies. Our tests reveal five critical results: 1) the ability to confidently constrain metallicity, population ages, and Av all require photometric accuracy better than what is currently achievable (i.e. less than a few percent); 2) the ability to confidently constrain stellar masses (within a factor of two) can be achieved without the need for high-precision photometry; 3) the addition of IRAC photometry does not guarantee that tighter constraints of the stellar masses and ages can be defined; 4) different assumptions about the star formation history can lead to significant biases in mass and age estimates; and 5) we are able to constrain stellar age and Av of objects that are both young and relatively dust free. In the second part of the paper we apply PiMC^2 to 17 4<z<8 objects, including the GRAPES Ly alpha sample (4<z<6), supplemented by HST/WFC3 near-IR observations, and several broad band selected z>6 galaxies. Using PiMC^2, we are able to constrain the stellar mass of these objects and in some cases their stellar age and find no evidence that any of these sources formed at a redshift much larger than z_f=8, a time when the Universe was ~ 0.6 Gyr old.Comment: Submitted to ApJ (Full abstract, 47 pages, 17 figures, 7 tables

Prospects for CDM sub-halo detection using high angular resolution observations

In the CDM scenario, dark matter halos are assembled hierarchically from smaller subunits. A long-standing problem with this picture is that the number of sub-halos predicted by CDM simulations is orders of magnitudes higher than the known number of satellite galaxies in the vicinity of the Milky Way. A plausible way out of this problem could be that the majority of these sub-halos somehow have so far evaded detection. If such "dark galaxies" do indeed exist, gravitational lensing may offer one of the most promising ways to detect them. Dark matter sub-halos in the 1e6 - 1e10 solar mass range should cause strong gravitational lensing on (sub)milliarcsecond scales. We study the feasibility of a strong lensing detection of dark sub-halos by deriving the image separations expected for density profiles favoured by recent simulations and comparing these to the angular resolution of both existing and upcoming observational facilities. We find that there is a reasonable probability to detect sub-halo lensing effects in high resolution observations at radio wavelengths, such as produced by the upcoming VSOP-2 satellite, and thereby test the existence of dark galaxies.Comment: 9 pages, 5 figures, Proceedings for "The Universe under the Microscope" (AHAR 2008), held in Bad Honnef (Germany) in April 2008, to be published in Journal of Physics: Conference Series by Institute of Physics Publishing, R. Schoedel, A. Eckart, S. Pfalzner, and E. Ros (eds.

The properties of SN Ib/c locations

Aims: To gain better insight on the physics of stripped-envelope core-collapse supernovae through studying their environments. Methods: We obtained low-resolution optical spectroscopy with the New Technology Telescope (+ EFOSC2) at the locations of 20 Type Ib/c supernovae. We measure the flux of emission lines in the stellar-continuum-subtracted spectra from which local metallicities are computed. For the supernova regions we estimate both the mean stellar age, interpreting the stellar absorption with population synthesis models, and the age of the youngest stellar populations using the H-alpha equivalent width as an age indicator. These estimates are compared with the lifetimes of single massive stars. Results: Based on our sample, we detect a tentative indication that Type Ic supernovae might explode in environments that are more metal-rich than those of Type Ib supernovae (average difference of 0.08 dex), but this is not a statistically significant result. The lower limits placed on the ages of the supernova birthplaces are overall young, although there are several cases where these appear older than what is expected for the evolution of single stars more massive than 25-30 M_{sun}. This is only true, however, assuming that the supernova progenitors were born during an instantaneous (not continuous) episode of star formation. Conclusions: These results do not conclusively favor any of the two evolutionary paths (single or binary) leading to stripped supernovae. We do note a fraction of events for which binary evolution is more likely, due to their associated age limits. The fact, however, that the supernova environments contain areas of recent (< 15 Myr) star formation and that the environmental metallicities are, at least, not against the single evolutionary scenario, suggest that this channel is also broadly consistent with the observations.Comment: Matches published version (after proofs

21-cm cosmology

Imaging the Universe during the first hundreds of millions of years remains one of the exciting challenges facing modern cosmology. Observations of the redshifted 21 cm line of atomic hydrogen offer the potential of opening a new window into this epoch. This would transform our understanding of the formation of the first stars and galaxies and of the thermal history of the Universe. A new generation of radio telescopes is being constructed for this purpose with the first results starting to trickle in. In this review, we detail the physics that governs the 21 cm signal and describe what might be learnt from upcoming observations. We also generalize our discussion to intensity mapping of other atomic and molecular lines.Comment: 64 pages, 20 figures, submitted to Reports on Progress in Physics, comments welcom

The HDUV survey: six Lyman continuum emitter candidates at z ∼ 2 revealed by HST UV imaging

We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$\beta-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope.Comment: 11 pages, 6 figures, submitted to Ap

Dark Stars and Boosted Dark Matter Annihilation Rates

Dark Stars (DS) may constitute the first phase of stellar evolution, powered by dark matter (DM) annihilation. We will investigate here the properties of DS assuming the DM particle has the required properties to explain the excess positron and elec- tron signals in the cosmic rays detected by the PAMELA and FERMI satellites. Any possible DM interpretation of these signals requires exotic DM candidates, with an- nihilation cross sections a few orders of magnitude higher than the canonical value required for correct thermal relic abundance for Weakly Interacting Dark Matter can- didates; additionally in most models the annihilation must be preferentially to lep- tons. Secondly, we study the dependence of DS properties on the concentration pa- rameter of the initial DM density profile of the halos where the first stars are formed. We restrict our study to the DM in the star due to simple (vs. extended) adiabatic contraction and minimal (vs. extended) capture; this simple study is sufficient to illustrate dependence on the cross section and concentration parameter. Our basic results are that the final stellar properties, once the star enters the main sequence, are always roughly the same, regardless of the value of boosted annihilation or concentration parameter in the range between c=2 and c=5: stellar mass ~ 1000M\odot, luminosity ~ 10^7 L\odot, lifetime ~ 10^6 yrs (for the minimal DM models considered here; additional DM would lead to more massive dark stars). However, the lifetime, final mass, and final luminosity of the DS show some dependence on boost factor and concentration parameter as discussed in the paper.Comment: 37 pages, 11 figure