We derive stellar masses from SED fitting to rest-frame optical and UV fluxes for 401 star-forming galaxies at z 4, 5, and 6 from Hubble-WFC3/IR observations of the ERS combined with the deep GOODS-S Spitzer/IRAC data (and include a previously-publis
hed z 7 sample). A mass-luminosity relation with strongly luminosity-dependent M/Luv ratios is found for the largest sample (299 galaxies) at z 4. The relation M propto L_{UV,1500}^(1.7+/-0.2) has a well-determined intrinsic sample variance of 0.5 dex. This relation is also consistent with the more limited samples at z 5-7. This z 4 mass-luminosity relation, and the well-established faint UV luminosity functions at z 4-7, are used to derive galaxy mass functions (MF) to masses M~10^8 at z 4-7. A bootstap approach is used to derive the MFs to account for the large scatter in the M--Luv relation and the luminosity function uncertainties, along with an analytical crosscheck. The MFs are also corrected for the effects of incompleteness. The incompleteness-corrected MFs are steeper than previously found, with slopes alpha_M-1.4 to -1.6 at low masses. These slopes are, however, still substantially flatter than the MFs obtained from recent hydrodynamical simulations. We use these MFs to estimate the stellar mass density (SMD) of the universe to a fixed M_{UV,AB}<-18 as a function of redshift and find a SMD growth propto(1+z)^{-3.4 +/-0.8} from z 7 to z 4. We also derive the SMD from the completeness-corrected MFs to a mass limit M~10^{8} Msun. Such completeness-corrected MFs and the derived SMDs will be particularly important for model comparisons as future MFs reach to lower masses.
We use ~88 arcmin**2 of deep (>~26.5 mag at 5 sigma) NICMOS data over the two GOODS fields and the HDF South to conduct a search for bright z>~7 galaxy candidates. This search takes advantage of an efficient preselection over 58 arcmin**2 of NICMOS H
-band data where only plausible z>~7 candidates are followed up with NICMOS J-band observations. ~248 arcmin**2 of deep ground-based near-infrared data (>~25.5 mag, 5 sigma) is also considered in the search. In total, we report 15 z-dropout candidates over this area -- 7 of which are new to these search fields. Two possible z~9 J-dropout candidates are also found, but seem unlikely to correspond to z~9 galaxies. The present z~9 search is used to set upper limits on the prevalence of such sources. Rigorous testing is undertaken to establish the level of contamination of our selections by photometric scatter, low mass stars, supernovae (SNe), and spurious sources. The estimated contamination rate of our z~7 selection is ~24%. Through careful simulations, the effective volume available to our z>~7 selections is estimated and used to establish constraints on the volume density of luminous (L*(z=3), or -21 mag) galaxies from these searches. We find that the volume density of luminous star-forming galaxies at z~7 is 13_{-5}^{+8}x lower than at z~4 and >25x lower (1 sigma) at z~9 than at z~4. This is the most stringent constraint yet available on the volume density of >~L* galaxies at z~9. The present wide-area, multi-field search limits cosmic variance to <20%. The evolution we find at the bright end of the UV LF is similar to that found from recent Subaru Suprime-Cam, HAWK-I or ERS WFC3/IR searches. The present paper also includes a complete summary of our final z~7 z-dropout sample (18 candidates) identified from all NICMOS observations to date (over the two GOODS fields, the HUDF, galaxy clusters).
We use a robust sample of 11 z~7 galaxies (z-dropouts) to estimate the stellar mass density of the universe when it was only ~750 Myr old. We combine the very deep optical to near-Infrared photometry from the HST ACS and NICMOS cameras with mid-Infra
red Spitzer IRAC imaging available through the GOODS program. After carefully removing the flux from contaminating foreground sources we have obtained reliable photometry in the 3.6 and 4.5 micron IRAC channels. The spectral shapes of these sources, including their rest frame optical colors, strongly support their being at z~7 with a mean photometric redshift of <z>=7.2+/-0.5. We use Bruzual & Charlot (2003) synthetic stellar population models to constrain their stellar masses and star formation histories. We find stellar masses that range over 0.1 -12x10^9 M_sol and average ages from 20 Myr to up to 425 Myr with a mean of ~300 Myr, suggesting that in some of these galaxies most of the stars were formed at z>8 (and probably at z>~10). The best fits to the observed SEDs are consistent with little or no dust extinction, in agreement with recent results at z~4-8. The star formation rates (SFR) are in the range from 5-20 M_sol/yr. From this sample we measure a stellar mass density of 6.6_{-3.3}^{+5.4}x10^5 M_sol/Mpc^3 to a limit of M_{UV,AB}<-20 (or 0.4 L*(z=3)). Combined with a fiducial lower limit for their ages (80 Myr) this implies a maximum SFR density of 0.008 M_sol/yr/Mpc^3. This is well below the critical level needed to reionize the universe at z~8 using standard assumptions. However, this result is based on luminous sources (>L*) and does not include the dominant contribution of the fainter galaxies. Strikingly, we find that the specific SFR is constant from z~7 to z~2 but drops substantially at more recent times.
