Using new ultradeep Spitzer/IRAC photometry from the IRAC Ultradeep Field program (IUDF), we investigate the stellar populations of a sample of 63 Y-dropout galaxy candidates at z~8, only 650Myr after the Big Bang. The sources are selected from HST/A
CS+WFC3/IR data over the Hubble Ultra Deep Field (HUDF), two HUDF parallel fields, and wide area data over the CANDELS/GOODS-South. The new Spitzer/IRAC data increase the coverage at 3.6 micron and 4.5 micron to ~120h over the HUDF reaching depths of ~28 (AB,1 sigma). The improved depth and inclusion of brighter candidates result in direct >3 sigma IRAC detections of 20/63 sources, of which 11/63 are detected at > 5 sigma. The average [3.6]-[4.5] colors of IRAC detected galaxies at z~8 are markedly redder than those at z~7, observed only 130Myr later. The simplest explanation is that we witness strong rest-frame optical emission lines (in particular [OIII]4959,5007+Hbeta) moving through the IRAC bandpasses with redshift. Assuming that the average rest-frame spectrum is the same at both z~7 and z~8 we estimate a rest-frame equivalent width of W([OIII]4959,5007+Hbeta) = 670 (+260,-170) Angstrom contributing 0.56 (+0.16,-0.11) mag to the 4.5 micron filter at z~8. The corresponding W(Halpha) = 430 (+160,-110) Angstrom implies an average specific star formation rate of sSFR = 11 (+11,-5) Gyr^-1 and a stellar population age of 100 (+100,-50) Myr. Correcting the spectral energy distribution for the contribution of emission lines lowers the average best-fit stellar masses and mass-to-light ratios by x3, decreasing the integrated stellar mass density to rho*(z=8,MUV<-18)=0.6 (+0.4,-0.3) x 10^6 Msun Mpc^-3.
Using deep NIR VLT/ISAAC and optical HST/WFPC2 imaging in the fields of the HDFS and MS1054-03, we study the rest-frame UV-to-optical colors and magnitudes of galaxies to z~3. While there is no evidence for a red sequence at z~3, there does appear to
be a well-defined color-magnitude relation (CMR) for blue galaxies at all redshifts, with more luminous galaxies having redder U-V colors. The slope of the blue CMR is independent of redshift d(U-V)/dMV = -0.09 (0.01) and can be explained by a correlation of dust-reddening with luminosity. The average color at fixed luminosity reddens strongly Delta(U-V) = 0.75 from z~3 to z=0, much of which can be attributed to aging of the stars. The color scatter of the blue sequence is relatively small sigma(U-V) = 0.25 (0.03) and constant to z~3, but notably asymmetrical with a sharp blue ridge and a wing towards redder colors. We explore sets of star formation histories to study the constraints placed by the shape of the scatter at z=2-3. One particular set of models, episodic star formation, reproduces the detailed properties very well. For a two-state model with high and low star formation, the duty cycle is constrained to be > 40% and the contrast between the states must be a factor > 5 (or a scatter in log(SFR) of > 0.35 dex around the mean). However, episodic models do not explain the observed tail of very red galaxies, primarily Distant Red Galaxies (DRGs), which may have ceased star formation altogether or are more heavily obscured. Finally, the relative number density of red, luminous MV < -20.5 galaxies increases by a factor of ~ 6 from z = 2.7 to z = 0.5, as does their contribution to the total rest-frame V-band luminosity density. We are likely viewing the progressive formation of red, passively evolving galaxies.
