We study cosmic variance in deep high redshift surveys and its influence on the determination of the luminosity function for high redshift galaxies. For several survey geometries relevant for HST and JWST instruments, we characterize the distribution
of the galaxy number counts. This is obtained by means of analytic estimates via the two point correlation function in extended Press-Schechter theory as well as by using synthetic catalogs extracted from N-body cosmological simulations of structure formation. We adopt a simple luminosity - dark halo mass relation to investigate the environment effects on the fitting of the luminosity function. We show that in addition to variations of the normalization of the luminosity function, a steepening of its slope is also expected in underdense fields, similarly to what is observed within voids in the local universe. Therefore, to avoid introducing artificial biases, caution must be taken when attempting to correct for field underdensity, such as in the case of HST UDF i-dropout sample, which exhibits a deficit of bright counts with respect to the average counts in GOODS. A public version of the cosmic variance calculator based on the two point correlation function integration is made available on the web.
Lyman-Break Galaxy (LBG) samples observed during reionization ($zgtrsim6$) with Hubble Space Telescopes Wide Field Camera 3 are reaching sizes sufficient to characterize their clustering properties. Using a combined catalog from the Hubble eXtreme De
ep Field and CANDELS surveys, containing $N=743$ LBG candidates at z>6.5 at a mean redshift of $z=7.2$, we detect a clear clustering signal in the angular correlation function (ACF) at $sim4sigma$, corresponding to a real-space correlation length $r_{0}=6.7^{+0.9}_{-1.0}h^{-1}$cMpc. The derived galaxy bias $b=8.6^{+0.9}_{-1.0}$ is that of dark-matter halos of $M=10^{11.1^{+0.2}_{-0.3}}$M$_{odot}$ at $z=7.2$, and highlights that galaxies below the current detection limit ($M_{AB}sim-17.7$) are expected in lower-mass halos ($Msim10^{8}-10^{10.5}$M$_{odot}$). We compute the ACF of LBGs at $zsim3.8-zsim5.9$ in the same surveys. A trend of increasing bias is found from $z=3.8$ ($bsim3.0$) to $z=7.2$ ($bsim8.6$), broadly consistent with galaxies at fixed luminosity being hosted in dark-matter halos of similar mass at $4<z<6$, followed by a slight rise in halo masses at $zsim7$ ($sim2sigma$ confidence). Separating the data at the median luminosity of the $z=7.2$ sample ($M_{UV}=-19.4$) shows higher clustering at $z=5.9$ for bright galaxies ($r_{0}=5.5^{+1.4}_{-1.5}h^{-1}$cMpc, $b=6.2^{+1.2}_{-1.5}$) compared to faint galaxies ($r_{0}=1.9^{+1.1}_{-1.0}h^{-1}$cMpc, $b=2.7pm1.2$) implying a constant mass-to-light ratio $frac{dlogM}{dlogL}sim1.2^{+1.8}_{-0.8}$. A similar trend is present in the $z=7.2$ sample with larger uncertainty. Finally, our bias measurements allow us to investigate the fraction of dark-matter halos hosting UV-bright galaxies (the duty-cycle, $epsilon_{DC}$). At $z=7.2$ values near unity are preferred, which may be explained by the shortened halo assembly time at high-redshift.
We present a tally of Milky Way late-type dwarf stars in 68 WFC3 pure-parallel fields (227 arcmin^2) from the Brightest of Reionizing Galaxies (BoRG) survey for high-redshift galaxies. Using spectroscopically identified M-dwarfs in two public surveys
, the CANDELS and the ERS mosaics, we identify a morphological selection criterion using the half-light radius (r50), a near-infrared J-H, G-J color region where M-dwarfs are found, and a V-J relation with M-dwarf subtype. We apply this morphological selection of stellar objects, color-color selection of M-dwarfs and optical-near-infrared color subtyping to compile a catalog of 274 M-dwarfs belonging to the disk of the Milky Way with a limiting magnitude of m_F125W < 24. Based on the M-dwarfs statistics, we conclude that (a) the previously identified North/South discrepancy in M-dwarf numbers persists in our sample; there are more M-dwarfs in the Northern fields on average than in Southern ones, (b) the Milky Ways single disk scale-height for M-dwarfs is 0.3-4 kpc, depending on sub-type, (c) {bf ERRATUM:} we present corrected coordinates (AstroPy) and distances and find a constant $z_0$=600 pc for all types. (d) a second component is visible in the vertical distribution, with a different, much higher scale-height. We report the M-dwarf component of the Sagittarius stream in one of our fields with 11 confirmed M-dwarfs, 7 of which are at the streams distance. The dwarf scale-height and the relative low incidence in our fields of L- and T-dwarfs in these fields makes it unlikely that these stars will be interlopers in great numbers in color-selected samples of high-redshift galaxies. The relative ubiquity of M-dwarfs however will make them ideal tracers of Galactic Halo substructure with EUCLID and reference stars for JWST observations.
Gamma Ray Bursts (GRBs) and galaxies at high redshift represent complementary probes of the star formation history of the Universe. In fact, both the GRB rate and the galaxy luminosity density are connected to the underlying star formation. Here, we
combine a star formation model for the evolution of the galaxy luminosity function from z=0 to z=10 with a metallicity-dependent efficiency for GRB formation to simultaneously predict the comoving GRB rate. Our model sheds light on the physical origin of the empirical relation often assumed between GRB rate and luminosity density-derived star formation rate: Rgrb(z) = epsilon(z)*SFR_{obs}(z), with epsilon(z) (1+z)^{1.2}. At z<4, epsilon(z) is dominated by the effects of metallicity evolution in the GRB efficiency. Our best-fitting model only requires a moderate preference for low-metallicity, that is a GRB rate per unit stellar mass about four times higher for log(Z/Zsun)<-3 compared to log(Z/Zsun)>0. Models with total suppression of GRB formation at log(Z/Zsun)>0 are disfavored. At z>4, most of the star formation happens in low-metallicity hosts with nearly saturated efficiency of GRB production per unit stellar mass. However at the same epoch, galaxy surveys miss an increasing fraction of the predicted luminosity density because of flux limits, driving an accelerated evolution of epsilon(z) compared to the empirical power-law fit from lower z. Our findings are consistent with the non-detections of GRB hosts in ultradeep imaging at z>5, and point toward current galaxy surveys at z>8 only observing the top 15-20 % of the total luminosity density.
It is widely believed that globular clusters evolve over many two-body relaxation times toward a state of energy equipartition, so that velocity dispersion scales with stellar mass as sigma ~ m^{-eta} with eta = 0.5. We show that this is incorrect, u
sing direct N-body simulations with a variety of realistic IMFs and initial conditions. No simulated system ever reaches a state close to equipartition. Near the center, the luminous main-sequence stars reach a maximum eta_{max} ~ 0.15 pm 0.03. At large times, all radial bins convergence on an asymptotic value eta_{infty} ~ 0.08 pm 0.02. The development of this partial equipartition is strikingly similar across our simulations, despite the range of initial conditions employed. Compact remnants tend to have higher eta than main-sequence stars (but still eta < 0.5), due to their steeper (evolved) mass function. The presence of an intermediate-mass black hole (IMBH) decreases eta, consistent with our previous findings of a quenching of mass segregation under these conditions. All these results can be understood as a consequence of the Spitzer instability for two-component systems, extended by Vishniac to a continuous mass spectrum. Mass segregation (the tendency of heavier stars to sink toward the core) has often been studied observationally, but energy equipartition has not. Due to the advent of high-quality proper motion datasets from the Hubble Space Telescope, it is now possible to measure eta. Detailed data-model comparisons open up a new observational window on globular cluster dynamics and evolution. Comparison of our simulations to Omega Cen observations yields good agreement, confirming that globular clusters are not generally in energy equipartition. Modeling techniques that assume equipartition by construction (e.g., multi-mass Michie-King models) are approximate at best.
Until now, investigating the early stages of galaxy formation has been primarily the realm of theoretical modeling and computer simulations, which require many physical ingredients and are challenging to test observationally. However, the latest Hubb
le Space Telescope observations in the near infrared are shedding new light on the properties of galaxies within the first billion years after the Big Bang, including our recent discovery of the most distant proto-cluster of galaxies at redshift z~8. Here, I compare predictions from models of primordial and metal-enriched star formation during the dark ages with the latest Hubble observations of galaxies during the epoch of reionization. I focus in particular on the luminosity function and on galaxy clustering as measured from our Hubble Space Telescope Brightest of Reionizing Galaxies (BoRG) survey. BoRG has the largest area coverage to find luminous and rare z~8 sources that are among the first galaxies to have formed in the Universe.
Small galaxies consisting entirely of population III (pop III) stars may form at high redshifts, and could constitute one of the best probes of such stars. Here, we explore the prospects of detecting gravitationally lensed pop III galaxies behind the
galaxy cluster J0717.5+3745 (J0717) with both the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). By projecting simulated catalogs of pop III galaxies at z~7-15 through the J0717 magnification maps, we estimate the lensed number counts as a function of flux detection threshold. We find that the ongoing HST survey CLASH, targeting a total of 25 galaxy clusters including J0717, potentially could detect a small number of pop III galaxies if ~1% of the baryons in these systems have been converted into pop III stars. Using JWST exposures of J0717, this limit can be pushed to ~0.1% of the baryons. Ultra-deep JWST observations of unlensed fields are predicted to do somewhat worse, but will be able to probe pop III galaxies with luminosities intermediate between those detectable in HST/CLASH and in JWST observations of J0717. We also explain how current measurements of the galaxy luminosity function at z=7-10 can be used to constrain pop III galaxy models with very high star formation efficiencies (~10% of the baryons converted into pop III stars).
Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest o
f Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_ABsim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright zsim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter zsim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_hsim(4-7)x10^11Msun (sim5sigma density peak) and is surrounded by several M_hsim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the zsim8 nature of the overdensity, discover new members, and measure their precise redshift.
We present the first results on the search for very bright (M_AB -21) galaxies at redshift z~8 from the Brightest of Reionizing Galaxies (BoRG) survey. BoRG is a Hubble Space Telescope Wide Field Camera 3 pure-parallel survey that is obtaining images
on random lines of sight at high Galactic latitudes in four filters (F606W, F098M, F125W, F160W), with integration times optimized to identify galaxies at z>7.5 as F098M-dropouts. We discuss here results from a search area of approximately 130 arcmin^2 over 23 BoRG fields, complemented by six other pure-parallel WFC3 fields with similar filters. This new search area is more than two times wider than previous WFC3 observations at z~8. We identify four F098M-dropout candidates with high statistical confidence (detected at greater than 8sigma confidence in F125W). These sources are among the brightest candidates currently known at z~8 and approximately ten times brighter than the z=8.56 galaxy UDFy-38135539. They thus represent ideal targets for spectroscopic followup observations and could potentially lead to a redshift record, as our color selection includes objects up to z~9. However, the expected contamination rate of our sample is about 30% higher than typical searches for dropout galaxies in legacy fields, such as the GOODS and HUDF, where deeper data and additional optical filters are available to reject contaminants.
Cosmological simulations of galaxy formation often rely on prescriptions for star formation and feedback that depend on halo properties such as halo mass, central over-density, and virial temperature. In this paper we address the convergence of indiv
idual halo properties, based on their number of particles N, focusing in particular on the mass of halos near the resolution limit of a simulation. While it has been established that the halo mass function is sampled on average down to N~30 particles, we show that individual halo properties exhibit significant scatter, and some systematic biases, as one approaches the resolution limit. We carry out a series of cosmological simulations using the Gadget2 and Enzo codes with N_p=64^3 to N_p=1024^3 total particles, keeping the same large-scale structure in the simulation box. We consider boxes from l_{box} = 8 Mpc/h to l_{box} = 512 Mpc/h to probe different halo masses and formation redshifts. We cross-identify dark matter halos in boxes at different resolutions and measure the scatter in their properties. The uncertainty in the mass of single halos depends on the number of particles (scaling approximately as N^{-1/3}), but the rarer the density peak, the more robust its identification. The virial radius of halos is very stable and can be measured without bias for halos with N>30. In contrast, the average density within a sphere containing 25% of the total halo mass is severely underestimated (by more than a factor 2) and the halo spin is moderately overestimated for N<100. If sub-grid physics is implemented upon a cosmological simulation, we recommend that rare halos (~3sigma peaks) be resolved with N>100 particles and common halos (~1sigma peaks) with N>400 particles to avoid excessive numerical noise and possible systematic biases in the results.
