Simulating Galaxy Evolution

The forwards approach to galaxy formation and evolution is extremely powerful but leaves several questions unanswered. Foremost among these is the origin of disks. A backwards approach is able to provide a more realistic treatment of star formation and feedback and provides a practical guide to eventually complement galaxy formation ab initio.Comment: 11 pages with 2 figures, to appear in "After the Dark Ages: When Galaxies were Young", proceedings of the 9th annual October Astrophysics Conference, ed. S. Holt and E. Smith, simulated images available at http://astro.berkeley.edu/~bouwens/simulation.htm

Expanded Search for z~10 Galaxies from HUDF09, ERS, and CANDELS Data: Evidence for Accelerated Evolution at z>8?

We search for z~10 galaxies over ~160 arcmin^2 of WFC3/IR data in the Chandra Deep Field South, using the public HUDF09, ERS, and CANDELS surveys, that reach to 5sigma depths ranging from 26.9 to 29.4 in H_160 AB mag. z>~9.5 galaxy candidates are identified via J_125-H_160>1.2 colors and non-detections in any band blueward of J_125. Spitzer IRAC photometry is key for separating the genuine high-z candidates from intermediate redshift (z~2-4) galaxies with evolved or heavily dust obscured stellar populations. After removing 16 sources of intermediate brightness (H_160~24-26 mag) with strong IRAC detections, we only find one plausible z~10 galaxy candidate in the whole data set, previously reported in Bouwens et al. (2011). The newer data cover a 3x larger area and provide much stronger constraints on the evolution of the UV luminosity function (LF). If the evolution of the z~4-8 LFs is extrapolated to z~10, six z~10 galaxies are expected in our data. The detection of only one source suggests that the UV LF evolves at an accelerated rate before z~8. The luminosity density is found to increase by more than an order of magnitude in only 170 Myr from z~10 to z~8. This increase is >=4x larger than expected from the lower redshift extrapolation of the UV LF. We are thus likely witnessing the first rapid build-up of galaxies in the heart of cosmic reionization. Future deep HST WFC3/IR data, reaching to well beyond 29 mag, can enable a more robust quantification of the accelerated evolution around z~10.Comment: 13 pages, 11 figures, ApJ resubmitted after referee repor

Inside-Out Infall Formation of Disk Galaxies: Do Predictions Differ from Models without Size Evolution?

We develop an idealized inside-out formation model for disk galaxies to include a realistic mix of galaxy types and luminosities that provides a fair match to the traditional observables. The predictions of our infall models are compared against identical models with no-size evolution by generating fully realistic simulations of the HDF, from which we recover the angular size distributions. We find that our infall models produce nearly identical angular size distributions to those of our no-size evolution models in the case of a Omega = 0 geometry but produce slightly smaller sizes in the case of a Omega = 1 geometry, a difference we associate with the fact that there is a different amount of cosmic time in our two models for evolving to relatively low redshifts (z \approx 1-2). Our infall models also predict a slightly smaller (11% - 29%) number of large (disk scale lengths > 4 h_{50} ^{-1} kpc) galaxies at z \approx 0.7 for the CFRS as well as different increases in the central surface brightness of the disks for early-type spirals, the infall model predicting an increase by 1.2 magnitudes out to z \approx 2 (Omega = 0), 1 (Omega = 1), while our no-size evolution models predict an increase of only 0.5 magnitude. This result suggests that infall models could be important for explaining the 1.2-1.6 magnitude increase in surface brightness reported by Schade et al. (1995, 1996a, 1996b).Comment: 12 pages, LaTeX (aaspp4.sty), accepted by ApJ Letter

Cloning Dropouts: Implications for Galaxy Evolution at High Redshift

The evolution of high redshift galaxies in the two Hubble Deep Fields, HDF-N and HDF-S, is investigated using a cloning technique that replicates z~ 2-3 U dropouts to higher redshifts, allowing a comparison with the observed B and V dropouts at higher redshifts (z ~ 4-5). We treat each galaxy selected for replication as a set of pixels that are k-corrected to higher redshift, accounting for resampling, shot-noise, surface-brightness dimming, and the cosmological model. We find evidence for size evolution (a 1.7x increase) from z ~ 5 to z ~ 2.7 for flat geometries (Omega_M+Omega_LAMBDA=1.0). Simple scaling laws for this cosmology predict that size evolution goes as (1+z)^{-1}, consistent with our result. The UV luminosity density shows a similar increase (1.85x) from z ~ 5 to z ~ 2.7, with minimal evolution in the distribution of intrinsic colors for the dropout population. In general, these results indicate less evolution than was previously reported, and therefore a higher luminosity density at z ~ 4-5 (~ 50% higher) than other estimates. We argue the present technique is the preferred way to understand evolution across samples with differing selection functions, the most relevant differences here being the color cuts and surface brightness thresholds (e.g., due to the (1+z)^4 cosmic surface brightness dimming effect).Comment: 56 pages, 22 figures, accepted for publication in Ap

z~7 galaxy candidates from NICMOS observations over the HDF South and the CDF-S and HDF-N GOODS fields

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).Comment: 13 pages, 6 figures, 6 tables, accepted for publication in the Astrophysical Journal, replaced to match accepted version, see http://firstgalaxies.org/astronomers-area/ for a link to a complete reduction of the NICMOS observations over the two GOODS field

Very blue UV-continuum slopes of low luminosity z~7 galaxies from WFC3/IR: Evidence for extremely low metallicities?

We use the ultra-deep WFC3/IR data over the HUDF and the Early Release Science WFC3/IR data over the CDF-South GOODS field to quantify the broadband spectral properties of candidate star-forming galaxies at z~7. We determine the UV-continuum slope beta in these galaxies, and compare the slopes with galaxies at later times to measure the evolution in beta. For luminous L*(z=3) galaxies, we measure a mean UV-continuum slope beta of -2.0+/-0.2, which is comparable to the beta~-2 derived at similar luminosities at z~5-6. However, for the lower luminosity 0.1L*(z=3) galaxies, we measure a mean beta of -3.0+/-0.2. This is substantially bluer than is found for similar luminosity galaxies at z~4, just 800 Myr later, and even at z~5-6. In principle, the observed beta of -3.0 can be matched by a very young, dust-free stellar population, but when nebular emission is included the expected beta becomes >~-2.7. To produce these very blue beta's (i.e., beta~-3), extremely low metallicities and mechanisms to reduce the red nebular emission are likely required. For example, a large escape fraction (i.e., f_{esc}>~0.3) could minimize the contribution from this red nebular emission. If this is correct and the escape fraction in faint z~7 galaxies is >~0.3, it may help to explain how galaxies reionize the universe.Comment: 5 pages, 5 figures, accepted for publication in Astrophysical Journal Letter

The Rest Frame UV to Optical Colors and SEDs of z~4-7 Galaxies

We use the ultra-deep HUDF09 and the deep ERS data from the HST WFC3/IR camera, along with the wide area Spitzer/IRAC data from GOODS-S to derive SEDs of star-forming galaxies from the rest-frame UV to the optical over a wide luminosity range (M_1500 ~ -21 to M_1500 ~ -18) from z ~ 7 to z ~ 4. The sample contains ~ 400 z ~ 4, ~ 120 z ~ 5, ~ 60 z ~ 6, and 36 prior z ~ 7 galaxies. Median stacking enables the first comprehensive study of very faint high-z galaxies at multiple redshifts (e.g., [3.6] = 27.4 +/- 0.1 AB mag for the M_1500 ~ -18 sources at z ~ 4). At z ~ 4 our faint median-stacked SEDs reach to ~ 0.06 L*(z=4) and are combined with recently published results at high luminosity L > L* that extend to M_1500 ~ -23. We use the observed SEDs and template fits to derive rest frame UV-to-optical colors (U - V) at all redshifts and luminosities. We find that this color does not vary significantly with redshift at a fixed luminosity. The UV-to-optical color does show a weak trend with luminosity, becoming redder at higher luminosities. This is most likely due to dust. At z >~ 5 we find blue colors [3.6]-[4.5] ~ -0.3 mag that are most likely due to rest-frame optical emission lines contributing to the flux in the IRAC filter bandpasses. The scatter across our derived SEDs remains substantial, but the results are most consistent with a lack of any evolution in the SEDs with redshift at a given luminosity. The similarity of the SEDs suggests a self-similar mode of evolution over a timespan from 0.7 Gyr to 1.5 Gyr that encompasses very substantial growth in the stellar mass density in the universe (from ~ 4x10^6 to ~ 2x10^7 Msun Mpc^-3).Comment: 15 pages, 11 figures, 3 tables, submitted to Ap